Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,109,551	C→T	P240S (CCA→TCA)	sbp →	sulfate transporter subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,109,551	0	C	T	100.0%	236.9 / NA	62	P240S (CCA→TCA)	sbp	sulfate transporter subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (26/36);  total (26/36)

GCGCGTGGCTCCACTAACACTTTTGTCGAGCGCGGAATTGGCGATGTACTGATTGCCTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGCCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGCCAGAAGGTCAGGAAATTGCCGCGAAAAACTACTACCGTCCGCGCGACGCTGAGGTGGCGA  >  NC_000913/4109407‑4109695                                                                                                                                                 |                                                                                                                                                 gcgcgTGGCTCCACTAACACTTTTGTCGAGCGCGGAATTGGCGATGTACTGATTGCCTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCg                                                                                                                                            <  1:626814/151‑1 (MQ=255)  cgcgTGGCTCCACTAACACTTTTGTCGAGCGCGGAATTGGCGATGTACTGATTGCCTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCgt                                                                                                                                           >  1:1006689/1‑151 (MQ=255)   gcgTGGCTCCACTAACACTTTTGTCGAGCGCGGAATTGGCGATGTACTGATTGCCTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCgtg                                                                                                                                          >  2:765091/1‑151 (MQ=255)      tGGCTCCACTAACACTTTTGTCGAGCGCGGAATTGGCGATGTACTGATTGCCTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCg                                                                                                                                       >  2:1085121/1‑151 (MQ=255)        gCTCCACTAACACTTTTGTCGAGCGCGGAATTGGCGATGTACTGATTGCCTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCggt                                                                                                                                     <  1:922828/151‑1 (MQ=255)           ccACTAACACTTTTGTCGAGCGCGGAATTGGCGATGTACTGATTGCCTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCggtg                                                                                                                                    <  2:850816/149‑1 (MQ=255)              cTAACACTTTTGTCGAGCGCGGAATTGGCGATGTACTGATTGCCTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCg                                                                                                                                <  1:1137235/150‑1 (MQ=255)                aaCACTTTTGTCGAGCGCGGAATTGGCGATGTACTGATTGCCTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATa                                                                                                                             <  1:735035/151‑1 (MQ=255)                 acttttttGTCGAGCGCGGAATTGGCGCTGTACTGATTGCCCGGGAAAACGAAGCCCTGCTGGCAGCGAATGCGCTGGGGAAAGATAAATTCTAATTCGTCCCGCCGAGTGAGTCTATCCTCGCCGGGTCAAGCGTGTCGGTGGTCGATaa                                                                                                                            <  2:357351/147‑1 (MQ=255)                     ttttGTCGAGCGCGGAATTGGCGATGTACTGATTGCCTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTg                                                                                                                        >  2:1239386/1‑151 (MQ=255)                      tttGTCGAGCGCGGAATTGGCGATGTACTGATTGCCTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGt                                                                                                                         <  1:378334/149‑1 (MQ=255)                      tttGTCGAGCGCGGAATTGGCGATGTACTGATTGCCTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTgg                                                                                                                       >  1:314753/1‑151 (MQ=255)                            gAGCGCGGAATTGGCGATGTACTGATTGCCTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCg                                                                                                                    >  2:614964/1‑148 (MQ=255)                               cgcgGAATTGGCGATGTACTGATTGCCTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGaaa                                                                                                               >  2:442985/1‑150 (MQ=255)                                 cgGAATTGGCGATGTACTGATTGCCTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGaaaaaa                                                                                                            <  1:728892/151‑1 (MQ=255)                                     aTTGGCGATGTACTGATTGCCTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGt                                                                                                         >  2:1160387/1‑150 (MQ=255)                                        ggCGATGTACTGATTGCCTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACt                                                                                                      <  1:375472/150‑1 (MQ=255)                                        ggCGATGTACTGATTGCCTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTa                                                                                                     <  1:1216479/151‑1 (MQ=255)                                         gCGATGTACTGATTGCCTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTa                                                                                                     >  2:45827/1‑150 (MQ=255)                                          cGATGTACTGATTGCCTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAGAGGTACTaa                                                                                                    >  2:1203073/1‑150 (MQ=255)                                          cGATGTACTGATTGCCTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTaa                                                                                                    <  2:452499/150‑1 (MQ=255)                                            aTGTACTGATTGCCTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTa                                                                                                     >  2:1241022/1‑147 (MQ=255)                                                aCTGATTGCCTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGt                                                                                              <  1:274430/150‑1 (MQ=255)                                                aCTGATTGCCTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTg                                                                                             <  2:226642/151‑1 (MQ=255)                                                       gCCTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGaa                                                                                       <  1:146772/150‑1 (MQ=255)                                                          tGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGcc                                                                                    >  1:11923/1‑150 (MQ=255)                                                              aaaaCGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACCAAAGAGGTGGCGGAAGTCTACCt                                                                               <  2:696318/151‑1 (MQ=255)                                                                    aaGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAat                                                                          >  2:707426/1‑150 (MQ=255)                                                                      gCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATAtc                                                                       <  2:889825/151‑1 (MQ=255)                                                                         ctgctgGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATAtct                                                                      <  1:32635/149‑1 (MQ=255)                                                                          tgctgGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATctct                                                                    >  1:587119/1‑150 (MQ=255)                                                                              ggCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACt                                                                 <  1:1022029/149‑1 (MQ=255)                                                                              ggCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACt                                                                 <  1:1210785/149‑1 (MQ=255)                                                                              ggCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACt                                                                 <  1:702764/149‑1 (MQ=255)                                                                                 aGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGcc                                                             >  2:566080/1‑150 (MQ=255)                                                                                 aGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGCCa                                                            <  2:839293/151‑1 (MQ=255)                                                                                  gCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGcc                                                             >  1:176763/1‑149 (MQ=255)                                                                                   cGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGcc                                                             >  2:211143/1‑148 (MQ=255)                                                                                         aaCTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGCCAGAAGGTc                                                     >  1:37054/1‑150 (MQ=255)                                                                                           cTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTAGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAACCTACCTGAAATATCTCTACTCGCCAGAAGGTCAg                                                   <  2:102279/150‑1 (MQ=255)                                                                                           cTGGGGAAAGATAAATTCAAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGCCAGAAGGTc                                                     <  1:1203059/148‑1 (MQ=255)                                                                                             ggggAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGCCAGAAGGTCAGGaa                                                <  1:614987/151‑1 (MQ=255)                                                                                             ggggAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGCCAGAAGGTCAGGa                                                 <  1:302323/150‑1 (MQ=255)                                                                                               ggAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGCCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCGACCCGAAATATCTCTACTCGCCAGAAGGTCAGGAAAt                                              <  1:243086/151‑1 (MQ=255)                                                                                                    gATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGCCAGAAGGTCAGGAAATTGCcg                                         <  1:313532/151‑1 (MQ=255)                                                                                                    gATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGCCAGAAGGTCAGGAAATTGCcg                                         <  1:270868/151‑1 (MQ=255)                                                                                                      tAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGCCAGAAGGTCAGGAAATTGCcgcg                                       >  1:108908/1‑151 (MQ=255)                                                                                                       aaaTTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGCCAGAAGGTCAGGAAATTGCcgcg                                       <  1:633225/150‑1 (MQ=255)                                                                                                             gAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGCCAGAAGGTCAGGAAATTGCCGCGAAAAact                                <  2:35688/151‑1 (MQ=255)                                                                                                               aaTCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGCCAGAAGGTCAGGAAATTGCCGCGAAAAactac                              >  1:866387/1‑151 (MQ=255)                                                                                                                     cACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGCCAGAAGGTCAGGAAATTGCCGCGAAAAACTACTACCg                         <  1:847731/150‑1 (MQ=255)                                                                                                                       cGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGCCAGAAGGTCAGGAAATTGCCGCGAAAAACTACTACCGTcc                      >  2:965764/1‑151 (MQ=255)                                                                                                                       cGCAGAGTGAGTCTATCCTCGCAGAATCAACCGTGACGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTGCTAAAGTGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGCCAGAAGGTCAGGAAATTGCCGCGAAAAACTACTACCGTcc                      <  2:639028/151‑1 (MQ=255)                                                                                                                           gagtgagtCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGCCAGAAGGTCAGGAAATTGCCGCGAAAAACTACTACCGTCcgcgc                  <  1:765101/151‑1 (MQ=255)                                                                                                                            agtgagtCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGCCAGAAGGTCAGGAAATTGCCGCGAAAAACTACTACCGTCcgcgcg                 >  2:121278/1‑151 (MQ=255)                                                                                                                              tgagtCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGCCAGAAGGTCAGGAAATTGCCGCGAAAAACTACTACCGTCCGCGCGAc               >  2:195174/1‑151 (MQ=255)                                                                                                                                agtCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGCCAGAAGGTCAGGAAATTGCCGCGAAAAACTACTACCGTCCGCGCGACGc             >  2:90517/1‑151 (MQ=255)                                                                                                                                      tCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGCCAGAAGGTCAGGAAATTGCCGCGAAAAACTACTACCGTCCGCGCGACGCTGAgg        >  1:980241/1‑150 (MQ=255)                                                                                                                                       ccTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGCCAGAAGGTCAGGAAATTGCCGCGAAAAACTACTACCGTCCGCGCGACGCTGAGGt       <  2:1035641/150‑1 (MQ=255)                                                                                                                                          cGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGCCAGAAGGTCAGGAAATTGCCGCGAAAAACTACTACCGTCCGCGCGACGCTGAGGTGGc    <  2:1201343/150‑1 (MQ=255)                                                                                                                                            cAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGCCAGAAGGTCAGGAAATTGCCGCGAAAAACTACTACCGTCCGCGCGACGCTGAGGTGGCGa  >  1:921963/1‑150 (MQ=255)                                                                                                                                             agagTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGCCAGAAGGTCAGGAAATTGCCGCGAAAAACTACTACCGTCCGCGCGACGCTGAGGTGGCGa  <  2:631084/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                 GCGCGTGGCTCCACTAACACTTTTGTCGAGCGCGGAATTGGCGATGTACTGATTGCCTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGCCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGCCAGAAGGTCAGGAAATTGCCGCGAAAAACTACTACCGTCCGCGCGACGCTGAGGTGGCGA  >  NC_000913/4109407‑4109695

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