Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,361,3370AC48.6% 171.2 / 19.8 110G26G (GGT→GGGrhmDL‑rhamnonate dehydratase
Reads supporting (aligned to +/- strand):  ref base A (33/23);  new base C (0/53);  total (33/77)
Fisher's exact test for biased strand distribution p-value = 3.82e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 112 total aligned reads displayed.

TGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATCCCTGG  >  NC_000913/2361188‑2361487
                                                                                                                                                      |                                                                                                                                                       
tGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCa                                                                                                                                                         >  1:207816/1‑150 (MQ=255)
    cGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCAccgcc                                                                                                                                                    <  1:1188018/151‑1 (MQ=255)
      ttttCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCAccgccg                                                                                                                                                   <  3:262272/150‑1 (MQ=255)
      ttttCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCACcgccgc                                                                                                                                                  >  3:264050/1‑151 (MQ=255)
       tttCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGACCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCACcgccgc                                                                                                                                                  <  1:1287885/150‑1 (MQ=255)
         tCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCa                                                                                                                                                <  2:264058/150‑1 (MQ=255)
           tGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCa                                                                                                                                                >  2:537231/1‑148 (MQ=255)
           tGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCa                                                                                                                                                <  3:537210/148‑1 (MQ=255)
            gCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAg                                                                                                                                               <  1:55277/148‑1 (MQ=255)
                  aCTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTAATCATCGGGGTGGCAATATGAACGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCttttt                                                                                                                                      >  1:922949/1‑151 (MQ=255)
                     tCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTct                                                                                                                                    <  4:774854/150‑1 (MQ=255)
                     tCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTct                                                                                                                                    <  1:611163/150‑1 (MQ=255)
                      caacaaCCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTctgc                                                                                                                                  >  1:648606/1‑151 (MQ=255)
                       aacaacCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTctgct                                                                                                                                 <  2:45372/151‑1 (MQ=255)
                       aacaacCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTctgc                                                                                                                                  <  4:137612/150‑1 (MQ=255)
                        acaacCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTctgctg                                                                                                                                <  3:256191/151‑1 (MQ=255)
                           accaacGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCg                                                                                                                            >  4:881761/1‑151 (MQ=255)
                           accaacGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCg                                                                                                                            <  1:935713/151‑1 (MQ=255)
                              aacGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCAc                                                                                                                         <  4:1495314/151‑1 (MQ=255)
                               acGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCAcc                                                                                                                        <  1:198329/151‑1 (MQ=255)
                                   gCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCa                                                                                                                    <  4:747793/151‑1 (MQ=255)
                                      aaGAACGTTAATGCCAAAGGACGGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGGGGCAATATGATCGTCTATCCAGTGATTCGCCCCCGGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTa                                                                                                                 <  4:685081/151‑1 (MQ=255)
                                         aaCGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAAc                                                                                                              >  3:323973/1‑151 (MQ=255)
                                          aCGTAAATGCCAAATGCCTGGCGTGACTGCTCGTAATCGCGGTATTTACTCAGCGGGGTGGCAATATGATCGTCTATCCAGTTATTCGCCCCCTGGTCGTTATAATCCCCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAAcc                                                                                                             <  4:1104853/151‑1 (MQ=255)
                                             ttAATGCCAAATGACTGGCGTGACTGCTCGTATCCGGGGTTTTTACTCGTCGGGGGGGCAATATGTCCGTCTACCCAGTGTTTCGCCCCCTGGTCGTAATAATCCCCGCCGCCAGCGCCTTTTTCTGCGGTCGCACCGCCAGTAAACCAg                                                                                                           <  1:1232532/150‑1 (MQ=255)
                                             ttAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAgg                                                                                                          >  4:1230441/1‑151 (MQ=255)
                                             ttAATGCAAAATGACTGGCGTTACTGCTCTTAATTGCGGTATTTTCTCATCGGGTTGGCAATATGACCTTCTATCCTTTATTTCGCCCCCTGTTCGAGATACTCCCCCCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAAcc                                                                                                             <  1:533000/148‑1 (MQ=255)
                                             ttAATCCCAAATGACTGGCGTGACTGCTCGTAATCGCGTTTTTTACTCATCGGGGTGGCATTATGATCGTCTATCCAGTAATTCGCCCCCTGGCCGTAATAATCCCCGCCGCCACCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCa                                                                                                            <  1:304178/149‑1 (MQ=255)
                                             taaaggcaaaaTGCCGGGGGTGCCGGCTCGAAACCGCGTTTTTTCCTCATCGGGGTGGAATTATGTTCGTCTTTCCATTGATTCCCCCCCTGGTCGTGAAAATCCCCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAg                                                                                                           <  1:1732367/142‑1 (MQ=255)
                                              tAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGgc                                                                                                         >  3:424283/1‑151 (MQ=255)
                                               aaggCCAAAGGACGGGCGTGACTGCTCGTAATCGCGTTATTTCCCCATGGGGGGGGAAATATGATCGTCTACCCAGTGATTCCCCCCCTGGTCGTAAAAATCCCCGCCGCCAGCGCTTTTTTCTGCTGTCGCACCGCCAGTAAACCAgg                                                                                                          <  1:979925/146‑1 (MQ=255)
                                               aaTGCAAAATGACTGGCGTGACGGCTCGAAATCGCGGTATTTCCCCATCGGGGTGGAAATATGATCGTCTACCCAGTGATTCGCCCCCGGGTCGTAATAACCCCCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGgcg                                                                                                        <  1:1613731/151‑1 (MQ=255)
                                               aaTGCAAAATGACGGGCGTGCCTGCTCGTAATCGCGGTATTTACTCATGGGGGTGGAAATATGATCGTCTACCCAGTGTTTCGCCCCCTGGTCGTAATAATCCCCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAgg                                                                                                          <  1:1523309/149‑1 (MQ=255)
                                                    cAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAAc                                                                                                   >  1:1560128/1‑151 (MQ=255)
                                                     aaaTGACTGGCGTGCCTGCTCGTAACCGCGTTATTTCCTCAGCGGGGGGGAAATATGATCGTCTACCCAGTGATTCCCCCCCTGGCCGGAATAACCCCCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAAc                                                                                                   <  1:1246929/150‑1 (MQ=255)
                                                     aaaTGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAAcc                                                                                                  >  4:1078782/1‑151 (MQ=255)
                                                     aaaTGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAAcc                                                                                                  >  2:373971/1‑151 (MQ=255)
                                                     aaaTGACGGGCGAGACTGCTCGTAAGCGGGTTATTTACTCTAGGGGGTGTAAATATGTTGGTCTATCCTTTGATTCGCCCCCTGGTCTGGAAAATCCCCGCCCCCAGCGCCTTTTTCTGCTGCCGCCCCGCCAGTAAACCAGgcgc                                                                                                       <  3:281699/146‑1 (MQ=255)
                                                      aaTGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCt                                                                                                 >  4:1126242/1‑151 (MQ=255)
                                                       aTGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTg                                                                                                >  3:367558/1‑151 (MQ=255)
                                                         gccTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTAATAACCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGtt                                                                                              <  4:1130238/149‑1 (MQ=255)
                                                         gACTGGCGTGACTGCTTGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCCCCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGtt                                                                                              <  4:1895893/151‑1 (MQ=255)
                                                         gACTGGCGTGACTGCTCTTAATCGCGGTTTTTACTCAGCGGGGTGCCAATATGCCCGTCTACCCAGTGTTTCGCCCCCTGGTGGTGATACTCCCCGCCGCCAGCGCCTTTTTCTGCTGTCGCGCCGCCAGTAAACCAGGCGCGAACCTGt                                                                                               <  3:62257/150‑1 (MQ=255)
                                                         gACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGtt                                                                                              >  1:1622740/1‑151 (MQ=255)
                                                           cTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGttt                                                                                             >  1:1679660/1‑150 (MQ=255)
                                                            tGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGt                                                                                                                  >  2:87640/1‑128 (MQ=255)
                                                             ggCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGTTCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCCCCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGt                                                                                                                  <  3:87640/127‑1 (MQ=255)
                                                              gCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAAtt                                                                                         >  4:811582/1‑151 (MQ=255)
                                                                gTGACTGCTCGTAATCGCGGTATTTCCTCATCGGGGTGGAAATATGATCGTCTATCCAGTGTTTCGCCCCCTGGTCGTGATAATCCCCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAAtt                                                                                         <  4:1669311/149‑1 (MQ=255)
                                                                 tgccTGCTCGTAATCGCGGTATTTCCCCATCGGGGTGGAAATTTGATCGTCTACCCAGTGATTCCCCCCCTGTGCGAGATAAGCCCCGCCGCCAGCGCTTTTTTCGGCGGTCGCACCGCCAGTAAACCAGGCGCGAACCTGttt                                                                                             <  3:59460/141‑1 (MQ=255)
                                                                 tGACTGCTCGAAACCGCGGTTTTTACCCAGCGGGGGGGCAATATGATCGTCTACCCAGTGATTCCCCCCCGGGTCTTGATAATCCCCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAAtttt                                                                                       <  1:1176961/150‑1 (MQ=255)
                                                                  gACTGCTCGTAATCGCGTTTTTTCCTCAGCGGGGTGCAAATATGCTCGTCTTTCCAGTGATTCGCCCCCGGGTCGTGAAAATCCCCGCCGCCAGCGCTTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAAttttt                                                                                      <  1:44599/150‑1 (MQ=255)
                                                                  gACTGCCCGTAATCGCGGTATTTACCGACGGGGGCGACAAAATGTTCGTCTACCCAGTGATCCGCCCCTGGGTCTTAATAATCGCCGCCGCCAGCGCCTTTTTCTTCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAAttttt                                                                                      <  1:897285/150‑1 (MQ=255)
                                                                     tGCTCGTAATCGCGGTATTTCCTCATCGGGGTGGCAATATGATCGTCTATCCGGTGATTCGCCCCCTGTTCGTGATAATCCCCGCCGCCAGCGCCTTTTTCTGCTGTCGCGCCGCCAGTAAACCAGGCGCGAACCTGTTTAAttttt                                                                                      <  4:720484/147‑1 (MQ=255)
                                                                      gCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGt                                                                                   >  4:450671/1‑149 (MQ=255)
                                                                       ctctTAATCGCGTTTTTTACTCATGGGGGTGGCAATATGATCGTCTACCCAGTGTTTCGCCCCCTGGTCGTGATAATCCCCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAg                                                                                 <  1:1027338/150‑1 (MQ=255)
                                                                       cTCGTAATCGCGGTATTTCCTCATCGGGGTGGAAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCCCCGCCGCCAGCGCTTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAAtttt                                                                                       <  1:650077/144‑1 (MQ=255)
                                                                         cGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGt                                                                              >  1:1924670/1‑151 (MQ=255)
                                                                          gTAATCGCGGTATTTACTCATGGGGGTGGAAATATGATCGTCTATCCAGTGTTTCGCCCCCTGGTCGTGATAATCCCCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTc                                                                             <  4:1675559/151‑1 (MQ=255)
                                                                            actcgggtttttttCTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCGGTTCGTTATAACCCCCGCCGCCAGCCCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTc                                                                             <  1:1229921/139‑1 (MQ=255)
                                                                            aaTCGGGGTATTTCCCCATCGGGGTGGCAATTTGTTCGTCTCCCCAGTGATTCCCCCCTTGGTCGTGATAACCCCCGCCCCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCat                                                                           <  3:81702/151‑1 (MQ=255)
                                                                            aaTCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCat                                                                           >  4:1447828/1‑151 (MQ=255)
                                                                               cgcgGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCatgat                                                                        >  1:1122450/1‑151 (MQ=255)
                                                                                gcgGTTTTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCCCCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCatgatg                                                                       <  4:1573579/151‑1 (MQ=255)
                                                                                gcgGTATTTACTCATCGGGGTGGAAATATGATCGTCTATCCAGTGATTCCCCCCCTGGTCGTGATAATCCCCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCatga                                                                         <  4:589500/149‑1 (MQ=255)
                                                                                  gttttttACTCATGGGGGGGGAAATATGATCGTCTTCCCAGTGATCCCCCCCCTGGTCGTGAAAACCCCCGCCGCCAGCGCTTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCatg                                                                          <  1:1799867/142‑1 (MQ=255)
                                                                                  ggTATTTACTCATGGGGGTGGAAATATGATCGTCTATCCAGTGATTCGCCCCCGGGTCGTGAAAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGt                                                                      <  2:277583/150‑1 (MQ=255)
                                                                                          cTCATGGGGGTGGAAATATGTCCTTCTATCCAGTGTTTCCCCCCCTGGTCGTGATAACCCCCGCCGCCAGCGCTTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCa                                                                <  1:1921277/148‑1 (MQ=255)
                                                                                           tCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGt                                                            >  4:369771/1‑151 (MQ=255)
                                                                                            cATCGGGGTGGCAATATGATCGTCTTCCCAGTGATTCGCCCCCTGGTCGTAATAACCCCCGCCGCCAGCGCCTTTTTCTGCTGTCGCGCCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGt                                                            <  1:322133/150‑1 (MQ=255)
                                                                                                gggtTGGAAATATGTTCGTCTACCAGTTGTTCCCCCCCTTGGTGGTAAAAACCCCCGCCGCCAGCGCCTTTTTCTGTTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGa                                                       <  1:541736/151‑1 (MQ=255)
                                                                                                ggggggGCAATATGCTCGTTTATCCAGTGATCCGCCCCCTGGGGGTAATAATCCCCGCCGCCAGCGTTTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATg                                                        <  3:92867/150‑1 (MQ=255)
                                                                                                 gggTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATg                                                        >  3:355057/1‑149 (MQ=255)
                                                                                                     ggaaaTATGATCGTCTACCCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCTTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCtt                                                  <  3:516917/148‑1 (MQ=255)
                                                                                                        aaTATGTTCTTCTCTCCAGTTTTTCCCCCCCTGGTCTTAATAATCCCCGCCGCCAGCGCTTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTa                                               <  3:64138/151‑1 (MQ=255)
                                                                                                           atGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGt                                            >  1:1624903/1‑151 (MQ=255)
                                                                                                             gATCGTCTATCCAGTGATTCCCCCCCGGGTCTTGATAACCCCCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTc                                           <  4:824409/150‑1 (MQ=255)
                                                                                                             gATCGTCTATCAAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGt                                            <  4:1746847/149‑1 (MQ=255)
                                                                                                                agTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCaa                                         >  4:850707/2‑149 (MQ=255)
                                                                                                                    taccctgtgtTTCGCCCCGTGGGCGTGAAAATCACCGCCACGAGGCAATTTTTCGGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGa                                    <  4:1623376/144‑1 (MQ=255)
                                                                                                                      tCCAGTGATCCGCCCCCGGGTCGTGAAAATCCCCGCCGCCAGCCCTTTTTTCTGCTGCCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAggg                                     <  2:264844/147‑1 (MQ=255)
                                                                                                                       ccggTGATTCGCCCCCTGGTCGTGATAATCCCCGCCGCCAGCGCTTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACttt                                <  1:431286/148‑1 (MQ=255)
                                                                                                                       ccAGTGTTTCGCCCCCTGGTCGTGATAATCCCCGCCGCCAGCGCTTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACt                                  <  1:805363/149‑1 (MQ=255)
                                                                                                                       ccAGTGATTCGCCCCCTGGTCGTGATAATCCCCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACttt                                <  4:283536/151‑1 (MQ=255)
                                                                                                                          gtgtTTCCCCCCCGGGTCGTGATAATCCCCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACt                                  <  4:1695723/146‑1 (MQ=255)
                                                                                                                            gATTCGCCCCCTGGTCTTAATAACCCCCGCCGCCAGCCCTTTTTTCGGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGAc                                   <  1:1213929/143‑1 (MQ=255)
                                                                                                                               tCGCCCCCTGGTCGTGATAATCCCCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATcc                         <  4:34326/150‑1 (MQ=255)
                                                                                                                                  cccctTGTTCGTGATAATCCCCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCa                       <  4:443542/149‑1 (MQ=255)
                                                                                                                                    cccTGGCCGTGATAATCCCCGCCGCCAGCGCTTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAg                   <  1:543468/151‑1 (MQ=255)
                                                                                                                                     ccTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAAc                     >  4:1635708/1‑148 (MQ=255)
                                                                                                                                      cTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGAc                 >  4:377705/1‑151 (MQ=255)
                                                                                                                                         ggTCGTGATAATCCCCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGAc                 <  2:276642/149‑1 (MQ=255)
                                                                                                                                         ggTCGTGAAAATCCCCGCCGCCAGCGCTTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACg                <  1:1455836/150‑1 (MQ=255)
                                                                                                                                         ggTCGTAATAACCCCCGCCGCCAGCGCTTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGaa                 <  1:686448/149‑2 (MQ=255)
                                                                                                                                          gTCGTTAAAATCCCCGCCGCAAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGAc                 <  4:1543191/148‑1 (MQ=255)
                                                                                                                                           tCGTGATAACCCCCGCCGCCAGCGCTTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACg                <  1:1802035/148‑1 (MQ=255)
                                                                                                                                             gTGATAATCCCCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGa             <  2:425543/149‑1 (MQ=255)
                                                                                                                                              taaaaaTCCCCCCCCCCAGCGCCTTTCTCGGTGTGCGCCCCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGaaa           <  1:1362569/146‑1 (MQ=255)
                                                                                                                                                aTAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAatat        >  1:1273369/1‑151 (MQ=255)
                                                                                                                                                      aCCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATCCCTgg  >  3:279832/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                       
TGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATCCCTGG  >  NC_000913/2361188‑2361487

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