Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,536,2070GA41.3% 106.9 / 36.1 75intergenic (+176/+407)yjhR/nanSpseudogene, helicase family; putative frameshift suppressor/9‑O‑acetyl N‑acetylneuraminic acid esterase
Reads supporting (aligned to +/- strand):  ref base G (11/33);  new base A (31/0);  total (42/33)
Fisher's exact test for biased strand distribution p-value = 6.20e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGCCCAAAACGTTCCA  >  NC_000913/4536064‑4536357
                                                                                                                                                 |                                                                                                                                                        
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                                                                                                               <  1:467958/151‑1 (MQ=255)
                                   cGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGaa                                                                                                               <  2:449993/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:912210/151‑1 (MQ=255)
                                                                                             ccTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCa                                                    <  2:989081/151‑1 (MQ=255)
                                                                                             ccTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCa                                                    <  1:57945/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)
                                                                                                                                          aaTGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGCCCaaaa         <  2:740810/151‑1 (MQ=255)
                                                                                                                                           aTGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGCCCAAAAc        <  1:546833/151‑1 (MQ=255)
                                                                                                                                           aTGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGCCCAAAAc        >  1:1126184/1‑151 (MQ=255)
                                                                                                                                               aagaaaATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGCCCAAAACGtt     >  2:327581/1‑150 (MQ=255)
                                                                                                                                                 gaaaATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGCCCAAAACGTTCCa  <  2:336918/151‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                        
AAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGCCCAAAACGTTCCA  >  NC_000913/4536064‑4536357

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: