Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,536,1970AT42.2% 84.1 / 62.0 71intergenic (+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/35);  new base T (30/0);  total (36/35)
Fisher's exact test for biased strand distribution p-value = 2.37e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATGATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTG  >  NC_000913/4536047‑4536342
                                                                                                                                                       |                                                                                                                                                 
cATGATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTa                                                                                                                                                   <  2:261048/151‑1 (MQ=255)
 aTGATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAg                                                                                                                                                  <  2:1058803/151‑1 (MQ=255)
   gATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTgg                                                                                                                                                 >  2:776914/1‑150 (MQ=255)
      aCAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAAt                                                                                                                                              <  1:772742/150‑1 (MQ=255)
          gAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTTGGATgaaa                                                                                                                                          >  2:429313/1‑150 (MQ=255)
              agtaAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATgaaagaaa                                                                                                                                      <  1:772972/150‑1 (MQ=255)
               ggaaaGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGTTTTTATTAATtaaaaa                                                                                                                                        >  1:1068046/3‑145 (MQ=255)
                 aaaGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATgc                                                                                                                                  <  2:976946/151‑1 (MQ=255)
                   aGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATgcaa                                                                                                                                <  2:90386/151‑1 (MQ=255)
                   aGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTCATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTTGGATGaaaaaaaaagaa                                                                                                                                 >  2:36525/1‑146 (MQ=255)
                         tatTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTTGGATTTAAAAAAAAgcaaacaa                                                                                                                            >  2:840091/1‑149 (MQ=255)
                           ttttACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGTAGCTTTTTGGACTGAAAGaaaaaaaaagaaata                                                                                                                          >  1:87096/1‑138 (MQ=255)
                              ttccAATGTTAAGGTGACTTTCCGTTTTATTTCGAGTCGCCTTACCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCAAGGATTAATGGCGGGGGCTTTTAGGAATGAAAGAAAATGCAATCAataatat                                                                                                                      <  1:144573/148‑1 (MQ=255)
                              tACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAataatat                                                                                                                      <  1:358170/150‑1 (MQ=255)
                              tACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAataatat                                                                                                                      <  1:558976/150‑1 (MQ=255)
                               aCCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAataatataa                                                                                                                    <  2:1009353/151‑1 (MQ=255)
                               aCCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAataatata                                                                                                                     <  1:429310/150‑1 (MQ=255)
                                ccAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAataatataat                                                                                                                   <  2:793583/151‑1 (MQ=255)
                                   aTGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGTTTTTTGGATGAAAAAAAAAAAAAACAataata                                                                                                                       >  1:489564/1‑144 (MQ=255)
                                     gTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATg                                                                                                               <  2:478571/150‑1 (MQ=255)
                                       tAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTTGAATGAAAGAAAAAACAAGCAAAAATAtaataaaagcc                                                                                                            >  2:179927/1‑146 (MQ=255)
                                       tAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGc                                                                                                             <  1:739668/150‑1 (MQ=255)
                                        aaGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTGTTTGGAATGGAAAAAAAATCAAGGGATAAAAtaaaaa                                                                                                                 >  2:279234/1‑145 (MQ=255)
                                             gACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTTGAATGAAAAAAAAAACAAGCAATAATATAATAATAGCCTAAttt                                                                                                      >  2:544253/1‑151 (MQ=255)
                                              aCTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGTTTTTTGAATGAAAAAAAAAACAAACAAAAATAtaaaaaaaggcgaatttt                                                                                                     >  1:156884/1‑139 (MQ=255)
                                              aCTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTc                                                                                                     <  2:926088/151‑1 (MQ=255)
                                               cTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTGGAATGAAAAAAAAAGGAAACAAAAATATAATAAAGGCCTAAttttt                                                                                                    >  2:468900/1‑149 (MQ=255)
                                                tttCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTTTGATTAAAAAAAAAGAAAACAATAATATAATAAAGGGCTAAtttttt                                                                                                   >  1:171840/1‑151 (MQ=255)
                                                 ttCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCtt                                                                                                   <  1:894381/150‑1 (MQ=255)
                                                 ttCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCtt                                                                                                   <  1:1107899/150‑1 (MQ=255)
                                                 ttCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCt                                                                                                    <  2:4324/149‑1 (MQ=255)
                                                  tCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTTGGATTAAAAAAAAAAAAAACAATAATATAATAAAAGGCTTATttttttt                                                                                                 >  1:698916/1‑150 (MQ=255)
                                                      ttttATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTGTGATGAAAAAAAATGCACGAAAAAAAATAATAAAAGACTCATtttttt                                                                                                  >  2:1174767/1‑146 (MQ=255)
                                                       tttATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTTGAATGAAAAAAAAAACAAGGAAAAAAATAATAAATGGCTAATTTCttt                                                                                                  >  2:45466/1‑145 (MQ=255)
                                                                 gTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTTGTATTAAAAAAAAAGCAAGCAATAAAATAAAAAAGGACTAATTTTTTTTTAAATCAtttaaca                                                                                   >  1:201098/1‑145 (MQ=255)
                                                                   cGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTGGGATGAAAAAAAAAGCAAGCAATAATATAATAAAAGGCTAATTTTTTTTTACATCATTAAACata                                                                                 >  2:1173421/1‑150 (MQ=255)
                                                                      cTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGTTTTTTGTATTAAAAAAAAAACAAACAATAAAAAAATAAAAGCCCTATTTTTTTTTAAATCATTTAGCa                                                                                   >  1:559301/1‑145 (MQ=255)
                                                                       tttCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGGTTTTTGAATGAAAGAAAAAAAAAGCAATAAAATAAAAATTGGCTAATTTTTTTTTAGATTATTTAAAATatcat                                                                             >  2:17073/1‑150 (MQ=255)
                                                                        ttCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATatcatc                                                                            <  2:1091908/150‑1 (MQ=255)
                                                                         tCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATatcatca                                                                           <  2:483215/150‑1 (MQ=255)
                                                                         tCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATatcatca                                                                           <  2:968584/150‑1 (MQ=255)
                                                                            tGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTTGAATTAAAAAAAAAGAAAACAATAATATAAAAAAAACCTAATTTTTTTTTAAATCATTAAACATATCATCATCGt                                                                       >  1:139846/1‑151 (MQ=255)
                                                                              aTTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCt                                                                     <  1:337680/151‑1 (MQ=255)
                                                                               ttATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAAca                                                                                             <  1:428031/126‑1 (MQ=255)
                                                                                 ataAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTTTCATTAAAAAAAAAACAAAAAAAAATATAAAAATGGGTTCATTTTTTTTAACATTATTTAACATATCATCATCCTCTAtt                                                                  >  1:488423/1‑151 (MQ=255)
                                                                                 ataAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCg                                                                        <  1:424513/145‑1 (MQ=255)
                                                                                 ataAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTAt                                                                   <  2:256445/150‑1 (MQ=255)
                                                                                 ataAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTTGAATGAAAAAAAAAGCAAGCAATAATATAAAAAAAGCCTAAATTTTTTTTACATCATTAAACa                                                                                   >  2:424516/1‑134 (MQ=255)
                                                                                   aaaTCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTTTAATTAAAAAAAAAGCAAGAAAAAAAATAAAAAAAGCCTTATTTCTTTTTAAATCAATAAAAATATAATTATAGTCTAttt                                                                 >  1:625873/1‑150 (MQ=255)
                                                                                   aaaTCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTAttt                                                                 <  1:216423/150‑1 (MQ=255)
                                                                                     aTCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTTGGATTAAAAAAAAAAAAAACCAAAAAATAAAAATAGGCTTATTTTTTTTTAAATAATTAAACATATAATCATCGTCTAttttt                                                               >  1:1077916/1‑148 (MQ=255)
                                                                                         tATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCa                                                             <  2:139848/148‑1 (MQ=255)
                                                                                         tATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTGGGATGAAAAAAAAAGCAAACAAAAATATAATAAAAGCCTAATTTTTTTTAAAATCATTAAACATATCATCATCGTCTATTTTTCAGa                                                           >  2:818013/1‑150 (MQ=255)
                                                                                          aTTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAAt                                                         <  1:524465/151‑1 (MQ=255)
                                                                                               ttACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGt                                                     <  2:746366/150‑1 (MQ=255)
                                                                                                 aCCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGa                                                  <  1:53103/151‑1 (MQ=255)
                                                                                                  ccTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTTGAATGAAAAAAAAAACAAGAAATAAAATAATAATAGGCTAATTTTTTTTTACATCATTAAGAATATCATCATCGTCTATTTTCCAGAATCCGTCGAc                                                 >  1:190520/1‑151 (MQ=255)
                                                                                                     aaTCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTta                                              <  1:592765/151‑1 (MQ=255)
                                                                                                      aTCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGTTTTTTTTGTTGATAAAAAAAGTAAACAAAAAAATAAAAATAAAGTCATTTTTTTGAATAAAAATAAGAATATCATCCTCGTCTATTTTCCAGAATACCTCGCCTtat                                             >  1:167391/1‑151 (MQ=255)
                                                                                                            ttCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCg                                        <  2:268269/150‑1 (MQ=255)
                                                                                                                tttGCATAATGCGCCTGGGATTAATGGCGGGGGGCTTTTTGGATGAAAAAAAAAGCAAGCAAAAAAAAAATAAAAGGCTTATTTTTTTTTAAAACATTAAACATATAATCATCGTCTATTTTTCAGAATACGTCGCCTTATaaactc                                       >  1:92907/1‑143 (MQ=255)
                                                                                                                     aTAATGCGCCTGGGATTAATGGCGGGGGCCTTTTGGGATGAAAAAAAAAACAAGAAATAATATAAAAAAAGGCCAATTTTTTTTAAAATCATTAAACATATTATCATCCTCTATTTTTCAGAATCCGTCGACTTTTAATAACCAGAACgg                               >  2:405585/1‑150 (MQ=255)
                                                                                                                        aTGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCata                           <  1:949875/151‑1 (MQ=255)
                                                                                                                        aTGCGCCTGGGATTAATGGCGGGGGGTTTTTGGAAAGAAAAAAAAATCAAACAAAAAAATAAAAAATGGCTAATTTTTTTTAAAAACATTAAACATATCATTATCCTCTATTTTCCAGAATCCCCCGACTTATAATACCCACAACGGCa                             >  2:395536/1‑149 (MQ=255)
                                                                                                                           cgcCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCa                        <  2:1026685/151‑1 (MQ=255)
                                                                                                                                gggATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGa                     <  2:772572/149‑1 (MQ=255)
                                                                                                                                   aTTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGt                 <  1:160424/150‑1 (MQ=255)
                                                                                                                                        tGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCATTGAACAACATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGtgctg             >  2:792292/1‑149 (MQ=255)
                                                                                                                                         ggCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCg           >  1:889984/1‑150 (MQ=255)
                                                                                                                                                 gCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCtgt   >  2:999080/1‑151 (MQ=255)
                                                                                                                                                  cTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCtgtg  >  2:616707/1‑151 (MQ=255)
                                                                                                                                                  cTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCtgtg  >  1:399781/1‑151 (MQ=255)
                                                                                                                                                       |                                                                                                                                                 
CATGATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTG  >  NC_000913/4536047‑4536342

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: