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

GATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGCC  >  NC_000913/4536050‑4536345
                                                                                                                                                   |                                                                                                                                                    
gATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGa                                                                                                                                                   <  2:890366/151‑1 (MQ=255)
   aCAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAAt                                                                                                                                                 <  2:595870/150‑1 (MQ=255)
     aGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATgaa                                                                                                                                              <  1:1040762/151‑1 (MQ=255)
       gAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATgaaag                                                                                                                                            <  1:406083/151‑1 (MQ=255)
         gtagtaAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATgaaagaa                                                                                                                                          <  2:156606/151‑1 (MQ=255)
           agtaAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATgaaagaaa                                                                                                                                         <  2:530374/150‑1 (MQ=255)
             taAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTTGGATGaaaaaaaa                                                                                                                                        >  1:688543/1‑149 (MQ=255)
              aaaGGGCATATTTTACCAATGTTGAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGTTTTTTGGATGAAAGaaaa                                                                                                                                        >  2:723495/1‑148 (MQ=255)
                aGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTGGAATTaaaaaaaaagaa                                                                                                                                    >  2:715127/1‑146 (MQ=255)
                 gggCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTTGAAAGAAAGaaaaagga                                                                                                                                    >  1:1125154/1‑145 (MQ=255)
                 gggCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTGTAATTaaaaaaaaaaaaa                                                                                                                                   >  2:1170935/1‑145 (MQ=255)
                      tatTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGTTTTTTATTATGAAAAAAAAAacaagaca                                                                                                                               >  2:1095944/1‑146 (MQ=255)
                         tttACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTTGAATGAAAGAAAAAGAAAACaaaaa                                                                                                                            >  1:966389/1‑149 (MQ=255)
                               aaTGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATAtaataa                                                                                                                    <  1:687872/151‑1 (MQ=255)
                                    tAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGc                                                                                                                <  2:161847/150‑1 (MQ=255)
                                     aaGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCt                                                                                                              <  2:641647/151‑1 (MQ=255)
                                       ggCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTaa                                                                                                            <  2:726831/151‑1 (MQ=255)
                                       ggCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTaa                                                                                                            <  1:613840/151‑1 (MQ=255)
                                         cGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAAtt                                                                                                          <  1:261475/151‑1 (MQ=255)
                                          gACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTGGAATGAAAAAAAAAGGAAAAAATAATATAATAATGGccaa                                                                                                             >  1:883203/1‑145 (MQ=255)
                                             tttCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTGGGATGCAAGAAAAAGAAAGGAATAAAAAAATAAATGCCTAAttttt                                                                                                       >  1:874694/1‑148 (MQ=255)
                                               tCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCttt                                                                                                     <  2:356870/150‑1 (MQ=255)
                                               tCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCttt                                                                                                     <  1:154900/150‑1 (MQ=255)
                                                 cGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGaa                                                                                                  <  1:346479/151‑1 (MQ=255)
                                                      tATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGTTTTTTGTATGAAAGAAAAAGCAAACAATAATATAATAAATGGCTAATttttttttaaa                                                                                                >  2:1137377/1‑143 (MQ=255)
                                                      tATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAAcatc                                                                                              <  1:188233/150‑1 (MQ=255)
                                                      tATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTTGAATGAAAAAAAAAGCAAACAAAAATATAAAAATAGGCCAATttttttttaaa                                                                                                >  1:531121/1‑143 (MQ=255)
                                                          gCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTGGAATGGAAAAAAAAGAAAGCAAAAATATAAAAAAAGGCTAATTTTTTTTTACATCAtt                                                                                           >  2:895785/1‑149 (MQ=255)
                                                           cGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGTTTTTTGAATGAAAAAAAAAAAAAACAAAAAAATAAAAAAAGCCTAATTTTTTTTAAaatta                                                                                             >  2:779685/1‑144 (MQ=255)
                                                           cGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTag                                                                                         <  1:674145/150‑2 (MQ=255)
                                                                cGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCatat                                                                                   <  2:69637/151‑1 (MQ=255)
                                                                cGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCatat                                                                                   <  1:423671/151‑1 (MQ=255)
                                                                 gCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATatc                                                                                  <  2:314556/151‑1 (MQ=255)
                                                                 gCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTGGAATGTAAGAAAAAGCAAACAAAAATATAAAAAAAGCCTAATTTTTTTTAACATCCTTTAGCa                                                                                      >  1:961692/1‑147 (MQ=255)
                                                                      tCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTTTGATGAAAAAAAAAACAAGACAAAATATAATAAAGGCCTAATTTTTTTTAAAAATATTAATCATatcatca                                                                              >  2:358691/1‑150 (MQ=255)
                                                                              ataAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTAtt                                                                     <  1:766667/151‑1 (MQ=255)
                                                                               taAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGTTTTTTGTATTAAAAAAAAAGCAAAAAAAAAAATAAAAATAGGCTAATTTTTTTTTACATCATTAATCATATTATCATCATCTAttt                                                                    >  2:860119/1‑151 (MQ=255)
                                                                                  aTCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGCTTTTTGAATGAAAAAAAAAGGAAAGAACAAAATAAAAAAGACCTAATTTTTTTTAACAACAATTAAAACATCATTATAGTTTATTTTcc                                                                 >  1:44865/1‑151 (MQ=255)
                                                                                   tCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTGGGATGAAAGAAAAAGCAAGCAATAATATAATAATAGGCTAATTTTTTTTTACATCATTAAACATATCATCATCGTCTtttt                                                                    >  1:945479/1‑147 (MQ=255)
                                                                                        ttGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTTTTATGAAAAAAAAAACAAGCAAAAAAAAAATAAATGCCTAAATTCTTTGTAAATAATTAAACAAATTATAATCCTTTTTTTTCCATaa                                                             >  2:383683/1‑149 (MQ=255)
                                                                                            ttACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTGGGATGAAAAAAAAAGCAAGCAATAATATAATAAAAGCCTAATTTTTTTTAAAATCATTTAACATATTATTATCCTCTATTTTTCAGAAACCGTc                                                       >  1:231677/1‑151 (MQ=255)
                                                                                            ttACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGt                                                        <  2:310145/150‑1 (MQ=255)
                                                                                            ttACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTGGGATGAAAAAAAAAGCAAACAAAAATATAAAAAAGGGCTAATTTTTTTTTACATCATTAAACAAATCATCATCCTCTTTTTTCCAGAATcc                                                          >  1:36562/1‑148 (MQ=255)
                                                                                             tACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTc                                                       <  2:398134/150‑1 (MQ=255)
                                                                                                 tAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTTGGATGAAAAAAAAAAGAAACAAAAATATAAAAAAAGCCCAATTTTTTTTTACATCATTAAAAATATCATCATCGTCTATTTTCCCGAATCCCTCGCCtt                                                  >  1:559045/1‑151 (MQ=255)
                                                                                                 tAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTGTTATGAAAAAAAAAGCAAGAAAAAAAAAAATAATAGGCTAATTTTTTTTTAAAAAATTAAAAAAAAAATCATCCTCTTTTTTTCAGGATCCGTCGACtt                                                  >  2:285493/1‑151 (MQ=255)
                                                                                                     cATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTtata                                               <  1:1147094/150‑1 (MQ=255)
                                                                                                     cATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTTGGGTGAAAAAAAAAACAAACAAAAATATAAAAATGGTCTTATTTTTTTTAAAATCATTTAGCATATCATCATCGTCTATTTTTCAGAAACCGTTGACTtat                                                >  2:457946/1‑149 (MQ=255)
                                                                                                          tCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTGGTATGAAAAAAAAAACAAGAAAAAATATAATAAAAGGCTAAATTCTTTTTTAATCATTAAACAAATAAAAATCGTCTATTTTTTTAAATTCATCGACTTATaaatgcc                                         >  2:253678/1‑146 (MQ=255)
                                                                                                                gCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTTGGATGAAAAAAAAAGCAAGGAAAAAAATAAAAATAGCCTAATTTTTTTTTAAATTATTAAGCATATCAACATCGTCTTTTTTTCCGAATCCGTCGACTTTTAAATACCAGaa                                     >  2:911248/1‑149 (MQ=255)
                                                                                                                 cATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACgg                                  <  1:280505/151‑1 (MQ=255)
                                                                                                                 cATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACgg                                  <  1:48228/151‑1 (MQ=255)
                                                                                                                  aTAATGCGCCTGGGATTAATGGCGGGTGCTTTTTGGAATGAAAGAAAAAGGAAGCAATAAAAAAAAAAATGCCTAATTTTTTTTAACAACATTTAACATATCATCATCGTCTCTTTTCCAGAATACCTCGGCTTTTAATCGCCAAAACgg                                  >  1:637051/1‑150 (MQ=255)
                                                                                                                  aTAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGc                                 <  1:723590/151‑1 (MQ=255)
                                                                                                                  aTAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGc                                 <  1:733768/151‑1 (MQ=255)
                                                                                                                    aaTGCGCCTGGGATTAATGGCGGGGGGTTTTTTGGATGAAAAAAAAAACAAGCAAAAAAATAAAAAATGGCTAATTTTTTTTTACATCATTAAGAATATCATCATCGTCTATTTTTCAGAATACCTCGACTTATAAAATCCAGAAAgg                                  >  2:133295/1‑148 (MQ=255)
                                                                                                                     aTGCGCCTGGGATTAATGGCGGGTGGTTTTTGGGAAGAAAAAAAAAGCAAGGAATAATATAAAAAAAGGCTAATTTTTTTTTAAATCATTAAACAAATTATAATTCTCTATTTTTCCGAATACGTCGACTTATAAtttcccgaa                                     >  2:774168/1‑136 (MQ=255)
                                                                                                                       gcgcCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATc                            <  2:44859/151‑1 (MQ=255)
                                                                                                                           cTGGGATTAATGGCGGGTGCTTTTTGGGATGGAAGAAAAAACAAGGAAAAATATAATAATAGCCTTATTTTTTTTTACATCATTAAGCAAATCATCATCTTCTTTTTTCCAGAATACGTCGACTTATAAATACCACAACGGCCTATCa                           >  2:1099442/1‑148 (MQ=255)
                                                                                                                                aTTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACgg                     <  2:652039/149‑1 (MQ=255)
                                                                                                                                  tAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGtgc                  <  1:715218/150‑1 (MQ=255)
                                                                                                                                   aaTGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGtgctg                <  2:64114/151‑1 (MQ=255)
                                                                                                                                   aaTGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGtgctg                <  1:53171/151‑1 (MQ=255)
                                                                                                                                   aaTGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGtgct                 <  2:1170207/150‑1 (MQ=255)
                                                                                                                                        cGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAAc           <  1:280583/151‑1 (MQ=255)
                                                                                                                                         gggTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACg          >  1:658049/1‑151 (MQ=255)
                                                                                                                                         gggTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACg          >  2:949923/1‑151 (MQ=255)
                                                                                                                                            tGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCt        <  2:127083/150‑1 (MQ=255)
                                                                                                                                                 ttAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGcc  >  1:999622/1‑151 (MQ=255)
                                                                                                                                                  tAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGcc  <  2:1064296/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGCC  >  NC_000913/4536050‑4536345

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