Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,409,402	Δ1 bp	11.7%	coding (128/732 nt)	rlmB →	23S rRNA mG2251 2'‑O‑ribose methyltransferase, SAM‑dependent

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,409,402	0	T	.	11.7%	192.9 / 19.8	60	coding (128/732 nt)	rlmB	23S rRNA mG2251 2'‑O‑ribose methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base T (30/23);  new base . (4/3);  total (34/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.32e-01

ACGTAAACAACGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCG  >  NC_000913/4409254‑4409547                                                                                                                                                     |                                                                                                                                                   aCGTAAACAACGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCt                                                                                                                                                    >  1:282430/1‑149 (MQ=255)    tAAACAACGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTg                                                                                                                                                  >  4:20192/1‑148 (MQ=255)     aaaCAACGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAg                                                                                                                                                >  2:116933/1‑149 (MQ=255)      aacaacGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGt                                                                                                                                               <  1:365758/149‑1 (MQ=255)        caacGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGtct                                                                                                                                             >  1:320906/1‑149 (MQ=255)          acGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑a                                                                                                                                          <  2:282430/149‑1 (MQ=255)          acGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑a                                                                                                                                          <  1:51321/149‑1 (MQ=255)              gTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGc                                                                                                                                      >  2:20383/1‑149 (MQ=255)                aCATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCC‑TGAGACTCAAGGGCGt                                                                                                                                    <  2:77151/149‑1 (MQ=255)                 cATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTg                                                                                                                                   <  1:20383/149‑1 (MQ=255)                 cATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCC‑TGAGACTCAAGGGCGt                                                                                                                                    <  3:185914/148‑1 (MQ=255)                       tGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGt                                                                                                                                    >  4:319356/1‑142 (MQ=255)                       tGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGt                                                                                                                                    <  3:319356/142‑1 (MQ=255)                            gAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGtt                                                                                                                        >  1:352846/1‑149 (MQ=255)                            gAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGtt                                                                                                                        <  1:172481/149‑1 (MQ=255)                                tGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCa                                                                                                                    <  1:187549/149‑1 (MQ=255)                                   tttACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCaaa                                                                                                                  >  4:68390/1‑148 (MQ=255)                                        ggCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCa                                                                                                            <  1:337282/149‑1 (MQ=255)                                         gCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCaa                                                                                                           <  1:310817/149‑1 (MQ=255)                                            tCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAATCCGCCAata                                                                                                         >  3:73439/1‑148 (MQ=255)                                               aCGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTgg                                                                                                                                  <  1:9493/120‑1 (MQ=255)                                               aCGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTgg                                                                                                                                  >  2:9493/1‑120 (MQ=255)                                                cGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATAtctc                                                                                                     >  3:219448/1‑148 (MQ=255)                                                  cAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTcga                                                                                                   <  4:219448/148‑1 (MQ=255)                                                    gTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTcgacg                                                                                                 <  4:30146/148‑1 (MQ=255)                                                        aGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGaaa                                                                                            >  1:278695/1‑149 (MQ=255)                                                             ctgctgGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAggg                                                                                        >  4:281583/1‑146 (MQ=255)                                                                  ggAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGt                                                                                   >  3:265844/1‑148 (MQ=255)                                                                     gcgcgcCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCg                                                                               >  2:381955/1‑149 (MQ=255)                                                                      cgcgcCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCg                                                                               <  3:167585/148‑1 (MQ=255)                                                                            ccTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCGGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATc                                                                         >  3:126640/1‑148 (MQ=255)                                                                             cTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAg                                                                       >  2:107399/1‑149 (MQ=255)                                                                              tGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAgg                                                                      <  2:17577/149‑1 (MQ=255)                                                                                  cGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCa                                                                   >  4:325831/1‑148 (MQ=255)                                                                                    ttttCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGcc                                                                                >  2:94374/1‑133 (MQ=255)                                                                                    ttttCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGcc                                                                                <  1:94374/133‑1 (MQ=255)                                                                                                cTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTg                                                     >  3:33372/1‑148 (MQ=255)                                                                                                 tttATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGaa                                                   <  1:107399/149‑1 (MQ=255)                                                                                                      tttAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCa                                               >  4:207414/1‑148 (MQ=255)                                                                                                        tAAAAGGCCGTGAAGATAAACGTCTGTTACCTCTGATTCACGCCCTTGAGTCTC‑ATGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTTCATCAGGGCATTATCGCCCGCGTGAAGGCAGGa                                            <  1:116933/149‑1 (MQ=255)                                                                                                             ggCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTc                                        <  4:265844/148‑1 (MQ=255)                                                                                                             ggCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCa                                       >  2:137912/1‑149 (MQ=255)                                                                                                                  tGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATc                                  >  1:392713/1‑149 (MQ=255)                                                                                                                    aaGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCa                                 <  4:54151/148‑1 (MQ=255)                                                                                                                             gTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGAt                       >  1:246886/1‑149 (MQ=255)                                                                                                                                 gTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGc                   >  2:68421/1‑149 (MQ=255)                                                                                                                                     ccGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGa               <  2:246886/149‑1 (MQ=255)                                                                                                                                      cGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGAt              >  1:96892/1‑149 (MQ=255)                                                                                                                                      cGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGAt              <  2:320906/149‑1 (MQ=255)                                                                                                                                        cTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATc             >  3:156011/1‑148 (MQ=255)                                                                                                                                           aTTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGATAACGATCTGCCGGATCTGa          <  4:73439/148‑1 (MQ=255)                                                                                                                                           aTTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGa          >  4:18539/1‑148 (MQ=255)                                                                                                                                              cACGCCC‑TGAGACTCAAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGAtcgc      >  2:240710/1‑149 (MQ=255)                                                                                                                                              cACGCCC‑TGAGACTCAAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGAtcg       >  4:189033/1‑148 (MQ=255)                                                                                                                                               aCGCCC‑TGAGACTCAAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGGAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGt                                     >  4:82728/1‑117 (MQ=255)                                                                                                                                               aCGCCC‑TGAGACTCAAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGGAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGt                                     <  3:82728/117‑1 (MQ=255)                                                                                                                                               aCGCCC‑TGAGACTCAAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGAtcgct     >  2:77411/1‑149 (MQ=255)                                                                                                                                                  cccTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGa               >  4:380928/1‑136 (MQ=255)                                                                                                                                                  cccTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGa               <  3:380928/136‑1 (MQ=255)                                                                                                                                                  cccTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGAtcgcttc   >  3:345098/1‑148 (MQ=255)                                                                                                                                                   ccTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGAtcgcttcg  >  4:22685/1‑148 (MQ=255)                                                                                                                                                     |                                                                                                                                                   ACGTAAACAACGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCG  >  NC_000913/4409254‑4409547

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