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.8%	57.4 / 54.1	62	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/30);  new base T (29/0);  total (32/30)
Fisher's exact test for biased strand distribution p-value = 1.22e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTG‑CTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGCCC  >  NC_000913/4536048‑4536346                                                                                                                                                       |                                                                                                                                                      aTGATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTG‑CTTTTAg                                                                                                                                                      <  1:876190/151‑1 (MQ=255) aTGATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTG‑CTTATAg                                                                                                                                                      <  1:876165/151‑1 (MQ=255)   gATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTG‑CTTTTAGGa                                                                                                                                                    <  1:192925/151‑1 (MQ=255)   gATACAGGAGTAGTAAAGGGAATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTTCCTAATCATTTTCCGTTGCATAATGTGCCTGGGATTAATGGCGGGTG‑CTTTTAgg                                                                                                                                                     <  1:254765/150‑1 (MQ=255)           agtagtAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTG‑CTTTTAGGAATgaaaga                                                                                                                                            <  1:214451/151‑1 (MQ=255)            gtagtaAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTG‑GTTTTTGGAATgaaaaa                                                                                                                                            >  2:949875/1‑148 (MQ=255)            gtagtaAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGG‑GTTTTTTGAATtaaaaa                                                                                                                                            >  2:610448/1‑148 (MQ=255)             tagtaAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTG‑CTTTTAGGAATgaaagaa                                                                                                                                           <  2:793179/150‑1 (MQ=255)               gtaAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTG‑CTTTTAGGAATGAAAGAAAAt                                                                                                                                        <  2:248550/151‑1 (MQ=255)                 aaaGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTG‑GTTTTTGGAATGAAAGaaaa                                                                                                                                         >  1:738739/1‑148 (MQ=255)                  aaGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGG‑GTTTTTTGGATGGAaaaaaaaac                                                                                                                                      >  2:776433/1‑147 (MQ=255)                  aaGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGG‑CTTTTTTGAAAGaaaaaaaaa                                                                                                                                        >  1:505215/1‑147 (MQ=255)                            tttACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTG‑CTTTTAGGAATGAAAGAAAATGCAAGCaataat                                                                                                                            <  1:850276/150‑1 (MQ=255)                             ttACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTG‑CTTTTAGGAATGAAAGAAAATGCAAGCAataata                                                                                                                           <  1:746465/150‑1 (MQ=255)                                ccAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGG‑CTTTTTTGAATGAAAAAAAAAGCAAGCAaaaatataa                                                                                                                        >  2:651885/1‑150 (MQ=255)                                    tGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGG‑CTTTTTGGGATGAAAAAAAAAGCAAGCAATAATAtaataa                                                                                                                     >  1:293891/1‑149 (MQ=255)                                        aaGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTG‑CTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCt                                                                                                               <  1:541422/151‑1 (MQ=255)                                         aGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTG‑CTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTa                                                                                                              <  2:167682/151‑1 (MQ=255)                                            cGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTG‑GTTTTTTGGATGTAAGAAAAAACAAACAAAAAAATAATAATAGCCTAAtt                                                                                                           >  2:785108/1‑151 (MQ=255)                                                tttCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGAGATTAATGGCGGGGGTGTTTTAGCAAGAAAAAAAAATCCAAG‑AAAAATATAAAAAAGGCCCAAtttttt                                                                                                       >  1:129305/1‑151 (MQ=255)                                                  tCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTG‑GTTTTTGGAATGAAAGAAAAAGCAAGAAAAAAAATAATAAAAGGCTAATttttttt                                                                                                     >  1:695728/1‑150 (MQ=255)                                                   ccGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTG‑CTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTg                                                                                                     <  2:289447/150‑1 (MQ=255)                                                   ccGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTG‑CCTTTTTGAATGAAAAAAAAATCAAGCAAAAAAATAATAATGGCCTCATtttttttt                                                                                                    >  2:381584/1‑149 (MQ=255)                                                     gTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTG‑CTTTTTTGACTTAAAGAAAAAGCGAGCGATAAAATAATAAAAGGCTGTTTTCttt                                                                                                      >  2:581970/1‑147 (MQ=255)                                                     gTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTG‑CTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAAc                                                                                                  <  1:783403/151‑1 (MQ=255)                                                       tttATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTG‑CTTTTTTGGATGAAAAAAAAAGCAAGCAATAATATAATAATTGCCTTATtttttt                                                                                                      >  1:50522/1‑145 (MQ=255)                                                        ttATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTG‑CTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAAcat                                                                                                <  2:283483/150‑1 (MQ=255)                                                              cGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTG‑CTTTTTTGGATGAAAAAAAAAGAAAGCAATAATATAATAATGGGCTAATTTTTTTTTACATCAtt                                                                                            >  1:896441/1‑148 (MQ=255)                                                              cGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGG‑GTTTTTTGTATGAAAGAAAAAGCAAGAAAAAATATAATAAATGGCTAATTTTTTTTaa                                                                                                   >  2:406760/1‑141 (MQ=255)                                                                aGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTG‑CTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCa                                                                                       <  1:343739/151‑1 (MQ=255)                                                                    gCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGG‑GTTTTTGGAATGGAAAAAAAAGCAAACAAAAAAATAATAATGGCCTCATTTTTTTTAAAATCATTaaaaaaa                                                                                     >  1:809628/1‑144 (MQ=255)                                                                      cTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGG‑CTTTTTGGGATGAAAAAAAAAGCAAACAATAATATAAAAAAAGACTAATTTTTTTTAAAATCATTaa                                                                                          >  1:308302/1‑142 (MQ=255)                                                                        ttCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGG‑GTTTTTGGAATGAAAAAAAAAGCAAACAAAAAAATAATAAATGGCTAATTTTTTTTAACAACAtttaacaaat                                                                                    >  1:910146/1‑138 (MQ=255)                                                                         tCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTG‑CTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATatcatca                                                                               <  1:776508/150‑1 (MQ=255)                                                                             gATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGG‑GTTTTTGGAATGAAAAAAAAAGCAAACAATAATATAATAATAGCCTAATTTTTTTTAACATCATTAAACATatcaacatcct                                                                           >  1:383181/1‑148 (MQ=255)                                                                                  taAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTG‑CTTTTTTTTATGAAAAAAAAATCAAGCAAAAATATAAAAATTGGCTAATTTTTTTGAAAATAATTAAACATATTATAAATCTCTAttt                                                                     >  2:641906/1‑151 (MQ=255)                                                                                   aaaTCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTG‑CTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTg                                                                    <  1:940308/151‑1 (MQ=255)                                                                                      tCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTG‑CTTTTTTGAATGAAAAAAAAAGCAAGCAATAATATAAAAAATGCCTTATTTTTTTTTACATAATTAAGCAAATCATCATCTTCTtttttt                                                                   >  2:704453/1‑147 (MQ=255)                                                                                      tCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTG‑CTTTTTGGAATAAAAAAAAAAACAAACAATAATTTAAAAAATGCCTAATTTTTTTTAAAAACATTAAACATATCATCATCGTCTAtttt                                                                    >  1:706757/1‑147 (MQ=255)                                                                                      tCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGG‑CTTTTTTGGAAGAAAAAAAAAGCAAGCAATAAAATAATAAAGGGCTAATTTTTTTTAAAAAAATTAAGCATATCATTATTGTCTAttt                                                                     >  1:687695/1‑147 (MQ=255)                                                                                         tATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTG‑CTTTTTTGAATGAAAGAAAAAGCAAACAATAAAATAAAAAAAGCCTAATTTTTTTTAACATCATTAAGCATATTATCATCGTCTttttttccga                                                               >  1:626524/1‑144 (MQ=255)                                                                                         tATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTG‑CTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGa                                                               <  1:203514/150‑1 (MQ=255)                                                                                              tttACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTG‑CTTTTTTGTATGAAAAAAAAAGCAAGCAAAAAAATAATAAAAGGCTAATTTTTTTTTACATCATTTAGAATATCATCATCATCTATTTTTCAGAATCCg                                                          >  2:522817/1‑150 (MQ=255)                                                                                                tACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTG‑CTTTTTGGAATGAAAAAAAAAACAAACAAAAATATAATAATAGGCTAATTTTTTTTAACATCATTAAACATATCATCATCGTCTATTTTCCAGAAACCGTc                                                        >  1:828470/1‑150 (MQ=255)                                                                                                tACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTG‑CTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTc                                                        <  2:10788/150‑1 (MQ=255)                                                                                                      aTCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTG‑CTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTtat                                                 <  1:24950/151‑1 (MQ=255)                                                                                                        cATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTG‑CTTTTTGGAATGAAAGAAAAAGCAAGCAATAATATAATAAATGCCTAATTTCTTTGAAAATCATTAAAAATATCATCATCGTCTTTTTTCCAGAATCCGTCGACTtat                                                 >  1:708297/1‑149 (MQ=255)                                                                                                        cATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTG‑CTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATaa                                               <  1:644968/151‑1 (MQ=255)                                                                                                           tttCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTG‑CTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATc                                             <  2:738668/150‑1 (MQ=255)                                                                                                             tCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTG‑CTTTTTGGGATGAAAAAAAAAGCAAGCAATAATATAAAAAATGCCTTATTTTTTTTTAAATCATTAAACATATAATCATCGTCTATTTTCCACAATACGTCGACCTATAATctc                                           >  2:352567/1‑148 (MQ=255)                                                                                                              ccTTTGCATAATGCGCCTGGGATTAATGGCGGGTG‑CTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCa                                         <  2:729205/151‑1 (MQ=255)                                                                                                                  tGCATAATGCGCCTGGGATTAATGGCGGGTG‑CTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGaa                                      <  2:681415/150‑1 (MQ=255)                                                                                                                    cATAATGCGCCTGGGATTAATGGCGGGTG‑CTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACgg                                   <  1:384853/151‑1 (MQ=255)                                                                                                                     aTAATGCGCCTGGGATTAATGGCGGGGG‑GTTTTTGGAATGAAAAAAAAAGGAAGCAATAATATAATAATGGGCTAATTTTTTTTAAAATCATTAAACATATTATCATCGTCTATTTTCCAGAATCCGTCGACTTTTAATAGCCAGAAAgg                                   >  1:456251/1‑150 (MQ=255)                                                                                                                      tAATGCGCCTGGGATTAATGGCGGGGG‑CTTTTTTGGATGAAAAAAAAAGCAAAGAAAAAAATAAAAATAGCCTAATTTCTTTTTAAATCATTAAACATATTAATATCGTCTATTTTCCAGAAACCGTCGGCTTTTAAtttccacaa                                      >  1:280606/1‑138 (MQ=255)                                                                                                                          gcgcCTGGGATTAATGGCGGGTG‑CTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATc                             <  1:561186/151‑1 (MQ=255)                                                                                                                                gggATTAATGGCGGGTG‑CTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACg                       <  2:562439/151‑1 (MQ=255)                                                                                                                                           cGGGTG‑CTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAAc            >  2:656207/1‑151 (MQ=255)                                                                                                                                             ggTG‑CTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGc          <  1:664325/151‑1 (MQ=255)                                                                                                                                              gTG‑CTTTTAGGAATGAAAGGAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCt         <  1:248899/151‑1 (MQ=255)                                                                                                                                                    tttAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCtgtg      <  2:153051/149‑1 (MQ=255)                                                                                                                                                     ttAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGcc   >  1:67531/1‑151 (MQ=255)                                                                                                                                                       aGGAAAGAAAGAAAATGCAAGCAGTAATATAATAATGGCATAATTTCTTTGAACACCAGAAAGCATATCATCATCGTCTATTTGCCGGAATCCGTCGACTTATAATCGCCAGAACGGCATATCATGTGGGTGGTGCGAACGCTGTGGtcc  <  2:811529/150‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                      ATGATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTG‑CTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGCCC  >  NC_000913/4536048‑4536346

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: ‑