Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,386,8620GA41.9% 63.1 / 51.4 63intergenic (+31/‑72)yfbL/yfbMputative M28A family peptidase/DUF1877 family protein
Reads supporting (aligned to +/- strand):  ref base G (4/32);  new base A (26/0);  total (30/33)
Fisher's exact test for biased strand distribution p-value = 1.31e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTTTTTATCGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGAT‑GC‑G‑AAAC‑ATGAAATTCAACGGA‑GGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTG  >  NC_000913/2386715‑2387011
                                                                                                                                                   |                                                                                                                                                          
gCTTTTTATCGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGa                                                                                                                                                         <  2:238599/151‑1 (MQ=255)
     ttATCGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGAAATGc                                                                                                                                                    >  2:217936/1‑151 (MQ=255)
      tATCGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGAGGGTCCAGAAAAAGCGACATGc                                                                                                                                                    >  2:1332042/1‑150 (MQ=255)
         cGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACat                                                                                                                                                <  2:565569/151‑1 (MQ=255)
            aataaACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACattt                                                                                                                                              <  1:1116242/150‑1 (MQ=255)
            aataaACAATACCACTTGCACGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCGGGAAAAAAGACaatggcatttt                                                                                                                                             >  2:1012974/1‑141 (MQ=255)
             ataaACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGGCAAGAACatttt                                                                                                                                             >  1:79726/1‑150 (MQ=255)
                 aCAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAAAGACATGCACaatttttttt                                                                                                                                        >  1:930097/1‑151 (MQ=255)
                 aCAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAACGACAAGCACattttttttt                                                                                                                                        >  2:219735/1‑151 (MQ=255)
                      aCCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCt                                                                                                                                    <  1:1123408/150‑1 (MQ=255)
                              cccGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGAAAAGGACAATTTTTTTTTCCCCATCTtttt                                                                                                                           >  2:473879/1‑151 (MQ=255)
                              cccGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAAAAAAAAAAACCCGCACACTTTTTTTTTCCTCATTTCttt                                                                                                                           >  2:542537/1‑151 (MQ=255)
                                    gATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGGCACAAAAAAAGACATGCACATTTTATTTTTTCTCTTCTGtttt                                                                                                                          >  2:1086216/1‑146 (MQ=255)
                                         cGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGAT‑GC‑G‑aa                                                                                                             <  2:1198011/151‑1 (MQ=255)
                                         cGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGGCCGAAAAAAAAACAAGCACATTTTTTTTTTTCTCCTCttttttttgatacg                                                                                                                  >  1:27309/1‑141 (MQ=255)
                                             gacagaTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGAT‑GC‑G‑AAAC‑At                                                                                                        <  1:817817/151‑1 (MQ=255)
                                                 gaTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGAT‑GC‑G‑AAAC‑ATGaaa                                                                                                    <  1:308322/151‑1 (MQ=255)
                                                     gAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGAT‑GC‑G‑AAAC‑ATGAAATTc                                                                                                 <  2:1052698/150‑1 (MQ=255)
                                                     gAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGAT‑GC‑G‑AAAC‑ATGAAATTc                                                                                                 <  1:367563/150‑1 (MQ=255)
                                                          aTCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGAT‑GC‑G‑AAAC‑ATGAAATTCAACgg                                                                                            <  2:1105766/150‑1 (MQ=255)
                                                            cAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTGCCGAAAAAAAGAAAAGCACATTTTTTTTTTTCTCTTCTGTTTTTAGAT‑GC‑G‑AAAA‑AAGaaaaataaaggag                                                                                          >  2:847987/1‑140 (MQ=255)
                                                                aaaTGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGAT‑GC‑G‑AAAC‑ATGAAATTCAACGGA‑GGaaaa                                                                                    <  1:137184/151‑1 (MQ=255)
                                                                    ggCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAACGACAAGCACATTTTTTTTTTCCTCATCTGTTTTTAAAT‑GC‑GAAAAA‑AAAAAATTAAAGGGG‑GGaaaaaa                                                                                  >  2:748291/1‑149 (MQ=255)
                                                                     gCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGGCCGAAAAAAGGACAAGGACATTTTTTTTTTCCTCTTCTGTTTTTAGAA‑GCGG‑AAAA‑ATAAATTTCACCGAG‑GGaaaaaa                                                                                  >  2:1163428/1‑148 (MQ=255)
                                                                        cAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGAT‑GC‑G‑AAAC‑ATGAAATTCAACGGA‑GGAAAAATAAATg                                                                             <  1:274155/150‑1 (MQ=255)
                                                                        cAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGCCCCGCAAAAACGAAAAGCCCAATTTATTTTTCCCCTTCTGTTTTTAGAT‑GC‑G‑GAAA‑ATGAAAATCAACAGG‑GGAaaaaaaaa                                                                               >  2:718850/1‑148 (MQ=255)
                                                                         aGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGAT‑GC‑G‑AAAC‑ATGAAATTCAACGGA‑GGAAAAATAAATggg                                                                           <  2:784248/151‑1 (MQ=255)
                                                                            tAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAAAAAAAAAAACATGGACAATTTATTTTTTCTTATCTTTTTTTAGAT‑GC‑G‑CAAA‑ATGAAAATTAAAGGA‑GGAAAAATAAATGGGGAt                                                                        >  1:1337271/1‑151 (MQ=255)
                                                                             aGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAAAGAAATGCCCAATTTTTTTTTCCTCCTCTTTTTTTAAAT‑GG‑G‑GAAA‑AAGAAATTTAACGGA‑GGGAAAAAAAAAggg                                                                           >  2:865343/1‑147 (MQ=255)
                                                                             aGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGAT‑GC‑G‑AAAC‑ATGAAATTCAACGGA‑GGAAAAATAAATGGGAAt                                                                        <  2:195122/150‑1 (MQ=255)
                                                                                ggATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGAT‑GC‑G‑AAAC‑ATGAAATTCAACGGA‑GGAAAAATAAATGGGAATGAt                                                                     <  1:571849/150‑1 (MQ=255)
                                                                                 gATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGAT‑GC‑G‑AAAC‑ATGAAATTCAACGGA‑GGAAAAATAAATGGGAATGATTg                                                                   <  1:700579/151‑1 (MQ=255)
                                                                                 gATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGAT‑GC‑G‑AAAC‑ATGAAATTCAACGGA‑GGAAAAATAAATGGGAATGATTg                                                                   <  1:1337247/151‑1 (MQ=255)
                                                                                     gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGCTGTAGAAAAAAAGAAAAGCACAATTTATTTTTTCTCTTCTGTTTTTAGAA‑TC‑G‑AAAC‑AAGAAAATTAAAGGG‑GGAAAAAAAAAAGGGAATGAGTgggt                                                                >  1:329188/1‑148 (MQ=255)
                                                                                       gtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGAT‑GC‑G‑AAAC‑ATGAAATTCAACGGA‑GGAAAAATAAATGGGAATGATTGGCTAtt                                                             <  1:863471/151‑1 (MQ=255)
                                                                                       gtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGATCCAGAAAAAAAGACAAAAACAATTTTTTTTTTCCCATCTTTTTTTTAATAAC‑G‑AAA‑‑ATAAAAATCTACGGA‑GGAAAAAAAAATAGGAAAGATTGGCTAtt                                                             >  1:598757/1‑151 (MQ=255)
                                                                                        ttATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAACGACATGCACCTTTTTTTTTTCCTCCTCTGTTTTTAGAA‑GC‑G‑AAAA‑ATGAAAATCAACGGA‑GGAAAAATAAATGGGAATGATTGGCTAttt                                                            >  2:1027976/1‑151 (MQ=255)
                                                                                          ataACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGAT‑GC‑G‑AAAC‑ATGAAATTCAACGGA‑GGAAAAATAAATGGGAATGATTGGCTAtttt                                                           <  2:108066/150‑1 (MQ=255)
                                                                                          ataACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGAT‑GC‑G‑AAAC‑ATGAAATTCAACGGA‑GGAAAAATAAATGGGAATGATTGGCTAtttt                                                           <  1:461231/150‑1 (MQ=255)
                                                                                                ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAACGAAAAGCACAATTTTTTTTTCCTCCTCCGTTTTTATAT‑GG‑G‑GAAA‑AAGAAATTTAAAGGA‑GGGAAAATAAAAGGGGATGATTGGCTAttttt                                                          >  1:1266221/1‑144 (MQ=255)
                                                                                                ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGAT‑GC‑G‑AAAC‑ATGAAATTCAACGGA‑GGAAAAATAAATGGGAATGATTGGCTATTTTGCAGa                                                      <  2:323804/149‑1 (MQ=255)
                                                                                                 tGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGACAAGCCCCTTTTTTTTTTCCTCCTCTTTTTTTAAAA‑GC‑GAAAAA‑ATAAAATTAAAGGGG‑GGAAAAA‑AAATAGGGATGATTTGGTATTTTTTAGaaa                                                    >  2:792772/1‑150 (MQ=255)
                                                                                                 tGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGACAAGAACAATTTTTTTTTTCTCCTCTTTTTTTAAAA‑GC‑G‑AAAA‑AAGAAATTTAACGGA‑GGGAAAAAAAATGGGGATGATTTGCTATTTTTTAGaaa                                                    >  1:299006/1‑150 (MQ=255)
                                                                                                  gTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGAT‑GC‑G‑AAAC‑ATGAAATTCAACGGA‑GGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAAtt                                                  <  1:334959/151‑1 (MQ=255)
                                                                                                    tataCAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGAT‑GC‑G‑AAAC‑ATGAAATTCAACGGA‑GGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTg                                                 <  2:766285/150‑1 (MQ=255)
                                                                                                     ataCAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGAT‑GC‑G‑AAAC‑ATGAAATTCAACGGA‑GGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGa                                                <  2:1295125/150‑1 (MQ=255)
                                                                                                       acaacaGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGAT‑GC‑G‑AAAC‑ATGAAATTCAACGGA‑GGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTc                                             <  1:1025189/151‑1 (MQ=255)
                                                                                                       acaacaGTAAATAATAAAAATTTCTCGAGGCGGGGGCCGAAAAAACGGAAAACACCATTTTTTTTTCCTCCTCTGTTTTTAAAA‑GA‑G‑AAACAATAAAATTAAAGGAG‑GGAAAAAAAAAGGGGAAGGTTTGG‑TATTTTTCGGaaaattgatc                                             >  2:840050/1‑144 (MQ=255)
                                                                                                         aacaGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGAT‑GC‑G‑AAAC‑ATGAAATTCAACGGA‑GGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAg                                           <  1:1315496/151‑1 (MQ=255)
                                                                                                         aacaGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGACATGCACATTTTTTTTTTCCTCCTCTTTTTTAAGAT‑GC‑G‑AAAC‑ATGAAATTTAACGGA‑GGGAAAATAAATGGGGAAGATTGGGTATTTTTCAGaaaat                                                  >  2:804010/1‑142 (MQ=255)
                                                                                                          acaGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGAT‑GC‑G‑AAAC‑ATGAAATTCAACGGA‑GGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGa                                          <  1:990007/151‑1 (MQ=255)
                                                                                                            aGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAACCACATGCACAATTTATTTTTCCTCCTCTGTTTTTAAAT‑GC‑G‑AAAA‑ATTAAAATAAACGGAGGGAAAAAAAAAAGGGAA‑GATTTGGTATTTTTCAGAAAATGACTCAGaa                                         >  2:1210105/1‑150 (MQ=255)
                                                                                                                  taataaAAATTTCTCGAGGCGGGGCCCGAAAAAAAGGCCAGCACATTTTATTTTTCCTCCTTCGTTTTTAGAT‑GG‑G‑AAAC‑AAGAAAATCAAAGGA‑GGAAAAAAAAATGGGGATGATTGGCTATTTTTTAGAAAATTACTTAGaaaa                                       >  1:1097593/1‑146 (MQ=255)
                                                                                                                    ataaAAATTTCTCGAGGCGGGTCCAGAAAAAAAGACATGCACAATTTATTTTTTCTCCTCTGTTTTTAGAT‑GC‑G‑AAAA‑ATGAAAATTAAACGA‑GGAAAAAAAAATGGGGATGATTGGGTTTTTTTTAGAAAATGACTCCGaaaa                                       >  1:610168/1‑144 (MQ=255)
                                                                                                                      aaaaaTTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGAT‑GC‑G‑AAAC‑ATGAAATTCAACGGA‑GGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAAcc                               <  2:549386/150‑1 (MQ=255)
                                                                                                                                cGAGGCGGGTCCATAAAAATCGACATGCACATTTTATTGTCCCCCATCTGTTTTAAGAT‑GC‑G‑AAAC‑ATGAAATTCAACGGA‑GGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTg                      <  1:175539/149‑1 (MQ=255)
                                                                                                                                cGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGAT‑GC‑G‑AAAC‑ATGAAATTCAACGGA‑GGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGa                     <  1:144050/150‑1 (MQ=255)
                                                                                                                                       ggTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGAT‑GC‑G‑AAAC‑ATGAAATTCAACGGA‑GGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTg             <  2:969485/151‑1 (MQ=255)
                                                                                                                                          ccAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGAT‑GC‑G‑AAAC‑ATGAAATTCAACGGA‑GGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGaa           >  1:1194090/1‑150 (MQ=255)
                                                                                                                                            aGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGAT‑GC‑G‑AAAC‑ATGAAATTCAACGGA‑GGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGaaaa         >  2:998748/1‑150 (MQ=255)
                                                                                                                                                aaaGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGAT‑GC‑G‑AAAC‑ATGAAATTCAACGGA‑GGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCtt    <  2:276591/151‑1 (MQ=255)
                                                                                                                                                aaaGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGAT‑GC‑G‑AAAC‑ATGAAATTCAACGGA‑GGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCtt    <  2:947871/151‑1 (MQ=255)
                                                                                                                                                  aGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGAT‑GC‑G‑AAAC‑ATGAAATTCAACGGA‑GGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTg  >  2:512504/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                          
GCTTTTTATCGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGAT‑GC‑G‑AAAC‑ATGAAATTCAACGGA‑GGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTG  >  NC_000913/2386715‑2387011

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