Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,062,814 A→C 8.1% N190H (AAT→CAT) ‡ yihN → MFS transporter family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,062,8140AC 8.1% 398.9 / 21.5 135N190H (AAT→CAT) ‡yihNMFS transporter family protein
Reads supporting (aligned to +/- strand):  ref base A (64/60);  new base C (4/7);  total (68/67)
Fisher's exact test for biased strand distribution p-value = 3.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.92e-01

Only 100 of 136 total aligned reads displayed.

CAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGGCTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACTGTGGCCTGAC  >  NC_000913/4062665‑4062962
                                                                                                                                                     |                                                                                                                                                     
     cGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGAc                                                                                                                                                  >  2:1342688/1‑150 (MQ=255)
     cGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTAAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGAc                                                                                                                                                  >  1:526625/1‑150 (MQ=255)
      gTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACa                                                                                                                                                 >  2:1717430/1‑150 (MQ=255)
       ttgccttGTCGCTCCCGGGGTGGCATTTTCTGGCTTGGCAGTCTGTCCCTGGTGTTTCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATACCCTTATTGTGATTGCCGTAGGCATTATTAATATCTTTGTCATGAATGACaa                                                                                                                                                <  1:280571/145‑1 (MQ=255)
         gCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCATTCAGGACAaag                                                                                                                                              <  1:1211295/150‑1 (MQ=255)
          cATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTTCCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAaaga                                                                                                                                             >  2:1001037/1‑150 (MQ=255)
          cATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAaaga                                                                                                                                             <  2:1447001/150‑1 (MQ=255)
          cATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAaaga                                                                                                                                             >  1:1505722/1‑150 (MQ=255)
            ttGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAaagaag                                                                                                                                           <  2:647424/150‑1 (MQ=255)
            ttGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCATTCAGGACAaagaaa                                                                                                                                           <  1:439119/150‑2 (MQ=255)
              ggcGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCATTCAGGACAaagaaaat                                                                                                                                         <  2:876206/148‑4 (MQ=255)
              gTCGATACCGTGGTGGCATTTTCTGCGTTGACAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCATAGGCATTATTATTTTCTTTGTCATTCAGGACAaagaaaat                                                                                                                                         <  1:592606/150‑4 (MQ=255)
                 gATACCGTGGGGGGATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCa                                                                                                                                      <  2:281165/150‑1 (MQ=255)
                  aTACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCAc                                                                                                                                     >  2:1359719/1‑150 (MQ=255)
                   tACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCAcc                                                                                                                                    >  2:866769/1‑150 (MQ=255)
                   tACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCAcc                                                                                                                                    <  2:1386476/150‑1 (MQ=255)
                     ccGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCAccgt                                                                                                                                  >  2:1712992/1‑150 (MQ=255)
                      cGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCAccgtc                                                                                                                                 >  1:19443/1‑150 (MQ=255)
                          gtggCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGtt                                                                                                                             >  2:1202140/1‑150 (MQ=255)
                            ggCTTTTTCTGAGTTGGCAGTATTTAGCTGGTTTGGGTGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCGCTTATTGTGATTGCCGTAGGCTTTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGttta                                                                                                                           <  2:428655/150‑3 (MQ=255)
                             gCATTTTCTGCGTTGGCAGTATTTACCAGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAg                                                                                                                          <  2:870555/150‑1 (MQ=255)
                              cATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAgg                                                                                                                         <  2:1509170/150‑1 (MQ=255)
                              cATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCATTCAGGACAaagaa                                                                                                                                            <  1:1517276/131‑1 (MQ=255)
                              cATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCATTCAGGACAaagaa                                                                                                                                            >  2:1517276/1‑131 (MQ=255)
                                 tttCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑a                                                                                                                     <  2:215999/150‑1 (MQ=255)
                                    cTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCAccgtc                                                                                                                                 <  1:463900/136‑1 (MQ=255)
                                    cTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑aaaa                                                                                                                  >  1:1446461/1‑150 (MQ=255)
                                     tGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑aaaaa                                                                                                                 >  1:680980/1‑150 (MQ=255)
                                       cGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTCTGTTTTAAAGCAGGCATCTGGTTCTATTACCTTATTGTGATTGCCGTAGGAATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAaag                                                                                                               <  2:661015/150‑1 (MQ=255)
                                        gTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAaaga                                                                                                              <  2:1290228/150‑1 (MQ=255)
                                        gTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAaaga                                                                                                              <  1:1235678/150‑1 (MQ=255)
                                           ggCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAagaaga                                                                                                           >  2:549584/1‑150 (MQ=255)
                                            gCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGAc                                                                                                          <  1:1342688/150‑1 (MQ=255)
                                            gCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGAc                                                                                                          <  2:1419008/150‑1 (MQ=255)
                                            gCAGGATTTACCTGGGTTGGCAGTGGCTTATTAGGGTTTAAAGCAGGCAGCTGGTTCTATTCCC‑GATTGTGGTTTGCCTCCGCATTTTTTTTTTCTTTGTCGTGAATGCCCAAGAAGAGGCAACCGCCGGTGGGGGG‑AAATatgaattcg                                                                                                         >  2:212464/1‑146 (MQ=255)
                                                tATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAg                                                                                                      <  2:969935/150‑1 (MQ=255)
                                                 aTTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGc                                                                                                     >  2:781240/1‑150 (MQ=255)
                                                  tttACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTAGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGCa                                                                                                    >  2:1069927/1‑150 (MQ=255)
                                                     aCCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGCATCg                                                                                                 >  1:625291/1‑150 (MQ=255)
                                                       cTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGCATCGaa                                                                                               >  2:944970/1‑150 (MQ=255)
                                                        tGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGCATCGaaa                                                                                              >  1:958479/1‑150 (MQ=255)
                                                        tGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGAACCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGCATCGaaa                                                                                              >  2:824321/1‑150 (MQ=255)
                                                        tGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGACGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAACAAGCCGGAGCTTCGaaa                                                                                              >  2:240358/1‑150 (MQ=255)
                                                         ggTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGCATCGaaaa                                                                                             >  2:786538/1‑150 (MQ=255)
                                                          gTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGCATCGaaaaa                                                                                            >  2:943536/1‑150 (MQ=255)
                                                           tttGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGCATCGAAAAac                                                                                           <  1:549584/150‑1 (MQ=255)
                                                           tttGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGCATCGAAAAac                                                                                           <  1:723128/150‑1 (MQ=255)
                                                           tttGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGCATCGAAAAac                                                                                           <  2:625291/150‑1 (MQ=255)
                                                           tttGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGCATCGAAAAac                                                                                           >  2:1135974/1‑150 (MQ=255)
                                                           tttGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGCATCGAAAAac                                                                                           >  2:339416/1‑150 (MQ=255)
                                                           tttGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGGG‑AAAAAAGAAGACGGAGCATCGAAAAac                                                                                           <  2:1505722/150‑1 (MQ=255)
                                                           tttGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACAGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGCATCGAAAAac                                                                                           >  1:846929/1‑150 (MQ=255)
                                                              ggCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGCATCGAAAAACAcc                                                                                        >  1:792146/1‑150 (MQ=255)
                                                                cAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGCATCGAAAAACACCAg                                                                                      <  1:1295374/150‑1 (MQ=255)
                                                                  gTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGCATCGAAAAACACCAGTa                                                                                    <  2:758403/150‑1 (MQ=255)
                                                                     gCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGa                                                                                 <  2:1383192/150‑1 (MQ=255)
                                                                      cTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGCATCGAAAAACACCAGTGTGAc                                                                                <  2:711001/150‑1 (MQ=255)
                                                                       ttattaGGTTTTAAAGCTTGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCTCCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGAcc                                                                               <  2:1577617/150‑1 (MQ=255)
                                                                       ttattaGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGAcc                                                                               >  1:1570447/1‑150 (MQ=255)
                                                                         attaGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTc                                                                             <  2:846929/150‑1 (MQ=255)
                                                                           tacgTTTTAAAGCAGGCATCTGGTTCTATACCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCgg                                                                           >  2:795594/4‑150 (MQ=255)
                                                                           taGGTTTTTAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCgg                                                                           >  2:1607572/1‑150 (MQ=255)
                                                                           taGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCgg                                                                           >  1:1388464/1‑150 (MQ=255)
                                                                           taGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCgg                                                                           <  1:1001037/150‑1 (MQ=255)
                                                                             gggtttAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTg                                                                         <  2:1225633/147‑1 (MQ=255)
                                                                             ggTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTg                                                                         >  1:368112/1‑150 (MQ=255)
                                                                              gTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAATGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGc                                                                        >  1:1097888/1‑150 (MQ=255)
                                                                                tttAAAGCAGGCATCTGGTTCTATTCTCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTg                                                                      >  2:1448180/1‑150 (MQ=255)
                                                                                 ttAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGa                                                                     <  2:19443/150‑1 (MQ=255)
                                                                                  tAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGa                                                                     <  2:673548/150‑1 (MQ=255)
                                                                                  tAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGaa                                                                    >  2:1267432/1‑150 (MQ=255)
                                                                                    aaGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAg                                                                  <  2:973584/150‑1 (MQ=255)
                                                                                         ggCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACaaa                                                             <  2:463959/150‑1 (MQ=255)
                                                                                           cATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAAc                                                           >  1:426498/1‑150 (MQ=255)
                                                                                            aTCTGGTTCTATTACCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACATCGGTGCTGAAAGACAAAACt                                                          >  2:1081234/1‑150 (MQ=255)
                                                                                               tGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTATCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGATAGACAAAACTAta                                                       >  2:705174/1‑149 (MQ=255)
                                                                                                 gTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTg                                                     <  1:1369465/150‑1 (MQ=255)
                                                                                                 gTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTg                                                     >  1:1588387/1‑150 (MQ=255)
                                                                                                 gTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTg                                                     >  1:1667167/1‑150 (MQ=255)
                                                                                                 gTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTg                                                     <  1:1016819/150‑1 (MQ=255)
                                                                                                        tCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGGCTTATc                                              <  1:1615805/150‑1 (MQ=255)
                                                                                                         cccTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGGCTTATCg                                             >  2:1686363/1‑150 (MQ=255)
                                                                                                             tATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGGCTTATCGCttt                                         <  1:1359719/150‑1 (MQ=255)
                                                                                                              aTTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGGCTTATCGCtttt                                        >  2:1631137/1‑150 (MQ=255)
                                                                                                                   gATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGTCAAAACTATCTGGCTTATCGCTTTTAACGt                                   <  2:1364806/150‑1 (MQ=255)
                                                                                                                    aTTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGATGTG‑AAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGGCTTATCGCTTTTAACGtc                                  >  1:224182/1‑150 (MQ=255)
                                                                                                                          gTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGGCTTATCGCTTTTAACGTCttcttc                            <  2:1627563/150‑1 (MQ=255)
                                                                                                                           tAGTCATTTTTTTTTTATTTGTCCTGAATGACATAGAAGTGGCACCGTCGTTTGAGGTG‑AAAATAGTAGACGGAGCATCGTAAAACACCAGTATGTCCTCGGTGCTGATAGACAAAACTATCTGGCTTATCTCTTTTAACGTCTTCTTCg                           <  2:908732/150‑1 (MQ=255)
                                                                                                                             ggCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGGCTTATCGCTTTTAACGTCTTCTTCGtt                         <  1:795594/150‑1 (MQ=255)
                                                                                                                               cATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGGCTTATCGCTTTTAACGTCTTCTTCGTTTa                       >  1:893116/1‑150 (MQ=255)
                                                                                                                                attattattTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGGCTTATCGCTTTTAACGTCTTCTTCGTTTAc                      >  1:1697546/1‑150 (MQ=255)
                                                                                                                                  tattattTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGGCTTATCGCTTTTAACGTCTTCTTCGTTTAcgc                    >  1:823778/1‑150 (MQ=255)
                                                                                                                                     tattTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGGCTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGt                 <  1:1202140/150‑1 (MQ=255)
                                                                                                                                     cctcttCTTTGTCATTCAGGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGGCTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGt                 >  1:1715750/5‑150 (MQ=255)
                                                                                                                                            tttGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGGCTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACtga          >  1:9244/1‑149 (MQ=255)
                                                                                                                                              tGTCCTGAATGACAAAGAAGATGCACCTTCCTTTGAGGTG‑AAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGGCTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACTGTgg        <  1:705067/150‑1 (MQ=255)
                                                                                                                                                tgcTGAATGACAAAGAAGAGGATCCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGGCTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACTGTGGcc      <  2:1655881/148‑1 (MQ=255)
                                                                                                                                                 ccTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGGCTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACTGTGGCCt     <  2:1446461/150‑1 (MQ=255)
                                                                                                                                                    gAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGGCTTATCGCTTTTAACGTCTTCTTCGTTTCCGCGGTTTACTGTGGCCTGAc  >  1:79048/1‑150 (MQ=255)
                                                                                                                                                    gAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGGCTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACTGTGGCCTGAc  >  2:459332/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGGCTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACTGTGGCCTGAC  >  NC_000913/4062665‑4062962

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑