Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,410,395	Δ1 bp	69.7%	coding (885/1374 nt)	dbpA →	ATP‑dependent RNA helicase, specific for 23S rRNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,410,395	0	G	.	69.7%	106.6 / 91.1	66	coding (885/1374 nt)	dbpA	ATP‑dependent RNA helicase, specific for 23S rRNA
Reads supporting (aligned to +/- strand):  ref base G (8/11);  new base . (24/22);  total (33/33)
Fisher's exact test for biased strand distribution p-value = 5.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.43e-01

TTTTGCAATACCAAAAAAGATTGCCAGGCTGTCTGCGACGCGCTGAATGAAGTAGGGCAAAGTGCATTGTCATTACACGGCGATTTGGAGCAACGCGATCGCGATCAGACCCTGGTACGTTTTGCTAACGGTAGCGCCCGTGTACTGGTCGCGACTGATGTTGCTGCGCGTGGTCTGGATATTAAATCGCTTGAGCTGGTGGTGAACTTTGAGCTGGCGTGGGACCCTGAAGTTCATGTACATCGCATCGGTCGTACAGCTCGTGCAGGAAATAGCGGTCTGGCG  >  NC_000913/1410249‑1410533                                                                                                                                                   |                                                                                                                                           ttttGCAATACCAAAAAAGATTGCCAGGCTGTCTGCGACGCGCTGAATGAAGTAGGGCAAAGTGGATTGTCATTACACGGCGATTTGGAGCAACGCGATCGCGATCAGACCCTGGTACGTTTTGCTAACGGTAGCGGCCGTGTACTgtc                                                                                                                                          >  2:86091/1‑147 (MQ=255)  tttGCAATACGAACAAAGATTGTCAGGCTGTCTGCGACGCGCTGAATGAAGTAGGGCAAAGTGCATTGTCATTACACGGCGATTTGGAGCAACGCGTTCGCGATCAGACCCTGGTACGTTTTGCTAACGGTAGCGCCCGTGTACTgtcg                                                                                                                                         <  2:591376/149‑4 (MQ=255)      cAATACCAAAAAAGATTGCCAGGCTGTCTGCGACGCGCTGAATGAAGTAGGGCAAAGTGCATTGTCATTACACGGCGATTTGGAGCAACGCGATCGCGATCAGACCCTGGTACGTTTTGCTAACGGTAGCGCCCGTGTACTGGTCGCGa                                                                                                                                     <  1:534682/149‑1 (MQ=255)       aaTACCAAAAAAGATTGCCAGGCTGTCTGCGACGCGCTGAATGAAGTAGGGCAAAGTGCATTGTCATTACACGGCGATTTGGAGCAACGCGATCGCGATCAGACCCTGGTACGTTTTGCTAACGGTAGCGACCGTGTCCT‑GTCGCGACt                                                                                                                                   >  1:332209/1‑149 (MQ=255)       aaTACCAAAAAAGATTGCCAGGCTGTCTGCGACGCGCTGAATGAAGTAGGGCAAAGTGCATTGTCATTACAAGGCGATTTGGAGCCACGCGAAAGCGATCAGACCCTGGTACGTTTTGCTAACGGAAGCGACCGGGTACT‑GTTGCGACt                                                                                                                                   >  2:204533/1‑149 (MQ=255)        aTACCAAAAAAGATTGCCAGGCTGTCTGCGACGCGCTGAATGAAGTAGGGCAAAGTGCATTGTCATTACACGGCGATTTGGAGCAACGCGATCGCGATCAGACCCTGGTACGTTTTGCTAACGGTAGCGCCCGTGTACT‑GTCGCGACTg                                                                                                                                  >  2:503052/1‑149 (MQ=255)            cAAAAACGATTGCCAGGCTGGCTGCGACGCGCTGAAGGAAGCAGGGCAAAGCGCATTGTCATTACACGGCGATCTGGAGCAACGCGATCGCGATCAGACCCTGGTACGTTTTGCTAACGGTAGCGCCCGTGTACT‑GTCGCGACTGATGt                                                                                                                              <  2:599146/149‑1 (MQ=255)            cAAAAAAGATTGCCAGGCTGTCTGCGACGCGCTGAATGAAGTAGGGCAAAGTGCATTGTCATTACACGGCGATTTGGAGCAACGCGATCGCGATCAGACCCTGGTACGTTTTGCTAACGGTAGCGCCCGTGTACT‑GTCGCGACTGATGt                                                                                                                              >  2:595660/1‑149 (MQ=255)              aaaaaGATTGCCAGGCTGTATGCGACGCGCTGAATGAAGTGGGGCAAAGTGCATTGTCATTACACGGCGATTTGGAGCAACGCGATCGCGATCAGACCCTGGTACGTTTTGCTAACGATAGCGCCCGGGTACT‑GTCGCCACTGATGTtg                                                                                                                            >  1:56464/1‑149 (MQ=255)                 aaGATTGCCAGGCTGTCTGCGACGCGCTGAATGAAGTAGGGCAAAGTGCATTGTCATTACACGGCGATTTGGAGCAACGCGATCGCGATCAGACCCTGGTACGTTTTGCTAACGGTAGCGCCCGTGTACT‑GTCGCGACGGATGTtgccg                                                                                                                         >  2:114827/1‑147 (MQ=255)                       gCCAGGCTGTCTGCGACGCGCTGAATGAAGTAGGGCAAAGTGCATTGTCATTACACGGCGATTTGGAGCAACGCGATCGCGATCAGACCCTGGTACGTTTTGCTAACGGTAGCGCCCGTGTACTGGTCGCGACTGATGTTGCTGCGCGt                                                                                                                    >  1:362169/1‑149 (MQ=255)                        ccAGGCTGTCTGCGACGCGCTGAATGAAGTAGGGCAAAGTGCATTGTCATTACACGGCGATTTGGAGCAACGCGATCGCGATCAGACCCTGGTACGTTTTGCTAACGGTAGCGCCCGTGTACT‑GTCGCGACTGATGTTGCTGCGCGTgg                                                                                                                  >  1:638203/1‑149 (MQ=255)                         cAGGCTGTCTGCGACGCGCTGAATGAAGTAGGGCAAAGTGCATTGTCATTACACGGCGATTTGGAGCAACGCGATCGCGATCAGACCCTGGTACGTTTTGCTAACGGTAGCGCCCGTGTACTGGTCGCGACTGATGTTGCTGCGCGTgg                                                                                                                  >  1:46729/1‑149 (MQ=255)                           gggtGTCTGCTACGCGCTGAATGAAGGAGGGCAAAGTGCATTGTCATTACACGGCGCGTTGGAGCAAGGCGAGCGCGATCAGACCCTGGTACGTTTTGCTAACGGTAGCGCCCGTGTACT‑GTCGCGACTGATGTTGCTGCGCGTGGTCt                                                                                                               <  2:680293/146‑1 (MQ=255)                           ggCTGTCTGCGACGCGCTGAATGAAGTAGGGCAAAGTGCATTGTCATTACACGGCGATTTGGAGCAACGCGATCGCGATCAGACCCTGGTACGTTTTGCTAACGGTAGCGCCCGTGTACT‑GTCGCGACTGATGTTGCTGCGCGTGGTCt                                                                                                               >  1:162918/1‑149 (MQ=255)                            gCTGTCTGCGACGCGCTGAATGAAGTAGGGCAAAGTGCATTGTCATTACACGGCGATTTGGAGCAACGCGATCGCGATCAGACCCTGGTACGTTTTGCTAACGGTAGCGCCCGTGTACT‑GTCGCGACTGATGTTGCTGCGCGTGGTCTg                                                                                                              <  1:62451/149‑1 (MQ=255)                                tCTGCGACGCGCTGAATGAAGTAGGGCAAAGTGCATTGTCATTACACGGCGATTTGGAGCAACGCGATCGCGATCAGACCCTGGTACGTTTTGCTAACGGTAGCGCCCGTGTACT‑GTCGCGACTGATGTTGCTGCGCGTGGTCTGGata                                                                                                          <  2:227345/149‑1 (MQ=255)                                  tGCGACGCGCTGAATGAAGTAGGGCAAAGTGCATTGTCATTACACGGCGATTTGGAGCAACGCGATCGCGATCAGACCCTTGGACGTTTTGCTAACGGTAGCGCCCGTGTACT‑GTCGCGACTGATGTTGCTGCGCGTGGTCTGGATAtt                                                                                                        <  2:405078/149‑1 (MQ=255)                                  tGCGACGCGCTGAATGAAGTAGGGCAAAGTGCATTGTCATTACACGGCGATTTGGAGCAACGCGATCGCGATCAGACCCTGGTACGTTTTGCTAACGGTAGCGCCCGTGTACTGGTCGCGACTGATGTTGCTGCGCGTGGTCTGGatat                                                                                                         >  1:394277/1‑149 (MQ=255)                                  tGCGACGCGCTGAATGAAGTAGGGCAAAGTGAATTGTCATTACACGGCGATTTGGAGCAACGCAATCGCGATCAGACCCGGGTACGTTTTGCTAACGGTAGCGCCCGGGTACT‑GTCGCGAATGGAGTTGCTCCGCGTGGGCTGGATAtt                                                                                                        >  2:348404/1‑149 (MQ=255)                                   gcgcCGCGCTGAATGAAGTAGGGCAAAGTGCATTGTCATTACACGGCGATTTGGAGCAACGCGATCGCGATCAGACCCTGGTACGTTTTGCTAACGGTAGCGCCCGTGTACT‑GTCGCGACTGATGTTGCTGCGCGTGGTCTGGATATTa                                                                                                       <  1:418494/149‑1 (MQ=255)                                    cGACGCGCTGAATGAAGTAGGGCAAAGTGCATTGTCATTACACGGCGATTTGGAGCAACGCGATCGCGATCAGACCCTGGTACGTTTTGCTAACGGTAGCGCCCGTGTACT‑GTCGCGACTGATGTTGCTGCGCGTGGTCTGGGTGTTaa                                                                                                      >  2:371204/1‑149 (MQ=255)                                    cGACGCGCTGAATGAAGTAGGGCAAAGTGCATTGTCATTACACGGCGATTTGGAGCAACGCGATCGCGATCAGACCCTGGTACGTTTTGCTAACGGTAGCGCCCGTGTACT‑GTCGCGACTGATGTTGCTGCGCGTGGTCTGGATATTaa                                                                                                      >  1:675510/1‑149 (MQ=255)                                     gACGCGCTGAATGAAGTAGGGCAAAGTGCATTGTCATTACACGGCGATTTGGAGCAACGCGATCGCGATCAGACCCTGGTACGTTTTGCTAACGGTAGCGCCCGTGTACT‑ATCGCGACTGATGTTGCTGCGCGTGGTCTGGATATTaaa                                                                                                     <  2:116573/149‑1 (MQ=255)                                      aCGCGCTGACTTAAGTAGGGCCAAGTGCATTGTCATTCCACGGCGATTTGGAGCAACGCGATCGCGATCAGACCCTGGTACGTTTTGCTAACGGGAGCGCCCGTGTACTGGTTGCGACTGATGTTGCTGCGCGTGGTCTGGATATTaaa                                                                                                     <  1:32962/149‑1 (MQ=255)                                        gcgcTGAATGAAGTAGGGCAAAGTGCATTGTCATTACACGGCGATTTGGAGCAACGCGATCGCGATCAGACCCTGGTACGTTTTGCTAACGGTAGCGCCCGTGTACT‑GTCGCGACTGATGTTGCTGCGCGCGGTCTGGATAATAAATCg                                                                                                  >  2:50871/1‑149 (MQ=255)                                           cTGAAGGAAGTAGGGCAAAGTGCAGTGTCATTACACGGCGATTTGGAGCAACGCGATCGCGATCAGACCCTGGTACGTTTTGCTAACGGTAGCGCCCGTGTACT‑GTCGCGACTGATGTTGCTGCGCGTGGTCTGGATATTAAATCGCtt                                                                                               <  2:282718/149‑1 (MQ=255)                                             gaatgaaGTAGGGCAAAGGGCATTGTCATTACCCGGCGATTTGGAGCAACGCGATCGCGATCAGACCCTGGTACGTTTTGCTAACGGTAGCGCCCGTGTACTGGTCGCGACTGATGTTGCTGCGCGTGGTCTGGATATTAAATCGCTTg                                                                                              <  1:319742/149‑1 (MQ=255)                                              aatgaaGTAGGGCAAAGTGCATTGTCATTACACGGCGATTTGGAGCAACGCGATCGCGATCAGACCCTGGTACGTTTTGCTAACGGTAGCGCCCGTGTACT‑GTCGCGACTGATGTTGCTGCGCGTGGTCTGGATATTAAATCGCTTGAg                                                                                            <  2:712185/149‑1 (MQ=255)                                              aatgaaGTAGGGCAAAGTGCATTGTCATTACACGGCGATTTGGAGCAACGCGATCGCGATCAGACCCTGGTACGTTTTGCTAACGGTAGCGCCCGTGTACT‑GTCGCGACTGATGTTGCTGCGCGTGGTCTGGATATTAAATCGCTTGAg                                                                                            <  2:712194/149‑1 (MQ=255)                                                tgaaGCAGGGCAAAGTGCATTGTCATTACACGGCGATTTGGAGCAACGCGATCGCGATCAGACCCTGGTACGTTTTGCTAACGGTAGCGCCCGTGTACT‑GTCGCGACTGATGTTGCTGCGCGTGGTCTGGATATTAAATCGCTTGAGCt                                                                                          <  2:695682/149‑1 (MQ=255)                                                     tAGGGCAAAGTGCATTGTCATTACACGACGATTTGGACCAACGCGCTCGGGATCAGACCCTGGTACGTTTTGCTAACGGTAGCGCCCGTGTACT‑GTCGCGACTGATGTTGCTGCGCGTGGTCTGGATATTAAATCGCTTGAGCtggtgg                                                                                     <  2:335029/149‑1 (MQ=255)                                                        ggCAAAGTGCATTGTCATTACACGGCGATTTGGAGCAACGCGATCGCGATCAGACCCTGGTACGTTTTGCTAACGGTAGCGCCCGTGTACT‑GTCGCGACTGATGTTGCTGCGCGTGGTCTGGATATTAAATCGCTTGAGCTGGTGGTGa                                                                                  >  1:249433/1‑149 (MQ=255)                                                          cAAAGTGCATTGTCATTACACGGCGATTTGGAGCAACGCGATCGCGATCAGACCCTGGTACGTTTTGCTAACGGTAGCGCCCGTGTACT‑GTCGCGACTGATGTTGCTGCGCGTGGTCTGGATATTAAATCGCTTGAGCTGGTGGTGAAc                                                                                >  2:731065/1‑149 (MQ=255)                                                            aaGCGCATTGTCATTGCACGGGGAGTTGGGGCAACGCGATCGCGATTAGACCCTCGTACGTTTTGATAACAGTAGCCTCCGTGTAC‑GGCTGTGACTGATGTTGCTGCGCGTGGTCTGGATATTAAACCGCTTGAGCTGGTGGTGAACtt                                                                              <  2:185032/149‑1 (MQ=255)                                                              gTGCATTGTCATTACACGGCGATTTGGAGCAACGCGATCGCGATCAGACCCTGGTACGTTTTGCTAACGGTAGCGCCCGTGTACT‑GTCGCGACTGATGTTGCTGCGCGTGGTCTGGATATTAAATCGCTTGAGCTGGTGGTGAACTTTg                                                                            <  1:704921/149‑1 (MQ=255)                                                                 cATTGTCATTACACTGCGATTTGGAGCAACGCGATCGCGATCAGACCCTGGGACGTTTTGCTAACGGTAGCGCCCGTGTACTGGTCGCGACTGATGTTGCTGCGCGTGGTCTGGATATTAAATCGCTTGAGCTGGTGGTGAACTTTGAg                                                                          <  1:609070/149‑1 (MQ=255)                                                                 cATTGTCATTACACGGCGATTTGGAGCAACGCGATCGCGATCAGACCCTGGTACGTTTTGCTAACGGTAGCGCCCGTGTACTGGTCGCGACTGATGTTGCTGCGCGTGGTCTGGATATTAAATCGCTTGAGCTGGTGGTGGACTTTGAg                                                                          >  1:456878/1‑149 (MQ=255)                                                                   ttGTCATTACACGGCGATTTGGAGCAACGCGATCGCGATCAGACCCTGGTACGTTTTGCTAACGGTAGCGCCCGTGTACT‑GTCGCGACTGATGTTGCTGCGCGTGGTCTGGATATTAAATCGCTTGAGCTGGTGGTGAACTTTGAGCTg                                                                       >  1:172693/1‑149 (MQ=255)                                                                       cATTACCCGGCGAGTTGGAGCCACGCGATCGCGATCACACCCTGGTACGTTTTGCTAACGGTAGCGCCCGTGTACT‑GTCGCGACTGATGTTGCTGCGCGTGGTCTGGATATTAAATCGCTTGAGCTGGTGGTGAACTTTGAGCTGGCGt                                                                   <  1:593680/149‑1 (MQ=255)                                                                        aTTACGTGGAGATTTGGAGCAAGGCGATCGGGATCAGACCCTGGTACTTTTTGCTAACGGTAGCGCCCGTGTACT‑GTCGCGACTGATGTTGCTGCGCGTGGTCTGGATATTAAATCGCTTGAGCTGGTGGTGAACTTTGAGCTGGCGTg                                                                  <  2:688278/149‑1 (MQ=255)                                                                         ttACACGGCGATTTGGAGCAACGCGATCGCGATCAGACCCTGGTACGTTTTGCTAACGGTAGCGCCCGTGTACT‑GTCGCGACTGATGTTGCTGCGCGTGGTCTGGATATTAAATCGCTTGAGCTGGTGGTGAACTTTGAGCTGGCGTgg                                                                 >  1:721741/1‑149 (MQ=255)                                                                         ttACACGGCGATTTGGAGCAACGCGATCGACATCAGGCCCTGGTACGTTTTGCTAACGGTAGCGCCCGTGTACTGGTCGCGACTGATGTTGCTGCGCGTGGTCTGGATATTAAATCGCTTGAGCTGGTGGTGAACTTTGAGCTGGCGTg                                                                  <  1:136620/149‑1 (MQ=255)                                                                                   aTTTGGAGCAACGCGATCGCGATCAGCCCCTGGTACGCTTTGCTAACGGTAGCGCCCGTGTACT‑GTCGCGACTGACGTTGCTGCGCGTGGTCTGGATATTAAATCGCTTGAGCTGGTGGTGAACTTTGAGCTGGCGTGGGACCCTGAAg                                                       <  2:204111/149‑1 (MQ=255)                                                                                        gAGCAACGCGATCGCGATCAGCCCCTGGTACGTTTTGCTAACGGTAGCGCCCGTGTACT‑GTCGCGACTGATGTTGCTGCGCGTGGTCTGGATATTAAATCGCTGGAGCTGGTGGTGAACTTTGAGCTGGCGTGGGACCCTGAAGTTCAt                                                  <  2:574486/149‑1 (MQ=255)                                                                                        gAGCAACGCGATCGCGATCAGACCCTGGTACGTTTTGCTAACGGTAGCGCCCGTGTACT‑GTCGCGACTGATGTTGCTGCGCGATGTGAATATTTTAAATCGCTTGTGGTGGTGGTGAACTTTGAGCTCGCGGGTGACCCTGAAGTTCAt                                                  <  2:701762/149‑1 (MQ=255)                                                                                         aGCAACGCGATCGCGATCAGACCCTGGTACGTTTTGCTAACGGTAGCGCCCGTGTACT‑GTCGCGACTGATGTTGCTGCGCGTGGTCTGGATATTAAATCGCTTGAGCTGGTGGTTAACTTTGAGCTGGCGTGGGCACCTGAAGTTCATg                                                 >  2:255804/1‑149 (MQ=255)                                                                                          gCAACGCGATCGCGCTCAGACCCTGGTACGTTTTGCTAATGGTAGCGCCCGTGTACT‑GTCGCGACTGATGTTGCTGCGCGTGGTCTGGATATTAAATCGCTTGAGCTGGTGGTGAACTTTGAGCTGGCGTGGGACCCTGAAGTTCATGt                                                <  2:376413/149‑1 (MQ=255)                                                                                          gCAACGCGATCGCGATCAGACCCTGGTACGTTTTGCTAACGGTAGCGCCCGTGTACTGGTCGCCACTGATGGTGCTGCGCGTGGTCTGGATATTAAATCGCTTGAGCGGGGGGTGAACTGTGAGCTGGCGTGGGACCCTGAAGTTCCTg                                                 >  2:118815/1‑149 (MQ=255)                                                                                           cAACGCGATCGCGATCAGACCCTGGTACGTTTTGCTAACGGTAGCGCCCGTGTACT‑GTCGCGACTGATGTTGCTGCGCGTGGTCTGGATATTAAATCGCTTGAGCTGGTGGTGAACTTTGAGCTGGCGTGGGACCCTGAAGTTCATGTa                                               <  1:315481/149‑1 (MQ=255)                                                                                            aaCGCGATCGCGATCAGACCCTGGTACGTTTTGCTAACGGTAGCGCCAGTGTACTGGTAGCGCCTGATGTTGCTGCGCGTGGTCTGGATATTAAATCGCTTGAGCTGGTGGTGAACTTTGACCTGGCGTGGGAGCCTGAAGTTCATGTa                                               >  1:339004/1‑149 (MQ=255)                                                                                                  atcgcgatcAGACCCTGGTACGTTTTGCTAACGGTAGCGCCCGTGTACTGGTCGCGACTGATGTTGCTGCGCGTGGTCTGGATATTAAATCGCTGGAGCTGGTGGTGAACTTTGAGCTGGCGTGGGGCCCTGAAGGTCCTGTCcatcgc                                         >  1:695284/1‑149 (MQ=255)                                                                                                    cgcgatcAGTCCCTGGTACGTTTTGCTAACGGTACCGCCCGTGTACGTGTCGCGACTGATGTTGCTGCGCGTGGGCTGGACATTGAATCGCTTGAGCTGGTGGCGCACTTTGAGCTGGCTTGGGACCCTGAACTTCATGTAcatcccct                                       >  2:290112/1‑145 (MQ=255)                                                                                                     gcgatcAGACCCTGGTACGTTTTGCTAACGGTAGCGCCCGTGTACT‑GTCGCGACTGATGTTGCTGCGCGTGGTCTGGATATTAAATCGCTTGAGCTGGTGGTGAACTTTGAGCTGGCGTGGGACCCTGAAGTTCATGTAcatcgcatcg                                     >  2:446581/1‑149 (MQ=255)                                                                                                     gcgatcAGACCCTGGTACGTTTTGCTAACGGTAGCGCCCGTGTACT‑GTCGCGACTGATGTTGCTGCGCGTGGTCTGGATATTAAATCGCTTGAGCTGGTGGTGAACTTTGAGCTGGCGTGGGACCCTGAAGTTCATGTAcatcgcatcg                                     >  1:693831/1‑149 (MQ=255)                                                                                                      cgatcAGACCCTGGTACGTTTTGCTAACGGTAGCGCCCGTGTACT‑GTCGCGACTGATGTTGCTGCGCGTGGTCTGGATATTAAATCGCTTGAGCTGGTGGTGAACTTTGAGCTGGCGTGGGACCCTGAAGTTCATGTACATCGCATCgg                                    <  1:524140/149‑1 (MQ=255)                                                                                                       gatcAGACCCTGGTACGTTTTGCTAACGGTAGCGCCCGGGTACT‑GTCGCGACTGATGTTGCTGCGCGTGGTCTGGATATTAAATCGCTTGAGCTGGTGGTGAACTTTGAGCTGGCGTGGGACCCTGACGTTCATGTACGTCGCTTCggc                                   >  1:541300/1‑148 (MQ=255)                                                                                                         tcAGACCCTGGTACGTTTTGCAAACGGTACCGCCCGAGTACTGGTCGCGACTGATGTTGCTGCGCGTGGTCTGGCTAGTAAATCGCTTGAGCTGGTGGTGAACTTCGCGCTGGCGTGGGACCCTGCAGGTCATGTACAACGGATCgggc                                  >  2:435472/1‑147 (MQ=255)                                                                                                              cccTGGTACGTTTTGCGAACGGTAGCGCGCTTGTACTGGTCGCGACTGATGTTGCTGCGCGTGGTCTGGATATTAAATCGCTTGAGCTGGTGGTGAACTTTGAGCTGGCGTGGGACCCTGAAGTTCATGTACATCGCATCGGTCGTACa                             <  2:118280/149‑1 (MQ=255)                                                                                                                cTGGTACGTTTTGCTAACGGTAGCGCCCGTGTACTGGTCGCGACTGATGTTGCTGCGCGTGGTCTGGATATTAAATCGCTTGAGCTGGTGGTGAACTTTGAGCTGGCGTGGGACCCTGAAGTTCATGTACATCGCATCGGTCGTACAGc                           <  1:384685/149‑1 (MQ=255)                                                                                                                  ggTACGTTTTGCTAACGGTAGCGCCCGTGTACT‑GTCGCGACTGATGTTGCTGCGCGTGGTCTGGATATTAAATCGCTTGAGCTGGTGGGGAACTTTGAGCTGGCGTGGGACCCTGAAGTTCATGTACATCGCATCGGTCGTACAGCTCg                        <  1:608663/149‑1 (MQ=255)                                                                                                                    tACGTTTTGCTAACGGTAGCGCCCGGGTGCTGGTCGCGACTGATGTTGCTGCGCGTGGTCTGGATATTAAATCGCTTGAGCTGGTGGTGAACTTTGAGCTGGCGTGGGACCCTGAAGTTCATGTACATCGCATCGGTCGTACAGCTCGt                       <  2:595504/149‑1 (MQ=255)                                                                                                                        ttttGCTAACGGTAGCGCCCGTGTACT‑GTCGCGACTGATGTTGCTGCGCGTGGTCTGGATATTAAATCGCTTGAGCTGGTGGTGAACTTTGAGCTGGCGTGGGACCCTGAAGTTCATGTACATCGCATCGGTCGTAAAGCTCGTGCAgg                  >  1:41660/1‑149 (MQ=255)                                                                                                                           tGCTAACGGTAGCGCCCGTGTACTGGTCGCGACTGATGTTGCTGCGCGTGGTCTGGATATTAAATCGCTTGAGCTGGTGGTGAACTTTGAGCTGGCGTGGGACCCTGAAGTTCATGTACATCGCATCGGTCGTACAGCTCGTGCAGGaa                <  2:362169/149‑1 (MQ=255)                                                                                                                                  ggTAGCGCCCGTGTACTGGTCGCGACGGGTGTTGCTGCGCGTGGTCTGGATATTAAATCGCTTGAGCTGGTGGTGAACTTTGAGCTGGCGTGGGACCCTGAAGTTCCTGTACATCGCATCGGCCGTACAGCTCCTGCAGGACATAGCgc         <  2:643956/149‑2 (MQ=255)                                                                                                                                   gTAGCGCCCGTGTACT‑GTCGCGACTGATGTTGCTGCGCGTGGTCTGGATATTAAATCGCTTGAGCTGGTGGTGAACTTTGAGCTGGCGTGGGACCCTGAAGTTCATGTACATCGCATCGGTCGTACAGCTCGTGCAGGAAATAGCGGTc       >  1:619191/1‑149 (MQ=255)                                                                                                                                        gcCCGTGTACT‑GTCGCGACTGATGTTGCTGCGCGTGGTCTGGATATTAAATCGCTTGAGCTGGTGGTGAACTTTGAGCTGGCGTGGGACCCTGAAGTTCATGTACATCGCATCGGTCGTACAGCTCGTGCAGGAAAGAGCGGTCTGGCg  >  2:329660/1‑149 (MQ=255)                                                                                                                                        gcCCGGGTACT‑GTCGCGACTGATGTTGCTGCGCGTGGTCTGGATATTAAATCGCTTGAGCTGGTGGTGAACTTTGAGCTGGCGTGGGACCCTGAAGGTCATGTACATCGCATCGGTCGGACGGCTCGTGCAGGAAATAGCGGTCTgggg  >  2:718050/1‑147 (MQ=255)                                                                                                                                                   |                                                                                                                                           TTTTGCAATACCAAAAAAGATTGCCAGGCTGTCTGCGACGCGCTGAATGAAGTAGGGCAAAGTGCATTGTCATTACACGGCGATTTGGAGCAACGCGATCGCGATCAGACCCTGGTACGTTTTGCTAACGGTAGCGCCCGTGTACTGGTCGCGACTGATGTTGCTGCGCGTGGTCTGGATATTAAATCGCTTGAGCTGGTGGTGAACTTTGAGCTGGCGTGGGACCCTGAAGTTCATGTACATCGCATCGGTCGTACAGCTCGTGCAGGAAATAGCGGTCTGGCG  >  NC_000913/1410249‑1410533

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