Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,386,8830CT43.9% 45.7 / 51.9 68intergenic (+52/‑51)yfbL/yfbMputative M28A family peptidase/DUF1877 family protein
Reads supporting (aligned to +/- strand):  ref base C (10/27);  new base T (29/0);  total (40/28)
Fisher's exact test for biased strand distribution p-value = 1.49e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATACATGACACG  >  NC_000913/2386735‑2387032
                                                                                                                                                       |                                                                                                                                                             
aTACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTccc                                                                                                                                                             <  2:692773/151‑1 (MQ=255)
    cACTAGCCCGGTGATACCGCAGACAGATTGAATATAAGAAAAAGGATCAGGTAATCGATGGTGTTATAACTTTGTTAGACAACATTACATAATAAAAAGTTCTCGAGGCGCTTCCAGGAAAAGCGCCATGCATATTTTAATTGGCCTCATc                                                                                                                                                        <  2:386328/151‑1 (MQ=255)
          cccGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCCAAAAAACAACATGCCCATTTTTTTTTTCCCCCTCTtttt                                                                                                                                                   >  1:607228/1‑151 (MQ=255)
              gTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAAAAATAAAATTTTTCGGGGGCGGGCCCCAAAAAAAGAGAATACAAATTTTTTTTTtccctttttttttttaa                                                                                                                                               >  1:449572/1‑137 (MQ=255)
               tGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGa                                                                                                                                             <  2:725486/151‑1 (MQ=255)
               tGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGa                                                                                                                                             <  2:400212/151‑1 (MQ=255)
                   aCCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCAGAAAAAAAGACAAGCACATTTTTTTTTTTCTCATCTTTTTTTAGATGCg                                                                                                                                        >  2:733790/1‑151 (MQ=255)
                   aCCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGACACGCACAATTTTTTTTTCCTCCTCTGtttttaaaaggg                                                                                                                                          >  1:142143/1‑143 (MQ=255)
                   aCCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCACAAAAAACAACAATCACCATTTTTTTTTTCCCCTCTGtttttttcagc                                                                                                                                           >  2:285707/1‑143 (MQ=255)
                         gacagaTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAAAGACATGCCCATTTTTTTTTTCCTCTTCTGtttttaaaatggaaaaat                                                                                                                                    >  1:1075203/1‑137 (MQ=255)
                             gaTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAACGAAAAGGACATTTTTTTTTTCCTCACCTGTTTTTAGAAGGGaaaaaatgaa                                                                                                                             >  1:971512/1‑144 (MQ=255)
                                tGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTc                                                                                                                      <  1:39959/151‑1 (MQ=255)
                                           aaaaTGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACggagg                                                                                                             <  1:881955/148‑1 (MQ=255)
                                                ggCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATGAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGACAAGCACAATTTATTTTTTCTCCTCTGTTTTTAAATGCGGAACATGAAAATTAAACGAGGaaaaaaaa                                                                                                     >  2:457832/1‑151 (MQ=255)
                                                ggCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGACAAGGCCAATTTTTTTTTCCTCCTCTTTTTTTAGATGCGGAAAAAGAAAATTAAAGGGGGGaaaaaaa                                                                                                      >  1:529184/1‑149 (MQ=255)
                                                 gCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGGCAAAAAAAACCAAAAGCCCAATTTTTTTTTTCTCATCTGTTTTTAGAAGGGGAAAAATGAAATTAAAGGAGGGAaaataaa                                                                                                    >  2:798359/1‑151 (MQ=255)
                                                      ggTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAACGACAAGCACATTTTTTTTTTCCCCCTCTTTTTTTTAATGCGAAAAAATAGATTTAACCGAGGGAAAAAAAAAGGGa                                                                                               >  1:730843/1‑151 (MQ=255)
                                                      ggTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGGAAAAAAGGCAAGCACATTTTTTTTTTTCTCCTCTTTTTTTAAAAGGGGAACAAAAAAATTAACGGGGGGAAAAAAAATgggg                                                                                               >  1:370490/1‑150 (MQ=255)
                                                      ggTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAAAAAAAAAGACAAGCACATTTTTTTTTTCCTCCTCTTTTTTTAAATGCGGAAAAAGAAAATCAACGGGGGGaaaaaaaaaagga                                                                                               >  1:113103/1‑146 (MQ=255)
                                                      ggTAGTGGATGGTGTTATAACTTTGTTATACAACAGGAAATAATAAAAATTTCTCGAGGCGGGTCCGGGAAAAACCAAAAACACCATTTATTTTTTCTCTTCTGTTTTTAAATTCGGAAACACAAATTTTAAAGAGGGAAAATAAATgggg                                                                                               >  2:670129/1‑150 (MQ=255)
                                                            ggATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCCGAAAAAAACGAAAGCAAATTTTATTTTTTCTTCTCTGGTTTTAGAAGGGAAAAATAGAAATTAAAGCGGGGAAAATAAAATGGGatggtt                                                                                          >  2:290570/1‑146 (MQ=255)
                                                            ggATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCTAAAAAAAAACATGCAAATTATTATTTTCCTCCTCTGTTTTTAAATGCGAAAACATAAAATCAAAAGGGGGAAAAAAAAAAGGGagtaat                                                                                          >  2:814978/1‑146 (MQ=255)
                                                                 gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCGGAAAAAAAGGCATGCACCTTTTTTTTTTCCCCTTCTGTTTTTAGATGGGTAAAATGAAAATTCAAGGGGGGAAAAAAAATGGGGATGaattggta                                                                                    >  1:58548/1‑144 (MQ=255)
                                                                 gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAAAGACATGAACAATTTTTTTTTCCCCATCTGTTTTTAGATGCGGAAAAAGAAATTTAACCGAGGGAAAATAAATTGGAATGAGTGGGTa                                                                                    >  1:966249/1‑151 (MQ=255)
                                                                 gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGACATGCACATTTTTTTTTTCCTCCTCTGTTTTTAGATGGGAAAAATGAAAATTAAAGGGGGAAAAAAAAATGGGGAAGATTgggga                                                                                    >  1:964180/1‑148 (MQ=255)
                                                                 gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGGCAATACCATTTTTTTTTTCCTCCTCTGTTTTTAAAAGCGAAAAAAAAAAATTAAAGGAGGAAAAAAAAATAGGGAAGATTgggtt                                                                                    >  1:1121599/1‑148 (MQ=255)
                                                                 gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGAAATGCACATTTTTTTTTTCCTCATCTGTTTTTAAATGCGAAAAAAGAAAATTAAAGGAGGGAAAAAAAAAGGGGAAGATTGGGTa                                                                                    >  1:1103162/1‑151 (MQ=255)
                                                                 gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAAACAAGCAAATTTTATTTTTTCTCATCTGTTTTTAGATGCGAAACATGAAATTTAAAGGAGGGAAAATAAATGGGAATGAATgggtt                                                                                    >  2:1081407/1‑148 (MQ=255)
                                                                 gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCGAAAAAAACGGCAAGGACAATTTTTTTTTCCTCATCTGTTTTTTGAAGCGCAACATTAAATTCAACTGGGGGAAAAAAAATGGGAATGATTgggtt                                                                                    >  1:961643/1‑148 (MQ=255)
                                                                  tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCCGGAAAAACGACATGCCAATTTTTTTTTTTCTCCTCTGTTTTTAGATGCGAAACATGAAAATAAAAGGAGGAAAAAAAAATAGGAATGGTTGGCTAt                                                                                   >  2:112470/1‑151 (MQ=255)
                                                                   gtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGGCACAAAAAAAAACATGCACCAATTTTTTTCCCTCATCTTTTTTTAAAAACGAAACAACAAAATTAAAAGAGGGCAAATAAATTGGAATGATTGGTTAtt                                                                                  >  2:184781/1‑151 (MQ=255)
                                                                            ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGACAAGCACATTTTTTTTTTCCCCATCTGTTTTTAAATGCGAAACAATAAAATTAAACGAGGAAAAAAAAAAGGGGATGATTTGCTTTTTTTCAGaaa                                                                         >  1:932329/1‑151 (MQ=255)
                                                                              gTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAACCCATTCACCTTTTTTTTTTCCTCATCTGTTTTTAAATGCGAAAAATGAAAATTAACCGGGGGAAAATAAAAGGGGATGTTTTGCTTTTTTTTAGaaaat                                                                       >  1:116980/1‑149 (MQ=255)
                                                                              gTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAACGACACACAAACTTTATTTTTCCCCATCTGTTTTTAGATGCGTAAAAAAAAATTTAAATGAGGTAAAAAAAAAGGGGAAGAATGGGTGTTTTTCAGaaaat                                                                       >  1:952147/1‑149 (MQ=255)
                                                                              gTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAACGAAATGCACAATTTTTTTTTCCCCATCTGTTTTAAAATGGGAAAAATGAAAATTAAACGAGGGAAAAAAAAATGGGAAGGTTGGGTATTTTTTAGaaaat                                                                       >  1:788413/1‑149 (MQ=255)
                                                                               ttATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAACGGCAAGCACATTTTTTTTTTCCTCATCTTTTTTAAGATGCGAAAAATTAAATTCAACGGGGGGAAAAAAAATGGGGATGATTTGCTATTTTTCAGAAAATg                                                                      >  1:668887/1‑151 (MQ=255)
                                                                                    caacaGTAAATAATAAAAATTTCTCGAGGCGGGGTCACAAAAAAAGAAATGAACATTTTTTTTTTTCTCATCTGTTTTTAGATGCGAAAAATAAAAATCTAACGGGGAAAAAAAAAAGGGAAAGATTGTGTATTTTTTAGAAAATTACTCa                                                                 >  2:658278/1‑151 (MQ=255)
                                                                                      acaGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAAAGACATGCACATTTTTTTTTTCCTCCTTTTTTTTTAAATGGGAAAAATTAAAATTAAAGGGGGGAAAAAAAAATGGGAAGATTTGGTATTTTTTAGAAAATTACttaaa                                                               >  1:1150928/1‑147 (MQ=255)
                                                                                          tAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGaaa                                                             <  1:600813/149‑1 (MQ=255)
                                                                                                  aaaaaTTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCa                                                   <  2:788063/151‑1 (MQ=255)
                                                                                                    aaaTTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAAc                                                 <  1:113521/151‑1 (MQ=255)
                                                                                                    aaaTTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAAc                                                 <  1:986581/151‑1 (MQ=255)
                                                                                                                gCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTc                                       <  2:67765/149‑1 (MQ=255)
                                                                                                                gCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTc                                       <  2:1102662/149‑1 (MQ=255)
                                                                                                                gCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTc                                       <  1:931922/149‑1 (MQ=255)
                                                                                                                gCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTc                                       <  1:343474/149‑1 (MQ=255)
                                                                                                                gCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTc                                       <  1:1108341/149‑1 (MQ=255)
                                                                                                                 cGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTAc                                    <  2:657426/151‑1 (MQ=255)
                                                                                                                   ggTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTa                                     >  2:1070839/1‑148 (MQ=255)
                                                                                                                   ggTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTa                                     >  2:788529/1‑148 (MQ=255)
                                                                                                                          gAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGaa                                                                                              <  1:578949/85‑1 (MQ=255)
                                                                                                                          gAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGaa                                                                                              >  2:578661/1‑85 (MQ=255)
                                                                                                                                 cGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAACCCTTTGAt                     <  2:189262/151‑1 (MQ=255)
                                                                                                                                 cGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGAt                     >  2:978438/1‑151 (MQ=255)
                                                                                                                                   acttGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATgg                   >  2:1034832/4‑151 (MQ=255)
                                                                                                                                    cATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGTTTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGa                  >  1:247642/1‑151 (MQ=255)
                                                                                                                                           caTTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAatat            <  1:222683/151‑1 (MQ=255)
                                                                                                                                           caTTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAatat            >  2:832143/1‑151 (MQ=255)
                                                                                                                                           caTTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAatat            <  1:55913/151‑1 (MQ=255)
                                                                                                                                           caTTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAatat            <  2:78118/151‑1 (MQ=255)
                                                                                                                                           caTTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAatat            >  2:289402/1‑151 (MQ=255)
                                                                                                                                           caTTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAatat            <  2:675035/151‑1 (MQ=255)
                                                                                                                                             ttttattttCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTAGTGAAAAACCTTTGATGGATAATATAc          <  2:532880/151‑1 (MQ=255)
                                                                                                                                             ttttattttCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATAc          <  2:248649/151‑1 (MQ=255)
                                                                                                                                             ttttattttCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATAc          <  1:80018/151‑1 (MQ=255)
                                                                                                                                             ttttattttCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATAc          <  1:444139/151‑1 (MQ=255)
                                                                                                                                             ttttattttCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATAc          <  2:425540/151‑1 (MQ=255)
                                                                                                                                             ttttattttCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATAc          <  2:457786/151‑1 (MQ=255)
                                                                                                                                                   atttCCCTCATCTGTTTTAAGATGCGAAACATGAAATGCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATACATGac     <  2:35612/150‑1 (MQ=255)
                                                                                                                                                      tCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATACATGACACg  >  2:1122643/1‑151 (MQ=255)
                                                                                                                                                       cccTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGAAGAAATTGACTCAGAAAAGATTAAGCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATACATGACACg  >  2:902469/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                             
ATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATACATGACACG  >  NC_000913/2386735‑2387032

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