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

ATGATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGCCCA  >  NC_000913/4536048‑4536347
                                                                                                                                                      |                                                                                                                                                       
aTGATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAg                                                                                                                                                        <  2:1031581/151‑1 (MQ=255)
       aGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATgaa                                                                                                                                                 <  1:947675/151‑1 (MQ=255)
       aGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTGGaa                                                                                                                                                     >  1:23037/1‑147 (MQ=255)
       aGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTCATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATgaa                                                                                                                                                 <  2:239749/151‑1 (MQ=255)
         gAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATgaaag                                                                                                                                               <  2:49655/151‑1 (MQ=255)
           gtagtaAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTGGGATgaaaaa                                                                                                                                              >  1:571656/1‑148 (MQ=255)
             agtaAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATgaaagaaa                                                                                                                                            <  1:212644/150‑1 (MQ=255)
              gtaAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAAt                                                                                                                                          <  2:4294/151‑1 (MQ=255)
              gtaAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAAt                                                                                                                                          <  1:407407/151‑1 (MQ=255)
              gtaAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAAt                                                                                                                                          <  2:244389/151‑1 (MQ=255)
                     gCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATgcaagca                                                                                                                                   <  1:321715/151‑1 (MQ=255)
                       atatTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTGGGATGAAAGAAAAAgcaa                                                                                                                                      >  2:860511/1‑146 (MQ=255)
                        tatTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTGGAAAGAAAAAAAAAgcaagcaa                                                                                                                                  >  1:203242/1‑149 (MQ=255)
                           tttACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTTGGATGAAAAAAAAAgcaaaaaaaaa                                                                                                                               >  1:922109/1‑142 (MQ=255)
                             tACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGTTTTTGGAATGaaaaaaaaagacaaaaaaaatat                                                                                                                            >  1:604503/1‑135 (MQ=255)
                             tACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAataata                                                                                                                             <  1:226768/149‑1 (MQ=255)
                              aCCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAataatataa                                                                                                                          <  2:521228/151‑1 (MQ=255)
                              aCCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAataatata                                                                                                                           <  1:278204/150‑1 (MQ=255)
                                   tGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTTGAAAGAAAGAAAAAAAAAACAATAAAAtaaaaa                                                                                                                      >  1:653607/1‑149 (MQ=255)
                                      tAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGc                                                                                                                  <  1:437928/150‑1 (MQ=255)
                                       aaGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCt                                                                                                                <  1:485598/151‑1 (MQ=255)
                                       aaGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCt                                                                                                                <  1:557132/151‑1 (MQ=255)
                                       aaGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCt                                                                                                                <  1:64590/151‑1 (MQ=255)
                                        aGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTa                                                                                                               <  1:535825/151‑1 (MQ=255)
                                              cTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCACCTGGGATTAATGGCGGGTGGTTTTTTGGATGAAAAAAAAAGAAAGGAAAAAAAAAATAATTGGCTAAttt                                                                                                           >  1:212236/1‑149 (MQ=255)
                                                ttCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCt                                                                                                         <  2:9156/149‑1 (MQ=255)
                                                 tCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCttt                                                                                                       <  2:370457/150‑1 (MQ=255)
                                                      tttATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAAca                                                                                                  <  1:483964/150‑1 (MQ=255)
                                                      tttATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAAca                                                                                                  <  2:1064136/150‑1 (MQ=255)
                                                         aTTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAAcatca                                                                                               <  2:982485/150‑1 (MQ=255)
                                                            gCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGTTTTTGGGATGAAAAAAAAAACAAGCAATAATATAATAAAGGCCTAATTTTTTTTTAAATCAtt                                                                                             >  1:437287/1‑149 (MQ=255)
                                                              gAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTTGAATGGAAAAAAAAGAAAGGTATAAAAAAATAAAAGACTAATttttttttatataactata                                                                                          >  2:384242/1‑137 (MQ=255)
                                                                  cGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCatat                                                                                     <  1:1037973/151‑1 (MQ=255)
                                                                  cGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCata                                                                                      <  1:1025573/150‑1 (MQ=255)
                                                                  cGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCata                                                                                      <  1:16313/150‑1 (MQ=255)
                                                                   gCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTTGGATTAAAAAAAAAAAAAACAAAAATATAAAAATAGTCTTATTTTTTTTTTAATCATTAAACata                                                                                      >  2:975109/1‑149 (MQ=255)
                                                                   gCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTTGAATAAAAAAAAAAGCAAGCAATAATATAAAAATAGCCTAATTTTTTTTTACATCAtttaaaaaa                                                                                      >  2:797677/1‑142 (MQ=255)
                                                                      tttCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGCGGGTTTTTGGAATGAAAGAAAAAGCAAGCAAAAATATAAAAATGGCCCAATTTTTTTTAAAAACATTAAGCa                                                                                        >  1:290192/1‑144 (MQ=255)
                                                                       ttCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGTTTTTTGGAAGAAAAAAAAAACAAACAAAAAAATAATAAAAGGCTAATTTTTTTTAACAACATTTAACata                                                                                      >  1:58689/1‑145 (MQ=255)
                                                                        tCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTGGAATTAAAGAAAAAGAAAGCAAAAATATAATAAAAGCCTAATTTTTTTTAACAACATTAAACATatcatca                                                                                >  2:1038566/1‑150 (MQ=255)
                                                                         ccTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATatcatcatc                                                                              <  2:515401/151‑1 (MQ=255)
                                                                           tGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGCGTTTTTTTGGAAGGAAAAAAAATCAAGCAAAAATATAATAAAAGGCCTATTTTTTTTAGACTCACTAATCATatc                                                                                    >  2:368062/1‑143 (MQ=255)
                                                                            gATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCAAATCATCATCg                                                                             <  2:827053/149‑1 (MQ=255)
                                                                                 taAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTGGGATGAAAAAAAAAGAAAGAAAAAAAATAAAAAAAGCCTAATTTTTTTTTAAAACATTTAGAATATTATCATCGTCTAttt                                                                      >  1:1046172/1‑151 (MQ=255)
                                                                                  aaaTCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTAttt                                                                      <  2:1046048/150‑1 (MQ=255)
                                                                                        tATTGTTTACCTAATGATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGTTTTTTGGATGGAAGAAAAAACAAGCAAAAATATAAAAAAGAACTAATTTTTTTTTAAAAAATTAAAAAAATCATCATCGGTTAtttt                                                                     >  2:761400/1‑144 (MQ=255)
                                                                                         aTTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTTGAATGAAAAAAAAAGCAAGCAAAAAAATAAAAAAAGCCTAATTTTTTTTAAAATCATTAAACAAATCATTATCGTCTATTTTCCAGaa                                                               >  1:801593/1‑150 (MQ=255)
                                                                                           tGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGTTTTTTGGAATGAAAAAAAAAGGAAACAAAAAAATAAAAAATGCCTAATTTTTTTTAACAACATTAAACATATCATCATCGTCTTTTTTCCAGaa                                                               >  1:111641/1‑148 (MQ=255)
                                                                                               tACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGTTTTTTTTATGAAAGAAAAAGCAAGCAATAATATAATAATGGGCTTATTTTTTTTTACATCATTAAGCATATAATCATATTCTttttttcaaaa                                                               >  2:160612/1‑137 (MQ=255)
                                                                                               tACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTTGGATGAAAGAAAAAGCAAGCAAAAATATAAAAAATGGCTAAATTTTTTTAAAATCATTTAACATATTATAAACCTCTATTTTTCAGAAACCCTc                                                         >  1:905890/1‑150 (MQ=255)
                                                                                                 ccTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGa                                                       <  2:980684/150‑1 (MQ=255)
                                                                                                    aaTCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACtt                                                    <  2:617123/150‑1 (MQ=255)
                                                                                                    aaTCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGTTTTTTGGGAAGGAAAAAAAAAAAAACAAAAAAATAATAAAAGCCTAATTTTTTTTTACATCATTAAACAAAATAACATTTTCTTTTTTTCACAATACGTCGGCtt                                                    >  1:498427/1‑150 (MQ=255)
                                                                                                       cATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTGGAATGAAAAAAAAAGCAAGCAAAAAAATAAAAAAGGGCTAATTTTTTTTTACAACATTAAACATATCATCATCGTCTTTTTTTCAGAAACCGTCGGCTtat                                                  >  1:266558/1‑149 (MQ=255)
                                                                                                          tttCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATc                                              <  2:35048/150‑1 (MQ=255)
                                                                                                            tCCTTTGCATAATGCGCCTGGGATTAATGGCGGGCGCTTTTTGGGGAGAAAGAAAATTCAAGCTAAAATATAAAAATAGCCTACTTTCTTTGGACAAACTTTAGCCAATTATTCAGGTCTATTTTACAGAATTCGTTGACCTATaaaaacc                                           >  1:1041369/1‑146 (MQ=255)
                                                                                                                ttGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGa                                        <  2:571552/150‑1 (MQ=255)
                                                                                                                  gCATAATGCGCCTGGGATTAATGGCGGGTGGGTTTTTGGATGAAAAAAAAAGCAAGAAATAAAATAATAAAAGGCTAATTTTTTTTAAAATCATTAAAAATATCATCATCGTCTATTTTTCCGAATACGCCCGCTTTTAATCACCAGaa                                       >  2:727491/1‑149 (MQ=255)
                                                                                                                          cgcCTGGGATTAATGGCGGGTGCTTTTTGGAATGAAAAAAAAAGCAAGAAAAAAAATAATAATAGGCTAATTTTTTTTTAAATCATTAAGAATATCAACATCGTCTATTTTTCAGAAAACGTCGACTTTTAATTCCCAGAACAGCat                                 >  1:534038/1‑147 (MQ=255)
                                                                                                                                        ggCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCg                >  2:939738/1‑150 (MQ=255)
                                                                                                                                        ggCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGa               <  1:880504/151‑1 (MQ=255)
                                                                                                                                          cGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAAc             <  2:543754/151‑1 (MQ=255)
                                                                                                                                                     tAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGccc   >  2:953734/1‑151 (MQ=255)
                                                                                                                                                      aGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGCCCa  <  2:203196/151‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                       
ATGATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGCCCA  >  NC_000913/4536048‑4536347

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