Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	440,662	C→A	6.2%	S147I (AGC→ATC) ‡	ispA ←	geranyltranstransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	440,662	0	C	A	6.2%	287.4 / 7.2	81	S147I (AGC→ATC) ‡	ispA	geranyltranstransferase
Reads supporting (aligned to +/- strand):  ref base C (35/41);  new base A (3/2);  total (38/43)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

GCACCAAGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTCACATGGCAGGTTGGC  >  NC_000913/440521‑440808                                                                                                                                              |                                                                                                                                                   gCACCAAGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAg                                                                                                                                             <  3:231454/149‑1 (MQ=255)   aCCAAGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGtt                                                                                                                                           >  1:297625/1‑149 (MQ=255)   aCCAAGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGAGCGCCAGtg                                                                                                                                           <  3:659892/149‑2 (MQ=255)   aCCAAGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGAGCGCCAGtg                                                                                                                                           >  4:659892/1‑148 (MQ=255)    ccAAGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGAGCGCCAGtgg                                                                                                                                          <  1:664398/149‑3 (MQ=255)     cAAGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCa                                                                                                                                         >  3:585939/1‑149 (MQ=255)     cAAGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCa                                                                                                                                         <  4:23034/149‑1 (MQ=255)     cAAGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCa                                                                                                                                         >  2:558914/1‑149 (MQ=255)     cAAGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCa                                                                                                                                         >  1:455689/1‑149 (MQ=255)     cAAGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCa                                                                                                                                         <  2:610462/149‑1 (MQ=255)       aGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGa                                                                                                                                       >  2:463375/1‑149 (MQ=255)       aGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGa                                                                                                                                       <  3:345556/149‑1 (MQ=255)            aaCGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAatca                                                                                                                                  <  2:223891/149‑1 (MQ=255)             aCGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAatcat                                                                                                                                 <  1:12590/149‑1 (MQ=255)             aCGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAatcat                                                                                                                                 <  3:515550/149‑1 (MQ=255)              cggcggCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAatcatc                                                                                                                                <  2:297625/149‑1 (MQ=255)                gcggcgCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGa                                                                                                                              <  1:559299/149‑1 (MQ=255)                gcggcgCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGa                                                                                                                              >  2:213989/1‑149 (MQ=255)                  ggcgCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAatcatc                                                                                                                                <  1:181220/145‑1 (MQ=255)                  ggcgCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAatcatc                                                                                                                                >  2:181220/1‑145 (MQ=255)                  ggcgCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGaaa                                                                                                                            <  1:617073/149‑1 (MQ=255)                  ggcgCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGaaa                                                                                                                            <  4:131295/149‑1 (MQ=255)                  ggcgCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGaaa                                                                                                                            <  2:549113/149‑1 (MQ=255)                        aaTCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGt                                                                                                                      >  1:514453/1‑149 (MQ=255)                        aaTCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGt                                                                                                                      >  3:42969/1‑149 (MQ=255)                                  ccGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGa                                                                                                            <  2:167964/149‑1 (MQ=255)                                   cGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGac                                                                                                           >  3:682438/1‑149 (MQ=255)                                   cGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGac                                                                                                           <  2:232/149‑1 (MQ=255)                                     gTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTc                                                                                                                     >  2:310648/1‑137 (MQ=255)                                     gTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTc                                                                                                                     <  1:310648/137‑1 (MQ=255)                                     gTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGacac                                                                                                         <  2:685917/149‑1 (MQ=255)                                     gTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGacac                                                                                                         >  4:402158/1‑149 (MQ=255)                                         tATGACGATGAATACGCTCAATCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTATTGTCTGACCACCGCACATTTCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTACGACACTTcc                                                                                                     <  3:425960/149‑1 (MQ=255)                                          aTGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAGATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCg                                                                                                    <  4:318064/149‑1 (MQ=255)                                           tGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTc                                                                                                                                          >  4:304774/1‑110 (MQ=255)                                           tGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTc                                                                                                                                          <  3:304774/110‑1 (MQ=255)                                           tGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCgg                                                                                                   >  3:651851/1‑149 (MQ=255)                                                 tGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATc                                                                                             <  3:298849/149‑1 (MQ=255)                                                  gAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCg                                                                                            <  2:57603/149‑1 (MQ=255)                                                      aCGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCAt                                                                                        <  2:588985/149‑1 (MQ=255)                                                          tCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCt                                                                                    >  3:8955/1‑149 (MQ=255)                                                              gcgcGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTaaa                                                                                >  4:263938/1‑149 (MQ=255)                                                                  gTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCg                                                                            <  3:2404/149‑1 (MQ=255)                                                                  gTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCg                                                                            >  4:287156/1‑149 (MQ=255)                                                                   tCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCga                                                                           <  2:410478/149‑1 (MQ=255)                                                                   tCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCga                                                                           <  1:463375/149‑1 (MQ=255)                                                                   tCCAGAGGTACATGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCCGCATCGCTTAAAATCga                                                                           <  4:578424/149‑1 (MQ=255)                                                                              gtgtTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCg                                                                >  3:353662/1‑149 (MQ=255)                                                                                 tttGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGttt                                                             <  3:439405/149‑1 (MQ=255)                                                                                 tttGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGttt                                                             <  4:353662/149‑1 (MQ=255)                                                                                  ttGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAAAACGCTAGCGTTTg                                                            >  3:578495/1‑149 (MQ=255)                                                                                      cTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTaaa                                                        >  4:582845/1‑149 (MQ=255)                                                                                      cTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTaaa                                                        >  4:210821/1‑149 (MQ=255)                                                                                       ttCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAg                                                       >  1:558783/1‑149 (MQ=255)                                                                                                    tCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAatcatc                                                                                                                                >  2:162493/1‑63 (MQ=255)                                                                                                    tCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAatcatc                                                                                                                                <  1:162493/63‑1 (MQ=255)                                                                                                      tAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGaa                                        >  4:265477/1‑149 (MQ=255)                                                                                                             tGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGtt                                 <  2:455689/149‑1 (MQ=255)                                                                                                               accaccGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTg                               >  4:216968/1‑149 (MQ=255)                                                                                                                 caccGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCt                             >  4:19474/1‑149 (MQ=255)                                                                                                                  accGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCtt                            <  1:213989/149‑1 (MQ=255)                                                                                                                     gCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGt                                                                                                                      <  1:312645/56‑1 (MQ=255)                                                                                                                     gCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGt                                                                                                                      >  2:312645/1‑56 (MQ=255)                                                                                                                         aTTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCaaa                     >  1:125942/1‑149 (MQ=255)                                                                                                                         aTTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCaaa                     >  1:293856/1‑149 (MQ=255)                                                                                                                            ccGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACtt                  >  2:71682/1‑149 (MQ=255)                                                                                                                            ccGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACtt                  >  2:440528/1‑149 (MQ=255)                                                                                                                              ggCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTca                >  1:516664/1‑149 (MQ=255)                                                                                                                               gCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTcac               <  4:682438/149‑1 (MQ=255)                                                                                                                               gCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTcac               <  2:558783/149‑1 (MQ=255)                                                                                                                               gCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTcac               >  2:541425/1‑149 (MQ=255)                                                                                                                               gCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTcac               <  3:582845/149‑1 (MQ=255)                                                                                                                                 aaTACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTCACAt             <  4:8955/149‑1 (MQ=255)                                                                                                                                  tctgaaCTGGCGAGCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGCGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTCACATg            >  2:241128/6‑149 (MQ=255)                                                                                                                                  aTACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTCACATg            <  2:125942/149‑1 (MQ=255)                                                                                                                                  aTACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTCACATg            <  2:293856/149‑1 (MQ=255)                                                                                                                                    tgaaCTGGCGAGCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTCACATGGc          >  4:636995/4‑149 (MQ=255)                                                                                                                                     ccACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATc                                      <  1:675038/120‑1 (MQ=255)                                                                                                                                     ccACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATc                                      >  2:675038/1‑120 (MQ=255)                                                                                                                                       aCTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTCACATGGCAgg       <  3:402158/149‑1 (MQ=255)                                                                                                                                        cTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTCACATGGCAGGt      >  3:100029/1‑149 (MQ=255)                                                                                                                                           gcgcTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTCACATGGCAGGTTgg   <  1:71682/149‑1 (MQ=255)                                                                                                                                           gcgcTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTCACATGGCAGGTTgg   <  1:440528/149‑1 (MQ=255)                                                                                                                                            cgcTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTCACATGGCAGGTTGGc  <  4:585939/149‑1 (MQ=255)                                                                                                                                              |                                                                                                                                                   GCACCAAGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTCACATGGCAGGTTGGC  >  NC_000913/440521‑440808

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