Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,536,1970AT50.0% 46.8 / 88.2 74intergenic (+166/+417)yjhR/nanSpseudogene, helicase family; putative frameshift suppressor/9‑O‑acetyl N‑acetylneuraminic acid esterase
Reads supporting (aligned to +/- strand):  ref base A (6/31);  new base T (37/0);  total (43/31)
Fisher's exact test for biased strand distribution p-value = 6.98e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGCCC  >  NC_000913/4536054‑4536346
                                                                                                                                                 |                                                                                                                                                       
cAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTggga                                                                                                                                                     >  2:244481/1‑146 (MQ=255)
  ggAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATgaaa                                                                                                                                                <  2:413892/151‑1 (MQ=255)
          aaaGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGTTTTTGGAATGaaaaaaaa                                                                                                                                           >  2:1017128/1‑148 (MQ=255)
          aaaGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATgc                                                                                                                                        <  2:708892/151‑1 (MQ=255)
            aGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATgcaa                                                                                                                                      <  1:1023498/151‑1 (MQ=255)
               gCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATgcaagca                                                                                                                                   <  1:524617/151‑1 (MQ=255)
                 atatTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCaat                                                                                                                                 <  2:684146/151‑1 (MQ=255)
                  tatTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGTTTTTTGAATTAAAAAAAAAaaaagcaa                                                                                                                                  >  1:345692/1‑149 (MQ=255)
                      ttACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTGGGATGAAAAAAAAAGAAAGCaataaaa                                                                                                                             >  1:898202/1‑148 (MQ=255)
                       tACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTGGGATGAAAAAAAAAGAAAACaaaaaaa                                                                                                                             >  2:722310/1‑147 (MQ=255)
                       tACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAataatat                                                                                                                            <  2:686305/150‑1 (MQ=255)
                          cAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGTTTTTTGGATGAAAAAAAAAACAAACAaaaaaataa                                                                                                                         >  2:701737/1‑149 (MQ=255)
                             tGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGCTTTTTGAATTAAAAAAAAAGCAAGCAATAAAAaaataa                                                                                                                      >  2:1126495/1‑149 (MQ=255)
                                tAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGc                                                                                                                  <  2:430570/150‑1 (MQ=255)
                                tAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTGGGATGAAAAAAAAAAAAAGAAAAAAAAtaataa                                                                                                                      >  1:952508/1‑146 (MQ=255)
                                 aaGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGTTTTTTGGATGAAAAAAAAAAGAAAGAATAATATAAAAAAgggc                                                                                                                 >  1:125485/1‑148 (MQ=255)
                                 aaGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCt                                                                                                                <  1:520334/151‑1 (MQ=255)
                                     cGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTTGAATGAAAGAAAAAGCAAGCAATAATATAAAAAAAGGCTAAtt                                                                                                            >  2:578447/1‑151 (MQ=255)
                                      gACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAAttt                                                                                                           <  2:107618/151‑1 (MQ=255)
                                         tttCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTTTAATTTAAAAAAAATCAAGCAAAAATATAAAAAATGCCTAAtttttt                                                                                                        >  1:324711/1‑151 (MQ=255)
                                          ttCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCt                                                                                                         <  1:315845/149‑1 (MQ=255)
                                           tCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCttt                                                                                                       <  1:900610/150‑1 (MQ=255)
                                             cGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGaa                                                                                                    <  2:449993/151‑1 (MQ=255)
                                             cGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGaa                                                                                                    <  1:467958/151‑1 (MQ=255)
                                             cGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTGGAATGAAAAAAAAAAAAAACAATAATAAAATAAATGCCTAAAtttttt                                                                                                       >  2:604525/1‑148 (MQ=255)
                                               ttttATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTTGGAAGAAAGAAAAAACAAGCAAAAAAATAAAAAAAGACTAATttttttt                                                                                                      >  2:372298/1‑146 (MQ=255)
                                               ttttATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGATGTTTTTTTTGAATAAAAAAAAAGGAAACAAAAAAATAATAATATTCTAATtttttttt                                                                                                     >  1:795492/1‑146 (MQ=255)
                                                       cGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAg                                                                                          <  1:695017/151‑1 (MQ=255)
                                                             gCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATGATATAATAATGTCCTAATTTCTTTGAACATCATTAAGCATatc                                                                                   <  2:1177986/151‑1 (MQ=255)
                                                               cTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGTTTTTGGGATGAAAGAAAAAAGAAACAAAAAAATAAAAAATGGCTAATTTTTTTTTAAATCATTaaaaa                                                                                        >  2:597044/1‑142 (MQ=255)
                                                               cTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTTGGATTAAAAAAAAAACAAGCAAAAAAATAAAAAAAGCCTAATTTTTTTTTAAAACATTAAACa                                                                                        >  1:170899/1‑145 (MQ=255)
                                                                  tCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGTTTTTGGAATGAAAAAAAAAGCAAGCAAAAAAATAAAAAAAGGCTAATTTTTTTTTAAAACATTAAAAAAatcatca                                                                               >  2:553887/1‑150 (MQ=255)
                                                                   ccTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATatcatcatc                                                                             <  1:485116/151‑1 (MQ=255)
                                                                   ccTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATatcatcat                                                                              <  1:501113/150‑1 (MQ=255)
                                                                      gATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTc                                                                          <  2:814892/151‑1 (MQ=255)
                                                                        ttATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTGGGATGAAAAAAAAAGCAAACAATAATATAAAAAAAGCCTAATTTTTTTTAACAACATTAAACAAATCATCATCGTCt                                                                         >  2:838860/1‑150 (MQ=255)
                                                                          ataAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTGGAATTAAAAAAAAAGCAAACAATAATATAAAAATGGGCTAATTTTTTTTTAAATCATTAAAAAAATTATCAACGTCTAtt                                                                      >  2:1022032/1‑151 (MQ=255)
                                                                            aaaTCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTTGAATGAAAAAAAAAACAAAGAATAAAATAATAAAAGGCTAATTTTTTTTTACAACATTTAGAAAATCATCCTCGTCTAtttt                                                                    >  1:959321/1‑150 (MQ=255)
                                                                            aaaTCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTg                                                                    <  1:163208/151‑1 (MQ=255)
                                                                             aaTCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGCTTTTTGGAACAAAAAAAAAATCAGCAACAATATAAAAAATGGCCAATTTCTTTTAAAATACTTAAAAATattatcaaactccattt                                                                     >  1:1159499/1‑139 (MQ=255)
                                                                             aaTCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTTGTATGAAAAAAAAAAAAAGCAATAAAATAAAAAATGGCTTATTTTTTTTTACATTATTAAACATATTATTATCCTCTtttttt                                                                   >  1:275674/1‑149 (MQ=255)
                                                                             aaTCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTTGTATTAAAAAAAAAAAAAACAAAAAAATAAAAAAAGGCTAAATTTTTTTTAAATCATTAAGCAAATAATAATCATCTAttt                                                                     >  1:5907/1‑149 (MQ=255)
                                                                               tCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTTGGAAGAAAAAAAAAGAAAGCAAAAAAATAAAAAAAGCCTAATTTTTTTTAAAATCATTTAGCATATTATCATCGTCTAtttt                                                                    >  2:522936/1‑147 (MQ=255)
                                                                                 ttATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCGTTTTGTTGTGAAAAGAAATTCAAAGAATAAAAAAATAAAGGCATTATTTCTTTGAACATAATTAATAATATTATTATCATCTtttt                                                                     >  1:725808/1‑145 (MQ=255)
                                                                                   aTTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTGGGATGAAAAAAAAAACAAACAAAAATATAATAAAAGCCTTATTTTTTTTTAAATTATTAAAAATATCATCATCGTCTAtttttccgaa                                                              >  2:245941/1‑143 (MQ=255)
                                                                                     tGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGTTTTTTGGATGAAAAAAAAAAAAAACAATAAAATAATAAAAGGCTAATTTTTTTTTAAATCATTAAACAAATTATTCTCTTCTTTTTTTCAGa                                                               >  1:1133771/1‑147 (MQ=255)
                                                                                     tGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTGGTATTAAAAAAAAAACAAACAATAATATAATAAAAGGCTTTTTTTTTTTTAAAACATTAATAATATTATCATCGTTTAtttttcacaa                                                              >  1:425493/1‑141 (MQ=255)
                                                                                          aCCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGTTTTTTGAATTAAAAAAAAAACAAGCAAAAAAATAAAAAATGGCTTATTTTTTTTAACAACATTTAACATAATATCATCGTCTATTTTTCAGAAACCGTc                                                        >  1:444851/1‑149 (MQ=255)
                                                                                           ccTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTTGAATTAAAAAAAAAAAAAACAAAAATATAATAATAGTCTTATTTTTTTTTACATCATTAAACATATCATCATATTCTATTTTTCAGAATCCCTCGAc                                                     >  1:47137/1‑151 (MQ=255)
                                                                                           ccTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGAc                                                     <  1:614006/151‑1 (MQ=255)
                                                                                            cTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTGGAATGAAAAAAAATGCAAGCAATAATATAATAATGGGCTAATTTTTTTTTACATCATTAAGCATATCATCATCGTCTATTTTCCACAATCCCTCGAc                                                     >  2:1121767/1‑150 (MQ=255)
                                                                                             tAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGTTTTTTTGGAAGAAAGAAAATATACACAAAAAAATAATAAAAGCCACATTTTTTTGAACAAAAATAATAATACAATAATCGTTTATTTTCCCGCATcactcggact                                                   >  2:1143988/1‑142 (MQ=255)
                                                                                              aaTCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTta                                                  <  1:317754/151‑1 (MQ=255)
                                                                                                    tttCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTGGGATGAAAAAAAAAACAAGAAATAATATAATAATAGCCTAATTTTTTTTAACATCATTTAAAATATTATCATCGTCTATTTTTCACAATCCGTCGACTTATaa                                               >  2:363300/1‑148 (MQ=255)
                                                                                                       ccTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCa                                         <  2:989081/151‑1 (MQ=255)
                                                                                                       ccTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCa                                         <  1:57945/151‑1 (MQ=255)
                                                                                                       ccTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCa                                         <  2:912210/151‑1 (MQ=255)
                                                                                                         tttGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTGGGATGAAAAAAAAAAAAAGCAATAAAATAATAATAGCCTAATTTCTTTTTAAATCATTAAGAAAATTATTATCATCTATTTTTCAGAATACCTCGACTTATAATAGcc                                          >  1:497704/1‑148 (MQ=255)
                                                                                                         tttGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTGGGAAGAAAAAAAAAGCAAAAAATAAAATAAAAAATGGCTAATTTTTTTTAAAAACATTAAAAATATCATCATCGTCTATTTTCCAGAAACCGTCGACTTATAATAGcc                                          >  2:767542/1‑148 (MQ=255)
                                                                                                           tGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGaa                                      <  1:932170/150‑1 (MQ=255)
                                                                                                             cATAATGCGCCTGGGATTAATGGCGGGGGCTTTTAGGGATGGAAAAAAAAGAAAACAATAAATTAAAAAAGGGCCAATTTTTTTTTACAACATTAAACATAATCACATTGTCTTTTTTTCAGAAATCGTCGACTTTTaaattccacaaagg                                   >  1:300603/1‑139 (MQ=255)
                                                                                                              aTAATGCGCCTGGGATTAATGGCGGGTGCTTTTTTGGATGAAAAAAAAAGCAAACAAAAATATAAAAAAAGGCTAATTTTTTTTAAAATCATTAAACATATCATCATCGTCTATTTTTCAGAATCCGCCGACTTATAATTACCAGaa                                      >  2:340346/1‑147 (MQ=255)
                                                                                                               tAATGCGCCTGGGATTAATGGCGGGGGGTTTTTGGGATGAAAGAAAAAGAAAGGAAAAAAATAAAAAAAGGCTAATTTTTTTTTAAATAATTAAGAATATAATCATCGTCTTTTTTTCAGAAAACGTCGGCTTATAAAATCCAGaa                                      >  1:827798/1‑146 (MQ=255)
                                                                                                                 aTGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCata                               <  1:527063/151‑1 (MQ=255)
                                                                                                                    cgcCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCa                            <  1:682344/151‑1 (MQ=255)
                                                                                                                         gggATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGa                         <  1:1059118/149‑1 (MQ=255)
                                                                                                                               aaTGGAGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGtgctg                 >  1:778982/1‑151 (MQ=255)
                                                                                                                                 tGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCg               >  1:1072072/1‑151 (MQ=255)
                                                                                                                                   gCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGaa             <  1:372293/151‑1 (MQ=255)
                                                                                                                                    cGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAAc            <  1:805919/151‑1 (MQ=255)
                                                                                                                                    cGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAAc            >  1:328410/1‑151 (MQ=255)
                                                                                                                                            cTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCtgtg      >  2:210516/1‑151 (MQ=255)
                                                                                                                                               ttAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGcc   >  1:195329/1‑151 (MQ=255)
                                                                                                                                                tAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGcc   <  1:30733/150‑1 (MQ=255)
                                                                                                                                                 aGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGccc  <  1:121553/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                       
CAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGCCC  >  NC_000913/4536054‑4536346

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