Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,144,160	C→G	9.4%	S2T (AGT→ACT)  E327D (GAG→GAC)	hybB ← hybA ←	hydrogenase 2 membrane subunit hydrogenase 2 iron‑sulfur protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,144,160	0	C	G	9.4%	229.6 / 11.1	73	S2T (AGT→ACT)  E327D (GAG→GAC)	hybB hybA	hydrogenase 2 membrane subunit hydrogenase 2 iron‑sulfur protein
Reads supporting (aligned to +/- strand):  ref base C (32/34);  new base G (3/4);  total (35/38)
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.14e-01

CCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAGACGCTTCACAATCAGGAGCATACAGATGACGATTAACGGTCCAAAAATCATGACCGGTTTACTGATGATTTTGCCGCCCAGCGGTTGTGGATCATGACTCATGATCGTCTCCTCCGTCGTGATGGTCGTTTTTGGTGTTGCGACGAACCAGCACGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTC  >  NC_000913/3144011‑3144309                                                                                                                                                      |                                                                                                                                                      cccACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAGACGCTTCACAATCAGGAGCATACAGATGACGATTAACGGTCCAAAAATCATGACCGGTTTACTGATGATTTTGCCGCCCAGCGGTTGTGGATCATGAc                                                                                                                                                       >  1:1268131/1‑150 (MQ=255)            ccgccgTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAGACGCTTCACAATCAGGAGCATACAGATGACGATTAACGGTCCAAAAATCATGACCGGTTTACTGATGATTTTGCCGCCCAGCGGTTGTGGATCATGAGTCATGATccac                                                                                                                                            <  1:436941/150‑4 (MQ=255)             cgccgTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAGACGCTTCACAATCAGGAGCATACAGATGACGATTAACGGTCCAAAAATCATGACCGGTTTACTGATGATTTTGCCGCCCAGCGGTTGTGGATCATGACTCATGATCGtct                                                                                                                                           <  2:669522/150‑1 (MQ=255)              gccgTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAGACGCTTCACAATCAGGAGCATACAGATGACGATTAACGGTCCAAAAATCATGACCGGTTTACTGATGATTTTGCCGCCCAGCGGTTGTGGATCATGAGTCATGATccacaa                                                                                                                                          <  1:166276/150‑6 (MQ=255)              gccgTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAGACGCTTCACAATCAGGAGCATACAGATGACGATTAACGGTCCAAAAATCATGACCGGTTTACTGATGATTTTGCCGCCCAGCGGTTGTGGATCATGACTCATGATCGtctc                                                                                                                                          >  2:1017469/1‑150 (MQ=255)              gccgTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAGACGCTTCACAATCAGGAGCATACAGATGACGATTAACGGTCCAAAAATCATGACCGGTTTACTGATGATTTTGCCGCCCAGCGGTTGTGGATCATGACTCATGATCGtctc                                                                                                                                          <  2:190672/150‑1 (MQ=255)                 gTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAGACGCTTCACAATCAGGAGCATACAGATGACGATTAACGGTCCAAAAATCATGACCGGTTTACTGATGATTTTGCCGCCCAGCGGTTGTGGATCATGACTCATGATCGTctcctc                                                                                                                                       <  2:1398431/150‑1 (MQ=255)                    cAGGTCAGAGACAGAGCCCAGACCGAACACCAGACGCTTCACAATCAGGAGCATACAGATGACGATTAACGGTCCAAAAATCATGACCGGTTTACTGATGATTTTGCCGCCCAGCGGTTGTGGATCATGACTCATGATCGTCTCCTCcgt                                                                                                                                    <  1:623572/150‑1 (MQ=255)                      ggTCAGAGACAGAGCCCAGACCGAACACCAGACGCTTCACAATCAGGAGCATACAGATGACGATTAACGGTCCAAAAATCATGACCGGTTTACTGATGATTTTGCCGCCCAGCGGTTGTGGATCATGACTCATGATCGTCTCCTCcgtcg                                                                                                                                  <  1:31620/150‑1 (MQ=255)                       gTCAGAGACAGAGCCCAGACCGAACACCAGACGCTTCACAATCAGGAGCATACAGATGACGATTAACGGTCCAAAAATCATGACCGGTTTACTGATGATTTTGCCGCCCAGCGGTTGTGGATCATGACTCATGATCGTCTCCTCcgtcgt                                                                                                                                 <  1:369470/150‑1 (MQ=255)                         cAGAGACAGAGCCCAGACCGAACACCAGACGCTTCACAATCAGGAGCATACAGATGACGATTAACGGTCCAAAAATCATGACCGGTTTACTGATGATTTTGCCGCCCAGCGGTTGTGGATCATGACTCATGATCGTCTCCTCCGTCGTGa                                                                                                                               <  2:1473293/150‑1 (MQ=255)                         cAGAGACAGAGCCCAGACCGAACACCAGACGCTTCACAATCAGGAGCATACAGATGACGATTAACGGTCCAAAAATCATGACCGGTTTACTGATGATTTTGCCGCCCAGCGGTTGTGGATCATGACTCATGATCGTCTCCTCCGTCGTGa                                                                                                                               >  2:1346896/1‑150 (MQ=255)                         cAGAGACAGAGCCCAGACCGAACACCAGACGCGTCACAATCAGGAGCATACAGATGACGATTAACGGTCCAACAATCATGACCGGTGTACTGATGATGTTGCCGCCCAGCGGTTGTTGATCATGACTCATGATCGTCTCCTCCGTCGTGa                                                                                                                               <  1:610921/150‑1 (MQ=255)                          agagACATAGCACAGACAGAACACCAGACGCTTCACAATCAGGAGCATACAGATGACGATTAACGGTCCAAAAATCATGACCGGTTTACTGATGATTTTGCCGCCCAGCGGTTGTGGATCATGACTCATGATCGTCTCCTCCGTCGTGAt                                                                                                                              <  1:692992/150‑1 (MQ=255)                               cAGAGCCCAGACCGAACACCAGACGCTTCACAATCAGGAGCATACAGATGACGATTAACGGTCCAAAAATCATGACCGGTTTACTGATGATTTTGCCGCCCAGCGGTTGTGGATCATGACTCATGATCGTCTCCTCCGTCGTGATGGTCg                                                                                                                         >  1:339572/1‑150 (MQ=255)                                 gagCCCAGACCGAACACCAGACGCTTCACAATCAGGAGCATACAGATGACGATTAACGGTCCAAAAATCATGACCGGTTTACTGATGATTTTGCCGCCCAGCGGTTGTGGATCATGACTCATGATCGTCTCCTCCGTCGTGATGGTCGtt                                                                                                                       <  2:217217/150‑1 (MQ=255)                                     ccAGACCGAACACCAGACGCTTCACAATCAGGAGCATACAGATGACGATTAACGGTCCAAAAATCATGACCGGTTTACTGATGATTTTGCCGCCCAGCGGTTGTGGATCATGACTCATGATCGTCTCCTCCGTCGTGATGGTCGTTTTTg                                                                                                                   <  2:56971/150‑1 (MQ=255)                                      cAGACCGAACACCAGACGCTTCACAATCAGGAGCATACAGATGACGATTAACGGTCCACAAATCATGACCGGTTTACTGAGGATTTTGCCGCCCAGCGGTTGT‑GATCAAGCCTCATGATCGTCTCCTCCGTCGTGATGGTCGTTGTTGgt                                                                                                                 >  2:305418/1‑150 (MQ=255)                                      cAGACCGAACACCAGACGCTTCACAATCAGGAGCATACAGATGACGATTAACGGTCCAAAAATCATGACCGGTTTACTGATGATTTTGCCGCCCAGCGGTTGTGGATCATGACTCATGATCGTCTCCTCCGTCGTGATGGGCGTTTTTgg                                                                                                                  <  2:859380/150‑1 (MQ=255)                                       aGACCGAACACCAGACGCTTCACAATCAGGAGCATACAGATGACGATTAACGGTCCAAAAATCATGACCGGTTTACTGATGATTTTGCCGCCCAGCGGTTGTGGATCATGACTCATGATCGTCTCCTCCGTCGTGATGGTCGTTTTTGgt                                                                                                                 <  1:862571/150‑1 (MQ=255)                                            gAACACCAGACGCTTCACAATCAGGAGCATACAGATGACGATTAACGGTCCAAAAATCATGACCGGTTTACTGATGATTTTGCCGCCCAGCGGTTGTGGATCATGACTCATGATCGTCTCCTCCGTCGTGATGGTCGTTTTTGGTGTTGc                                                                                                            <  2:1195219/150‑1 (MQ=255)                                            gAACACCAGACGCTTCACAATCAGGAGCATACAGATGACGATTAACGGTCCAAAAATCATGACCGGTTTACTGATGATTTTGCCGCCCAGCGGTTGTGGATCATGACTCATGATCGTCTCCTCCGTCGTGATGGTCGTTTTTGGTGTTGc                                                                                                            >  1:144207/1‑150 (MQ=255)                                               cacCAGACGCTTCACAATCAGGAGCATACAGATGACGATTAACGGTCCAAAAATCATGACCGGTTTACTGATGATTTTGCCGCCCAGCGGTTGTGGATCATGACTCATGATCGTCTCCTCCGTCGTGATGGTCGTTTTTGGTGTTGcgac                                                                                                         <  1:653259/150‑1 (MQ=255)                                               cacCAGACGCTTCACAATCAGGAGCATACAGATGACGATTAACGGTCCAAAAATCATGACCGGTTTACTGATGATTTTGCCGCCCAGCGGTTGTGGATCATGACTCATGATCGTCTCCTCCGTCGTGATGGTCGTTTTTGGTGTTGcgac                                                                                                         >  2:926122/1‑150 (MQ=255)                                                 ccAGACGCTTCACAATCAGGAGCATACAGATGACGATTAACGGTCCAAAAATCATGACCGGTTTACTGATGATTTTGCCGCCCAGCGGTTGTGGATCATGACTCATGATCGTCTCCTCCGTCGTGATGGTCGTTTTTGGTGTTGcgacga                                                                                                       >  2:800685/1‑150 (MQ=255)                                                   aGACGCTTCACAATCAGGAGCATACAGATGACGATTAACGGTCCAAAAATCATGACCGGTTTACTGATGATTTTGCCGCCCAGCGGTTGTGGATCATGACTCATGATCGTCTCCTCCGTCGTGATGGTCGTTTTTGGTGTTGCGACGAAc                                                                                                     >  1:771832/1‑150 (MQ=255)                                                          tCACAATCAGGAGCATACAGATGACGATTAACGGTCCAAAAATCATGACCGGTTTACTGATGATTTTGCCGCCCAGCGGTTGTGGATCATGACTCATGATCGTCTCCTCCGTCGTGATGGTCGTTTTTGGTGTTGCGACGAACCAGCACg                                                                                              >  2:1507261/1‑150 (MQ=255)                                                           cacaATCAGGAGCATACAGATGACGATTAACGGTCCAAAAATCATGACCGGTTTACTGATGATTTTGCCGCCCAGCGGTTGTGGATCATGACTCATGATCGTCTCCTCCGTCGTGATGGTCGTTTTTGGTGTTGCGACGAACCAGCACgg                                                                                             >  2:463157/1‑150 (MQ=255)                                                           cacaATCAGGAGCATACAGATGACGATTAACGGTCCAAAAATCATGACCGGTTTACTGATGATATTGCCGCCCAGCGGTTGTGGATCATGACTCATGATCCTCTCCTCCGTCCTGATGGTCCTATTAGGTCTTGCGACGATCCAGCACgg                                                                                             >  1:996579/1‑150 (MQ=255)                                                             caATCAGGAGCATACAGATGACGATTAACGGTCCAAAAATCATGACCGGTTTACTGATGATTTTGCCGCCCAGCGGTTGTGGATCATGACTCATGATCGTCTCCTCCGTCGTGATGGTCGTTTTTGGTGTTGCGACGAACCAGCACGGtt                                                                                           >  2:69012/1‑150 (MQ=255)                                                                 cAGGAGCATACAGATGACGATTAACGGTCCAAAAATCATGACCGGTTTACTGATGATTTTGCCGCCCAGCGGTTGTGGATCATGACTCATGATCGTCTCCTCCGTCGTGATGGTCGTTTTTGGTGTTGCGACGAACCAGCACGGTTAAGc                                                                                       <  2:771832/150‑1 (MQ=255)                                                                  aGGAGCATACAGATGACGATTAACGGTCCAAAAATCATGACCGGTTTACTGATGATTTTGCCGCCCAGCGGTTGTGGATCATGACTCATGATCGTCTCCTCCGTCGTGATGGTCGTTTTTGGTGTTGCGACGAACCAGCACGGTTAAGcc                                                                                      >  2:1078540/1‑150 (MQ=255)                                                                   ggAGCATACAGATGACGATTAACGGTCCAAAAATCATGACCGGTTTACTGATGATTTTGCCGCCCAGCGGTTGTGGATCATGACTCATGATCGTCTCCTCCGTCGTGATGGTCGTTTTTGGTGTTGCGACGAACCAGCACGGTTAAGccc                                                                                     >  2:312978/1‑150 (MQ=255)                                                                        aTACAGATGACGATTAACGGTCCAAAAATCATGACCGGTTTACTGATGATTTTGCCGCCCAGCGGTTGTGGATCATGACTCATGATCGTCTCCTCCGTCGTGATGGTCGTTTTTGGTGTTGCGACGAACCAGCACGGTTAAGCCCGCCAg                                                                                >  1:44647/1‑150 (MQ=255)                                                                         tACAGATGACGATTAACGGTCCAAAAATCATGACCGGTTTACTGATGATTCTGCCGCCCAGCGGTTGTGGATCATGACTCATGATCGTCTCCTCCGTCGTGATGGTCGTTTTTGGTGTTGCGACGAACCAGCACGGTTAAGCCCGCCAGc                                                                               >  2:560644/1‑150 (MQ=255)                                                                          aCAGATGACGATTAACGGTCCAAAAATCATGACCGGTTTACTGATGATTTTGCCGCCCAGCGGTTGTGGATCATGACTCATGATCGTCTCCTCCGTCGTGATGGTCGTTTTTGGTGTTGCGACGAACCAGCACGGTTAAGCCCGCCAGca                                                                              <  2:339572/150‑1 (MQ=255)                                                                               tGACGATTAACGGTCCAAAAATCATGACCGGTTTACTGATGATTTTGCCGCCCAGCGGTTGTGGATCATGACTCATGATCGTCTCCTCCGTCGTGATGGTCGTTTTTGGTGTTGCGACGAACCAGCACGGTTAAGCCCGCCAGCACAGcc                                                                         <  1:926122/150‑1 (MQ=255)                                                                                taCGATTAACGGTCCAAAAATCATGACCGGTTTACTGATGATTTTGCCGCCCAGCGGTTGTGGATCATGACTCATGATCGTCTCCTCCGTCGTGATGGTCGTTTTTGGTGTTGCGACGAACCAGCACGGTTAAGCCCGCCAGCACAGCCa                                                                        <  1:445103/149‑1 (MQ=255)                                                                                  cGATTAACGGTCCAAAAATCATGGCCGGTTTACTGATGATTTTGCCGCCCAGCGGTTGTGGATCATGACTCATGATCGTCTCCTCCGTCGTGATGGTCGTTTTTGGTGTTGCGACGAACCAGCACGGTTAAGCCCGCCAGCCCAGCCAGt                                                                      >  2:440134/1‑150 (MQ=255)                                                                                  cGATTAACGGTCCAAAAATCATGACCGGTTTACTGATGATTTTGCCGCCCAGCGGTTGTGGATCATGACTCATGATCGTCTCCTCCGTCGTGATGGTCGTTTTTGGTGTTGCGACGAACCAGCACGGTTAAGCCCGCCAGCACAGCCAGt                                                                      >  1:343179/1‑150 (MQ=255)                                                                                   gATTAACGGTCCAAAAATCATGACCGGTTTACTGATGATTTTGCCGCCCAGCGGTTGTGGATCATGACTCATGATCGTCTCCTCCGTCGTGATGGTCGTTTTTGGTGTTGCGACGAACCAGCACGGTTAAGCCCGCCAGCACAGCCAGTg                                                                     >  2:898183/1‑150 (MQ=255)                                                                                     ttAACGGTCCAAAAATCATGACCGGTTTACTGATGATTTTGCCGCCCAGCGGTTGTGGATCATGACTCATGATCGTCTCCTCCGTCGTGATGGTCGTTTTTGGTGTTGCGACGAACCAGCACGGTTAAGCCCGCCAGCACAGCCAGTGGt                                                                   >  1:241892/1‑150 (MQ=255)                                                                                       aaCGGTCCAAAAATCATGACCGGTTTACTGATGATTTTGCCGCCCAGCGGTTGTGGATCATGACTCATGATCGTCTCCTCCGTCGTGATGGTCGTTTTTGGTGTTGCGACGAACCAGCACGGTTAAGCCCGCCAGCACAGCCAGTGGTAg                                                                 >  1:1006052/1‑150 (MQ=255)                                                                                           gTCCAAAAATCATGACCGGTTTACTGATGATTTTGCCGCCAAGCGGTTGTGGATCATGACTCATGATCGTCTCCTCCGTCGTGATGGTCGTTTTTGGTGTTGCGACGAACCAGCACGGTTAAGCCCGCCAGCACAGCCAGTGGTAGcatc                                                             >  2:628989/1‑150 (MQ=255)                                                                                            tCCAAAAATCATGACCGGTTTACTGATGATTTTGCCGCCCAGCGGTTGTGGATCATGACTCATGATCGTCTCCTCCGTCGTGATGGTCGTTTTTGGTGTTGCGACGAACCAGCACGGTTAAGCCCGCCAGCACAGCCAGTGGTAGcatca                                                            <  2:1345320/150‑1 (MQ=255)                                                                                                aaaaTCATGACCGGTTTACTGATGATTTTGCCGCCCAGCGGTTGTGGATCATGACTCATGATCGTCTCCTCCGTCGTGATGGTCGTTTTTGGTGTTGCGACGAACCAGCACGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGcc                                                        <  1:1507261/150‑1 (MQ=255)                                                                                                  aaTCATGACCGGTTTACTGATGATTTTGCCGCCCAGCGGTTGTGGATCATGACTCATGATCGTCTCCTCCGTCGTGATGGTCGTTTTTGGTGTTGCGACGAACCAGCACGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCtt                                                      <  2:224514/150‑1 (MQ=255)                                                                                                       tGACCGGTTTACTGATGATTTTGCCGCCCAGCGGTTGTGGATCATGACTCATGATCGTCTCCTCCGTCGTGATGGTCGTTTTTGGTGTTGCGACGAACCAGCACGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATAc                                                 <  2:1624883/150‑1 (MQ=255)                                                                                                        gACCGGTTTACTGATGATTTTGCCGCCCAGCGGTTGTGGATCATGACTCATGATCGTCTCCTCCGTCGTGATGGTCGTTTTTGGTGTTGCGACGAACCAGCACGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACa                                                >  1:1086119/1‑150 (MQ=255)                                                                                                        gACCGGTTTACTGATGATTTTGCCGCCCAGCGGTTGTGGATCATGACTCATGATCGTCTCCTCCGTCGTGATGGTCGTTTTTGGTGTTGCGACGAACCAGCACGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACa                                                >  2:280319/1‑150 (MQ=255)                                                                                                         aCCGGTTTACTGTTGATTTTGCCGCCCAGCGGTTGTGGTTCATGAGTCATGATCGTCTCCTACGTCGTGATGGTCGTTTTTGGTGTTGAGACGAAGCAGAAGGGTTTAGCCCGCCAGCACACCCAGTGGTAGCATCATGCCTTTATACAg                                               <  1:153451/150‑1 (MQ=255)                                                                                                           cGGTTTACTGATGATTTTGCCGCCCAGCGGTTGTGGATCATGACTCATGATCGTCTCCTCCGTCGTGATGGTCGTTTTTGGTGTTGCGACGAACCAGCACGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAggg                                             <  1:1346896/150‑1 (MQ=255)                                                                                                            ggTTTACTGATGATTTTGCCGCCCAGCGGTTGTGGATCATGACTCATGATCGTCTCCTCCGTCGTGATGGTCGTTTTTGGTGTTGCGACGAACCAGCACGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGgt                                            <  1:1017469/150‑1 (MQ=255)                                                                                                              tttACTGATGATTTTGCCGCCCAGCGGTTGTGGATCATGACTCATGATCGTCTCCTCCGTCGTGATGGTCGTTTTTGGTGTTGCGACGAACCAGCACGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGgtgt                                          <  2:1043265/150‑1 (MQ=255)                                                                                                               ttACTGATGATTTTGCCGCCCAGCGGTTGTGGATCATGACTCATGATCGTCTCCTCCGTCGTGATGGTCGTTTTTGGTGTTGCGACGAACCAGCACGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGgtgtg                                         <  1:645391/150‑1 (MQ=255)                                                                                                                 aCTGATGATTTTGCCGCCCAGCGGTTGTGGATCATGACTCATGATCGTCTCCTCCGTCGTGATGGTCGTTTTTGGTGTTGCGACGAACCAGCACGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGtt                                       <  1:463157/150‑1 (MQ=255)                                                                                                                  cTGATGATTTTGCCGCCCAGCGGTTGTGGATCATGACTCATGATCGTCTCCTCCGTCGTGATGGTCGTTTTTGGTGTTGCGACGAACCAGCACGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTg                                      >  1:622770/1‑150 (MQ=255)                                                                                                                  cTGATGATTTTGCCGCCCAGCGGTTGTGGATCATGACTCATGATCGTCTCCTCCGTCGTGATGGTCGTTTTTGGTGTTGCGACGAACCAGCACGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTg                                      <  2:44647/150‑1 (MQ=255)                                                                                                                   tgatgaTTTTGCCGCCCAGCGGTTGTGGATCATGACTCATGATCGTCTCCTCCGTCGTGATGGTCGTTTTTGGTGTTGCGACGAACCAGCACGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGa                                     >  1:1561925/1‑150 (MQ=255)                                                                                                                      tgatTTTGCCGCCCAGCGGTTGTGGATCATGACTCATGATCGTCTCCTCCGTCGTGATGGTCGTTTTTGGTGTTGCGACGAACCAGCACGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGAata                                  >  1:1159749/1‑150 (MQ=255)                                                                                                                           ttGCCGCCCAGCGGTTGTGGATCATGACTCATGATCGTCTCCTCCGTCGTGATGGTCGTTTTTGGTGTTGCGACGAACCAGCACGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGAATATTTTc                             <  2:807274/150‑1 (MQ=255)                                                                                                                             gccgccCAGCGGTTGTGGATCATGACTCATGATCGTCTCCTCCGTCGTGATGGTCGTTTTTGGTGTTGCGACGAACCAGCACGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGAATATTTTCgg                           <  1:117320/150‑1 (MQ=255)                                                                                                                                        gTTGTGGATCATGACTCATGATCGTCTCCTCCGTCGTGATGGTCGTTTTTGGTGTTGCGACGAACCAGCACGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGAATATTTTCGGAACGCGCACCg                <  1:1230834/150‑1 (MQ=255)                                                                                                                                          tgtgGATCATGACTCATGATCGTCTCCTCCGTCGTGATGGTCGTTTTTGGTGTTGCGACGAACCAGCACGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGAATATTTTCGGAACGCGCACCGGt              >  1:1420043/1‑150 (MQ=255)                                                                                                                                           gacgATCATGAGTCATGATCGTCTCCTCCGTCGTGATGGTCGTTTTTGGTGTTGCGACGAACCAGCACGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGAATATTTTCg                            <  2:1553200/132‑1 (MQ=255)                                                                                                                                           gacgATCATGAGTCATGATCGTCTCCTCCGTCGTGATGGTCGTTTTTGGTGTTGCGACGAACCAGCACGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGAATATTTTCg                            >  1:1553200/4‑135 (MQ=255)                                                                                                                                           gacgATCATGAGTCATGATCGTCTCCTCCGTCGTGATGGTCGTTTTTGGTGTTGCGACGAACCAGCACGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGAATATTTTCGGAACGCGCACCGGTa             >  2:1205685/4‑150 (MQ=255)                                                                                                                                            tgGATCATGACTCATGATCGTCTCCTCCGTCGTGATGGTCGTTTTTGGTGTTGCGACGAACCAGCACGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGAATATTTTCGGAACGCGCACCGGTAg            >  2:321427/1‑150 (MQ=255)                                                                                                                                             ggATCATGACTCATGATCGTCTCCTCCGTCGTGATGGTCGTTTTTGGTGTTGCGACGAACCAGCACGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGAATATTTTCGGAACGCGCACCGGTAGa           >  1:176692/1‑150 (MQ=255)                                                                                                                                             cgATCATGAGTCATGATCGTCTCCTCCGTCGTGATGGTCGTTTTTGGTGTTGCGACGAACCAGCACGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGAATATTTTCGGAACGCGCACCGGTAGa           >  2:1104704/2‑150 (MQ=255)                                                                                                                                               aTCATGACTCATGATCGTCTCCTCCGTCGTGATGGTCGTTTTTGGTGTTGCGACGAACCAGCACGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGAATATTTTCGGAACGCGCACCGGTAGaaa         <  1:132635/150‑1 (MQ=255)                                                                                                                                                  aTGACTCATGATCGTCTCCTCCGTCGTGATGGTCGTTTTTGGTGTTGCGACGAACCAGCACGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGAATATTTTCGGAACGCGCACCGGTAGAAAGAt      >  2:472020/1‑150 (MQ=255)                                                                                                                                                   tGACTCATGATCGTCTCCTCCGTCGTGATGGTCGTTTTTGGTGTTGCGACGAACCAGCACGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGAATATTTTCGGAACGCGCACCGGTa             <  2:642865/142‑1 (MQ=255)                                                                                                                                                   tGACTCATGATCGTCTCCTCCGTCGTGATGGTCGTTTTTGGTGTTGCGACGAACCAGCACGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGAATATTTTCGGAACGCGCACCGGTa             >  1:642865/1‑142 (MQ=255)                                                                                                                                                      cTCATGATCGTCTCCTCCGTCGTGATGGTCGTTTTTGGTGTTGCGACGAACCAGCACGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTc  >  2:1152569/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                      CCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAGACGCTTCACAATCAGGAGCATACAGATGACGATTAACGGTCCAAAAATCATGACCGGTTTACTGATGATTTTGCCGCCCAGCGGTTGTGGATCATGACTCATGATCGTCTCCTCCGTCGTGATGGTCGTTTTTGGTGTTGCGACGAACCAGCACGGTTAAGCCCGCCAGCACAGCCAGTGGTAGCATCATGCCTTTATACAGGGTGTGTTGAATATTTTCGGAACGCGCACCGGTAGAAAGATCGTC  >  NC_000913/3144011‑3144309

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