Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	659,940	0	T	A	43.9%	88.8 / 34.4	66	I93F (ATT→TTT)	lipA	lipoate synthase
Reads supporting (aligned to +/- strand):  ref base T (9/28);  new base A (29/0);  total (38/28)
Fisher's exact test for biased strand distribution p-value = 3.72e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCGCAGGTCATCACGGTCAACGGAGGTGATAACCACATAACGCAGCGCCATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGG‑GTACAAATAGCGCCGAGGATCATAAACGTTGCTGTGCCGTGGTTGAAGC‑ATTCCGCCAGGTTAGGGCAGGAGGCTTCCTCGCAGACAGAATGCAGGCCATTTTTGCGCATTGCGGCTTTGATGCCCTGGATACGTGTAGAGTCCGCTGGAAGCTTG  >  NC_000913/659790‑660088                                                                                                                                                        |                                                                                                                                                      cgcgCAGGTCATCACGTTCAACGGAGGTGATAACCACATAACGCAGCGCCATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGG‑GTACAAAt                                                                                                                                                       <  2:303471/151‑1 (MQ=255)    gcAGGTCATCACGGTCAACGGAGGTGATAACCACATAACGCAGCGCCATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGG‑GTACAAATAgc                                                                                                                                                    <  2:907082/151‑1 (MQ=255)     cAGGTCATCACGGTCAACGGAGGTGATAACCACATAACGCAGCGCCATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGG‑GTACAAATAgcg                                                                                                                                                   <  2:500919/151‑1 (MQ=255)     cAGGTCATCACGGTCAACGGAGGTGATAACCACATAACGCAGCGCCATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGG‑GCAGAAATa                                                                                                                                                      >  2:819939/1‑148 (MQ=255)        gTCATCACGGTCAACGGAGGTGATAACCACATAACGCAGCGCCATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGG‑GTACAAATAGCGCCg                                                                                                                                                <  2:311600/151‑1 (MQ=255)              aCGGTCAACGGAGGTGATAACCACATAACGCAGCGCCATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGG‑GTACAAATAGCGCCGAGGATc                                                                                                                                          <  2:233687/151‑1 (MQ=255)                  tCAACGGAGGTGATAACCACATAACGCAGCGCCATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGG‑GGACAAAAAACGCCCAGGGGCATa                                                                                                                                       >  1:117587/1‑150 (MQ=255)                   cAACGGAGGTGATAACCACATAACGCAGCGCCATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGG‑GTACAAATAGCGCCGAGGATCATaaa                                                                                                                                     <  2:425659/151‑1 (MQ=255)                   cAACGGAGGTGATAACCACATAACGCAGCGCCATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGG‑GTACAAATAGCGCCGAGGATCATaaa                                                                                                                                     <  1:298409/151‑1 (MQ=255)                    aaCGGAGGTGATAACCACATAACGCAGCGCCATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGG‑GTACAAATAGCGCCGAGGATCATAAAc                                                                                                                                    <  1:173500/151‑1 (MQ=255)                      cGGAGGTGATAACCACATAACGCAGCGCCATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGG‑GTACAAATAGCGCCGAGGATCATAAACGt                                                                                                                                  <  2:220231/151‑1 (MQ=255)                       ggAGGTGATAACCACATAACGCAGCGCCATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGG‑GTACAAATAGCGCCGAGGATCATAAACGtt                                                                                                                                 <  2:877665/151‑1 (MQ=255)                        gAGGTGATAACCACATAACGCAGCGCCATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGG‑GGACAAAAAAAGACGAGGGGCATAAACgat                                                                                                                                 >  2:576907/1‑148 (MQ=255)                          ggTGATAACCACATAACGCAGCGCCATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGG‑GTACAAATAGCGCCGAGGATCATAAACGTTGCt                                                                                                                              <  2:385970/151‑1 (MQ=255)                            tGATAACCACATAACGCAGCGCCATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGG‑GGACAAAAAACGCCCGGGGACATaaaaattgt                                                                                                                               >  2:346114/1‑142 (MQ=255)                             gATAACCACATAACGCAGCGCCATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGG‑GTACAAATAGCGCCGAGGATCATAAACGTTGCtgta                                                                                                                           <  2:384/151‑2 (MQ=255)                               tAACCACATAACGCAGCGCCATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGG‑GGACAAAAAAAGCCGAGGGGTATAAAAGATGCTGTGcc                                                                                                                         >  1:505311/1‑151 (MQ=255)                                aaCCACATAACGCAGCGCCATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGG‑GGACAAAAAGCGCCGGGGGGAATAAACGTTTCAGTGcc                                                                                                                         >  2:976132/1‑150 (MQ=255)                                 aCCACATAACGCAGCGCCATATCGGCAATGGTCTTCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGG‑GTACAAATAGCGCCGAGGATCATAAACGTTGCTGTGCCGt                                                                                                                       <  2:947564/151‑1 (MQ=255)                                     caTAACGCAGCGCCATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGG‑GTACAAATAGCGCCGAGGATCATAAACGTTGCTGTGCCGTGGtt                                                                                                                   <  2:370508/151‑1 (MQ=255)                                     caTAACGCAGCGCCATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGG‑GTACAAATAGCGCCGAGGATCATAAACGTTGCTGTGCCGTGGtt                                                                                                                   <  1:941326/151‑1 (MQ=255)                                        agcGCTGCACCATCGCGTCAATGGTCTGCGCCATCTTCACTGGTTCACTGGCATCAGGAGCTACCGGGCGACCGTGAGCAGCATCACAGAACGGACAGCGGCGG‑GTACATATAGCCCCGAGGATCATAAACTTTGCTGTGCCGTGTTTgc                                                                                                                 <  2:968513/148‑2 (MQ=255)                                           gCAGCGCCATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGG‑GTACAAATAGCGCCGAGGATCATAAACGTTGCTGTGCCGTGGTTGAAGC‑a                                                                                                            <  1:636663/151‑1 (MQ=255)                                            cAGCGCCATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGG‑GTACAAATAGCGCCGAGGATCATAAACGTTGCTGTGCCGTGGTTGAAGC‑At                                                                                                           <  2:883899/151‑1 (MQ=255)                                            cAGCGCCATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGG‑GGACAAAAAACGCCGGGGGGCATAAAAGGTGCTGTGCCGGGGTTGAAg                                                                                                               >  1:562834/1‑148 (MQ=255)                                            cAGCGCCATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGG‑GGACAAAAAACGCCCCGGAGCATAAACAATTCTGTGCCGCgggt                                                                                                                   >  2:236675/1‑142 (MQ=255)                                                 ccATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGG‑GGACAAAAAACGCCCCGGGGCATAAAACATGCTGTGCCGTGGGTGGAGC‑aattc                                                                                                        >  2:628515/1‑145 (MQ=255)                                                  cATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGG‑GTACAAATAGCGCCGAGGATCATAAACGTTGCTGTGCCGTGGTTGAAGC‑ATTCCGcc                                                                                                     <  1:583488/151‑1 (MQ=255)                                                     atCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGG‑GTACAAATAGCGCCGAGGATCATAAACGTTGCTGTGCCGTGGTTGAAGC‑ATTCCGCCAgg                                                                                                  <  2:529100/151‑1 (MQ=255)                                                      tCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGG‑GGACAAAAAACGCCGAGGATCATAAAAATTGCTGTGGCGCGGGTGGAGC‑ATTCCGcc                                                                                                     >  1:496688/1‑147 (MQ=255)                                                      tCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGG‑GGACAAAAAACGCCCAGGATCCTAAACACTGCTGTGCCGCGGGTGGAGC‑ATTCCCCCCggg                                                                                                 >  2:527175/1‑150 (MQ=255)                                                           aaTGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGG‑GGACAAAAAACGCCCGGGGACATAAAAGATGGTGTGCCGGGGGTGGAGG‑AATTCCCCCGGGTAggg                                                                                            >  2:384439/1‑151 (MQ=255)                                                                tCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGG‑GTACAAAAAGAGCCGAGGGTCATAAAAGTTGCTGTGGCGCGGTTGTAGC‑AGTCCGCCCGGTTAGGGCagg                                                                                        >  1:205354/1‑150 (MQ=255)                                                                   gcgcCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGG‑GTACAAATAGCGCCGAGGATCATAAACGTTGCTGTGCCGTGGTTGAAGC‑ATTCCGCCAGGTTAGGGCagga                                                                                       <  2:146151/148‑1 (MQ=255)                                                                       cAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGG‑GGACAAATAACGCCGGGGGGCATAAAAAATGCTGTGCCGGGGGTGGAGC‑ATTCCGCCCGGGTAGGGGAGGGGGGTTcc                                                                                >  2:56140/1‑151 (MQ=255)                                                                         gTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGG‑GGACAAAAAGCGCCGAGGGACATAAAAGGTGCTGTGCCGGGGGTGGAGC‑ATTCCCCCCGGGTTGGGGaggggg                                                                                     >  1:233708/1‑144 (MQ=255)                                                                           ttCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGG‑GGAAAAAAAGAGCCGCGGAACATAAAAGTTGCTGTGGCGGGGGTGGAGG‑ATACCCCCCGGGTTGGGGAGGGGGGGTCCTccc                                                                            >  1:420273/1‑149 (MQ=255)                                                                             cACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGG‑GTACAAATAGCGCCGAGGATCATAAACGTTGCTGTGCCGTGGTTGAAGC‑ATTCCGCCAGGTTAGGGCAGGAGGCTTCCTCGca                                                                           <  2:887297/150‑1 (MQ=255)                                                                             cACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGG‑GGACAAAAAACGCCCAGGGTCATAAAAGTTGCTGTGGCGGGGGTGAAGC‑ATTTCCCCCGGGTAGGGGagggggggtcc                                                                                >  1:907823/1‑140 (MQ=255)                                                                                  gTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGG‑GTACAAATAGCGCCGAGGATCATAAACGTTGCTGTGCCGTGGTTGAAGC‑ATTCCGCCAGGTTAGGGCAGGAGGCTTCCTCGcagacaga                                                                     <  2:599162/151‑1 (MQ=255)                                                                                       ttGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGG‑GTACAAATAGCGCCGAGGATCATAAACGTTGCTGTGCCGTGGTTGAAGC‑ATTCCGCCAGGTTAGGGCAGGAGGCTTCCTCGCAGACAGAATGc                                                                 <  1:110400/150‑1 (MQ=255)                                                                                         ggCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGG‑GGACAAAAAACGACGAGGAACATAAACCTTGCTGGGGCGGGGGTGAAAC‑AATCCGCCAGGTTAGGGGAGGGGGGCTCCTCGCAGACAGAAAGCAg                                                               >  2:868441/1‑150 (MQ=255)                                                                                          gCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGG‑GGAAAAAAAACGCCGAGGGTGATAAAACATGCTGTGCCGCGGGTGTAGCAATTCCCCCCGGTT‑GGGGAGGGGGGTTCCTCCCCGACAGAAAGCAg                                                               >  2:514817/1‑149 (MQ=255)                                                                                           cATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGG‑GTACAAATAGCGCCGAGGATCATAAACGTTGCTGTGCCGTGGTTGAAGC‑ATTCCGCCAGGTTAGGGCAGGAGGCTTCCTCGCAGACAGAATGCAGGcc                                                            <  1:427876/151‑1 (MQ=255)                                                                                           cATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGG‑GTACAAATAGCGCCGAGGATCATAAACGTTGCTGTGCCGTGGTTGAAGC‑ATTCCGCCAGGTTAGGGCAGGAGGCTTCCTCGCAGACAGAATGCAGGcc                                                            <  2:60822/151‑1 (MQ=255)                                                                                             tCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGG‑GTACAAAAAAAGCCGAGGGTCATAAACGGTGCTGTGCCGTGGGTGAAGC‑ATTCCCCCCGGGTAGGGGAGGGGGGTTCCTCGCCGACAGAAAGCAgg                                                              >  1:939758/1‑147 (MQ=255)                                                                                               aGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGG‑GTACAAAAAACGCCGGGGGGCATAAAACTTGCTTTGCCGGGGGTGGAGA‑AATTCGCCCGGGTAGGGGGGGGGGGTTCCTCCCCGGCAGAAAGCAGGCCCttt                                                        >  2:481863/1‑151 (MQ=255)                                                                                                     tACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGG‑GTACAAAAAACGCCCGGGATCATAAACGTTGCTGTGGCGGGGGTGAAAC‑ATTTCCCCCGGGTTGGGGGGGAGGGTTCCTCGCCGGCAGAATGCAGGCCACTTTTtcg                                                   >  1:28064/1‑150 (MQ=255)                                                                                                     tACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGG‑GGACAAAAAACGCCGAGGGTCATAAAAATTTCTGTGCCGGGGGTGGAGG‑ATTTCGCCCGGGTAGGGGAGGGGGGGTCCTCGCCGGCAGAAAGCAGGCCAtttttttg                                                   >  1:856151/1‑147 (MQ=255)                                                                                                      aCCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGG‑GTACAAATAGCGCCGAGGATCATAAACGTTGCTGTGCCGTGGTTGAAGC‑ATTCCGCCAGGTTAGGGCAGGAGGCTTCCTCGCAGACAGAATGCAGGCCATTTTTgcgc                                                  <  2:5188/150‑1 (MQ=255)                                                                                                      aCCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGG‑GTACAAATAGCGCCGAGGATCATAAACGTTGCTGTGCCGTGGTTGAAGC‑ATTCCGCCAGGTTAGGGCAGGAGGCTTCCTCGCAGACAGAATGCAGGCCATTTTTgcg                                                   <  1:823359/149‑1 (MQ=255)                                                                                                      aCCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGG‑GTACAAATAGCGCCGAGGATCATAAACGTTGCTGTGCCGTGGTTGAAGC‑ATTCCGCCAGGTTAGGGCAGGAGGCTTCCTCGCAGACAGAATGCAGGCCATTTTTgcg                                                   <  1:832154/149‑1 (MQ=255)                                                                                                       ccGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGG‑GGACAAAAAACGCCGCGGGACATAAAAATTGCTGTGCCGGGGGTGGAGC‑ATTCCGCCAGGGTAGGGCAGGGGGGCTCCTCCCAGACAGAAAGCAGGCCCttttt                                                      >  2:298713/1‑145 (MQ=255)                                                                                                        cGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGCGTA‑AAATAGCGCCGAGGATCATAAAGGAATGTTTGCCGGGGGTGAAAC‑ATTTCGCCACGGTCGGGCGGGTGGGGTCCTCGCAGGCAGAACGGAGGCCAATTTTGTGCCtt                                               >  1:661344/1‑151 (MQ=255)                                                                                                        cGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGG‑GGACAAAAAACGCCGAGGGACATAAAAGTTGCTGTGCCGTGGGTGAAGC‑ATTCCGCCAGGGTAGGGGAGGAGGGGTCCTCGCCGACAGAAAGCAGGCCATTTTTTTGCAtt                                               >  1:884668/1‑151 (MQ=255)                                                                                                        cGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGG‑GGACAAAAAACGCCCAGGGGAATAAAACATGCTGTGCCGTGGGTGGAGC‑AATCCGCCCGGGTAGGGCAGGGGGGGTCCTCCCAGACAGAAAGCAGGGCACTTTTGTGCAtt                                               >  2:965927/1‑151 (MQ=255)                                                                                                         gggCGACCGTGGGCAACATCACAGAACGGACAACGGCGG‑GGACAAAAAACGCCCAGGGACATAAAAATTGGTGTGCCGGGGGTGGAGCAATTCCCCCCGG‑TAGGGCGGGAGGGCTCCTCCCCGACAGAAAGGCGGGCCTTTTTtcgcct                                                >  2:855265/1‑147 (MQ=255)                                                                                                            cGACCGTGGGCAACATCACAGAACGGACAACGGCGG‑GGACAAAAAACGCCCGGGAACATAAAAGATGCTGGGCCGGGGGTGGAGC‑ATTTCGCCCGGGTAGGGGAGGGGGGTTTCTCGCCGGCAGAAAGAAGGCCCTTTTTTTGCATTGCgg                                           >  2:333266/1‑151 (MQ=255)                                                                                                             gACCGTGGGCAACATCACAGAACGGACAACGGCGG‑GTACAAAAAACGCCGAGGGTCATAAACATTGCTGTGCCGGGGGTTAAGC‑ATTCCGCCCGGGTTGGGGAGGGGGGTTTCTCCCCGACAGAAAGCAGGGCATTTTTTCGCCTTGCgg                                           >  1:736946/1‑150 (MQ=255)                                                                                                                 gTGGGCAACATCACAGAACGGACAACGGCGG‑GGACAAAAAACGCCCAGGAACATAAACGTTGCTGTGGCGTGGGTGGAGC‑ATTCCGCCCGGGTTGGGGAGGGGGGTTCCCCGCAGACAGAAAGCAGGCCCTTTTTTCGCATTTCGGGttt                                       >  1:770470/1‑150 (MQ=255)                                                                                                                  tGGGCAACATCACAGAACGGACAACGGCGG‑GGACAAAAAACGCCGGGGAACATAAAAGTTGCTGTGCCGTGGGTGGAGG‑ATTCCGCCCGGGTTGGGCAGGGGGGTTCCTCCCAGACAGAAAGCAGGGCCTTTTTTCGCATTGCGGCttt                                       >  1:41443/1‑149 (MQ=255)                                                                                                                                           ggcgg‑GTACAAATAGCGCCGAGGATCATAAACGTTGCTGTGCCGTGGTTGAAGC‑ATTCCGCCAGGTTAGGGCAGGAGGCTTCCTCGCAGACAGAATGCAGGCCATTTTTGCGCATTGCGGCTTTGATGCCCTGGATACGTGTAGAGTCCGc            >  2:192779/1‑151 (MQ=255)                                                                                                                                              gg‑gTACAAATAGCGCCGAGGATCATAAACGTTGCTGTGCCGTGGTTGAAGC‑ATTCCGCCAGGTTAGGGCAGGAGGCTTCCTCGCAGACAGAATGCAGGCCATTTTTGCGCATTGCGGCTTTGATGCCCTGGATACGTGTAGAGTCCGCTgg         >  1:686086/1‑151 (MQ=255)                                                                                                                                                  tACAAATAGCGCCGAGGATCATAAACGTTGCTGTGCCGTGGTTGAAGC‑ATTCCGCCAGGTTAGGGCAGGAGGCTTCCTCGCAGACAGAATGCAGGCCATTTTTGCGCATTGCGGCTTTGATGCCCTGGATACGTGTAGAGTCCGCTGGa        >  2:904400/1‑149 (MQ=255)                                                                                                                                                   aCAAATAGCGCCGAGGATCATAAACGTTGCTGTGCCGTGGTTGAAGC‑ATTCCGCCAGGTTAGGGCAGGAGGCTTCCTCGCAGACAGAATGCAGGCCATTTTTGCGCATTGCGGCTTTGATGCCCTGGATACGTGTAGAGTCCGCTGGAAGc     >  1:731459/1‑151 (MQ=255)                                                                                                                                                     aaaTAGCGCCGAGGATCATAAACGTTGCTGTGCCGTGGTTGAAGC‑ATTCCGCCAGGTTAGGGCAGGAGGCTTCCTCGCAGACAGAATGCAGGCCATTTTTGCGCATTGCGGCTTTGATGCCCTGGATACGTGTAGAGTCCGCTGGAAGCtt   >  2:428032/1‑151 (MQ=255)                                                                                                                                                      aaTAGCGCCGAGGATCATAAACGTTGCTGTGCCGTGGTTGAAGC‑ATTCCGCCAGGTTAGGGCAGGAGGCTTCCTCGCAGACAGAATGCAGGCCATTTTTGCGCATTGCGGCTTTGATGCCCTGGATACGTGTAGAGTCCGCTGGAAGCTTg  >  1:394873/1‑151 (MQ=255)                                                                                                                                                        |                                                                                                                                                      CGCGCAGGTCATCACGGTCAACGGAGGTGATAACCACATAACGCAGCGCCATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGG‑GTACAAATAGCGCCGAGGATCATAAACGTTGCTGTGCCGTGGTTGAAGC‑ATTCCGCCAGGTTAGGGCAGGAGGCTTCCTCGCAGACAGAATGCAGGCCATTTTTGCGCATTGCGGCTTTGATGCCCTGGATACGTGTAGAGTCCGCTGGAAGCTTG  >  NC_000913/659790‑660088

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