Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,536,197	0	A	T	46.0%	58.2 / 63.6	63	intergenic (+166/+417)	yjhR/nanS	pseudogene, helicase family; putative frameshift suppressor/9‑O‑acetyl N‑acetylneuraminic acid esterase
Reads supporting (aligned to +/- strand):  ref base A (3/31);  new base T (29/0);  total (32/31)
Fisher's exact test for biased strand distribution p-value = 7.14e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATGATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGCCCA  >  NC_000913/4536047‑4536347                                                                                                                                                       |                                                                                                                                                       cATGATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTa                                                                                                                                                        <  2:1163926/151‑1 (MQ=255)           gAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTAGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATgaaag                                                                                                                                              <  1:101532/151‑1 (MQ=255)               agtaAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATgaaagaaa                                                                                                                                           <  1:648937/150‑1 (MQ=255)               agtaAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATgaaagaaa                                                                                                                                           <  2:234956/150‑1 (MQ=255)                  aaaGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGTTTTTTGGATGaaaaaaaa                                                                                                                                          >  1:905367/1‑148 (MQ=255)                  aaaGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATg                                                                                                                                        <  2:1109897/150‑1 (MQ=255)                  aaaGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTGGGATGaaaaaaaaa                                                                                                                                         >  1:923253/1‑148 (MQ=255)                           atTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCaata                                                                                                                               <  1:564808/150‑1 (MQ=255)                            ttttACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGTTTTTTGAAATTAAAAAAAAACAAGGaataa                                                                                                                              >  2:280833/1‑150 (MQ=255)                             tttACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTTTGATGAAAAAAAAAGCAAAGaataa                                                                                                                              >  1:748227/1‑149 (MQ=255)                              ttACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTGGAATTAAAAAAAAAGCAAACaataa                                                                                                                              >  2:957346/1‑148 (MQ=255)                                 ccAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAataatataat                                                                                                                        <  2:200481/151‑1 (MQ=255)                                   aaTGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTTTAATGTAAGAAAAAACAAGCAaaaaaataaaa                                                                                                                       >  2:889817/1‑148 (MQ=255)                                      gTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATg                                                                                                                    <  2:42337/150‑1 (MQ=255)                                            gCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTTGGATGGAAAAAAAAAAAAGCAATAAaataataatatcc                                                                                                                 >  1:254407/1‑143 (MQ=255)                                              gACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGTTTTTGGTATTTAAAAAAAAAAAAACAAAAAAATAATAAATGGCTAAttt                                                                                                           >  2:103743/1‑151 (MQ=255)                                              gACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTTGTATTAAAAAAAAAAGAAAAAATAATAtaataaaaggctcattt                                                                                                           >  2:113798/1‑140 (MQ=255)                                                 tttCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTTGAATGAAAGAAAAAGCAAGCAAAAAAATAATAATTGCCTAAtttttt                                                                                                        >  1:110684/1‑151 (MQ=255)                                                  ttCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCt                                                                                                         <  1:51260/149‑1 (MQ=255)                                                   tCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCttt                                                                                                       <  2:826365/150‑1 (MQ=255)                                                    ccGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTTGGATTAAAAAAAATGTAAGCAAAAATATAATAATGGTCTTATttttttt                                                                                                      >  2:1183445/1‑149 (MQ=255)                                                     cGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGaa                                                                                                    <  2:1080117/151‑1 (MQ=255)                                                        tttATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGCGCGTTTTTGGGTGAAAAAAAATTCAAGAAATAATATAATAAAAGGCCAATTTCttt                                                                                                       >  1:582282/1‑145 (MQ=255)                                                           aTTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTTGAATTAAAGAAAAAGCAAGGAAAAATATAAAAAAAGGCTAAATTTTTTTaa                                                                                                    >  1:468451/1‑145 (MQ=255)                                                            ttGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGTTTTTTGAATGAAAAAAAAAACAAGCAATAATATAATAATAGCCTAATttttttttaaatca                                                                                               >  1:1313018/1‑141 (MQ=255)                                                                gAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGc                                                                                         <  2:644090/151‑1 (MQ=255)                                                                    cGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTTGAATTAAAAAAAATGCAAACAAAAATATAATAATAGCCTAATTTTTTTTAAAATCATTAAACata                                                                                      >  1:474957/1‑150 (MQ=255)                                                                    cGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTGGAATGAAAAAAAAAGCAAGCAATAATATAATAAAAGGCTAATTTTTTTTTACAACAtt                                                                                             >  1:1111095/1‑143 (MQ=255)                                                                        tttCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGTTTTTTTGTATGTAAAAAAAAGTAAACAATAAAATAAAAAATGGCTAATTTTTTTTATAATCATTaaaaa                                                                                        >  2:352597/1‑141 (MQ=255)                                                                           ccTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATatcatcatc                                                                              <  1:55239/151‑1 (MQ=255)                                                                              gATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTc                                                                           <  2:1000911/151‑1 (MQ=255)                                                                               aTTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGt                                                                            <  1:1116686/149‑1 (MQ=255)                                                                                ttATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTTGGATGAAAGAAAAAGAAAGCAAAAAAATAATAATGGCCTTATTTTTTTTAACATCATTAAAAATAATATCATTGTCt                                                                          >  1:639923/1‑150 (MQ=255)                                                                                 tataAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTTGGATGAAAAAAAAAAAAAAAAATAAAATAATAAAGGGCTAATttttttttaaaacatt                                                                                             >  1:613254/1‑120 (MQ=255)                                                                                  ataAATCTTATTGATTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTTGTATTGAAGAAAAAACAAGCAAAAATATAAAAATGAGCTTATTTTTTTTTAAATCAATAATCAAAATATCCACGTCTAtt                                                                       >  2:369070/1‑151 (MQ=255)                                                                                          tATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCa                                                                  <  1:1190310/148‑1 (MQ=255)                                                                                             tGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATc                                                             <  1:359461/150‑1 (MQ=255)                                                                                               tttACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGTTTTTGGAATGAAAGAAAAAACAACGAAAAAAATAAAAATAAACTAATTTTTTTTAAAATCATTAAACAAATCAACATCGCCTAtttttccgaa                                                               >  1:686000/1‑139 (MQ=255)                                                                                                ttACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGt                                                          <  2:517217/150‑1 (MQ=255)                                                                                                   ccTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTTGGATTTAAAAAAAAACAAGCAATAATATAATAATAGCCTTATTTTTTTTTACATCATTAATCCTATTATAATCGTTTTTTTTTCAGAATTCGTCTAc                                                      >  2:555019/1‑151 (MQ=255)                                                                                                   ccTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTTGAATGAAAAAAAAAACAAGCAAAAAAATAAAAATAGGCTAATTTTTTTTTAAATCATTAAACATATCATCATCGTCTATTTTCCCGAAACCCTc                                                         >  1:551915/1‑148 (MQ=255)                                                                                                   ccTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCCTTTTTGTATTAAAAAAAAAGAAAAAAATAATATAATAATAAACTAATTTTTTTTTTCATAATTAATAAAATCATCATTTTCTTTTTTCCCGAATTCCTCGa                                                       >  2:843428/1‑150 (MQ=255)                                                                                                    cTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTGGAAAGAAAAAAAAAACAAGCAATAATATAAAAATGGCCTAATTTTTTTTTAAATCATTAAACATATCATCATCCTCTATTTTCCCGAATCCCTCGAc                                                      >  1:256042/1‑150 (MQ=255)                                                                                                       aTCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGTTTTTGGAATGAAAAAAAAAACAAACAAAAATAAAAAAAAAGCCTAATTTTTTTTAAAATCATTAAACATAATATCATACTCTTTTTTTCCGAAACCGTCGCCtt                                                    >  1:343589/1‑149 (MQ=255)                                                                                                          aTTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATaa                                                <  2:529293/150‑1 (MQ=255)                                                                                                                   tGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTGGGATGAAAGAAAAAGCAAGGAATAATATAATAATGGGCTAATTTTTTTTTAAATCATTAAACATATCATCATCGTCTATTTTTCAGAATCCGTCGACTTATAATTGCCAGaa                                       >  1:38090/1‑150 (MQ=255)                                                                                                                     cATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACgg                                    <  1:889398/151‑1 (MQ=255)                                                                                                                         aTGCGCCTGGGATTAATGGCGGGTGGTTTTTTGGATGAAAGAAAATAAAAGCAAAAAAATAATAAAAGGCTAATTTTTTTTTAAAAAATTTAATATATCATCATCGTCTTTTTTTCAGAAACCGTCGGCCTTTAAAACCCAAAACGGCat                                 >  1:447729/1‑150 (MQ=255)                                                                                                                          tGCGCCTGGGATTAATGGCGGGTGCGTTTTTGAATGAAAGAAAAAGCAAGCAAAAAAATAAAAATGGCCTAATTTTTTTTTAAATAATTAAACAAAATATCATCGTCTATTTTTCAAAATCCTTCGACTTATAATAACCCGAACAGCata                                >  2:638014/1‑150 (MQ=255)                                                                                                                           gcgcCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCatat                               <  2:486747/150‑1 (MQ=255)                                                                                                                           gcgcCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATc                              <  1:939219/151‑1 (MQ=255)                                                                                                                             gcCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAg                            <  1:118805/151‑1 (MQ=255)                                                                                                                                 gtgATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCAGCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGAc                         <  2:960971/148‑1 (MQ=255)                                                                                                                                       aaTGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGtgct                   <  1:188111/150‑1 (MQ=255)                                                                                                                                       aaTGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGtgct                   <  1:98200/150‑1 (MQ=255)                                                                                                                                        aTGGCGAGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATAGCATCATCGTCGATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGtgctgc                 <  2:806661/151‑1 (MQ=255)                                                                                                                                          ggCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCg                <  1:1270997/150‑1 (MQ=255)                                                                                                                                          ggCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCg                >  1:825121/1‑150 (MQ=255)                                                                                                                                          ggCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGa               <  2:1233857/151‑1 (MQ=255)                                                                                                                                          ggCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGa               <  2:905793/151‑1 (MQ=255)                                                                                                                                                  cTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCtgtg       >  1:1050003/1‑151 (MQ=255)                                                                                                                                                    tttAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCtgtg       <  2:755681/149‑1 (MQ=255)                                                                                                                                                    tttAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCtgtg       <  2:305395/149‑1 (MQ=255)                                                                                                                                                      tAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGccc   >  2:1026344/1‑151 (MQ=255)                                                                                                                                                       aGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGccc   >  1:928005/1‑150 (MQ=255)                                                                                                                                                       aGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGCCCa  <  1:541576/151‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                       CATGATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGCCCA  >  NC_000913/4536047‑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: ‑