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

ATGATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGG  >  NC_000913/4536048‑4536343
                                                                                                                                                     |                                                                                                                                                  
aTGATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTa                                                                                                                                                    <  2:820640/150‑1 (MQ=255)
  gATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAgg                                                                                                                                                  <  1:996984/150‑1 (MQ=255)
     aCAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAAt                                                                                                                                               <  1:492277/150‑1 (MQ=255)
          agtagtAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATgaaag                                                                                                                                          <  2:901905/150‑1 (MQ=255)
           gtagtaAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATgaaaga                                                                                                                                         <  2:804003/150‑1 (MQ=255)
             agtaAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTTGAATgaaagaaa                                                                                                                                       >  1:602528/1‑150 (MQ=255)
               taAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTTGGATGaaaaaaaaa                                                                                                                                     >  1:519471/1‑149 (MQ=255)
                aaaGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATg                                                                                                                                    <  1:837364/150‑1 (MQ=255)
                 aaGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATgca                                                                                                                                  <  1:357157/151‑1 (MQ=255)
                      cATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTGGAATGAAAGAAAAagcaagca                                                                                                                              >  1:634583/1‑150 (MQ=255)
                      cATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTTGAATGaaaaaaaaa                                                                                                                                     >  2:631491/1‑142 (MQ=255)
                       atatTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCaat                                                                                                                            <  1:104187/151‑1 (MQ=255)
                       atatTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCaat                                                                                                                            <  2:1073498/151‑1 (MQ=255)
                       atatTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCaat                                                                                                                            <  1:1016806/151‑1 (MQ=255)
                         atTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGTTTTTTGAATGAAAGAAAAAaaaagcaaaa                                                                                                                           >  2:535190/1‑148 (MQ=255)
                         atTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCaataa                                                                                                                          <  2:647373/151‑1 (MQ=255)
                            ttACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTTGAATGAAAGaaaaaacaaaaaaaaata                                                                                                                        >  1:678046/1‑136 (MQ=255)
                             tACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTTGTATGaaaaaaaaaaaaaaaaaaaaaa                                                                                                                        >  2:676501/1‑135 (MQ=255)
                              aCCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAataatata                                                                                                                      <  2:728351/150‑1 (MQ=255)
                              aCCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTAGGTAAGaaaaaaaaaacaaaaaata                                                                                                                           >  1:618476/1‑134 (MQ=255)
                               ccAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTCTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTGGGATGAAAAAAAAAGCAAACAaaaatataa                                                                                                                     >  2:1042354/1‑150 (MQ=255)
                                        aGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTTGAATTAAAAAAAAAAAAAACAATAATAaaaaaaaagcc                                                                                                             >  2:272790/1‑144 (MQ=255)
                                             aCTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTTGAATGAAAAAAAAAAAAAAAAAAAATATAAAAAAGGCCTAAtttt                                                                                                      >  2:167618/1‑150 (MQ=255)
                                               tttCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTTGGATTAAAAAAAAAACAAAAAAAAATATAATAAATGGCTAAtttttt                                                                                                    >  2:61340/1‑151 (MQ=255)
                                                ttCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCt                                                                                                     <  2:189688/149‑1 (MQ=255)
                                                ttCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTGGGATGAAAAAAAAAACAAGCAATAATATAAAAAAGGCCTTATtttttt                                                                                                   >  1:981297/1‑151 (MQ=255)
                                                   cGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGaa                                                                                                <  2:55370/151‑1 (MQ=255)
                                                    gTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGTTTTTTGAATGAAAAAAAAAACAAAGAAAAAAATAAAAAAGGCCTTAAtttttt                                                                                                   >  2:656950/1‑147 (MQ=255)
                                                    gTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTGGGATGAAAAAAAAAAGAAACAAAAAAATAATAAAAGGCTAATtttttt                                                                                                   >  1:778777/1‑147 (MQ=255)
                                                     ttttATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTGGAATGAAAAAAAAAACAAACAAAAATAAAAAAAATGGCTAATtttttt                                                                                                   >  1:51786/1‑146 (MQ=255)
                                                        tATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAAcatc                                                                                            <  2:519456/150‑1 (MQ=255)
                                                        tATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAAcata                                                                                            <  2:882140/150‑2 (MQ=255)
                                                         aTTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAAcatcat                                                                                          <  1:495857/151‑1 (MQ=255)
                                                            gCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTaa                                                                                       <  1:365654/151‑1 (MQ=255)
                                                            gCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTaa                                                                                       <  2:903630/151‑1 (MQ=255)
                                                            gCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTa                                                                                        <  1:1191810/150‑1 (MQ=255)
                                                                   gCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATatc                                                                                <  1:829144/151‑1 (MQ=255)
                                                                    ccTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTTGTATTAAAAAAAAAACAAGAAAAAAAATAATAAAAGCCTAATTTTTTTTAAAATCATTAAACaaa                                                                                  >  2:322439/1‑146 (MQ=255)
                                                                      tttCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTTGGATGAAAGAAAAAGGAAACAAAAAAATAAAAAAAGGCCAATTTTTTTTAAAAACATTaaaaa                                                                                    >  1:54592/1‑141 (MQ=255)
                                                                       ttCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTTGTAAGAAAAAAAAATCAAGCAATAATATAATAATGGTCCAAATTTTTTGAACATCATTAAACATatcat                                                                              >  1:238579/1‑149 (MQ=255)
                                                                        tCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATatcatca                                                                            <  2:825204/150‑1 (MQ=255)
                                                                         ccTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATatcatcat                                                                           <  2:773770/150‑1 (MQ=255)
                                                                            gATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGt                                                                        <  2:983364/150‑1 (MQ=255)
                                                                             aTTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTTGAATGAAAAAAAAAGCAAGCAAAAATATAAAAATAGCCTTATTTTTTTTTACATTATTAAGCATATCATCATCATCt                                                                      >  2:276518/1‑151 (MQ=255)
                                                                             aTTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTGGAATGAAAAAAAATACAAGCAAAAATATAATAAAAGCCTAATTTTTTTTAAAATCATTAAACATATCATTATTGt                                                                        >  1:1195136/1‑149 (MQ=255)
                                                                              ttATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTTGGATTAAAAAAAAAGCAAGCAAAAAAATAATAAATGGCTAATTTTTTTTTACAACATTAAAAATATCATCATCGTCt                                                                      >  2:277810/1‑150 (MQ=255)
                                                                                ataAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTAtt                                                                   <  2:848304/151‑1 (MQ=255)
                                                                                  aaaTCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGTTTTTGGAATGAAAAAAAAAACAAGCAATAATATAATAAATGGCTAATTTTTTTTAAAATCATTAAACATATCATCATCGTCTAtttt                                                                 >  1:580309/1‑150 (MQ=255)
                                                                                  aaaTCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTg                                                                 <  1:448985/151‑1 (MQ=255)
                                                                                   aaTCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGTTTTTTGGATCTAAAAAAAAGGAAGCCAAAAAATTAAAATTGCCCCATTTTTTTGTAAATTTTTTAACATAATAACATCATCTAtttt                                                                 >  2:454001/1‑149 (MQ=255)
                                                                                   aaTCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGc                                                                <  1:806009/151‑1 (MQ=255)
                                                                                      cTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTTTAATGAAAAAAAATTAAAGCAATAATATAATAAAGGCCTAATTTTTTTTTACATTATTTAACAAATAAATATCGTCTAtttt                                                                 >  2:497184/1‑146 (MQ=255)
                                                                                         aTTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGaa                                                           <  1:1187571/150‑1 (MQ=255)
                                                                                                aCCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCg                                                    <  2:964812/150‑1 (MQ=255)
                                                                                                 ccTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGAc                                                  <  2:678030/151‑1 (MQ=255)
                                                                                                    aaTCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGTTTTTTGGATGAAAAAAAAAGCAAGCAAAAATATAATAAAAGGCTAATTTTTTTTAAAATCATTAAACATATCATCATCTTCTATTTTCCAGAATCCCTCGACtt                                                >  2:587277/1‑150 (MQ=255)
                                                                                                    aaTCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTta                                               <  2:392147/151‑1 (MQ=255)
                                                                                                      tCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGTTTTTGGAATGAAAAAAAAAGCAAGCAATAATAAAAAAAAAGCCTAATTTTTTTTAAAATCATTAAGCATATTATCATCGTCTATTTTTCAGAATACGTCGACtt                                                >  2:24292/1‑148 (MQ=255)
                                                                                                             ccTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGGTTTTGTAGTGAAAAAAAAAACAAGGCAAAAAATAAAAAAAGGCGAATTTTTTTTTAAAACATTAAGTAAATCTTCATAGTCTATTTTTCAGAATCCGTCGACTTATaa                                            >  2:256559/1‑145 (MQ=255)
                                                                                                               tttGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTGGAATGAAAGAAAAAGGACACAAAAATATAAAAAATGCCTAATTTCTTTGTACATCATTAAGCATATTATCATCGTCTATTTTCCAGAATCCGTCGACTTATAATCTcccg                                     >  1:82237/1‑148 (MQ=255)
                                                                                                                 tGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGaa                                   <  2:699638/150‑1 (MQ=255)
                                                                                                                  gCATAATGCGGCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAAc                                  <  2:286512/150‑1 (MQ=255)
                                                                                                                   cATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACgg                                <  2:756239/151‑1 (MQ=255)
                                                                                                                     tAATGCGCCTGGGATTAATGGCGGGGGCTTTTTTGGATTAAAAAAAAAAAAAAAAAAAATAAAAAAAAAGTCTAATTTTTTTTTACAACAATAAACATATTAACATTTTCTATTTTTCAGAATTCGTCGGCTTATAAttaccaaaa                                   >  2:715702/1‑138 (MQ=255)
                                                                                                                          cgcCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCa                         <  1:1139388/151‑1 (MQ=255)
                                                                                                                           gcCTGGGATTAATGGCGGGTGGTTTTTTGTATGAAAGAAAAAGCAAGCAATAATATAATAAAAGCCTATTTTTTTTTAAAAACATTAAAAAAATAATTATCCTTTATTTTTCCGAATTCGTCGCCTTATAAATTCCAGAAAGGCcta                            >  2:130931/1‑147 (MQ=255)
                                                                                                                           gcCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAg                        <  2:798241/151‑1 (MQ=255)
                                                                                                                               gggATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACg                    <  1:913230/151‑1 (MQ=255)
                                                                                                                                 gATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGt                  <  2:225532/151‑1 (MQ=255)
                                                                                                                                     aaTGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGtgct               <  1:564992/150‑1 (MQ=255)
                                                                                                                                     aaTGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGtgct               <  1:634341/150‑1 (MQ=255)
                                                                                                                                           gggTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAAc         >  2:40367/1‑150 (MQ=255)
                                                                                                                                                 ttttAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTgg  >  1:605626/1‑151 (MQ=255)
                                                                                                                                                 ttttAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTgg  >  1:291449/1‑151 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
ATGATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGG  >  NC_000913/4536048‑4536343

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: