Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	741,486	G→A	5.5%	A357A (GCC→GCT)	dtpD ←	dipeptide and tripeptide permease D

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	741,486	0	G	A	5.5%	226.0 / 2.6	72	A357A (GCC→GCT)	dtpD	dipeptide and tripeptide permease D
Reads supporting (aligned to +/- strand):  ref base G (36/32);  new base A (2/2);  total (38/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.04e-01

CGTAATTTGCGACATGGCAACCGGGTCGATAAACAGTTCCGCAAAGCCCATCACCGCCAGGCCTAATACCATCAGTGGCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCGACACTCTCTTTTACCACCCACGCCAGGAACACACCGCACAGCATAACTGCGAAGGCATTAATCGACTGGAACATCGCGGT  >  NC_000913/741348‑741629                                                                                                                                           |                                                                                                                                                cgTAATTTGCGACATGGCAACCGGGTCGATAAACAGTTCCGCAAAGCCCATCACCGCCAGGCCTAATACCATCAGTGGCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCtcatc                                                                                                                                        <  7:392565/148‑1 (MQ=255) cgTAATTTGCGACATGGCAACCGGGTCGATAAACAGTTCCGCAAAGCCCATCACCGCCAGGCCTAATACCATCAGTGGCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCtcat                                                                                                                                         >  7:332348/1‑147 (MQ=255) cgTAATTTGCGACATGGCAACCGGGTCGATAAACAGTTCCGCAAAGCCCATCACCGCCAGGCCTAATACCATCAGTGGCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGAGCTTAAGGTCAGAATGCAGAATCCGGCGCtcat                                                                                                                                         >  1:444009/1‑147 (MQ=255)     aTTTGCGACATGGCAACCGGGTCGATAAACAGTTCCGCAAAGCCCATCACCGCCAGGCCTAATACCATCAGTGGCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCaggc                                                                                                                                    >  1:370057/1‑148 (MQ=255)       ttGCGACATGGCAACCGGGTCGATAAACAGTTCCGCAAAGCCCATCACCTCCATGCCTAATACCATCAGTGGCAGAGAAGAGTGCCCATACATTGAGGACCAGCGGGCGCTTAAGGGCAGAATGCAGAATCCGGCGATCATCaggac                                                                                                                                   >  6:193771/1‑145 (MQ=255)        tGCTACATGGCAACCGGGTCGATAAACAGGTCCGCAAAGCCCATCACCGCCAGGCCTAATCCCATCAGTGGCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATCCAGAATCCGGCGCTCATCaggcca                                                                                                                                  <  7:154897/147‑1 (MQ=255)          cGACATGGCAACCGGGTCGATAAACAGTTCCGCAAAGCCCATCACCGCCAGGCCTAATACCATCAGTGGCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCaggccaag                                                                                                                                >  2:80207/1‑147 (MQ=255)          cGACATGGCAACCGGGTCGATAAACAGTTCCGCAAAGCCCATCACCGCCAGGCCTAATACCATCAGTGGCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCaggccaag                                                                                                                                <  1:80207/147‑1 (MQ=255)               tGGCAACCGGGTCGATAAACAGTTCCGCAAAGCCCATCACCGCCAGGCCTAATACCATCAGTGGCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGgccaaggccaag                                                                                                                          <  8:332348/148‑1 (MQ=255)                ggCAACCGGGTCGATAAACAGTTCCGCAAAGCCCATCACCGCCAGGCCTAATACCATCAGTGGCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCAGGCGCTCATCAggccaaggccaa                                                                                                                           >  7:397827/1‑146 (MQ=255)                      cGGGTCGATAAACAGTTCCGCAAAGCCCATCACCGCCAGGCCTAATACCCTCTGTGGCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAAt                                                                                                                   <  2:421348/148‑1 (MQ=255)                          tCGATAAACAGTTCCGCAAAGCCCATCACCGCCAGGCCTAATACCATCAGTGGCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTcc                                                                                                               >  3:383324/1‑148 (MQ=255)                          tCGATAAACAGTTCCGCAAAGCCCATCACCGCCAGGCCTAATACCATCAGTGGCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTcc                                                                                                               <  8:205623/148‑1 (MQ=255)                                        cGCAAAGCCCATCACCGCCCGGCCTTTTACCATCAGTGGCAGAGTAGAGTGACCATAAATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTGATCAGTCCAAGGCCAAGAGCAAATTTCCCCCAGAGgcgc                                                                                                    <  6:387654/145‑1 (MQ=255)                                              gCCCATCACCGCCAGGCCTAATACCATCAGTGGCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGa                                                                                           <  8:174360/148‑1 (MQ=255)                                                  aTCACCGCCAGGCCTAATACCATCAGTGGCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTAc                                                                                       >  1:23382/1‑148 (MQ=255)                                                   tCACCGCCAGGCCTAATACCATCAGTGGCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTAcc                                                                                      >  5:47955/1‑148 (MQ=255)                                                   tCACCGCCAGGCCTAATACCATCAGTGGCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTAcc                                                                                      >  5:380187/1‑148 (MQ=255)                                                   tCACCGCCAGGCCTAATACCATCAGTGGCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTAcc                                                                                      >  4:15617/1‑148 (MQ=255)                                                   tCACCGCCAGGCCTAATACCATCAGTGGCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAACAGCAAATTTCCCCCAGATGCGCACGGTACGATTAcc                                                                                      >  7:133198/1‑148 (MQ=255)                                                     aCCGCCAGGCCTAATACCATCAGTGGCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCcg                                                                                    <  2:370057/148‑1 (MQ=255)                                                      ccGCCAGGCCTAATACCATCAGTGGCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCcgc                                                                                   >  3:166279/1‑148 (MQ=255)                                                      ccGCCAGGCCTAATACCATCAGTGGCAGAGAAGAGTGACCATACATTCCGGACCAGCGGGCGCATAAGGTCAGAATGCACAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCAAGATGCGCACGCTACGATTACCcgc                                                                                   >  3:269209/1‑148 (MQ=255)                                                           aGGCCTAATACCATCAGTGGCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCGaca                                                                               <  1:134195/147‑1 (MQ=255)                                                           aGGCCTAATACCATCAGAGGCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCGaca                                                                               >  2:134195/1‑147 (MQ=255)                                                               cTAATACCATCAGTGGCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCCAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCGACActctc                                                                          <  8:156439/148‑1 (MQ=255)                                                                 aaTACCATCAGTGGCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTcc                                                                                                               >  6:335404/1‑109 (MQ=255)                                                                 aaTACCATCAGTGGCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTcc                                                                                                               <  5:335404/109‑1 (MQ=255)                                                                 aaTACCATCAGTGGCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCAGCGCtc                                                                                                                                           <  1:9036/81‑1 (MQ=255)                                                                 aaTACCATCAGTGGCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCAGCGCtc                                                                                                                                           >  2:9036/1‑81 (MQ=255)                                                                  aTACCATCAGTGGCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCGACActctct                                                                         <  8:72350/146‑1 (MQ=255)                                                                   tACCCTCAGTCGCAGAGAAGCGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCGACACTCTCtttt                                                                      <  5:193771/148‑1 (MQ=255)                                                                     ccATCAGTGGCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCGACACTCTCTTTTac                                                                    <  4:166279/148‑1 (MQ=255)                                                                        tCAGTGGCAGAGAAGAGTGACCATTCATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCAGCGCTCATCaggcc                                                                                                                                   >  6:185807/1‑82 (MQ=255)                                                                        tCAGTGGCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCAGCGCTCATCaggcc                                                                                                                                   <  5:185807/82‑1 (MQ=255)                                                                         cAGTGGCAGAGAAGAGTGACCATACATTGCGGACCATCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCGACACTCTCTTTTaccacc                                                                >  1:15524/1‑148 (MQ=255)                                                                         cAGTGGCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCcg                                                                                    <  8:110422/128‑1 (MQ=255)                                                                         cAGTGGCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCcg                                                                                    >  7:110422/1‑128 (MQ=255)                                                                            tGGCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCGACACTCTCTTTTAccacccac                                                             >  1:126958/1‑148 (MQ=255)                                                                             ggCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCACGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCGACACTCTCTGTTACCACCCACg                                                            >  6:335353/1‑148 (MQ=255)                                                                             ggCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAg                                                                                                          >  5:31502/1‑102 (MQ=255)                                                                             ggCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAg                                                                                                          <  6:31502/102‑1 (MQ=255)                                                                                agagAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCGACACTCTCTTTTACCACCCACGCCa                                                         <  2:23382/148‑1 (MQ=255)                                                                                  agaagaGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCaaa                                                                                                                    >  4:244068/1‑87 (MQ=255)                                                                                  agaagaGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCaaa                                                                                                                    <  3:244068/87‑1 (MQ=255)                                                                                   gaagaGTGACCATACATTGCGGACCATCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCGACACTCTCTTTTACCACCCACGCCAGGa                                                      <  2:15524/148‑1 (MQ=255)                                                                                   gaagaGTGACCATACATTGCGGACCAGCGTGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCGACACTCTCTTTTACCACCCACGCCAGGa                                                      <  3:247361/148‑1 (MQ=255)                                                                                   gaagaGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCGACACTCTCtt                                                                        >  7:162166/1‑130 (MQ=255)                                                                                   gaagaGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCGACACTCTCtt                                                                        <  8:162166/130‑1 (MQ=255)                                                                                    aagaGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCGACACTCTCTTTTACCACCCACGCCAGGaa                                                     >  8:182135/1‑148 (MQ=255)                                                                                       agTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCGACACTCTCTTTTACCACCCACGCCAGGAacac                                                  >  7:61765/1‑148 (MQ=255)                                                                                        gTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGTCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCGACACTCTCTTTTACCACCCACGCCAGGAacac                                                  <  2:126958/147‑1 (MQ=255)                                                                                             cATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCGACACTCTCTTTTACCACCCACGCCAGGAACACACCGca                                            <  4:383324/148‑1 (MQ=255)                                                                                                aCATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGAAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCGAAACTCTCTTTTACCACCCACGCCAGGAACACACCGCACAg                                         >  8:422145/1‑148 (MQ=255)                                                                                                 cATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCGACACTCTCTTTTACCACCCACGCCAGGAACACACCGCACAg                                         >  7:238047/1‑147 (MQ=255)                                                                                                    tGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGt                                                                                               >  2:230283/1‑90 (MQ=255)                                                                                                    tGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGt                                                                                               <  1:230283/90‑1 (MQ=255)                                                                                                     gCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCGACACTCTCTTTTACCACCCACGCCAGGAACACACCGCACAGCATa                                     <  6:380187/147‑1 (MQ=255)                                                                                                           cAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACTGTACGATTACCCGCGACACTCTCTTTTACCACCCACGCCAGGAACACACCGCACAGCATAACTGCGa                              <  7:422145/148‑1 (MQ=255)                                                                                                            aGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCGACACTCTCTTTTACCACCCACGCCAGGAACACACCGCACAGCATAACTGCGa                              <  5:385946/147‑1 (MQ=255)                                                                                                            aGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCGACACTCTCTTTTACCACCCACGCCAGGAACACACCGCACAGCATAACTGCGa                              >  6:385946/1‑147 (MQ=255)                                                                                                               gggCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCGACACTCTCttt                                                                       >  3:254692/1‑103 (MQ=255)                                                                                                               gggCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCGACACTCTCttt                                                                       <  4:254692/103‑1 (MQ=255)                                                                                                                      tAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCGACACTCTCTTTTACCACCCCCGCCCGGCACACACCGCACAGCATAACTGCGAAGGCCTTAATc                   >  8:379284/1‑148 (MQ=255)                                                                                                                                 tGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCGACACTCTCTTTTACCACCCACGCCAGGAACACACCGCACAGCATAACTGCGAAGGCATTAATCGACTGGAACAt        <  7:182135/148‑1 (MQ=255)                                                                                                                                     gAATCCGGCGCTCATCAGGCCACGGCCAAGAGAAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCCACACTCGCTTTTACCACCCACGCCAGTAACACACCGCACAGCATAACTGCGAAGGCATTAATCGACTGGAACATcgcg    >  2:167427/1‑148 (MQ=255)                                                                                                                                     gAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCGACACTCTCTTTTACCACCCACGCCAGGAACACACCGCACAGCATAACTGCGAAGGCATTAATCGACTGGAACATcgcg    >  8:250313/1‑148 (MQ=255)                                                                                                                                       aTCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCGACACTCTCTTTTACCACCCACGCCAGGAACACACCGCACAGCATAACTGCGa                              >  1:301375/1‑120 (MQ=255)                                                                                                                                       aTCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCGACACTCTCTTTTACCACCCACGCCAGGAACACACCGCACAGCATAACTGCGa                              <  2:301375/120‑1 (MQ=255)                                                                                                                                       aTCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCGACACTCTCTTTTACCACCCACGCCAGGAACACACCGCACAGCATAACTGCGAAGGCATTAATCGACTGGAACATCGCGGt  >  3:150985/1‑148 (MQ=255)                                                                                                                                         ccGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTcc                                                                                                               >  7:383910/1‑37 (MQ=255)                                                                                                                                         ccGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTcc                                                                                                               <  8:383910/37‑1 (MQ=255)                                                                                                                                           |                                                                                                                                                CGTAATTTGCGACATGGCAACCGGGTCGATAAACAGTTCCGCAAAGCCCATCACCGCCAGGCCTAATACCATCAGTGGCAGAGAAGAGTGACCATACATTGCGGACCAGCGGGCGCTTAAGGTCAGAATGCAGAATCCGGCGCTCATCAGGCCAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCGACACTCTCTTTTACCACCCACGCCAGGAACACACCGCACAGCATAACTGCGAAGGCATTAATCGACTGGAACATCGCGGT  >  NC_000913/741348‑741629

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