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

ATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATG  >  NC_000913/4456958‑4457251
                                                                                                                                                   |                                                                                                                                                   
atTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATcc                                                                                                                                                  <  1:193713/151‑1 (MQ=255)
       ggcaggtggcagAAGTGGCAGAAGCCTGCGTTCAGCCTGCCCACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCa                                                                                                                                                 <  1:323836/145‑1 (MQ=255)
           ggtggcagAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAActg                                                                                                                                      <  1:812517/151‑1 (MQ=255)
           ggtggcagAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAAct                                                                                                                                        <  1:420208/150‑1 (MQ=255)
            gtggcagaagtggcagaaGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTGGTGGGGGCCATCAGAAAatg                                                                                                                                      >  1:118938/1‑150 (MQ=255)
             tggcagaagtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTTGTGGGGTCCCTCAGaaa                                                                                                                                          >  1:1099796/1‑146 (MQ=255)
                cagaagtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCCTCCGAAACTGCTgg                                                                                                                                  >  2:389641/1‑150 (MQ=255)
                cagaagtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTGGTGGGGTCCATCAGAAAATGCTgg                                                                                                                                  >  1:811616/1‑150 (MQ=255)
                      tggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCCTTACAAACAGCTgggaggg                                                                                                                             >  2:1221762/1‑144 (MQ=255)
                      tggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGGCCATCAGAAACTGGTGGAAggg                                                                                                                             >  1:42154/1‑149 (MQ=255)
                        gcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGTCCATCAGAAAATGCTGGATGGGCTgg                                                                                                                         >  2:946726/1‑151 (MQ=255)
                        gcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTgg                                                                                                                         <  1:808973/151‑1 (MQ=255)
                        gcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTgg                                                                                                                         <  1:919285/151‑1 (MQ=255)
                          agaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGTCCATCAGAAAATGCTGGATGGGCTGGCg                                                                                                                       >  2:654388/1‑151 (MQ=255)
                              gCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTTGTGGGGTCCATAAGAAACAGCTGGATGGGGTgg                                                                                                                         >  1:511172/1‑145 (MQ=255)
                               ccTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCATCAGAAAATGCTGGGTGGGCTGGCg                                                                                                                       >  2:170214/1‑146 (MQ=255)
                                  gtgTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGaagaag                                                                                                                 <  2:1196693/147‑1 (MQ=255)
                                  gCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGaagaag                                                                                                                 <  1:158962/149‑1 (MQ=255)
                                         gCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGcc                                                                                                        <  2:1231924/151‑1 (MQ=255)
                                           cTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCATCAGAAACTGCTGGATGGGGTGGCGAAGAAAGCGGAAGcc                                                                                                        >  2:137540/1‑149 (MQ=255)
                                           cTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTTGTGGGGACCATCAGAAACTTGTGggtggggtgggggaaaaagggg                                                                                                             >  1:718000/1‑132 (MQ=255)
                                                   tGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAAATGCTGGATGGGGTGGCGGAGAAGGCGGGAGCCGCGCAGTaa                                                                                               >  2:1169273/1‑150 (MQ=255)
                                                     gAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTc                                                                                            <  1:1100340/151‑1 (MQ=255)
                                                        tGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTGGTGGGGTCCATCAGAAACTGGTTGGAGGGCAGGGGGAAGAAGCGGGAGCCGCGCAGGAATTctg                                                                                          >  1:452961/1‑148 (MQ=255)
                                                         ggCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGc                                                                                         <  2:381580/150‑1 (MQ=255)
                                                           cGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCATCCGAAACTGCTGGGTGGGCTGGCGGAGAAGGGGGGAGCCGCGCCGGAATTTGGGCtt                                                                                      >  1:636533/1‑151 (MQ=255)
                                                             aTGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAg                                                                                    <  1:504162/151‑1 (MQ=255)
                                                             aTGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAg                                                                                    <  1:819687/151‑1 (MQ=255)
                                                                   aTACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGCTCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGat                                                                              <  2:603323/151‑1 (MQ=255)
                                                                     aCGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGata                                                                             <  2:847395/150‑1 (MQ=255)
                                                                          ggCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTAt                                                                       <  2:1042395/151‑1 (MQ=255)
                                                                             aGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCAACAGAAACTGCTGGGTGGGCTGGCGAAGAAGGGGGGAGGCGGGCAGGAATTTGGCCTTAGCCAGATATAGGAAGCg                                                                    >  2:175838/1‑151 (MQ=255)
                                                                                      ggtgAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCATCAGAAAATGCTGGATGGGCTGGCGGAGAAGGCGGAAGGCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCgg                                                            >  2:1109461/1‑150 (MQ=255)
                                                                                         gAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCATTAGAAACTGCTGGGTGGGGTGGGGGAGAAAGGGGGAGCCGCGCCGGAATTTGGGCCTTGCCCAATATAGTATTGGTAAGcccggt                                                          >  2:443124/1‑145 (MQ=255)
                                                                                           aaaCCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAACTGCTGGATGGGCTGGCGGAGAAGGCGGGAGCCGCGCAGGAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGGTAAgg                                                      >  2:1195807/1‑151 (MQ=255)
                                                                                              ccGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGt                                                   <  1:1026166/151‑1 (MQ=255)
                                                                                                     gCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGATTCATCAGAAACTGCTGGGTGGGCTGGCGAAGAAGGCGGGAGGCGCGCAGTAATTCGGCCTTAACCAAATATAGTATGCGTTGGCCGGGTAAAGCGTTTAcgcc                                            >  2:685497/1‑151 (MQ=255)
                                                                                                     gCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCAcgcc                                            <  2:690060/151‑1 (MQ=255)
                                                                                                      ccTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCAACAGAAAATGCTGGGTTGGGGGGCGGAGAAGGCGGGAGACCCGCAGTAATTTCGCCTTATCCAGATATAAGATGCGGAGGGCGGAGAAGGGGTTTAcac                                             >  1:911188/1‑147 (MQ=255)
                                                                                                      ccTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCAcgccg                                           <  1:1168351/151‑1 (MQ=255)
                                                                                                      ccTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCATCAGAAAATGCTGGGTGGGCGGGGGAAGAAGGGGGAAGCCGCGCCGGAAATCGGGCCTAGCCCGATATAGTATGCCTAGGGCGGGTAAGGCGGTTAcgcc                                            >  1:514608/1‑150 (MQ=255)
                                                                                                        tGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAAATGCTGGGTGGGCTGGCGGAGAAAGCGGGAGCCGCGCAGTAATTCGGCCCTAGCCAGATATAGTATTGGTAAGCCGGGTAAGGCGTTCAcgcccc                                          >  2:964539/1‑148 (MQ=255)
                                                                                                         ggCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCGGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCa                                         <  2:941974/150‑1 (MQ=255)
                                                                                                           cGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCATCAGAAACTGCTGGGTGGGGTGGCGGAGAAAGGGGAAGCCGCGCAGGAAATTGGCCCTAGCCCGATATTGTAAGCGTAGGCCGGGTAAGGCGTTCCCCCCCCATcc                                      >  1:623608/1‑151 (MQ=255)
                                                                                                           cGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATcc                                      <  1:58692/151‑1 (MQ=255)
                                                                                                            gACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCATCAGAAAATGCTGGATGGGGTGGCGAAGAAAGCGGGAGCCGCGCAGTAAATCGGCCTTTGCCAGATATAGTATGCGTAGGCCGGGTAAGGCGGTCACGCCGCCTcc                                      >  2:592541/1‑150 (MQ=255)
                                                                                                             aCCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCATCAGAAACTGCTGGGTGGGCTGGCGAAGAAGGCGGGAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGGTAAAGCGGTCACGCCGCATCCgg                                    >  2:714281/1‑151 (MQ=255)
                                                                                                              ccATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGTCCATCAAAAACTGCTGGATGGGCTGGCGAAGAAAGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGGTAAGGCGTTCACCCCCCATcc                                      >  1:254616/1‑148 (MQ=255)
                                                                                                                atatTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATTAGAAACTGCTGGGTGGGGTGGGGGAAAAGGCGGGAGCCGCGCAGTAATTTCGCCTTAGCCAGAAATAGTATGCGTAGGCCGGGTAAGGCGTTCACGCCGCATCCCGc                                   >  1:1221766/1‑149 (MQ=255)
                                                                                                                atatTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAt                                 <  1:1035238/151‑1 (MQ=255)
                                                                                                                atatTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAt                                 <  2:376300/151‑1 (MQ=255)
                                                                                                                  atTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGTCCATCAGAAAATGCTGGATTGGCTGGCGAAGAAAGCGGAAGCCGCGCAGTAATTTGGCCTTAGCCAGATATAGTATTCGTAGGCCGGGTAAGGCGGTTACGCCGCATCCGGCAttt                               >  2:233268/1‑151 (MQ=255)
                                                                                                                  atTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAttt                               <  1:28282/151‑1 (MQ=255)
                                                                                                                   ttCTGGTGATGAGCAACGGCGGTTTTGGGGGGGGCCATCAAAAACTGCTGGATGGGGTGGCGAAGAAGGGGGGAGCCGCGCAGGAATTCGGCCTTTGCCCGATATAGGATGCGGAGGCCGGGTAAGGGGTTCACGCCGCCCCCCGCAttt                               >  2:1269449/1‑150 (MQ=255)
                                                                                                                    tCTGGTGATGAGCAACTGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTg                              <  1:875207/150‑1 (MQ=255)
                                                                                                                      tGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCCTCCGAAAATGGTGGAAGGGCTCGGGAAGAAGGCGGAAGCCGCGCCGGAAATCTGCCCTATCCAGATATAGTATTCGGAAGCCGGATAACGCGGTCACGCCCCATCCCGCAtttt                              >  1:753918/1‑147 (MQ=255)
                                                                                                                      tGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCCTCAGAAACTGCTGGGTGGGCTGGCGGAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGGTAAGGCGTTCACGCCGCCTCCGGCATTTGAg                            >  2:943805/1‑150 (MQ=255)
                                                                                                                      tGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCATCAGAAACTGCTGGGTGGGGGGGCGAAGAAAGCGGGAGCCGCGCAGGAATTCGGCCCTAGACCGATATAGAATGCGTAGGGCGGGTAAGGCGTTCCCGCCCCATCCCGCATTTTAg                            >  1:61361/1‑150 (MQ=255)
                                                                                                                       ggTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCa                          <  2:828670/151‑1 (MQ=255)
                                                                                                                           atgaGCAACGGCGGTTTTGGTGGGGTCCATCAGAAACTGCTTGGTGGGCTGGCGAAGAAAGCGGAAGCCGCGCAGTAAATCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGGTAAGGGGGTCACCCCGCCTCCCGCATTTTTGCATa                        >  1:353786/1‑149 (MQ=255)
                                                                                                                           atgaGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCGGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTGACGCCGCATCCGGCATTTGAGCAt                         <  1:276086/148‑1 (MQ=255)
                                                                                                                             gaGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTg                     <  1:732570/150‑1 (MQ=255)
                                                                                                                                    ggcggTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGc             >  2:244838/1‑151 (MQ=255)
                                                                                                                                     gcggTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCg            >  2:571921/1‑151 (MQ=255)
                                                                                                                                     gcggTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCg            <  2:61960/151‑1 (MQ=255)
                                                                                                                                      cggTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGa           <  2:15157/151‑1 (MQ=255)
                                                                                                                                           ttGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCt       <  1:662596/150‑1 (MQ=255)
                                                                                                                                             ggtggGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTg     <  2:261599/150‑1 (MQ=255)
                                                                                                                                             ggtggGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTg     <  2:1079536/150‑1 (MQ=255)
                                                                                                                                             ggtggGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGa    <  2:1099835/151‑1 (MQ=255)
                                                                                                                                                gggATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATg  >  1:788938/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATG  >  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: