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

GCTTTTTATCGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTG  >  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)
                                         cGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGaa                                                                                                             <  2:1198011/151‑1 (MQ=255)
                                         cGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGGCCGAAAAAAAAACAAGCACATTTTTTTTTTTCTCCTCttttttttgatacg                                                                                                                  >  1:27309/1‑141 (MQ=255)
                                             gacagaTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACAt                                                                                                        <  1:817817/151‑1 (MQ=255)
                                                 gaTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGaaa                                                                                                    <  1:308322/151‑1 (MQ=255)
                                                     gAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTc                                                                                                 <  2:1052698/150‑1 (MQ=255)
                                                     gAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTc                                                                                                 <  1:367563/150‑1 (MQ=255)
                                                          aTCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACgg                                                                                            <  2:1105766/150‑1 (MQ=255)
                                                            cAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTGCCGAAAAAAAGAAAAGCACATTTTTTTTTTTCTCTTCTGTTTTTAGATGCGAAAAAAGaaaaataaaggag                                                                                          >  2:847987/1‑140 (MQ=255)
                                                                aaaTGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGaaaa                                                                                    <  1:137184/151‑1 (MQ=255)
                                                                    ggCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAACGACAAGCACATTTTTTTTTTCCTCATCTGTTTTTAAATGCGAAAAAAAAAAATTAAAGGGGGGaaaaaa                                                                                  >  2:748291/1‑149 (MQ=255)
                                                                     gCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGGCCGAAAAAAGGACAAGGACATTTTTTTTTTCCTCTTCTGTTTTTAGAAGCGGAAAAATAAATTTCACCGAGGGaaaaaa                                                                                  >  2:1163428/1‑148 (MQ=255)
                                                                        cAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATg                                                                             <  1:274155/150‑1 (MQ=255)
                                                                        cAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGCCCCGCAAAAACGAAAAGCCCAATTTATTTTTCCCCTTCTGTTTTTAGATGCGGAAAATGAAAATCAACAGGGGAaaaaaaaa                                                                               >  2:718850/1‑148 (MQ=255)
                                                                         aGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATggg                                                                           <  2:784248/151‑1 (MQ=255)
                                                                            tAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAAAAAAAAAAACATGGACAATTTATTTTTTCTTATCTTTTTTTAGATGCGCAAAATGAAAATTAAAGGAGGAAAAATAAATGGGGAt                                                                        >  1:1337271/1‑151 (MQ=255)
                                                                             aGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAAAGAAATGCCCAATTTTTTTTTCCTCCTCTTTTTTTAAATGGGGAAAAAGAAATTTAACGGAGGGAAAAAAAAAggg                                                                           >  2:865343/1‑147 (MQ=255)
                                                                             aGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAAt                                                                        <  2:195122/150‑1 (MQ=255)
                                                                                ggATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGAt                                                                     <  1:571849/150‑1 (MQ=255)
                                                                                 gATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTg                                                                   <  1:700579/151‑1 (MQ=255)
                                                                                 gATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTg                                                                   <  1:1337247/151‑1 (MQ=255)
                                                                                     gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGCTGTAGAAAAAAAGAAAAGCACAATTTATTTTTTCTCTTCTGTTTTTAGAATCGAAACAAGAAAATTAAAGGGGGAAAAAAAAAAGGGAATGAGTgggt                                                                >  1:329188/1‑148 (MQ=255)
                                                                                       gtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTAtt                                                             <  1:863471/151‑1 (MQ=255)
                                                                                       gtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGATCCAGAAAAAAAGACAAAAACAATTTTTTTTTTCCCATCTTTTTTTTAATAACGAAAATAAAAATCTACGGAGGAAAAAAAAATAGGAAAGATTGGCTAtt                                                             >  1:598757/1‑151 (MQ=255)
                                                                                        ttATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAACGACATGCACCTTTTTTTTTTCCTCCTCTGTTTTTAGAAGCGAAAAATGAAAATCAACGGAGGAAAAATAAATGGGAATGATTGGCTAttt                                                            >  2:1027976/1‑151 (MQ=255)
                                                                                          ataACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTAtttt                                                           <  2:108066/150‑1 (MQ=255)
                                                                                          ataACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTAtttt                                                           <  1:461231/150‑1 (MQ=255)
                                                                                                ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAACGAAAAGCACAATTTTTTTTTCCTCCTCCGTTTTTATATGGGGAAAAAGAAATTTAAAGGAGGGAAAATAAAAGGGGATGATTGGCTAttttt                                                          >  1:1266221/1‑144 (MQ=255)
                                                                                                ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGa                                                      <  2:323804/149‑1 (MQ=255)
                                                                                                 tGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGACAAGCCCCTTTTTTTTTTCCTCCTCTTTTTTTAAAAGCGAAAAAATAAAATTAAAGGGGGGAAAAAAAATAGGGATGATTTGGTATTTTTTAGaaa                                                    >  2:792772/1‑150 (MQ=255)
                                                                                                 tGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGACAAGAACAATTTTTTTTTTCTCCTCTTTTTTTAAAAGCGAAAAAAGAAATTTAACGGAGGGAAAAAAAATGGGGATGATTTGCTATTTTTTAGaaa                                                    >  1:299006/1‑150 (MQ=255)
                                                                                                  gTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAAtt                                                  <  1:334959/151‑1 (MQ=255)
                                                                                                    tataCAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTg                                                 <  2:766285/150‑1 (MQ=255)
                                                                                                     ataCAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGa                                                <  2:1295125/150‑1 (MQ=255)
                                                                                                       acaacaGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTc                                             <  1:1025189/151‑1 (MQ=255)
                                                                                                       acaacaGTAAATAATAAAAATTTCTCGAGGCGGGGGCCGAAAAAACGGAAAACACCATTTTTTTTTCCTCCTCTGTTTTTAAAAGAGAAACAATAAAATTAAAGGAGGGAAAAAAAAAGGGGAAGGTTTGGTATTTTTCGGaaaattgatc                                             >  2:840050/1‑144 (MQ=255)
                                                                                                         aacaGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAg                                           <  1:1315496/151‑1 (MQ=255)
                                                                                                         aacaGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGACATGCACATTTTTTTTTTCCTCCTCTTTTTTAAGATGCGAAACATGAAATTTAACGGAGGGAAAATAAATGGGGAAGATTGGGTATTTTTCAGaaaat                                                  >  2:804010/1‑142 (MQ=255)
                                                                                                          acaGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGa                                          <  1:990007/151‑1 (MQ=255)
                                                                                                            aGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAACCACATGCACAATTTATTTTTCCTCCTCTGTTTTTAAATGCGAAAAATTAAAATAAACGGAGGGAAAAAAAAAAGGGAAGATTTGGTATTTTTCAGAAAATGACTCAGaa                                         >  2:1210105/1‑150 (MQ=255)
                                                                                                                  taataaAAATTTCTCGAGGCGGGGCCCGAAAAAAAGGCCAGCACATTTTATTTTTCCTCCTTCGTTTTTAGATGGGAAACAAGAAAATCAAAGGAGGAAAAAAAAATGGGGATGATTGGCTATTTTTTAGAAAATTACTTAGaaaa                                       >  1:1097593/1‑146 (MQ=255)
                                                                                                                    ataaAAATTTCTCGAGGCGGGTCCAGAAAAAAAGACATGCACAATTTATTTTTTCTCCTCTGTTTTTAGATGCGAAAAATGAAAATTAAACGAGGAAAAAAAAATGGGGATGATTGGGTTTTTTTTAGAAAATGACTCCGaaaa                                       >  1:610168/1‑144 (MQ=255)
                                                                                                                      aaaaaTTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAAcc                               <  2:549386/150‑1 (MQ=255)
                                                                                                                                cGAGGCGGGTCCATAAAAATCGACATGCACATTTTATTGTCCCCCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTg                      <  1:175539/149‑1 (MQ=255)
                                                                                                                                cGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGa                     <  1:144050/150‑1 (MQ=255)
                                                                                                                                       ggTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTg             <  2:969485/151‑1 (MQ=255)
                                                                                                                                          ccAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGaa           >  1:1194090/1‑150 (MQ=255)
                                                                                                                                            aGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGaaaa         >  2:998748/1‑150 (MQ=255)
                                                                                                                                                aaaGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCtt    <  2:276591/151‑1 (MQ=255)
                                                                                                                                                aaaGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCtt    <  2:947871/151‑1 (MQ=255)
                                                                                                                                                  aGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTg  >  2:512504/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                          
GCTTTTTATCGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTG  >  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: