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

ACCGCTTTTTATCGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTG  >  NC_000913/2386712‑2387011
                                                                                                                                                        |                                                                                                                                                        
acCGCTTTTTATCGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAg                                                                                                                                                          <  1:1202461/151‑1 (MQ=255)
      ttttATCGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACa                                                                                                                                                     <  1:528248/150‑1 (MQ=255)
           tCGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCACAAAAAACGACATGc                                                                                                                                                  >  1:247329/1‑148 (MQ=255)
             gTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATTGCTCGGGTCGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACatt                                                                                                                                            <  2:783622/151‑1 (MQ=255)
             gTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACatt                                                                                                                                            <  1:1032639/151‑1 (MQ=255)
                    aCAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAAAGGCAAGCCCattttttttt                                                                                                                                     >  2:62476/1‑151 (MQ=255)
                       aTACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGACAAGCACattttattttttc                                                                                                                                  >  1:137084/1‑148 (MQ=255)
                        tACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAGGAAATGCAAattttatttt                                                                                                                                     >  2:617017/1‑147 (MQ=255)
                        tACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGAGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGATCCATAAAAGCGACATTCAAattttttttt                                                                                                                                     >  2:1094783/1‑147 (MQ=255)
                           cACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGACAAGAAAatttttttttttc                                                                                                                                  >  1:799413/1‑144 (MQ=255)
                            aCTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCCCAAAAAAAGACAAAAACaattttttttt                                                                                                                                    >  2:214344/1‑143 (MQ=255)
                               tGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTg                                                                                                                           <  2:731085/150‑1 (MQ=255)
                                gCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAACAACAAGCAAAATTTATTTTTTCTCCTCtttt                                                                                                                         >  2:591780/1‑151 (MQ=255)
                                gCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGtt                                                                                                                         <  1:134384/151‑1 (MQ=255)
                                 cccGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGttt                                                                                                                        <  1:3636/151‑1 (MQ=255)
                                  ccGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGCAAAAACGAAAAGCACAATTTTTTTTTCCTCCTCtttttt                                                                                                                       >  2:417114/1‑151 (MQ=255)
                                  ccGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGtttt                                                                                                                       <  2:831234/151‑1 (MQ=255)
                                   cGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTa                                                                                                                      <  2:65292/151‑1 (MQ=255)
                                    ggTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTaa                                                                                                                     <  1:1210759/151‑1 (MQ=255)
                                    ggTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTaa                                                                                                                     <  2:523168/151‑1 (MQ=255)
                                     gTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAg                                                                                                                    <  1:411808/151‑1 (MQ=255)
                                        aTACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATg                                                                                                                 <  1:985220/151‑1 (MQ=255)
                                        aTACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATg                                                                                                                 <  1:441473/151‑1 (MQ=255)
                                         tACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGACAAGCCAATTTTTTTTTTCCTCCTCtttttt                                                                                                                       >  1:225101/1‑144 (MQ=255)
                                          aCCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCAAAAAAAACCACATGAAAATTTTTTTTTTCCTCTTCTGttttt                                                                                                                      >  1:346870/1‑143 (MQ=255)
                                          aCCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCg                                                                                                               <  1:123661/151‑1 (MQ=255)
                                          aCCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGACATGCACATTTTTTTTTTCCCCCTCTGTTTTTAGAAGCg                                                                                                               >  1:1145351/1‑151 (MQ=255)
                                            cGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGaa                                                                                                             <  1:571874/151‑1 (MQ=255)
                                             gCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGAAAAGAACCATTTATTTTTCCTCCTCTGTTTTTTGATGCGaa                                                                                                             >  2:1185253/1‑150 (MQ=255)
                                               agacagaTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCAGAAAAAAAGACAAGCACATTTTTTTTTTTCTTCTCTTTTTTTAAAAGCGaaaa                                                                                                           >  2:1128355/1‑149 (MQ=255)
                                                 acagaTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATg                                                                                                        <  2:442979/151‑1 (MQ=255)
                                                   agaTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAACGACATGCACATTTTATTTTTTCTCATCTGTTTTTAGATGCGaaaaa                                                                                                          >  1:507336/1‑145 (MQ=255)
                                                    gaTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAAAGACAAGCACATTTTTTTTTTCCTCCTCTGTTTTTAAATGCGGAACaaaaa                                                                                                      >  1:760740/1‑146 (MQ=255)
                                                    gaTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGaaa                                                                                                    <  2:465597/151‑1 (MQ=255)
                                                       tGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAAtt                                                                                                  <  1:385406/150‑1 (MQ=255)
                                                        gAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCa                                                                                                <  2:934618/151‑1 (MQ=255)
                                                         aaTTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCCGAAAAAACGACCTGGACAATTTTTTTTTCCCCATCTCTTTTTAGATGTGGAaaattaaatt                                                                                                  >  1:70480/1‑148 (MQ=255)
                                                         aaTTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAACGACATGCACATTTTTTTTTTCCTCCTCTGTTTTTAAAAGGGAAAAACGaaa                                                                                                    >  1:779497/1‑146 (MQ=255)
                                                         aaTTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCa                                                                                                <  1:846230/150‑1 (MQ=255)
                                                               cAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAACGAAAAGCACATTTTTTTTTTTCTCCTCTGTTTTTAGATACGAAAAATGAAAATTAACggggg                                                                                         >  1:927710/1‑151 (MQ=255)
                                                                 gAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGaa                                                                                       <  2:33511/151‑1 (MQ=255)
                                                                  aaaaTGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGTGCCCGAAAAAAAGAGATGTACAATTTATTTTCTCTCATCTGTTTTTAAAAGGGAAACAAAAAAAATAACggagg                                                                                         >  2:343470/1‑148 (MQ=255)
                                                                    aaTGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGaaaaa                                                                                    <  1:889499/151‑1 (MQ=255)
                                                                     aTGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAaaat                                                                                   <  1:1084781/151‑1 (MQ=255)
                                                                        gCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCAAAAAAAAGACATAGACCTTTTTTTTTTCCCCATCTGTTTTTTAATTGGCAACATGAAATTTAACGGAGGGAaaataa                                                                                 >  1:360416/1‑150 (MQ=255)
                                                                             ggTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGa                                                                          <  1:693774/151‑1 (MQ=255)
                                                                                aGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATg                                                                       <  2:160883/151‑1 (MQ=255)
                                                                                   ggATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGAt                                                                     <  1:644624/150‑1 (MQ=255)
                                                                                       ggTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTGCAGAAAAAAGGAAATGGACATTTTTTTTTTCCTCTTCTGTTTTTAAAAGCGAAAAATGAAATTTAACGGGGGAAAAAAAAAAGGGGATGGTTTGCt                                                                >  1:410170/1‑151 (MQ=255)
                                                                                        gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTa                                                               <  1:880449/151‑1 (MQ=255)
                                                                                         tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCCGAAAAAAAGACAAGAAAAATTTTTTTTTCCTCCTCTTTTTTTAAATGCGAAAAAAGAAAATTAAAGGAGGGAAAAAAAAAAGGGATGAATGGCTAt                                                              >  2:587029/1‑151 (MQ=255)
                                                                                         tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAAAGACATGCACATTTTTTTTTTCCTCCTCTGTTTTTAGAAGCGAAAAATGAAAATTAACGGAGGAAAAATAAATTGGGATGATTGGCTAt                                                              >  1:748856/1‑151 (MQ=255)
                                                                                           ttATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCAAAAAAAAGACAAGAAAAATTTTTTTTTTCTCATCTGTTTTTATATGGGAAAAAATAAATTAAACAGAGGGAAAAAAAAAGGGGATGATTTGATAttt                                                            >  2:1204989/1‑151 (MQ=255)
                                                                                                aCTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCCAAAAAAAGACATGCCCACTTTTTTTTTCCTCCTCTTTTTTTAAATTCGGAAAAAGAAAATCAACCGGGGGAAAAAAAATGGGGATGATTTGCTAttttt                                                          >  1:746537/1‑147 (MQ=255)
                                                                                                    tGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGaaa                                                    <  2:877081/150‑1 (MQ=255)
                                                                                                     gTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAAtt                                                  <  1:1123877/151‑1 (MQ=255)
                                                                                                     gTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGGAAATCACAATTTTTTTTTTCTCATCTGTTTTTAAATGCGAAAAATGAAATTTTACGGAGGAAAAAAAAAAAGGAATGATTGGCTATTTTTCAGaaa                                                    >  2:231730/1‑149 (MQ=255)
                                                                                                        ataCAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGa                                                <  2:247288/150‑1 (MQ=255)
                                                                                                        ataCAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGa                                                <  2:796982/150‑1 (MQ=255)
                                                                                                          acaacaGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACt                                              <  1:143847/150‑1 (MQ=255)
                                                                                                             acaGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGCAAAAACGGCAAGCACAATTTTTTTTTTCTCCTCTTTTTTTAAAAGGGAAAAAATAAAATTAACGGAGGGAAAATAAATAGGGAATATTGGTTATTTTTTAGaaaat                                                  >  2:646440/1‑141 (MQ=255)
                                                                                                             acaGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCa                                            <  2:174351/149‑1 (MQ=255)
                                                                                                                  aaaTAATAAAAATTTCTCGAGGCGGGTGCCACAAAAATGAAAAGCACAATTTTTTTTTCCCCCTACTTTTTTAAATGCGAAAAACGAAAATTTAACGAGGGAAAAAAAATTGGGAAGATGGGCTCTTTTTCAGAAAAGTAATCAGaaaa                                       >  2:863834/1‑149 (MQ=255)
                                                                                                                  aaaTAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGACAAGCACAATTTTTTTTTTCTCCTCTGTTTTTAAAAGCGGAACAAAAAAATTAAAGGGGGGAAAAAAAATGGGGAAGAATTGGTTTTTTTTAGAAAATAAATTAGaaaaa                                      >  2:999239/1‑149 (MQ=255)
                                                                                                                  aaaTAATAAAAATTTCTCGAGGCGGAGGCATAAAAAACGACCAGCACAATTTTTTTTTCCTCATCTGTTTTTTTATGGGGAACAATAAATTCAACGGAGGAAAAAAAAATGGGAATGGTTTGGTATTTTTTAGAAAATGACttcgca                                         >  2:1149378/1‑142 (MQ=255)
                                                                                                                       ataaAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTa                                  <  1:688591/149‑1 (MQ=255)
                                                                                                                          aaaaTTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCaa                             <  2:1179751/151‑1 (MQ=255)
                                                                                                                            aaTTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACt                           <  2:326781/151‑1 (MQ=255)
                                                                                                                             aTTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTg                          <  1:1095550/151‑1 (MQ=255)
                                                                                                                                   cGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCtt                       <  2:267483/148‑1 (MQ=255)
                                                                                                                                    gAGGCGGGTCCAGAAAAAACGACATGCACATTTTTTTTTTCCTCCTCTGTTTTTAGATGGGGAAAATGAAAATCAAAGGAGGAAAAAAAAATGGGAATGATTGGCTATTTTTGAGAAAATGACTCAGAAAAAATTAACCAActtctt                       >  2:939102/1‑147 (MQ=255)
                                                                                                                                       gCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTa                <  1:62482/151‑1 (MQ=255)
                                                                                                                                         gggTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACt              <  1:214860/151‑1 (MQ=255)
                                                                                                                                               ccAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGaa           <  1:204505/150‑1 (MQ=255)
                                                                                                                                                 aGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGaaaa         >  1:252670/1‑150 (MQ=255)
                                                                                                                                                        gCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTg  >  2:360480/1‑150 (MQ=255)
                                                                                                                                                        |                                                                                                                                                        
ACCGCTTTTTATCGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTG  >  NC_000913/2386712‑2387011

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