Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,390,661	T→A	6.2%	K456N (AAA→AAT)	tldD ←	putative peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,390,661	0	T	A	6.2%	184.8 / 2.7	74	K456N (AAA→AAT)	tldD	putative peptidase
Reads supporting (aligned to +/- strand):  ref base T (36/33);  new base A (3/2);  total (39/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.60e-01

CCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTTGGCTGGCCCACGCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGCCGTTATCCAGTTTCAGGTCGTTGCCAACCATCGAAATCTGCTGCATGGTTTCGATACCGGAACCAATCAACGTTGCGCCTTTCACCGGCTTCGTTACTTTACCGTTTTCAATCAGATATGCTTCTGAA  >  NC_000913/3390513‑3390808                                                                                                                                                     |                                                                                                                                                    ccccAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTTGGCTGGCCCACGCCAACCGGCAAACTTTGCCcttct                                                                                                                                                     <  1:575637/149‑1 (MQ=255)  cccAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTTTCGACTTTCAACGTTGGCTGGCCCACGCCAACCGGCAAACTTTGCCcttctt                                                                                                                                                    >  2:680043/1‑149 (MQ=255)  cccAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTTGGCTGGCCCACGCCAACCGCCAAACTTTTCCcttctt                                                                                                                                                    >  1:275368/1‑149 (MQ=255)        ggggCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTTGGCTGTCCCACGCCAACCGGCAAACTTTGCCCTTCTTTACCGc                                                                                                                                              <  2:734309/149‑1 (MQ=255)           gCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTTGGCTGGCCCACGCCAACCGGCAAACTTTGCCCTTCTTTACCGCAgc                                                                                                                                           >  1:403158/1‑148 (MQ=255)                      gTTCGTGCACGTAGAAAGATTAGTTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTTGGCTGGCCCACGCCAACCGTCAAACTTTGCCCTTCTTTACCGAAGACACCCccgccg                                                                                                                                >  2:55805/1‑149 (MQ=255)                      gTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTTGGCTGGCCCACGCCAACCGGCAAACTTTGCCCTTCTTTACCGCCGACACCCACGCCg                                                                                                                                >  2:299228/1‑149 (MQ=255)                      gTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTTGGCTGGCCCACGCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGCCg                                                                                                                                >  1:277577/1‑149 (MQ=255)                      gTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTTGGCTGGCCCACGCCAACCGGCAAACTTTGCCCTTCATTTCCGCAGACACCCACGCCg                                                                                                                                >  1:1038540/1‑149 (MQ=255)                      gTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTTGGCTGGCCCACGCCAACCGGCAAACTTTGCCCTTCATTACCGCAGACACCCACGCCg                                                                                                                                >  1:1042631/1‑149 (MQ=255)                          gTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTCCCGCCAACAGTCAGGTTATCGACTGTCAACGTTGGCTGGCCCAAGCCAACCGGCAAACTTTGCCCTTCTTTACAGCAGCCACCCACTCCGGTAt                                                                                                                            >  2:263714/1‑149 (MQ=255)                          gTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTAGGCTGGCCCACGCCAACCGGCAAACTGTGCCCTTCTTTACCGCAGCCACCCACGCCGttca                                                                                                                            >  2:259418/1‑147 (MQ=255)                           tGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTTGGCTGGCCCACGCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGCCGTTATc                                                                                                                           <  2:114321/149‑1 (MQ=255)                            gCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTTGGCTGGCCCACGCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGCCGTTATcc                                                                                                                          <  2:687426/149‑1 (MQ=255)                              aCGTAGAAAGATTAATTATCCTTCTTAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTTGGCTGGCCCACGCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACAACCACGCCGTTATCCAg                                                                                                                        <  1:294642/149‑1 (MQ=255)                                gTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTTGGCTGGCCCACGCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGCCGTTATCCAGtt                                                                                                                      <  2:158053/149‑1 (MQ=255)                                gTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTTGGCTGGCCCACGCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGCCGTTATCCAGtt                                                                                                                      <  2:735416/149‑1 (MQ=255)                                  aGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCACCGTTGGCTGGCCCACGCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACCCCGTTATCCAGttgc                                                                                                                    >  2:1024141/1‑147 (MQ=255)                                  aGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTTGGCTGGCCCACGCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGCCGTTATCCAGTTTc                                                                                                                    >  1:411788/1‑149 (MQ=255)                                  aGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTTGGCTGGCCCACGCCAACCGGCAAACTTTGCCCTTCTTAACCCCAGACACCCACGCCGTTAACCAGTTTc                                                                                                                    >  2:254301/1‑149 (MQ=255)                                    aaaGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTTGGCTGGCCCACGCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGCCGTTATCCAGTTTCAg                                                                                                                  <  1:442767/149‑1 (MQ=255)                                        attaattaTCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTTGGCTGGCCCACGCCATCCGGCAATCTTTGCCCTTCTTTACCGCAGACACCCACGCCGTTATCCAGTTTCAGGTCg                                                                                                              >  1:777954/1‑149 (MQ=255)                                        attaattaTCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTACAACGATGGCTGCCCCACGCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGCCACCCACGCCGTTTTCCCGTTGCGGGTCg                                                                                                              >  2:209136/1‑149 (MQ=255)                                        attaattaTCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGACAACATTCAGGTTATCGACTTTCAACTTTTGCTGGCCCACGCCAACCGGCAAACTTTGCCCTTATTTACCGCAGACACCCACGCCGTTATCCAGTTTCAGGTCg                                                                                                              >  1:27608/1‑149 (MQ=255)                                            attaTCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTTGGCTGGCCCACGCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGCCCCCCACGCCGTGATCCAGTTTCAGGTCGTTGc                                                                                                          >  2:967347/1‑149 (MQ=255)                                             ttaTCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTTGGCTGGCCCACGCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGCCGTTATCCAGTTTCAGGTCGTTGcc                                                                                                         <  1:70694/149‑1 (MQ=255)                                                tCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTTGGCTGGCCCACGCCAACCGGCATACTTTGCACTTCTTTACCGCAGACCCCCACGCGGTTATCCAGTTTCAGGTCGTTGCCAAc                                                                                                      >  2:1133303/1‑149 (MQ=255)                                                     cTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTTGGCTGGCCCACGCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGCCGTTATCCAGTTTCAGGTCGTTGCCAACCATCg                                                                                                 >  1:843968/1‑149 (MQ=255)                                                      tGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCTACTTTCAACGTTGGCTGGCCCACGCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACCCCCACGCCGTTATCCCGTTTCAGGTCGTTGCCAACCTTCGa                                                                                                >  1:537372/1‑149 (MQ=255)                                                      tGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTTGGCTGGCCCACGCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGCCGTTATCCAGTTTCAGGTCGTTGCCAACCATCGa                                                                                                >  1:493247/1‑149 (MQ=255)                                                       gAAAATAGTGAAATTATTACGCAGGACCGCCCACAGTCAGCTTATCGACTTTCAACGTTGGCTGGCCCACGCCAACCGGCAACATTTGCACATCTTTACCGCAGACCCCCCCGCCGTTATCCAGTTTCAGCTCGTGGCCACCCATCGaa                                                                                               >  1:130611/1‑149 (MQ=255)                                                         aaaTAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTTGGCTGGCCCACGCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGCCGTTATCCAGTTTCAGGTCGTTGCCAACCATCGAAAt                                                                                             <  1:1055242/149‑1 (MQ=255)                                                          actAGTGAAATTAATCCGCAGGACCGCCAACAGTCAGGTTATCGACTATCAACGTTGGCTGGCCCACGCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACTCCCACGCCGTTATCCAGTTTCAGGTCGTTGCCAACCATCGAAATc                                                                                            <  2:504252/147‑1 (MQ=255)                                                           aTAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGAATTTCAACGTTGGCTGGCCCACGCCAACCGGCAAACTTAGCCCTTCTTTACCTCAGACACCCACGCCGTTATCCAGTTTCAGGTCGTTGCCAACCATCGAAATct                                                                                           <  1:773232/149‑1 (MQ=255)                                                               tGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTTGGCTGGCCCACGCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGCCGTTATCCAGTTTCAGGTCGTTGCCAACCATCGAAATctgctg                                                                                       <  2:411788/149‑1 (MQ=255)                                                                 aaaTTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTTGGCTGGCCCACGCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGCCGTTATCCAGTTTCAGGTCGTTGCCACCCATCGAAATCTGCTGCa                                                                                     >  2:636482/1‑149 (MQ=255)                                                                 aaaTTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTTGGCTGGCCCACGCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGCCGTTATCCAGTTTCAGGTCGTTGCCAACCATCGAAATCTGCTGCa                                                                                     >  1:846119/1‑149 (MQ=255)                                                                      attaCGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTTGGCTGGCCCACGCCAACCGGCAAACTTTGCCCTTCATTACCGCAGACACCCACGCCGTTATCCAGTTTCAGGTCGTTGCCAACCATCGAAATCTGCTGCATGGtt                                                                                <  2:1042631/149‑1 (MQ=255)                                                                      attaCGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTTGGCTGGCCCACGCCAACCGGCAAACTTTGCCCTTCATTACCGCAGACACCCACGCCGTTATCCAGTTTCAGGTCGTTGCCAACCATCGAAATCTGCTGCATGGtt                                                                                <  2:1038540/149‑1 (MQ=255)                                                                       ttaCGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTTGGCTGGCCCACGCCACCCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGCCGTTATCCAGTTTCAGGTCGTTGCCAACCATCGAAATCTGCTGCATGGttt                                                                               >  1:157803/1‑149 (MQ=255)                                                                        taCGCCGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTTGGCTGGCCCACGCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGCCGTTATCCAGTTTCAGGTCGTTGCCAACCATCGAAATCTGCTGCATGGTTTc                                                                              <  1:330987/149‑1 (MQ=255)                                                                          cGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTTGGCTGGCCCACGCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGCCGTTATCCAGTTTCGGGTCGTTGCCAACCATCGAAATCTGCTGCATGGTTTCGa                                                                            <  2:268450/149‑1 (MQ=255)                                                                          cGCAGTACCGCCAACAGTCAGGATATCGACATTCCACGTTGGCTGGCCCACGCCAACCGGCAAACTTGGCCCTTCTTTACCGCAGACCCCCACGCCGTTATCCAGTGTCACGTCGTTGCCAACCATCGAAATCTGCTGCATGGTTTCGa                                                                            >  2:454144/1‑149 (MQ=255)                                                                              gTACCGCCAACAGTCAGGTTATCGACTTTCAACGTTGGCTGGCCCACGCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGCCGTTATCCAGTTTCAGGTCGTTGCCAACCATCGAAATCTGCTGCATGGTTTCGATAcc                                                                        <  2:290548/149‑1 (MQ=255)                                                                              gTACCGCCAACAGTCAGGTGATCGACTTTCAACGTTGGCTGGCCCACGCCCACCGGCAAACTTTGCCCAGCTTTTCCGCAGACACCCACGCCGGGATCCAGATTCAGGTCGTTGCGAGCGCTGGATCTCGGCTGCCTGGGTTCGATAcc                                                                        >  2:814475/1‑149 (MQ=255)                                                                               tACCGCCAACAGTCAGGTTATCGACTTTCAACGTTGGCTGGCCCCCGCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGCCGTTATCCAGTTTCAGGTCGTTGCCAACCATCGAAATCTGCTGCATGGTTTCGATACCg                                                                       >  1:400065/1‑149 (MQ=255)                                                                                  cGCCAACAGTCAGGTTATCGACTTTCAACGTTGGCTGGCCCACGCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGCCGTTATCCAGTTTCAGGTCGTTGCCAACCATCGAAATCTGCTGCATGGTTTCGATACCGGaa                                                                    >  2:867446/1‑149 (MQ=255)                                                                                       aCAGTCAGGTTATCGACTTTCAACGTTGGCTGGCCCACGCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGCCGTTATCCAGTTTCAGGTCGTTGCCAACCTTCGAAATCTGCTGCATGGTTTCGATACCGGAACCAAt                                                               <  2:171566/149‑1 (MQ=255)                                                                                       aCAGTCAGGTTATCGACTTTCAACGTTGGCTGGCCCACGCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGCCGTTATCCAGTTTCAGGTCGTTGCCAACCATCGAAATCTGCTGCATGGTTTCGATACCGGAACCAAt                                                               <  2:640532/149‑1 (MQ=255)                                                                                        cAGTCAGGTTATCGACTTTCAACGTTGGCTGGCCCACGCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGCCGTTATCCAGTTTCAGGTCGTTGCCAACCATCGAAATCTGCTGCATGGTTTCGATACCGGAACCAATc                                                              <  1:316474/149‑1 (MQ=255)                                                                                          ggccggTTATCGACTTTCAACGTAGGCTGGCCCTAGGCAACCGGCAAACTTTGCCCTGCTTTACCGCACCCACCCCCGCCTGTATCCACTTTCAGGGCGTTGCCAACCATCGAAATCTGCTGCATGGTTTCGATACCGGAACCAATCaa                                                            <  1:856427/145‑1 (MQ=255)                                                                                              ggTTATCGACTTTCAACGTTGGCTGGCCCACGCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGCCGTTATCCAGTTTCAGGTCGTTGCCAACCATCGAAATCTGCTGCATGGTTTCGATACCGGAACCAATCAACGtt                                                        >  1:186449/1‑149 (MQ=255)                                                                                               gTTATCGACTTTCAACGTTGGCTGGCCCACGCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGCCGTTATCCAGTTTCAGGTCGTTGCCAACCATCGAAATCTGCTGCATGGTTTCGATACCGGAACCAATCAACGTTg                                                       <  2:323640/149‑1 (MQ=255)                                                                                                  aTCGACTTTCAACGTTGGCTGGCCCACGCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGCCGTTATCCAGTTTCAGGTCGTGGCCAACCCTCGAAATCTGGTGAATTGTTTCGATACCGGAACCAATCAACGTTgcgc                                                    >  2:874194/1‑149 (MQ=255)                                                                                                      aCTTTCAACGTTGGCTGGCCCACTCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGCCGTTATCCAGTTTCAGGTCGTTGCCAACCATCGAAATCTGCTGCATGGTTTCGATACCGGAACCAATCAACGTTGCGCCttt                                                <  1:917387/149‑1 (MQ=255)                                                                                                      aCTTTCAACGTTGGCTGGCCCACGCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGCCGTTATCCAGTTTCAGGTCGTTGCCAACCATCGAAATCTGCTGCATGGTTTCGATACCGGAACCAATCAACGTTGCGCCttt                                                >  2:476946/1‑149 (MQ=255)                                                                                                       cTTTCAACGTTGGCTGGCCCACGCCAACCGGCAAACTTTGACCTTCTTTACCGCAGACACCCACGCCGTTATCCAGTTTCAGGTCGTTGCCAACCATCGAAATCTGCTGCATGGTTTCGATACCGGAACCAATCAACGTTGCGCCTTTc                                               <  2:59507/149‑1 (MQ=255)                                                                                                        tttCAACGTTGGCTGGCCCACGCCAACCGGCAAACTTTGCCCATCTTTACCGCAGACACCCAAGCCGTTATCCAGTTTCAGGTCGTTGCCAACCAACGAAATCTGCTGCATGGTTTCGCTAACGGAACCAATCAACGTTCCGCCTTTCa                                              >  2:826104/1‑149 (MQ=255)                                                                                                        tttCAACGTTGGCTGGCCCACGACAACCGGCAAACTTTGCCCATAATTACCGCCGACCCCCACGCCGTTATCCAGTTTCAGGACGTTGCCCACCATCGAAATCTGCTGCATGGTTTCGATACCGGACCCAATCAACCTTGCGCCTTTcc                                              >  1:256624/1‑148 (MQ=255)                                                                                                         ttCAACGTTGGCTGGCCCACGCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGCCGTTATCCAGTTTCAGGTCGTTGCCAACCATCGAAATCTGCTGCATGGTTTCGATACCGGAACCAATCAACGTTGCGCCTTTCAc                                             <  2:403158/149‑1 (MQ=255)                                                                                                             aCGTTGGCTGGCCCACGCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCCAGCCGTTATCCAGTTTCAGGTCGTTGCCAACCATCGAAATCTGCTGCATGGTTTCGATACCGGAACCAATCAACGTTGCGCCTTTCACCGGc                                         >  1:390300/1‑149 (MQ=255)                                                                                                                ttGGCTGGCCCACGCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGCCGTTATCCAGTTTCAGGTCGTTGCCAACCATCGAAATCTGCTGCATGGTTTCGATACCGGAACCAATCAACGTTGCGCCTTTCACCGGCTTc                                      <  1:1087258/149‑1 (MQ=255)                                                                                                                        cccACGCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGCCGTTATCCAGTTTCAGGTCGTTGCCAACCATCGAAATCTGCTGCATGGTTTCGATACCGGAACCAATCAACGTTGCGCCTTTCACCGGCTTCGTTACttt                              <  1:974901/149‑1 (MQ=255)                                                                                                                        cccACGCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGCCGTTATCCAGTTTCAGGTAGTTGCCAACCATCGAAATCTGCTGCATGGTTTCGATCCCGGAACCAAACAGCGTTGCGCCTTTCCCCGGCGTCGTTCCttt                              >  2:267490/1‑149 (MQ=255)                                                                                                                          cACGCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGCCGTTATCCAGTTTCAGGTCGTTGCCAACCATCGAAATCTGCTGCATGGTTTCGATACCGGAACCAATCAACGTTGCGCCTTTCACCGGCTTCGTTACTTTAc                            >  1:136908/1‑149 (MQ=255)                                                                                                                          cACGCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGCCGTTATCCAGTTTCACGTCGTTGCCAACCATCGAAATCTGCTGCAAGTTTTCGATACCGGAACCAATCAACGTTGAGCCTTTCGCCGACTACgttacgttac                            >  1:259190/1‑149 (MQ=255)                                                                                                                           aCGCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGCCGTTATCCAGTTTCAGGTCGTTGCCAACCATCGAAATCTGCTGCATGGTTTCGATACCGGAACCAATCAACGTTGCGCCTTTCACCGGCTTCGTTACTTTAcg                           <  1:252672/149‑2 (MQ=255)                                                                                                                           aCGCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGCCGTTATCCAGTTTCAGGGCGTTGCCAACCATCGAAATCTGCTGCATGGTTTCGATACCGGAACCAATCAACGTTGCGCCTTTCACCGGCTTCGTTACTTTAcc                           <  2:132157/149‑1 (MQ=255)                                                                                                                              ccAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGCCGTTATCCAGTTTCAGGTCGTTGCCAACCATCGAAATCTGCTGCATGGTTTCGATACCGGAACCAATCAACGTTGCGCCTTTCACCGGCTTCGTTACTTTACCGtt                        <  2:601758/149‑1 (MQ=255)                                                                                                                               cAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGCCGTTATCCAGTTTCAGGTGGTTGCCCACCATCGAAAACTGCTGCATGGTTTCGATACCGGATCCAATCTACGTTGCGCCTTTCACAGGCTTCGTTACTTTAACGttt                       >  2:70520/1‑149 (MQ=255)                                                                                                                               cAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGCCGTTATCCAGTTTCAGGTCGTTGCCAACCATCGAAATCTGCTGCATGGTTTCGATACCGGAACCAATCAACGTTGCGCCTTTCACCGGCTTCGTTACTTTACCGttt                       >  2:981577/1‑149 (MQ=255)                                                                                                                                    ggCTAACTTTGCACATCTTTATCTCAGACACCAACGCCGTTCTCCAGTTTAACGTCGTTGCCAACCATCTATATCTGCTGCATGTTTTCGATACCGGAACCAATCAACGTTGCGCCTTTCACCGGCTTCGTTACTTTACCGTTTTCAAt                  <  1:812093/149‑1 (MQ=255)                                                                                                                                    ggCAAACTTTGCCCTTCTTTACCGCAGACACCCACGCCGTTATCCAGTTTCAGGTCGTAGCCAACCATCGAAATCTGCTGCATGGTTTCGATACCGGAACCAATCAACGTTGCGCCTTTCACCGGCTTCGTTACTTTACCGTTTTCAAt                  <  2:564273/149‑1 (MQ=255)                                                                                                                                    ggCAAACTTGGCCCGTCTTTACCGCAGACACCCACGCCGTTATCCAGTTTCAGGTCGTTGCCAACCATCGAAATCTGCTGCATGGGTTCGATACCGGAACCAATCAACGTTGCGCCTTTCACCGGCTTCGTTACTTTACCGTTTTCAAt                  <  2:935091/149‑1 (MQ=255)                                                                                                                                     gCAAACTTTGCCCTTCTTTACCGCAGACACCCACGCCGTTATCCAGTTTCAGGTCGTTGCCAACCATCGAAATCTGCTGCATGGTTTCGATACCGGAACCAATCAACGTTGCGCCTTTCACCGGCTTCGTTACTTTACCGTTTTCAATc                 <  2:575023/149‑1 (MQ=255)                                                                                                                                        aaCTTTGCCCTTCTTTACCGCAGACACCCACGCCGTTATCCAGTTTCAGGTCGTTGCCAACCATCGAAATCTGCTGCATGGTTTCGATACCGGAACCAATCAACGTTGCGCCTTTCACCTGCTTCGTTACTTTACCGTTTTCAATCAGa              <  2:462772/149‑1 (MQ=255)                                                                                                                                              gCCCTTCTTTACCGCAGACACCCACGCCGTTTTCCAGTTTCAGGTCGTTGCCAACCATCGAAATCTGCTGCATGGTTGCGATACCGGAACCAATCAACGTTGCGCCTTTCACCGGCTTCGTTACTTTACCGTTATCAATCAGATATGCt        >  2:796226/1‑149 (MQ=255)                                                                                                                                                    cttTACCGCAGACACCCACGCCGTTATCCAGTTTCAGGTCGTTGCCAACCATCGAAATCTGCTGCATGGTTTCGATACCGGAACCAATCAACGTTGCGCCTTTCACCGGCTTCGTTACTTTACCGTTTTCAATCAGATATGCTTCTGaa  >  1:1135987/1‑149 (MQ=255)                                                                                                                                                    cttTACCGCAGACACCCACGCCGTTATCCAGTTTCAGGTCGTTGCCAACCATCGAAATCTGCTGCATGGTTTCGATACCGGAACCAATCAACGTTGCGCCTTTCACCGGCTTCGTTACTTTACCGTTTTCAATCAGATATGCTTCTGaa  <  2:186449/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    CCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTTGGCTGGCCCACGCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGCCGTTATCCAGTTTCAGGTCGTTGCCAACCATCGAAATCTGCTGCATGGTTTCGATACCGGAACCAATCAACGTTGCGCCTTTCACCGGCTTCGTTACTTTACCGTTTTCAATCAGATATGCTTCTGAA  >  NC_000913/3390513‑3390808

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: ‑