Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	942,588	T→A	8.4%	C544S (TGC→AGC)	dmsA →	dimethyl sulfoxide reductase subunit A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	942,588	0	T	A	8.4%	200.9 / 6.4	67	C544S (TGC→AGC)	dmsA	dimethyl sulfoxide reductase subunit A
Reads supporting (aligned to +/- strand):  ref base T (33/28);  new base A (3/3);  total (36/31)
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.00e+00

ACCGTACCCATGAAATCCTTCAGGATGATAAGAAGTGCGAGCTGATTGTGGTTATCGACTGCCACATGACCTCATCGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGCGAACTGGCAAAACGTCTTGGCGTTGAGCAACAGTTTACTGAAGGCCGTACCCAGGAAGAGTGGATG  >  NC_000913/942439‑942737                                                                                                                                                      |                                                                                                                                                      aCCGTACCCATGAAATCCTTCAGGATGATAAGAAGTGCGAGCTGATTGTGGTTATCGACTGCCACATGACCTCATCGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCt                                                                                                                                                       >  2:865627/1‑150 (MQ=255)  ccGTACCCATGAAATCCTTCAGGATGATAAGAAGTGCGAGCTGATTGTGGTTATCGACTGCCACATGACCTCATCGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTg                                                                                                                                                      <  1:401484/150‑1 (MQ=255)     tACCCATGAAATCCTTCAGGATGATAAGAAGTGCGAGCTGATTGTGGGTATCGACTGCCACATGACCTCATCGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCAGGATGCATCCAGCgc                                                                                                                                                   <  1:1536018/150‑2 (MQ=255)      aCCCATGAAATCCTTCAGGATGATAAGAAGTGCGAGCTGATTGTGGTTATCGACTGCCACATGACCTCATCGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCggg                                                                                                                                                  <  1:1226195/150‑1 (MQ=255)       cccATGAAATCCTTCAGGATGATAAGAAGTGCGAGCTGATTGTGGTTATCGACTGCCACATGACCTCATCGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGa                                                                                                                                                 >  2:115282/1‑150 (MQ=255)       cccATGAAATCCTTCAGGATGATAAGAAGTGCGAGCTGATTGTGGTTATCGACTGCCACATGACCTCATCGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGa                                                                                                                                                 <  1:1176731/150‑1 (MQ=255)         cATGAAATCCTTCAGGATGATAAGAAGTGCGAGCTGATTGTGGTTATCGACTGCCACATGACCTCATCGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAat                                                                                                                                               <  2:841798/150‑1 (MQ=255)         cATGAAATCCTTCAGGATGATAAGAAGTGCGAGCTGATTGTGGTTATCGACTGCCACATGACCTCATCGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCAGGATGCATCCAGCgcaaag                                                                                                                                               <  1:748412/150‑6 (MQ=255)              aaTCCTTCAGGATGATAAGAAGTGCGAGCTGATTGTGGTTATCGACTGCCACATGACCTCATCGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTc                                                                                                                                          <  1:865627/150‑1 (MQ=255)                tCCTTCAGGATGATAAGAAGTGCGAGCTGATTGTGGTTATCGACTGCCACATGACCTCATCGGCGAAATATGGTGACATCCTGCTGCCTGACTTGACCCCTTCCGAATAGATGTACTTTGCGCTGGATGCATCCTGCGGGAATATGTCtt                                                                                                                                        <  2:134818/150‑1 (MQ=255)                 ccTTCAGGATTATAAGAAGTGCGAGCTGATTGTGGTTATCGACTGCCACATGACCTCATCGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTa                                                                                                                                       <  1:1477710/150‑1 (MQ=255)                 ccTTCAGGATGATAAGAAGTGCGAGCTGATTGTGGTTATCGACTGCCACATGACCTCATCGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTa                                                                                                                                       <  1:548465/150‑1 (MQ=255)                  cTTCAGGATGATAAGAAGTGCGAGCTGATTGTGGTTATCGACTGCCACATGACCTCATCGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGAACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTAc                                                                                                                                      <  1:1175119/150‑1 (MQ=255)                     aaGGATGATAAGAAGTGCGAGCTGATTGTGGTTATCGACTGCCACATGACCTCATCGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTg                                                                                                                                   <  2:173277/149‑1 (MQ=255)                          tgatAAGAAGTGCGAGCTGATTGTGGTTATCGACTGCCACATGACCTCATCGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGAtttt                                                                                                                              >  1:110882/1‑150 (MQ=255)                             tAAGAAGTGCGAGCTGATTGTGGTTATCGACTGCCACATGACCTCATCGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCaa                                                                                                                           >  1:791883/1‑150 (MQ=255)                             tAAGAAGTGCGAGCTGATTGTGGTTATCGACTGCCACATGACCTCATCGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCaa                                                                                                                           <  1:1293054/150‑1 (MQ=255)                                gaagTGCGAGCTGATTGTGGTTATCGACTGCCACATGACCTCATCGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGa                                                                                                                        >  2:1176237/1‑150 (MQ=255)                                     gCGAGCTGATTGTGGTTATCGACTGCCACATGACCTCATCGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAgg                                                                                                                   >  1:859879/1‑150 (MQ=255)                                                   gTTATCGACTGCCACATGACCTCATCGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCgcgc                                                                                                     <  2:110882/150‑1 (MQ=255)                                                    ttATCGACTGCCACATGACCTCATCGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCt                                                                                                    >  2:485515/1‑150 (MQ=255)                                                           cTGCCACATGACCTCATCGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGGTTAAACCGCGCTTTGAATg                                                                                             >  1:729216/1‑150 (MQ=255)                                                             gCCACATGACCTCATCGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTa                                                                                           >  2:884213/1‑150 (MQ=255)                                                              ccACATGACCTCATCGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTaa                                                                                          >  1:1342070/1‑150 (MQ=255)                                                               cacaTGACCTCATCGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGa                                                                                                >  2:1371328/1‑143 (MQ=255)                                                               cacaTGACCTCATCGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGa                                                                                                <  1:1371328/143‑1 (MQ=255)                                                               cacaTGACCTCATCGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAg                                                                                         >  2:1153583/1‑150 (MQ=255)                                                                   tGACCTCATCGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAGCGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCa                                                                                     <  2:354562/150‑1 (MQ=255)                                                                      ccTCATCGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGACCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCt                                                                                  >  1:1237095/1‑150 (MQ=255)                                                                         cATCGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATg                                                                               <  2:1098440/150‑1 (MQ=255)                                                                         cATCGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGAGCAGGTGATTAAACCGCGCTTTGAATGTAAGCCCATATATg                                                                               >  1:124955/1‑150 (MQ=255)                                                                            cGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGAC‑TTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAAt                                                                           >  2:338190/1‑150 (MQ=255)                                                                              gCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATg                                                                          >  1:38235/1‑150 (MQ=255)                                                                              gCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATg                                                                          <  2:1146659/150‑1 (MQ=255)                                                                               cGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGa                                                                         >  1:452696/1‑150 (MQ=255)                                                                                 aaaTATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCCATGAAATGAcc                                                                       <  2:859879/150‑1 (MQ=255)                                                                                  aaTATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCAGAGGGAATATGTCTTACGTGATTTTCAACCATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGACCa                                                                      >  2:1461682/1‑150 (MQ=255)                                                                                      tGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGCACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAACCCGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGCgc                                                                  >  2:414002/1‑149 (MQ=255)                                                                                        cTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGCGAAc                                                                >  2:300428/1‑150 (MQ=255)                                                                                          gACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCAACTATGAAATGACCAGCGAACTg                                                              >  2:1109130/1‑150 (MQ=255)                                                                                          gACATCCTGCTGCCTGACTGCACAGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGATTTGAATGTAATACCATCTATGAAATGACCAGCGAACTg                                                              >  2:495612/1‑150 (MQ=255)                                                                                           aCATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGTGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGCGAACTgg                                                             <  1:726935/150‑1 (MQ=255)                                                                                              tCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGc                                                                    <  2:932316/140‑1 (MQ=255)                                                                                              tCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGc                                                                    >  1:932316/1‑140 (MQ=255)                                                                                              tCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGCGAACTGGc                                                            >  1:294520/1‑148 (MQ=255)                                                                                              tCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGCGAACTGGc                                                            <  2:294520/148‑1 (MQ=255)                                                                                               ccTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGCGAACTGGCaaa                                                         >  2:1398361/1‑150 (MQ=255)                                                                                                     gcCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGCGAACTGGCAAAACGTCt                                                   >  2:1080844/1‑150 (MQ=255)                                                                                                      ccTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGCGAACTGGCAAAACGTCtt                                                  >  1:200966/1‑150 (MQ=255)                                                                                                      ccTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGCGAACTGGCAAAACGTCtt                                                  >  2:40140/1‑150 (MQ=255)                                                                                                          aCTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGCGAACTGGCAAAACGTCTTGGCg                                              <  1:115282/150‑1 (MQ=255)                                                                                                           cTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGCGAACTGGCAAAACGTCTTGGCGt                                             >  1:730724/1‑150 (MQ=255)                                                                                                             gCACCGCTTCCGAACAGATGGAC‑TTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGCGAACTGGCAAAACGTCTTGGCGTTGa                                          <  1:338190/150‑1 (MQ=255)                                                                                                              cACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGCGAACTGGCAAAACGTCTTGGCGTTGa                                          <  1:1194381/150‑1 (MQ=255)                                                                                                              cACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGCGAACTGGCAAAACGTCTTGGCGTTGa                                          <  1:5330/150‑1 (MQ=255)                                                                                                              cACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGCGAACTGGCAAAACGTCTTGGCGTTGa                                          >  2:35173/1‑150 (MQ=255)                                                                                                                ccGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGCGAACTGGCAAAACGTCTTGGCGTTGAGc                                        >  1:1438758/1‑150 (MQ=255)                                                                                                                        gAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGCCCATCTATGAAATGACCAGCGAACTGGCAAACCGTCTTGGCGGTGAGCAACAGttt                                >  1:1532515/1‑150 (MQ=255)                                                                                                                        gAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGCGAACTGGCAAAACGTCtt                                                  >  1:1435794/1‑132 (MQ=255)                                                                                                                        gAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGCGAACTGGCAAAACGTCtt                                                  <  2:1435794/132‑1 (MQ=255)                                                                                                                               tGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGCGAACTGGCAAAACGTCTTGGCGTTGAGCAACAGTTTACTGAAg                         <  1:1428400/150‑1 (MQ=255)                                                                                                                                 gACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGCGAACTGGCAAAACGTCTTGGCGTTGAGCAACAGTTTACTGAAGGc                       <  1:1153583/150‑1 (MQ=255)                                                                                                                                 gACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGCGAACTGGCAAAACGTCTTGGCGTTGAGCAACAGTTTACTGAAGGc                       <  1:414002/150‑1 (MQ=255)                                                                                                                                  aCTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGCGAACTGGCAAAACGTCTTGGCGTTGAGCAACAGTTTACTGAAGGcc                      >  1:837405/1‑150 (MQ=255)                                                                                                                                    ttcccgcAGGATGCATCCAGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGCCCAGCGAACTGGCAAAACGTCTTGGCGTTGAGCAACAGTTTACTGAAGGCCGt                    >  2:540540/5‑150 (MQ=255)                                                                                                                                     tcccgcAGGATGCATCCAGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGCGAACTGGCAAAACGTCTTGGCGTTGAGCAACAGTTTACTGAAGGc                       >  2:357209/4‑146 (MQ=255)                                                                                                                                     tcccgcAGGATGCATCCAGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGCGAACTGGCAAAACGTCTTGGCGTTGAGCAACAGTTTACTGAAGGc                       <  1:357209/143‑1 (MQ=255)                                                                                                                                            ggATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGCGAACTGGCAAAACGTCTTGGCGTTGAGCAACAGTTTACTGAAGGCCGTACCCAGGa            >  1:1110806/1‑150 (MQ=255)                                                                                                                                                      tGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGCGAACTGGCAAAACGTCTTGGCGTTGAGCAACAGTTTACTGAAGGCCGTACCCAGGAAGAGTGGATg  >  2:606709/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                      ACCGTACCCATGAAATCCTTCAGGATGATAAGAAGTGCGAGCTGATTGTGGTTATCGACTGCCACATGACCTCATCGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGCGAACTGGCAAAACGTCTTGGCGTTGAGCAACAGTTTACTGAAGGCCGTACCCAGGAAGAGTGGATG  >  NC_000913/942439‑942737

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