Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,231,4440TG42.5% 178.4 / 12.8 94G142G (GGT→GGGyohKLrgB family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (12/42);  new base G (40/0);  total (52/42)
Fisher's exact test for biased strand distribution p-value = 3.90e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACCATCTGTTTTATCGGCAGCGTGGTTGCGATGGTGACGGGGACTTCCGTGGCATTGTTGATGGGCGCTTCACCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGTATTCCGGCAATCAGCGCAGTTTGCGTGATTTTCGTCGGCATCCTCGGCGCGGTATTTGGCCATACATTGCTTAATGCGATGCGTATTCGTACCAAAGCTGCGCGCGGTCTGGCAATGGGGACTGCCTCGCACGCCCTCGGTACGGCGCG  >  NC_000913/2231294‑2231593
                                                                                                                                                       |                                                                                                                                                          
caCCATCTGTTTTATCGGCAGCGTGGTTGCGATGGTGACGGGGACTTCCGTGGCATTGTTGATGGGCGCTTCACCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGt                                                                                                                                                            <  4:169864/151‑1 (MQ=255)
   cATCTGTTTTATCGGCAGCGTGGTTGCGATGGTGACGGGGACTTCCGTGGCATTGTTGATGGGCGCTTCACCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGTAtt                                                                                                                                                        <  4:126277/151‑1 (MQ=255)
       tGTTTTATCGGCAGCGTGGTTGCGATGGTGACGGGGACTTCCGTGGCATTGTTGATGGGCGCTTCACCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGGAtttccg                                                                                                                                                    >  1:1230702/1‑147 (MQ=255)
        gTTTTATCGGCAGCGTGGTTGCGATGGTGACGGGGACTTCCGTGGCATTGTTGATGGGCGCTTCACCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGTATTcc                                                                                                                                                      <  4:631086/148‑1 (MQ=255)
        gTTTTATCGGCAGCGTGGTTGCGATGGTGACGGGGACTTCCGTGGCATTGTTGATGGGCGCTTCACCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGTATTcc                                                                                                                                                      <  4:431826/148‑1 (MQ=255)
        gTTTTATCGGCAGCGTGGTTGCGATGGTGACGGGGACTTCCGTGGCATTGTTGATGGGCGCTTCACCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGGAtttccg                                                                                                                                                    >  1:224447/1‑146 (MQ=255)
          tttATCGGCAGCGTGGTTGCGATGGTGACGGGGACTTCCGTGGCATTGTTGATGGGCGCTTCACCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGTATTCCGGCa                                                                                                                                                  <  4:1347732/150‑1 (MQ=255)
           ttATCGGCAGCGTGGTTGCGATGGTGACGGGGACTTCCGTGGCATTGTTGATGGGCGCTTCACCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGGATTTCGGGaa                                                                                                                                                 >  1:667547/1‑150 (MQ=255)
               cGGCAGCGTGGTTGCGATGGTGACGGGGACTTCCGTGGCATTGTTGATGGGCGCTTCACCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGGAtttccggaa                                                                                                                                                 >  1:42491/1‑139 (MQ=255)
                 gCAGCGTGGTTGCGATGGTGACGGGGACTTCCGTGGCATTGTTGATGGGCGCTTCACCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGGATTGGCGGGAATTCGGcgaaacgcg                                                                                                                                           >  2:100805/1‑142 (MQ=255)
                    gCGTGGTTGCGATGGTGACGGGGACTTCCGTGGCATTGTTGATGGGCGCTTCACCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGTATTCCGGCAATCAGCGCAGt                                                                                                                                       <  4:1361529/151‑1 (MQ=255)
                      gTGGTTGCGATGGTGACGGGGACTTCCGTGGCATTGTTGATGGGCGCTTCACCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGGTTTTCCGCCATAAgcgccggtt                                                                                                                                     >  4:154645/1‑146 (MQ=255)
                         gTTGCGATGGTGACGGGGACTTCCGTGGCATTGTTGATGGGCGCTTCACCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGTATTCCGGCAATCAGCGCAGTTTGCg                                                                                                                                  <  1:502357/151‑1 (MQ=255)
                           tGCGATGGTGACGGGGACTTCCGTGGCATTGTTGATGGGCGCTTCACCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGTATTTCCGCAAACAGCGCAGGTTTCGTg                                                                                                                               >  1:697455/1‑151 (MQ=255)
                           tGCGATGGTGACGGGGACTTCCGTGGCATTGTTGATGGGCGCTTCACCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGTATTCCGGCAATCAGCGCAGTTTGCGt                                                                                                                                 <  4:169489/150‑1 (MQ=255)
                              gATGGTGACGGGGACTTCCGTGGCATTGTTGATGGGCGCTTCACCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGGATTTCCGGCAACAgcgccggtt                                                                                                                                     >  1:159121/1‑138 (MQ=255)
                              gATGGTGACGGGGACTTCCGTGGCATTGTTGATGGGCGCTTCACCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGGATTTCCGCCAACAGCGCAGGTTTCggggtt                                                                                                                             >  1:1130782/1‑146 (MQ=255)
                              gATGGTGACGGGGACTTCCGTGGCATTGTTGATGGGCGCTTCACCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGGATTTCCGCAAACAGCGCCGGTTGCggggtt                                                                                                                             >  1:1273402/1‑146 (MQ=255)
                               aTGGTGACGGGGACTTCCGTGGCATTGTTGATGGGCGCTTCACCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGGATTTCCGCAAACAGCGCCGGTTggggggttt                                                                                                                            >  1:20828/1‑143 (MQ=255)
                                  gTGACGGGGACTTCCGTGGCATTGTTGATGGGCGCTTCACCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGTATTCCGGCAATCAGCGCAGTTTGCGTGATTTtcg                                                                                                                        <  4:1205706/151‑1 (MQ=255)
                                    gACGGGGACTTCCGTGGCATTGTTGATGGGCGCTTCACCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGGATTGGCGGGATTTCCGCCAACAGCGCCGGTTGCGTGGTTTttgt                                                                                                                       >  4:159419/1‑150 (MQ=255)
                                     aCGGGGACTTCCGTGGCATTGTTGATGGGCGCTTCACCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGTATTCCGGCAATCAGCGCAGTTTGCGTGATTTtcgtcg                                                                                                                     <  4:1135996/151‑1 (MQ=255)
                                     aCGGGGACTTCCGTGGCATTGTTGATGGGCGCTTCACCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGGATTTCCGGCATCAGCGCCGGTTTCGGGGtttt                                                                                                                          >  1:33998/1‑146 (MQ=255)
                                          gACTTCCGTGGCATTGTTGATGGGCGCTTCACCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGTATTCCGGCAATCAGCGCAGTTTGCGTGATTTTCGTCGGCATc                                                                                                                <  1:728051/151‑1 (MQ=255)
                                           aCTTCCGTGGCATTGTTGATGGGCGCTTCACCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGGATTTCCGCAAACAGCGCCGTTTGCGGGATTTTTGTCCGCAAcc                                                                                                               >  1:983409/1‑151 (MQ=255)
                                             ttCCGTGGCATTGTTGATGGGCGCTTCACCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGGATTTCCGCCAACAACGCCGGTTTCGGGGTTTTTGTCGGCATcccc                                                                                                             >  2:49340/1‑149 (MQ=255)
                                              tCCGTGGCATTGTTGATGGGCGCTTCACCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGGATTGGCGGGATTTCGGCCAACAACGCCGGTTTCGGGGttttttttgggctcc                                                                                                               >  4:504907/1‑137 (MQ=255)
                                                cGTGGCATTGTTGATGGGCGCTTCACCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGTATTCCGGCAATCAGCGCAGTTTGCGTGATTTTCGTCGGCATCCTCGgc                                                                                                          <  4:1123154/151‑1 (MQ=255)
                                                 gTGGCATTGTTGATGGGCGCTTCACCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGGATTTCCGCAAACAACGGCGGTTGGGTGGTTTTTGTCGGGATCCTccggg                                                                                                         >  1:44164/1‑147 (MQ=255)
                                                   ggCATTGTTGATGGGCGCTTCACCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGTATTCCGGCAATCAGCGCAGTTTGCGTGATTTTCGTCGGCATCCTCGGcgcg                                                                                                     <  1:810837/151‑1 (MQ=255)
                                                      aTTGTTGATGGGCGCTTCACCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGTATTCCGGCAATCAGCGCAGTTTGCGTGATTTTCGTCGGCATCCTCGGCGCGGTa                                                                                                  <  1:1209982/151‑1 (MQ=255)
                                                       ttgttgATGGGCGCTTCACCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTTGCGGGATTTCCGCCAAACACGCAGGTTTCGTGGGTTTTCTTGGGCATCCCCgggcggg                                                                                                      >  4:542619/1‑144 (MQ=255)
                                                       ttgttgATGGGCGCTTCACCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGGATTTCCGCCATCAGCGCCGGTTGGGGGGTTTTTTTCGGGATCCTCGGCGCgg                                                                                                    >  1:128896/1‑148 (MQ=255)
                                                        tgttgATGGGCGCTTCACCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTTGCGGGATTCCCGCCATCAACGCCGTTTTCGTGGTTTTTTTCGGCGTCCCCGGCGCGGGAtt                                                                                                >  4:1291277/1‑151 (MQ=255)
                                                        tgttgATGGGCGCTTCACCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGTATTCCGGCAATCAGCGCAGTTTGCGTGATTTTCGTCGGCATCCTCGGCGCGGTAt                                                                                                 <  1:574515/150‑1 (MQ=255)
                                                        tgttgATGGGCGCTTCACCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGGATTTCCGCCAATAGCGCCGGTTTCGGGGTTTTTGTCGGCAACCTCCGCGCGGGAtt                                                                                                >  1:63851/1‑151 (MQ=255)
                                                         gttgATGGGCGCTTCACCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGTATTCCGGCAATCAGCGCAGTTTGCGTGATTTTCGTCGGCATCCTCGGCGCGGTAtt                                                                                                <  1:837540/150‑1 (MQ=255)
                                                           tgATGGGCGCTTCACCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGTATTCCGGCAATCAGCGCAGTTTGCGTGATTTTCGTCGGCATCCTCGGCGCGGTATTTg                                                                                              <  2:49748/150‑1 (MQ=255)
                                                            gATGGGCGCTTCACCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGGATTGGCGGGGTTCCCGCCCTCAGCGCCGGTTGCGTGGTTTTTGTTGGCCTCCTCCGCGCGGGAttttg                                                                                             >  4:125103/1‑148 (MQ=255)
                                                             aTGGGCGCTTCACCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGTATTCCGGCAATCAGCGCAGTTTGCGTGATTTTCGTCGGCATCCTCGGCGCGGTATTTGGcc                                                                                           <  1:1337834/151‑1 (MQ=255)
                                                               gggCGCTTCACCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGTATTCCGGCAATCAGCGCAGTTTGCGTGATTTTCGTCGGCATCCTCGGCGCGGTATTTGGCCa                                                                                          <  1:767679/150‑1 (MQ=255)
                                                                  cgcTTCACCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGTATTCCGGCAATCAGCGCAGTTTGCGTGATTTTCGTCGGCATCCTCGGCGCGGTATTTGGCCATACa                                                                                      <  1:1187373/151‑1 (MQ=255)
                                                                   gcTTCACCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGTATTCCGGCAATCAGCGCAGTTTGCGTGATTTTCGTCGGCATCCTCGGCGCGGTATTTGGCCATACAt                                                                                    <  1:1237779/151‑1 (MQ=255)
                                                                   gcTTCACCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGGATTCCCGCAAACAGCGCCGGTTGCGGGGTTTTTTTCGGCCTCCTCCGCGCGGGATTTTGcc                                                                                           >  1:867572/1‑145 (MQ=255)
                                                                   gcTTCACCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGGATTCCCGCAAACAACGCCGGTTTTGGGGTTTTTGTCCGCATCCCCCGGGCGGGAttttg                                                                                             >  1:424349/1‑141 (MQ=255)
                                                                      tCACCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGTATTCCGGCAATCAGCGCAGTTTGCGTGATTTTCGTCGGCATCCTCGGCGCGGTATTTGGCCATACAtt                                                                                   <  4:1145970/149‑1 (MQ=255)
                                                                       cACCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGTATTCCGGCAATCAGCGCAGTTTGCGTGATTTTCGTCGGCATCCTCGGCGCGGTATTTGGCCATACATTGc                                                                                 <  4:386182/150‑1 (MQ=255)
                                                                               aTCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCCGTATTGGCGGGATTTCCGCAAACAACGCCGGTTGGGGGGTTTTTTTCGGGATCCCCGGCGCGGGGTTTTGGCATACACTGCtttaatcg                                                                         >  4:777271/1‑144 (MQ=255)
                                                                               aTCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGGATTTCCGGAATACGCGCAGGTTTTGGGGTTTTTTTCGGGGTACTCGGCGGGGGGTTTTGGCATACACTTGTTTatgcgg                                                                        >  4:1213606/1‑150 (MQ=255)
                                                                                 cgccgcGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTTGCGGGATTTCCGCCCACAGCGCCGGTTTCGGGGTTTTTTTCGGCCTCCTCGGCGCGGGGTTTTGCCCCACATTTCtt                                                                               >  4:903213/1‑142 (MQ=255)
                                                                                 cgccgcGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGTATTCCGGCAATCAGCGCAGTTTGCGTGATTTTCGTCGGCATCCTCGGCGCGGTATTTGGCCATACATTGCTTAatgcgatg                                                                      <  2:109822/151‑1 (MQ=255)
                                                                                 cgccgcGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGGATTCCCGCAATCAGCGCAGGTTTGGTGGTTTTTTTCGGGCTCCTCCGCGCGGGGTTTTGGCCCACATTTCTTTatgcg                                                                         >  4:749940/1‑148 (MQ=255)
                                                                                  gccgcGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGGGTTTCCGCCAATCACCCCGCTTTTCGGGGTTTTTTCGGCATCCCCCGCGCCGGGTTTTGGcccaaaatttttta                                                                               >  3:100905/1‑131 (MQ=255)
                                                                                  gccgcGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGGATTTCGGCAATCCGCGCCGGTTGCGGGGTTTTTGTCGGGATCCTCGGCGCGGGATTTTGCCCTAAATTTCtt                                                                               >  1:733237/1‑141 (MQ=255)
                                                                                    cgcgTCAATCCCGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGGATTTCCGCCAACAGCGCAGGTTTCGGGGTTTTTGTCGGCATCCTCCGCGCGGGGTTTTGCCCTACATTTCTTAatgcgatg                                                                      >  1:310228/1‑148 (MQ=255)
                                                                                     gcgTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGTATTCCGGCAATCAGCGCAGTTTGCGTGATTTTCGTCGGCATCCTCGGCGCGGTATTTGGCCATACATTGCTTAATGCGATGCGTa                                                                  <  1:225788/151‑1 (MQ=255)
                                                                                      cgTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGTATTCCGGCAATCAGCGCAGTTTGCGTGATTTTCGTCGGCATCCTCGGCGCGGTATTTGGCCATACATTGCTTAATGCGATGCGTAt                                                                 <  4:1132071/151‑1 (MQ=255)
                                                                                           aTCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGTATTCCGGCAATCAGCGCAGTTTGCGTGATTTTCGTCGGCATCCTCGGCGCGGTATTTGGCCATACATTGCTTAATGCGATGCGTATTCGTa                                                            <  1:1043937/151‑1 (MQ=255)
                                                                                             ccTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGTATTCCGGCAATCAGCGCAGTTTGCGTGATTTTCGTCGGCATCCTCGGCGCGGTATTTGGCCATACATTGCTTAATGCGATGCGTATTCGTAc                                                           <  4:867512/150‑1 (MQ=255)
                                                                                             ccTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGTATTCCGGCAATCAGCGCAGTTTGCGTGATTTTCGTCGGCATCCTCGGCGCGGTATTTGGCCATACATTGCTTAATGCGATGCGTATTCGTAc                                                           <  1:724700/150‑1 (MQ=255)
                                                                                               tGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGGATTTCCGCCAACAGCGCCGGTTTCGGGGTTTTTGTCGGCATCCTCCGCGCGGGATTTTGCCCTACATTTGCTTAATCGATGGGGATTTGTAcc                                                          >  1:718399/1‑149 (MQ=255)
                                                                                                gCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGTATTCCGGCATTCAGCGCAGTTTGCGTGATTTTCGTCGGCATCCTCGGCGCGGTATTTGGCCATACATTGCTTAATGCGATGCGTATTCGTACCaa                                                        <  4:1139631/150‑1 (MQ=255)
                                                                                                gCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGTATTCCGGCAATCAGCGCAGTTTGCGTGATTTTCGTCGGCATCCTCGGCGCGGTATTTGGCCATACATTGCTTAATGCGATGCGTATTCGTACCa                                                         <  4:20813/149‑1 (MQ=255)
                                                                                                  cGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGGATTTCCGCAATCAGCGCAGGTTGCGGGGTTTTTGTCGGCCACCTCCGCGCGGGATTTTGGCATAAATTTCTTAATTGGAAGCGGATTTGTACCaaa                                                       >  1:1125476/1‑149 (MQ=255)
                                                                                                   gAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGTATTCCGGCAATCAGCGCAGTTTGCGTGATTTTCGTCGGCATCCTCGGCGCGGTATTTGGCCATACATTGCTTAATGCGATGCGTATTCGTACCAAAGCt                                                    <  1:1005565/151‑1 (MQ=255)
                                                                                                        cAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGGATTTCCGGCAACAACGCCGGTTTGGGGGTTTTTTTTGGCCTTCCCCGCGCGGGATTTTGGCAAACAATTTTTTATTCGGTGCGTATTccttcccaa                                                       >  4:941716/1‑135 (MQ=255)
                                                                                                        cAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGGATTTCCGCCAACCACGCCGGTTTCGGGGATTTTCTCGGCCTCCCCCGGGCGGGATTTTGGCCTACATTTCTTTATGCGGTGCGGATTTGTaacca                                                        >  1:1250801/1‑138 (MQ=255)
                                                                                                         aGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGTATTCCGGCAATCAGCGCAGTTTGCGTGATTTTCGTCGGCATCCTCGGCGCGGTATTTGGCCATACATTGCTTAATGCGATGCGTATTCGTACCAAAGCTgcgcgc                                              <  4:83221/151‑1 (MQ=255)
                                                                                                          gTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGTATTCCGGCAATCAGCGCAGTTTGCGTGATTTTCGTCGGCATCCTCGGCGCGGTATTTGGCCATACATTGCTTAATGCGATGCGTATTCGTACCAAAGCTGcgcgcg                                             <  1:667592/151‑1 (MQ=255)
                                                                                                           tCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTTGCGGGATTTCCGCCAAAAGCGCCGGTTTCGGGGTTTTTTTCCGCATCCCCCGCGCGGGGTTTTGGCCTACAATTCTTTATTCGCTGCGTATTTGTACCCAAACCGCGCGCgg                                            >  1:1458918/1‑151 (MQ=255)
                                                                                                                 cGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGGATTTCGGCAAACAGCGCCGGTTTCGGGGTTTTTGTCGGCATCCTCCGCGCGGGATTTTGGCCTACATTGCTTTATGCGATGCGTATTTGTACCAAAACCGCGCGCGGGCTgg                                       >  1:1225709/1‑150 (MQ=255)
                                                                                                                 cGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGGATTTCCGCCAACAACGCCGTTTTGGGGGTTTTTTTCGGGCTCCTCCGCGCGGGTTTTTGGCCTACATTTCTTTATTCGGGGGGTATTTGTACCCAAACCGCGCGCggg                                           >  4:1210318/1‑145 (MQ=255)
                                                                                                                  gCCTATTGCAATGGCGGTTGGCGGCCGTATTGGCGGGATTTCCGCAAACAGCGCAGTTTTCGGGGTTTTTGTCGGCAACCTCGGCGCGGGATTTTGGCCTACAATTCTTTATTGGATGCGTATTTGTAACCAAACTGCGCGCGGGCTgg                                       >  2:64138/1‑149 (MQ=255)
                                                                                                                  gCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGTATTTCGGCAATTCGCGCCGTTTTTGGGGTTTTTGGGGGCATCCTCTGGGCGGGATTTTGCCATACATTTCTTTATAGGGTTGGTATTCCGACCAAAACTGCGCGCGGGCTgg                                       >  1:1215657/1‑149 (MQ=255)
                                                                                                                  gCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGGATTTCCGCAATCAGCGCAGGTTTCGGGGTTTTTGTCGGCATCCTCGGGGCGGGATTTTGGCCTACATTTCTTTATGCGATGCGGATTTGTAACAAAACTGCGCGCGGGCTgg                                       >  1:699125/1‑149 (MQ=255)
                                                                                                                   ccTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGTATTCCGGCAATCAGCGCAGTTTGCGTGATTTTCGTCGGCATCCTCGGCGCGGTATTTGGCCATACATTGCTTAATGCGATGCGTATTCGTACCAAAGCTGCGCGCGGTCTGGCaa                                    <  1:76990/151‑1 (MQ=255)
                                                                                                                        tGCAATGGCGGTTGGCGGCAGTATTGGCGGTATTCCGGCAATCAGCGCAGTTTGCGTGATTTTCGTCGGCATCCTCGGCGCGGTATTTGGCCATACATTGCTTAATGCGATGCGTATTCGTACCAAAGCTGCGCGCGGTCTGGCAATggg                                <  1:1140318/150‑1 (MQ=255)
                                                                                                                         gCAATGGCGGTTGGCGGCAGTATTGGCGGTATTCCGGCAATCAGCGCAGTTTGCGTGATTTTCGTCGGCATCCTCGGCGCGGTATTTGGCCATACATTGCTTAATGCGATGCGTATTCGTACCAAAGCTGCGCGCGGTCTGGCAATgggg                               <  1:1252392/150‑1 (MQ=255)
                                                                                                                          cAATGGCGGTTGGCGGCAGTATTGGCGGTATTCCGGCAATCAGCGCAGTTTGCGTGATTTTCGTCGGCATCCTCGGCGCGGTATTTGGCCATACATTGCTTAATGCGATGCGTATTCGTACCAAAGCTGCGCGCGGTCTGGCAATGGGGAc                             <  1:1335846/151‑1 (MQ=255)
                                                                                                                            aTGGCGGTTGGCGGCAGTATTGGCGGTATTCCGGCAATCAGCGCAGTTTGCGTGATTTTCGTCGGCATCCTCGGCGCGGTATTTGGCCATACATTGCTTAATGCGATGCGTATTCGTACCAAAGCTGCGCGCGGTCTGGCAATGGGGACTg                           <  1:1302494/151‑1 (MQ=255)
                                                                                                                                cGGTTGGCGGCAGTATTGGCGGTATTCCGGCAATCAGCGCAGTTTGCGTGATTTTCGTCGGCATCCTCGGCGCGGGATTTTGCCATACATTGCTTAATGCGATGCGTATTCGTACCAAAGCTGCGCGCGGTCTGGCAATGGGGACTGCCTc                       >  4:870940/1‑151 (MQ=255)
                                                                                                                                    tGGCGGCAGTATTGGCGGTATTCCGGCAATCAGCGCAGTTTGCGTGATTTTCGTCGGCATCCTCGGCGCGGTATTTGGCCATACATTGCTTAATGCGATGCGTATTCGTACCAAAGCTGCGCGCGGTCTGGCAATGGGGACTGCCTCGCa                    <  4:302961/150‑1 (MQ=255)
                                                                                                                                     ggcggcAGTATTGGCGGGATTCCGGCAATCAGCGCAGTTTGCGTGGTTTTCGTCGGCATCCTCGGCGCGGGATTTGGCCATACATTGCTTAATGCGATGCGTATTCGTACCAAAGCTGCGCGCGGTCTGGCAATGGGGGCTGCCTCGc                     >  4:1153551/1‑148 (MQ=255)
                                                                                                                                      gcggcAGTATTGGCGGTATTCCGGCAATCAGCGCAGTTTGCGTGATTTTCGTCGGCATCCTCGGCGCGGTATTTGGCCCTACATTGCTTAATGCGATGCGTATTCGTACCAAAGCTGCGCGCGGTCTGGCAATGGGGACTGCCTCGCACGc                 >  2:191133/1‑151 (MQ=255)
                                                                                                                                         gcAGTATTGGCGGTATTCCGGCAATCAGCGCAGTTTGCGTGATTTTCGTCGGCATCCTCGGCGCGGTATTTGGCCATACATTGCTTAATGCGATGCGTATTCGTACCAAAGCTGCGCGCGGTCTGGCAATGGGGACTGCCTCGCACGCCCt              <  1:1390939/151‑1 (MQ=255)
                                                                                                                                           aGTATTGGCGGTATTCCGGCAATCAGCGTAGTTTGCGTGATTTTCGTCGGCATCCTCGGCGCGGGATTTGGCCATACATTGCTTAATGCGATGCGTATTCGTACCAAAGCTGCGCGCGGTCTGGCAATGGGGACTGCCTCGCACGCCCTCg            >  1:200685/1‑151 (MQ=255)
                                                                                                                                           aGTATTGGCGGTATTCCGGCAATCAGCGCAGTTTGCGTGATTTTCGTCGGCATCCTCGGCGCGGGATTTGGCCATACATTGCTTAATGCGATGCGTATTCGTACCAAAGCTGCGCGCGGTCTGGCAATGGGGACTGACTCGCACGCCCTCg            >  4:379414/1‑151 (MQ=255)
                                                                                                                                             tATTGGCGGTATTTCGGCAAACAGCGCAGTTTGCGGGATTTTTGTCGGCATCCTCGGCGCGGGATTTTGCCATACATTGCTTAATGCGATGCGTATTCGTACCAAAACTGCGCGCGGGCTGGCAATGGGGACTGCCTCGCACGCCCTCgg           >  4:1404722/1‑150 (MQ=255)
                                                                                                                                               ttGGCGGTATTCCGGCAATCAGCGCAGTTTGCGTGATTTTCGTCGGCATCCTCGGCGCGGTATTTGGCCATACATTGCTTAATGCGATGCGTATTCGTACCAAAGCTGCGCGCGGTCTGGCAATGGGGACTGCCTCGCACGCCCTCGGTa         <  2:190136/150‑1 (MQ=255)
                                                                                                                                                tGGCGGTATTCCGGCAATCAGCGCAGTTTGCGTGATTTTCGTCGGCATCCTCGGCGCGGTATTTGGCCATACATTGCTTAATGCGATGCGTATTCGTACCAAAGCTGCGCGCGGTCTGGCAATGGGGACTGCCTCGCACGCCCTCGGTACg       >  4:409828/1‑151 (MQ=255)
                                                                                                                                                tGGCGGTATTCCGGCAATCAGCGCAGTTTGCGTGATTTTCGTCGGCATCCTCGGCGCGGTATTTGGCCATACATTGCTTAATGCGATGCGTATTCGTACCAAAGCTGCGCGCGGTCTGGCAATGGGGACTGCCTCGCACGCCCTCGGTACg       >  1:848941/1‑151 (MQ=255)
                                                                                                                                                   cGGTATTCCGGCAATCAGCGCAGTTTGCGTGATTTTCGTCGGCATCCTCGGCGCGGTATTTGGCCATACATTGCTTAATGCGATGCGTATTCGTACCAAAGCTGCGCGCGGTCTGGCTATGGGGACTGCCTCGCACGCCCTCGGTACGgcg    <  1:948419/151‑1 (MQ=255)
                                                                                                                                                   cGGTATTCCGGCAATCAGCGCAGTTTGCGTGATTTTCGTCGGCATCCTCGGCGCGGTATTTGGCCATACATTGCTTAATGCGATGCGTATTCGTACCAAAGCTGCGCGCGGTCTGGCAATGGGGACTGCCTCGCACGCCCTCGGTACGgcg    >  1:201530/1‑151 (MQ=255)
                                                                                                                                                     ggTATTCCGGCAATCAGCGCAGTTTGCGTGATTTTCGTCGGCATCCTCGGCGCGGTATTTGGCCATACATTGCTTAATGCGATGCGTATTCGTACCAAAGCTGCGCGCGGTCTGGCAATGGGGACTGCCTCGCACGCCCTCGGTACGgcgc   >  1:827054/1‑151 (MQ=255)
                                                                                                                                                      gTATTCCGGCAATCAGCGCAGTTTGCGTGATTTTCGTCGGCATCCTCGGCGCGGTATTTGGCCATACATTGCTTAATGCGATGCGTATTCGTACCAAAGCTGCGAGCGGTCTGGCAATGGGGACTGCCTCGCACGCCCTCGGAACGgcgcg  >  1:490236/1‑151 (MQ=255)
                                                                                                                                                       |                                                                                                                                                          
CACCATCTGTTTTATCGGCAGCGTGGTTGCGATGGTGACGGGGACTTCCGTGGCATTGTTGATGGGCGCTTCACCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGTATTCCGGCAATCAGCGCAGTTTGCGTGATTTTCGTCGGCATCCTCGGCGCGGTATTTGGCCATACATTGCTTAATGCGATGCGTATTCGTACCAAAGCTGCGCGCGGTCTGGCAATGGGGACTGCCTCGCACGCCCTCGGTACGGCGCG  >  NC_000913/2231294‑2231593

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