Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,062,811 C→A 9.0% L189M (CTG→ATG) ‡ yihN → MFS transporter family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,062,8110CA 9.0% 363.6 / 20.1 132L189M (CTG→ATG) ‡yihNMFS transporter family protein
Reads supporting (aligned to +/- strand):  ref base C (62/57);  new base A (4/8);  total (67/65)
Fisher's exact test for biased strand distribution p-value = 2.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.63e-01

Only 100 of 136 total aligned reads displayed.

AAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGGCTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACTGTGGCCT  >  NC_000913/4062662‑4062959
                                                                                                                                                     |                                                                                                                                                     
aaaCAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTcc                                                                                                                                                       <  2:462121/150‑1 (MQ=255)
        cGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGAc                                                                                                                                               >  2:1342688/1‑150 (MQ=255)
         gTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACa                                                                                                                                              >  2:1717430/1‑150 (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:1484979/1‑149 (MQ=255)
               ttGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCATTCAGGACAaagaaa                                                                                                                                        <  1:439119/150‑2 (MQ=255)
                 ggcGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCATTCAGGACAaagaaaat                                                                                                                                      <  2:876206/148‑4 (MQ=255)
                 gTCGATACCGTGGTGGCATTTTCTGCGTTGACAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCATAGGCATTATTATTTTCTTTGTCATTCAGGACAaagaaaat                                                                                                                                      <  1:592606/150‑4 (MQ=255)
                    gATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCa                                                                                                                                   >  1:758403/1‑150 (MQ=255)
                    gATACCGTGGGGGGATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCa                                                                                                                                   <  2:281165/150‑1 (MQ=255)
                     aTACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCAc                                                                                                                                  >  2:1359719/1‑150 (MQ=255)
                      tACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCAcc                                                                                                                                 <  2:1386476/150‑1 (MQ=255)
                      tACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCAcc                                                                                                                                 >  2:866769/1‑150 (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                                                                                                                                         >  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)
                                          cGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGTAGAGGCACCGTCCGTTGAGGTG‑AAAAaag                                                                                                            <  1:428592/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:1717430/150‑1 (MQ=255)
                                               gCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGAc                                                                                                       <  2:1419008/150‑1 (MQ=255)
                                               gCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGAc                                                                                                       <  1:1342688/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)
                                                          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                                                                                        >  2:1135974/1‑150 (MQ=255)
                                                              tttGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGCATCGAAAAac                                                                                        <  1:549584/150‑1 (MQ=255)
                                                              tttGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGCATCGAAAAac                                                                                        >  2:339416/1‑150 (MQ=255)
                                                              tttGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGCATCGAAAAac                                                                                        <  2:625291/150‑1 (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)
                                                                      tGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGCATCGAAAAACACCAGTAt                                                                                >  1:1655881/1‑150 (MQ=255)
                                                                        gCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGa                                                                              <  2:1383192/150‑1 (MQ=255)
                                                                         cTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGCATCGAAAAACACCAGTGTGAc                                                                             <  2:711001/150‑1 (MQ=255)
                                                                         cTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGAc                                                                             <  1:54268/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)
                                                                                             gCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACaaaa                                                         <  1:281677/150‑1 (MQ=255)
                                                                                               aTCTGGTTCTATTACCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACATCGGTGCTGAAAGACAAAACt                                                       >  2:1081234/1‑150 (MQ=255)
                                                                                                  tGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTATCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGATAGACAAAACTAta                                                    >  2:705174/1‑149 (MQ=255)
                                                                                                    gTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTg                                                  <  1:1016819/150‑1 (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)
                                                                                                            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)
                                                                                                                                              cTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGGCTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACtg        >  1:669803/1‑150 (MQ=255)
                                                                                                                                              cTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGGCTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACtg        <  1:1568050/150‑1 (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)
                                                                                                                                                     |                                                                                                                                                     
AAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTG‑AAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGGCTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACTGTGGCCT  >  NC_000913/4062662‑4062959

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