Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	968,109	C→T	5.8%	R497C (CGT→TGT)	msbA →	lipid ABC transporter permease/ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	968,109	0	C	T	5.8%	222.5 / 2.4	69	R497C (CGT→TGT)	msbA	lipid ABC transporter permease/ATPase
Reads supporting (aligned to +/- strand):  ref base C (30/35);  new base T (2/2);  total (32/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.05e-01

AAATTGAAGAAGCGGCGCGTATGGCCTACGCCATGGACTTCATCAATAAGATGGATAACGGTCTCGATACAGTGATTGGTGAAAACGGCGTGCTGCTCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAA  >  NC_000913/967966‑968251                                                                                                                                                |                                                                                                                                               aaaTTGAAGAAGCGGCGCGTATGGCCTACGCCATGGACTTCATCAATAAGATGGATAACGGTCTCGATACAGTGATTGGTGAAAACGGCGTGCTGCTCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGAt                                                                                                                                           >  3:578456/1‑149 (MQ=255) aaaTTGAAGAAGCGGCGCGTATGGCCTACGCCATGGACTTCATCAATAAGATGGATAACGGTCTCGATACAGTGATTGGTGAAAACGGCGTGCTGCTCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGAt                                                                                                                                           >  4:428015/1‑149 (MQ=255)    ttGAAGAAGCGGCGCGTATGGCCTACGCCATGGACTTCATCAATAAGATGGATAACGGTCTCGATACAGTGATTGGTGAAAACGGCGTGCTGCTCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGc                                                                                                                                        >  2:429435/1‑149 (MQ=255)    ttGAAGAAGCGGCGCGTATGGCCTACGCCATGGACTTCATCAATAAGATGGATAACGGTCTCGATACAGTGATTGGTGAAAACGGCGTGCTGCTCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGc                                                                                                                                        >  4:28354/1‑149 (MQ=255)     tGAAGAAGCGGCGCGTATTTCCTACGCCATGGACTTCATCAATAAGATGGATATCGGTCTAGATACTGTGATTGGTGAAAACGGCGTGCTGCTCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCGCAATCCTTGTTGCGGGATAGcc                                                                                                                                       <  4:86936/149‑1 (MQ=255)     tGAAGAAGCGGCGCGTATGGCCTACGCCATGGACTTCATCAATAAGATGGATAACGGTCTCGATCCAGTGATTGGTGAAAACGGCGTGCTGCTCTCTGGAGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGcc                                                                                                                                       >  4:485728/1‑149 (MQ=255)         gaagCGGCGCGTATGGCCTACGCCATGGACTTCATCAATAAGATGGATAACGGTCTCGATACAGTGATTGGTGAAAACGGCGTGCTGCTCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCgat                                                                                                                                   >  2:34237/1‑149 (MQ=255)              ggcgCGTATGGCCTACGCCATGGACTTCATCAATAAGATGGATAACGGTCTCGATACAGTGATTGGTGAAAACGGCGTGCTGCTCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCgattctga                                                                                                                              <  4:578456/149‑1 (MQ=255)              ggcgCGTATGGCCTACGCCATGGACTTCATCAATAAGATGGATAACGGTCTCGATACAGTGATTGGTGAAAACGGCGTGCTGCTCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCgattctga                                                                                                                              <  1:469900/149‑1 (MQ=255)               gcgcGTATGGCCTACGCCATGGACTTCATCAATAAGATGGATAACGGTCTCGATACAGTGATTGGTGAAAACGGCGTGCTGCTCTCTGGTGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCgattctgat                                                                                                                             <  1:465722/149‑1 (MQ=255)                 gcgTATGGCCTACGCCATGGACTTCATCAATAAGATGGATAACGGTCTCGATACAGTGATTGGTGAAAACGGCGTGCTGCTCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCgattctgattc                                                                                                                           >  4:156634/1‑149 (MQ=255)                  cgTATGGCCTACGCCATGGACTTCATCAATAAGATGGATAACGGTCTCGATACAGTGATTGGTGAAAACGGCGTGCTGCTCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCgattctgattct                                                                                                                          <  2:544080/149‑1 (MQ=255)                    tATGGCCTACGCCATGGACTTCATCAATAAGATGGATAACGGTCTCGATACAGTGATTGGTGAAAACGGCGTGCTGCTCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTgg                                                                                                                        <  4:322012/149‑1 (MQ=255)                         ccTACGCCATGGACTTCATCAATAAGATGGATAACGGTCTCGATACAGTGATTGGTGAAAACGGCGTGCTGCTCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGaa                                                                                                                   <  3:28354/149‑1 (MQ=255)                         ccTACGCCATGGACTTCATCAATAAGATGGATAACGGTCTCGATACAGTGATTGGTGAAAACGGCGTGCTGCTCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGaa                                                                                                                   >  3:221777/1‑149 (MQ=255)                         ccTACGCCATGGACTTCATCAATAAGATGGATAACGGTCTCGATACAGTGATTGGTGAAAACGGCGTGCTGCTCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGaa                                                                                                                   >  2:402637/1‑149 (MQ=255)                          cTACGCCATGGACTTCATCAATAAGATGGATAACGGTCTCGATACAGTGATTGGTGAGAACGGCGTGCTGCTCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAg                                                                                                                  >  2:486520/1‑149 (MQ=255)                          cTACGCCATGGACTTCATCAATAAGATGGATAACGGTCTCGATACAGTGATTGGTGAAAACGGCGTGCTGCTCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAg                                                                                                                  <  2:290733/149‑1 (MQ=255)                                 aTGGACTTCATCAATAAGATGGATAACGGTCTCGATACAGTGATTGGTGAAAACGGCGTGCTGCTCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTc                                                                                                           <  2:388289/149‑1 (MQ=255)                                      cTTCATCAATAAGATGGATAACGGTCTCGATACAGTGATTGGTGAAAACGGCGTGCTGCTCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGctc                                                                                                      >  1:454281/1‑149 (MQ=255)                                       ttCATCAATAAGATGGATAACGGTCTCGATACAGTGATTGGTGAAAACGGCGTGCTGCTCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGctct                                                                                                     >  1:219377/1‑149 (MQ=255)                                        tcatcaATAAGATGGATAACGGTCTCGATACAGTGATTGGTGAAAACGGCGTGCTGCTCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTg                                                                                                    <  1:32743/149‑1 (MQ=255)                                             aaTAAGATGGATAACGGTCTCGATACAGTGATTGGTGAAAACGGCGTGCTGCTCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATAc                                                                                               <  1:204667/149‑1 (MQ=255)                                                 aGATGGATAACGGTCTCGATACAGTGATTGGTGAAAACGGCGTGCTGCTCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGaa                                                                                           <  3:263410/149‑1 (MQ=255)                                                     ggATAACGGTCTCGATACAGTGATTGGTGAAAACGGCGTGCTGCTCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCg                                                                                       >  1:476386/1‑149 (MQ=255)                                                        tAACGGTCTCGATACAGTGATTGGTGAAAACGGCGTGCTGCTCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAAc                                                                                    <  4:266425/149‑1 (MQ=255)                                                         aaCGGTCTCGATACAGTGATTGGTGAAAACGGCGTGCTGCTCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGTGTGATAGCCCgattctgtttct                                                                                                                          >  2:35292/1‑110 (MQ=255)                                                         aaCGGTCTCGATACAGTGATTGGTGAAAACGGCGTGCTGCTCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGTGTGATAGCCCgattctgattct                                                                                                                          <  1:35292/110‑1 (MQ=255)                                                             gTCTCGATACAGTGATTGGTGAAAACGGCGTGCTGCTCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCg                                                                               >  3:450612/1‑149 (MQ=255)                                                              tctcGATACAGTGATTGGTGAAAACGGCGTGCTGCTCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGa                                                                              <  2:412448/149‑1 (MQ=255)                                                              tctcGATACAGTGATTGGTGAAAACGGCGTGCTGCTCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGa                                                                              <  3:248945/149‑1 (MQ=255)                                                                tcGATACAGTGATTGGTGAAAACGGCGTGCTGCTCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGAtt                                                                            >  1:3583/1‑149 (MQ=255)                                                                     aCAGTGATTGGTGAAAACGGCGTGCTGCTCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAggc                                                                       <  2:26421/149‑1 (MQ=255)                                                                            ttGGTGAAAACGGCGTGCTGCTCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTg                                                                >  3:378301/1‑149 (MQ=255)                                                                              ggTGAAAACGGCGTGCTGCTCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGa                                                              <  4:498864/149‑1 (MQ=255)                                                                              ggTGAAAACGGCGTGCTGCTCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGa                                                              <  3:299288/149‑1 (MQ=255)                                                                               gTGAAAACGGCGTGCTGCTCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGAt                                                             <  2:346834/149‑1 (MQ=255)                                                                               gTGAAAACGGCGTGCTGCTCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGAt                                                             >  3:171924/1‑149 (MQ=255)                                                                               gTGAAAACGGCGTGCTGCTCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGAt                                                             <  2:256913/149‑1 (MQ=255)                                                                                  aaaaCGGCGTGCTGCTCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAg                                                          >  3:187690/1‑149 (MQ=255)                                                                                   aaaCGGCGTGCTGCTCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGt                                                         <  3:428015/149‑1 (MQ=255)                                                                                    aaCGGCGTGCTGCTCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGtt                                                        >  1:121823/1‑149 (MQ=255)                                                                                      cGGCGTGCTGCTCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACCGGATGAGTTGc                                                      <  3:342855/149‑1 (MQ=255)                                                                                              tgctCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACt                                                                 >  3:308851/1‑130 (MQ=255)                                                                                              tgctCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACt                                                                 <  4:308851/130‑1 (MQ=255)                                                                                              tgctCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACt                                                                 <  2:459079/130‑1 (MQ=255)                                                                                              tgctCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACt                                                                 >  1:459079/1‑130 (MQ=255)                                                                                              tgctCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAAc                                              >  4:112296/1‑149 (MQ=255)                                                                                               gctCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAAcc                                             <  4:221777/149‑1 (MQ=255)                                                                                               gctCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAAcc                                             <  1:56808/149‑1 (MQ=255)                                                                                                 tctctGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGt                                           >  3:270942/1‑149 (MQ=255)                                                                                                  ctctGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTa                                          <  3:165082/149‑1 (MQ=255)                                                                                                   tctGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTAc                                         >  2:161823/1‑149 (MQ=255)                                                                                                           tcagcgtcagcgTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTgg                                 <  4:548840/149‑1 (MQ=255)                                                                                                           tcagcgtcagcgTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTgg                                 <  4:378301/149‑1 (MQ=255)                                                                                                                      gtATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCga                                                                                                                                    <  3:466365/39‑1 (MQ=255)                                                                                                                      gtATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCga                                                                                                                                    >  4:466365/1‑39 (MQ=255)                                                                                                                      gtATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCAc                      >  4:397796/1‑149 (MQ=255)                                                                                                                              aTCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTc              >  4:437542/1‑149 (MQ=255)                                                                                                                               tCGTTCGAGCCTTGTTGCGTGATAGCCGGATTCTGATTCTGGACTAAGCTACCTCGGCTCTGGATACCGAAGCCGAACGTGCGATTCAGGGGGCACTGGATGAGTTGCAGAAAAACCGTACCGCTCTGGTGATTGCCCACCGCTTGTCt             <  3:485728/149‑1 (MQ=255)                                                                                                                                   tCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGa                                                           <  1:172932/99‑1 (MQ=255)                                                                                                                                   tCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGa                                                           >  2:172932/1‑99 (MQ=255)                                                                                                                                    cGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCAt        <  3:397796/149‑1 (MQ=255)                                                                                                                                     gAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCAtt       <  4:450612/149‑1 (MQ=255)                                                                                                                                       gCCTTGTTGTGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGaa                                                                                     >  4:132266/1‑69 (MQ=255)                                                                                                                                       gCCTTGTTGTGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGaa                                                                                     <  3:132266/69‑1 (MQ=255)                                                                                                                                          ttgttgCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGa                                                              <  2:206352/89‑1 (MQ=255)                                                                                                                                          ttgttgCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGa                                                              >  1:206352/1‑89 (MQ=255)                                                                                                                                          ttgttgCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGaaaa  >  1:199324/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                               AAATTGAAGAAGCGGCGCGTATGGCCTACGCCATGGACTTCATCAATAAGATGGATAACGGTCTCGATACAGTGATTGGTGAAAACGGCGTGCTGCTCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAA  >  NC_000913/967966‑968251

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