Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,457,1050AG45.7% 84.2 / 28.2 71I441V (ATC→GTC) mplUDP‑N‑acetylmuramate:L‑alanyl‑gamma‑D‑glutamyl‑ meso‑diaminopimelate ligase
Reads supporting (aligned to +/- strand):  ref base A (5/33);  new base G (32/0);  total (37/34)
Fisher's exact test for biased strand distribution p-value = 5.23e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACT‑GCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATG  >  NC_000913/4456958‑4457251
                                                                                                                                                   |                                                                                                                                                   
atTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATcc                                                                                                                                                  <  1:193713/151‑1 (MQ=255)
       ggcaggtggcagAAGTGGCAGAAGCCTGCGTTCAGCCTGCCCACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCa                                                                                                                                                 <  1:323836/145‑1 (MQ=255)
           ggtggcagAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAAct‑g                                                                                                                                      <  1:812517/151‑1 (MQ=255)
           ggtggcagAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAAct                                                                                                                                        <  1:420208/150‑1 (MQ=255)
            gtggcagaagtggcagaaGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTGGTGGGGGCCATCAGAAAat‑g                                                                                                                                      >  1:118938/1‑150 (MQ=255)
             tggcagaagtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTTGTGGGGTCCCTCAGaaa                                                                                                                                          >  1:1099796/1‑146 (MQ=255)
                cagaagtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCCTCCGAAACT‑GCTgg                                                                                                                                  >  2:389641/1‑150 (MQ=255)
                cagaagtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTGGTGGGGTCCATCAGAAAAT‑GCTgg                                                                                                                                  >  1:811616/1‑150 (MQ=255)
                      tggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCCTTACAAACA‑GCTgggaggg                                                                                                                             >  2:1221762/1‑144 (MQ=255)
                      tggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGGCCATCAGAAACT‑GGTGGAAggg                                                                                                                             >  1:42154/1‑149 (MQ=255)
                        gcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGTCCATCAGAAAAT‑GCTGGATGGGCTgg                                                                                                                         >  2:946726/1‑151 (MQ=255)
                        gcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACT‑GCTGGATGGGCTgg                                                                                                                         <  1:808973/151‑1 (MQ=255)
                        gcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACT‑GCTGGATGGGCTgg                                                                                                                         <  1:919285/151‑1 (MQ=255)
                          agaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGTCCATCAGAAAAT‑GCTGGATGGGCTGGCg                                                                                                                       >  2:654388/1‑151 (MQ=255)
                              gCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTTGTGGGGTCCATAAGAAACA‑GCTGGATGGGGTgg                                                                                                                         >  1:511172/1‑145 (MQ=255)
                               ccTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCATCAGAAAAT‑GCTGGGTGGGCTGGCg                                                                                                                       >  2:170214/1‑146 (MQ=255)
                                  gtgTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACT‑GCTGGATGGGCTGGCGaagaag                                                                                                                 <  2:1196693/147‑1 (MQ=255)
                                  gCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACT‑GCTGGATGGGCTGGCGaagaag                                                                                                                 <  1:158962/149‑1 (MQ=255)
                                         gCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACT‑GCTGGATGGGCTGGCGAAGAAGGCGGAAGcc                                                                                                        <  2:1231924/151‑1 (MQ=255)
                                           cTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCATCAGAAACT‑GCTGGATGGGGTGGCGAAGAAAGCGGAAGcc                                                                                                        >  2:137540/1‑149 (MQ=255)
                                           cTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTTGTGGGGACCATCAGAAACT‑TGTGggtggggtgggggaaaaagggg                                                                                                             >  1:718000/1‑132 (MQ=255)
                                                   tGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAAAT‑GCTGGATGGGGTGGCGGAGAAGGCGGGAGCCGCGCAGTaa                                                                                               >  2:1169273/1‑150 (MQ=255)
                                                     gAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACT‑GCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTc                                                                                            <  1:1100340/151‑1 (MQ=255)
                                                        tGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTGGTGGGGTCCATCAGAAACTGGTTGGAGGGCAGGGGGAAGAA‑GCGGGAGCCGCGCAGGAATTctg                                                                                          >  1:452961/1‑148 (MQ=255)
                                                         ggCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACT‑GCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGc                                                                                         <  2:381580/150‑1 (MQ=255)
                                                           cGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCATCCGAAACT‑GCTGGGTGGGCTGGCGGAGAAGGGGGGAGCCGCGCCGGAATTTGGGCtt                                                                                      >  1:636533/1‑151 (MQ=255)
                                                             aTGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACT‑GCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAg                                                                                    <  1:504162/151‑1 (MQ=255)
                                                             aTGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACT‑GCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAg                                                                                    <  1:819687/151‑1 (MQ=255)
                                                                   aTACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGCTCCATCAGAAACT‑GCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGat                                                                              <  2:603323/151‑1 (MQ=255)
                                                                     aCGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACT‑GCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGata                                                                             <  2:847395/150‑1 (MQ=255)
                                                                          ggCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACT‑GCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTAt                                                                       <  2:1042395/151‑1 (MQ=255)
                                                                             aGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCAACAGAAACT‑GCTGGGTGGGCTGGCGAAGAAGGGGGGAGGCGGGCAGGAATTTGGCCTTAGCCAGATATAGGAAGCg                                                                    >  2:175838/1‑151 (MQ=255)
                                                                                      ggtgAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCATCAGAAAAT‑GCTGGATGGGCTGGCGGAGAAGGCGGAAGGCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCgg                                                            >  2:1109461/1‑150 (MQ=255)
                                                                                         gAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCATTAGAAACT‑GCTGGGTGGGGTGGGGGAGAAAGGGGGAGCCGCGCCGGAATTTGGGCCTTGCCCAATATAGTATTGGTAAGcccggt                                                          >  2:443124/1‑145 (MQ=255)
                                                                                           aaaCCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAACT‑GCTGGATGGGCTGGCGGAGAAGGCGGGAGCCGCGCAGGAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGGTAAgg                                                      >  2:1195807/1‑151 (MQ=255)
                                                                                              ccGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACT‑GCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGt                                                   <  1:1026166/151‑1 (MQ=255)
                                                                                                     gCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGATTCATCAGAAACT‑GCTGGGTGGGCTGGCGAAGAAGGCGGGAGGCGCGCAGTAATTCGGCCTTAACCAAATATAGTATGCGTTGGCCGGGTAAAGCGTTTAcgcc                                            >  2:685497/1‑151 (MQ=255)
                                                                                                     gCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACT‑GCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCAcgcc                                            <  2:690060/151‑1 (MQ=255)
                                                                                                      ccTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCAACAGAAAAT‑GCTGGGTTGGGGGGCGGAGAAGGCGGGAGACCCGCAGTAATTTCGCCTTATCCAGATATAAGATGCGGAGGGCGGAGAAGGGGTTTAcac                                             >  1:911188/1‑147 (MQ=255)
                                                                                                      ccTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACT‑GCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCAcgccg                                           <  1:1168351/151‑1 (MQ=255)
                                                                                                      ccTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCATCAGAAAAT‑GCTGGGTGGGCGGGGGAAGAAGGGGGAAGCCGCGCCGGAAATCGGGCCTAGCCCGATATAGTATGCCTAGGGCGGGTAAGGCGGTTAcgcc                                            >  1:514608/1‑150 (MQ=255)
                                                                                                        tGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAAAT‑GCTGGGTGGGCTGGCGGAGAAAGCGGGAGCCGCGCAGTAATTCGGCCCTAGCCAGATATAGTATTGGTAAGCCGGGTAAGGCGTTCAcgcccc                                          >  2:964539/1‑148 (MQ=255)
                                                                                                         ggCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACT‑GCTGGATGGGCGGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCa                                         <  2:941974/150‑1 (MQ=255)
                                                                                                           cGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCATCAGAAACT‑GCTGGGTGGGGTGGCGGAGAAAGGGGAAGCCGCGCAGGAAATTGGCCCTAGCCCGATATTGTAAGCGTAGGCCGGGTAAGGCGTTCCCCCCCCATcc                                      >  1:623608/1‑151 (MQ=255)
                                                                                                           cGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACT‑GCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATcc                                      <  1:58692/151‑1 (MQ=255)
                                                                                                            gACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCATCAGAAAAT‑GCTGGATGGGGTGGCGAAGAAAGCGGGAGCCGCGCAGTAAATCGGCCTTTGCCAGATATAGTATGCGTAGGCCGGGTAAGGCGGTCACGCCGCCTcc                                      >  2:592541/1‑150 (MQ=255)
                                                                                                             aCCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCATCAGAAACT‑GCTGGGTGGGCTGGCGAAGAAGGCGGGAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGGTAAAGCGGTCACGCCGCATCCgg                                    >  2:714281/1‑151 (MQ=255)
                                                                                                              ccATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGTCCATCAAAAACT‑GCTGGATGGGCTGGCGAAGAAAGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGGTAAGGCGTTCACCCCCCATcc                                      >  1:254616/1‑148 (MQ=255)
                                                                                                                atatTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATTAGAAACT‑GCTGGGTGGGGTGGGGGAAAAGGCGGGAGCCGCGCAGTAATTTCGCCTTAGCCAGAAATAGTATGCGTAGGCCGGGTAAGGCGTTCACGCCGCATCCCGc                                   >  1:1221766/1‑149 (MQ=255)
                                                                                                                atatTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACT‑GCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAt                                 <  1:1035238/151‑1 (MQ=255)
                                                                                                                atatTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACT‑GCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAt                                 <  2:376300/151‑1 (MQ=255)
                                                                                                                  atTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGTCCATCAGAAAAT‑GCTGGATTGGCTGGCGAAGAAAGCGGAAGCCGCGCAGTAATTTGGCCTTAGCCAGATATAGTATTCGTAGGCCGGGTAAGGCGGTTACGCCGCATCCGGCAttt                               >  2:233268/1‑151 (MQ=255)
                                                                                                                  atTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACT‑GCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAttt                               <  1:28282/151‑1 (MQ=255)
                                                                                                                   ttCTGGTGATGAGCAACGGCGGTTTTGGGGGGGGCCATCAAAAACT‑GCTGGATGGGGTGGCGAAGAAGGGGGGAGCCGCGCAGGAATTCGGCCTTTGCCCGATATAGGATGCGGAGGCCGGGTAAGGGGTTCACGCCGCCCCCCGCAttt                               >  2:1269449/1‑150 (MQ=255)
                                                                                                                    tCTGGTGATGAGCAACTGCGGTTTTGGTGGGATCCATCAGAAACT‑GCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTg                              <  1:875207/150‑1 (MQ=255)
                                                                                                                      tGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCCTCCGAAAAT‑GGTGGAAGGGCTCGGGAAGAAGGCGGAAGCCGCGCCGGAAATCTGCCCTATCCAGATATAGTATTCGGAAGCCGGATAACGCGGTCACGCCCCATCCCGCAtttt                              >  1:753918/1‑147 (MQ=255)
                                                                                                                      tGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCCTCAGAAACT‑GCTGGGTGGGCTGGCGGAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGGTAAGGCGTTCACGCCGCCTCCGGCATTTGAg                            >  2:943805/1‑150 (MQ=255)
                                                                                                                      tGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCATCAGAAACT‑GCTGGGTGGGGGGGCGAAGAAAGCGGGAGCCGCGCAGGAATTCGGCCCTAGACCGATATAGAATGCGTAGGGCGGGTAAGGCGTTCCCGCCCCATCCCGCATTTTAg                            >  1:61361/1‑150 (MQ=255)
                                                                                                                       ggTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACT‑GCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCa                          <  2:828670/151‑1 (MQ=255)
                                                                                                                           atgaGCAACGGCGGTTTTGGTGGGGTCCATCAGAAACT‑GCTTGGTGGGCTGGCGAAGAAAGCGGAAGCCGCGCAGTAAATCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGGTAAGGGGGTCACCCCGCCTCCCGCATTTTTGCATa                        >  1:353786/1‑149 (MQ=255)
                                                                                                                           atgaGCAACGGCGGTTTTGGTGGGATCCATCAGAAACT‑GCGGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTGACGCCGCATCCGGCATTTGAGCAt                         <  1:276086/148‑1 (MQ=255)
                                                                                                                             gaGCAACGGCGGTTTTGGTGGGATCCATCAGAAACT‑GCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTg                     <  1:732570/150‑1 (MQ=255)
                                                                                                                                    ggcggTTTTGGTGGGATCCATCAGAAACT‑GCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGc             >  2:244838/1‑151 (MQ=255)
                                                                                                                                     gcggTTTTGGTGGGATCCATCAGAAACT‑GCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCg            >  2:571921/1‑151 (MQ=255)
                                                                                                                                     gcggTTTTGGTGGGATCCATCAGAAACT‑GCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCg            <  2:61960/151‑1 (MQ=255)
                                                                                                                                      cggTTTTGGTGGGATCCATCAGAAACT‑GCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGa           <  2:15157/151‑1 (MQ=255)
                                                                                                                                           ttGGTGGGATCCATCAGAAACT‑GCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCt       <  1:662596/150‑1 (MQ=255)
                                                                                                                                             ggtggGATCCATCAGAAACT‑GCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTg     <  2:261599/150‑1 (MQ=255)
                                                                                                                                             ggtggGATCCATCAGAAACT‑GCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTg     <  2:1079536/150‑1 (MQ=255)
                                                                                                                                             ggtggGATCCATCAGAAACT‑GCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGa    <  2:1099835/151‑1 (MQ=255)
                                                                                                                                                gggATCCATCAGAAACT‑GCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATg  >  1:788938/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACT‑GCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATG  >  NC_000913/4456958‑4457251

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