Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	220,849	(T)5→4	100%	coding (80/816 nt)	metQ ←	DL‑methionine transporter subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	220,845	0	T	.	54.7%	‑4.5 / 189.8	86	coding (84/816 nt)	metQ	DL‑methionine transporter subunit
Reads supporting (aligned to +/- strand):  ref base T (23/16);  new base . (24/23);  total (47/39)
Fisher's exact test for biased strand distribution p-value = 5.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.16e-01

AACGTCCAGGCCATATTTGTCTTTCGCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATCTTTTTCATCCTGACCGCAGCCTACCAGTGCCAGTGATCCGATCAGGGCTCCCACTGCCGCAAAGGTTTTGAATTTGAACGCCATACCTTATTCCTTCTTC  >  NC_000913/220746‑220944                                                                                                    |                                                                                                    aaCGTCCAGGCCATATTTGTCTTTCGCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATCt                                                                                                     >  2:1403272/1‑100 (MQ=255)     tCCAGGCCATGTTTGTCTTTCGCAACTTTCTGCTCTCCTTCTGCACGCTGTGGTTCGGCACCAATATTCACTCGGTCTTTAATGTGGTTTGGATCttttc                                                                                                 <  2:1903007/100‑2 (MQ=255)      ccAGGCCATATTTGTCTTTCGCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCGGCACAAACAATCACGCCGACTTTAATGTGGTTTGGATCttttt                                                                                                 >  2:50437/1‑99 (MQ=255)       cAGGCCATATTTGTCTTTCGCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATCttttcat                                                                                               >  2:1252742/1‑97 (MQ=255)       cAGGCCATATTTGTCTTTCGCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATCttttcat                                                                                               <  1:846441/100‑4 (MQ=255)       cAGGCCATATTTGTCTTTCGCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATCttttcat                                                                                               >  1:961735/1‑97 (MQ=255)       cAGGCCATATTTGTCTTTCGCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATCTTTTTCa                                                                                               >  1:119174/1‑100 (MQ=255)        aGGCCATATTTGTCTTTCGCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATCttttcatc                                                                                              <  2:1969792/100‑5 (MQ=255)        aGGCCATATTTGTCTTTCGCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATCttttcatc                                                                                              >  1:1589007/1‑96 (MQ=255)        aGGCCATATTTGTCTTTCGCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATCTTTTTCAt                                                                                              >  2:1909246/1‑100 (MQ=255)          gCCATATTTGTCTTTCGTAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATCTTTTTCATcc                                                                                            <  1:630525/100‑1 (MQ=255)          gCCATATTTGTCTTTCGCAACTTTCTGCGCGACTTCGGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATCTTTTTCATcc                                                                                            >  1:1578493/1‑100 (MQ=255)                 ttGTCTTTCGCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATC‑TTTTCATCCTGACCGCa                                                                                    <  1:1019524/100‑1 (MQ=255)                  tGTCTTTCGCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATCTTTTTCATCCTGACCGCa                                                                                    <  1:1862610/100‑1 (MQ=255)                    tCTTTCGCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATC‑TTTTCATCCTGACCGCAGcc                                                                                 <  2:1089424/100‑1 (MQ=255)                     cTTTCGCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATCTTTTTCATCCTGACCGCAGcc                                                                                 <  1:893186/100‑1 (MQ=255)                     cTTTCGCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATC‑TTTTCATCCTGACCGCAGCCt                                                                                >  1:473991/1‑100 (MQ=255)                      tttCGCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATC‑TTTTCATCCTGACCGCAGCCTa                                                                               <  2:1809244/100‑1 (MQ=255)                      tttCGCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATC‑TTTTCATCCTGACCGCAGCCTa                                                                               >  2:2201757/1‑100 (MQ=255)                       ttCGCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATC‑TTTTCATCCTGACCGCAGCCTAc                                                                              <  1:1616775/100‑1 (MQ=255)                         cGCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATC‑TTTTCATCCTGACCGCAGCCTAcca                                                                            >  1:1022859/1‑100 (MQ=255)                          gCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATCTTTTTCATCCTGACCGCAGCCTAcca                                                                            >  1:712687/1‑100 (MQ=255)                          gCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATC‑TTTTCATCCTGACCGCAGCCTAccgg                                                                           >  2:1355244/1‑98 (MQ=255)                          gCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATC‑TTTTCATCCTGACCGCAGCCTAcca                                                                            >  2:35366/1‑99 (MQ=255)                          gCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCGACACAAACAATCCCGCCGACTTTAATGTGGTTTGGATC‑TTTTCATCCTGACCGCAGGCTAccag                                                                           <  2:821184/100‑1 (MQ=255)                           cAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATC‑TTTTCATCCTGACCGCAGCCTAccagt                                                                          >  1:433984/1‑100 (MQ=255)                              cTTTCTGCGCGACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATCTTTTTCATCCTGACCGCAGCCTAccagtgc                                                                        <  2:1108530/100‑1 (MQ=255)                              cTTTCTGCGCGACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATCTTTTTCATCCTGACCGCAGCCTAccagtgc                                                                        >  1:1033276/1‑100 (MQ=255)                              cTTTCTGCGCGACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATC‑TTTTCATCCTGTCCGCAGCCTAccagtgcg                                                                       <  2:452173/100‑2 (MQ=255)                               tttCTGCGCGACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATC‑TTTTCATCCTGACCGCAGCCTAccagtgcca                                                                      <  2:54917/100‑1 (MQ=255)                                ttCTGCGCGACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATCTTTTTCATCCTGACCGCAGCCTAccagtgcca                                                                      >  2:665883/1‑100 (MQ=255)                                  cTGCGCGACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAACGTGGTTTGGATCTTTTTCATCCTGACCGCAGCCTAccagtgccagt                                                                    >  1:1299094/1‑100 (MQ=255)                                   tGCGCGACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATCTTTTTCATCCTGACCGCAGCCTAccagtgccag                                                                     >  1:1749222/1‑98 (MQ=255)                                    gcgcGACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATC‑TTTTCATCCTGACCGCAGCCTACCAGTGCCAGTGAt                                                                 >  1:121457/1‑100 (MQ=255)                                    gcgcGACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTGTAATGTGGTTTGGATC‑TTTTCATCCTGACCGCAGCCTACCAGTGCCAGTGAt                                                                 >  2:876887/1‑100 (MQ=255)                                     cgcgACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATCTTTTTCATCCTGACCGCAGCCTACCAGTGCCAGTGAt                                                                 >  1:479563/1‑100 (MQ=255)                                     cgcgACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATCTTTTTCATCCTGACCGCAGCCTACCAGTGCCAGTGAt                                                                 <  2:712687/100‑1 (MQ=255)                                     cgcgACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATCTTTTTCATCCTGACCGCAGCCTACCAGTGCCAGTGAt                                                                 >  1:1433727/1‑100 (MQ=255)                                      gcgACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATCTTTTTCATCCTGACCGCCGCCTACCAGTGCTAGTGATc                                                                >  2:1276140/1‑100 (MQ=255)                                      gcgACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATC‑TTTTCATCCTGACCGCAGCCTACCAGTGCCAGTGATcc                                                               >  1:153980/1‑100 (MQ=255)                                        gACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATC‑TTTTCATCCTGACCGCAGCCTACCAGTGCCAGTGATCCGa                                                             >  2:1755755/1‑100 (MQ=255)                                          cTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATC‑TTTTCATCCTGACCGCAGCCTACCAGTGCCAGTGATCCGATc                                                           <  2:570110/100‑1 (MQ=255)                                             cTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATC‑TTTTCATCCTGACCGCAGCCTACCAGTGCCAGTGATCCGATCAgg                                                        >  2:179011/1‑100 (MQ=255)                                             cTGCAACCTGCTGTTCAGCACCAACAATCACGCCGACTGTAATGTGGTTTGGATCTTTTTCATCCTGACCGCAGCCTACCAGTGCCAGTGATCCGATCAg                                                         >  2:1376709/1‑100 (MQ=255)                                              tGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATCTTTTTCATCCTGACCGCAGCCTACCAGTGCCAGTGATCCGATCAgg                                                        <  2:1678290/100‑1 (MQ=255)                                              tGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATC‑TTTTCATCCTGACCGCAGCCTACCAGTGCCAGTGATCCGATCAggg                                                       <  1:1058245/100‑1 (MQ=255)                                              tGCAACCTGCTGTTCGGCACCAACAATCACGCCGACGTTAATGTGCTTTGGATC‑TTTTAATCCACACCGCCGCCTACCAGTGCCAGTGATCCGATc                                                           >  2:1714261/1‑96 (MQ=255)                                               gCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATCTTTTTCATCCTGACCGCAGCCTACCAGTGCCAGTGATCCGATCAggg                                                       <  2:22985/100‑1 (MQ=255)                                               gCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATC‑TTTTCATCCTGACCGCAGCCTACCAGTGCCAGTGATCCGATCAgg                                                        >  2:1149006/1‑98 (MQ=255)                                               gCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATC‑TTTTCATCCTGACCGCAGCCTACCAGTGCCAGTGATCCGATCAGGGc                                                      >  1:223165/1‑100 (MQ=255)                                               gCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATC‑TTTTCATCCTGACCGCAGCCTACCAGTGCCAGTGATCCGATCAGGGc                                                      >  2:1117629/1‑100 (MQ=255)                                                  aCCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATC‑TTTTCATCCTGACCGCAGCCTACCAGTGCCAGTGATCCGATCAGGGCTcc                                                   <  1:1701592/100‑1 (MQ=255)                                                  aCCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATC‑TTTTCATCCTGACCGCAGCCTACCAGTGCCAGTGATCCAATCAGGGCTcc                                                   <  2:1439440/100‑1 (MQ=255)                                                     tgctgTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATC‑TTTTCATCCTGACCGCAGCCTACCAGTGCCAGTGATCCGATCAGGGCTCCCAc                                                <  2:1667284/100‑1 (MQ=255)                                                       ctgTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATCTTTTTCATCCTGACCGCAGCCTACCAGTGCCAGTGATCCGATCAGGGCTCCCACt                                               >  2:1262410/1‑100 (MQ=255)                                                       ctgTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATCTTTTTCATCCTGACCGCAGCCTACCAGTGCCAGTGATCCGATCAGGGCTCCCACt                                               <  2:119174/100‑1 (MQ=255)                                                         gTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATC‑TTTTCATCCTGACCGCAGCCTACCAGTGCCAGTGATCCGATCAGGGCTCCCACTGcc                                            <  2:1264054/100‑1 (MQ=255)                                                         gTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATC‑TTTTCATCCTGACCGCAGCCTACCAGTGCCAGTGATCCGATCAGGGCTCCCACTGcc                                            <  1:1453017/100‑1 (MQ=255)                                                         gTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATC‑TTTTCATCCTGACCGCAGCCTACCAGTGCCAGTGATCCGATCAGGGCTCCCACTGcc                                            >  2:1694202/1‑100 (MQ=255)                                                          ttCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATCTTTTTCATCCTGACCGCAGCCTACCAGTGCCAGTGATCCGATCAGGGCTCCCACTGcc                                            >  2:404852/1‑100 (MQ=255)                                                           tCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATC‑TTTTCATCCTGACCGCAGCCTACCAGTGCCAGTGATCCGATCAGGGCTCCCACTGCCGc                                          <  1:172324/100‑1 (MQ=255)                                                             ggCACCAACAATCACGCCGACTTTAATGTGGTTTGGATCTTTTTCATCCTGACCGCAGCCTACCAGTGCCAGTGATCCGATCAGGGCTCCCACTGCCGCa                                         >  1:1660844/1‑100 (MQ=255)                                                               cACCAACAATCACGCCGACTTTAATGTGGTTTGGATC‑TTTTCATCCTGACCGCAGCCTACCAGTGCCCGTGATCCGATCAGGGCTCCCACTGCCGCAAAg                                      >  2:47573/1‑100 (MQ=255)                                                                   aacaaTCACGCCGACTTTAATGTGGTTTGGATC‑TTTTCATCCTGACCGCAGCCTACCAGTGCCAGTGATCCGATCAGGGCTCCCACTGCCGCAAAGGttt                                  <  1:1580623/100‑1 (MQ=255)                                                                     caaTCACGCCGACTTTAATGTGGTTTGGATCTTTTTCATCCTGACCGCAGCCTACCAGTGCCAGTGATCCGATCAGGGCTCCCACTGCCGCAAAGGtttt                                 <  2:723824/100‑1 (MQ=255)                                                                     caaTCACGCCGACTTTAATGTGGTTTGGATC‑TTTTCATCCTGACCGCAGCCTACCAGTGCCAGTGATCCGATCAGGGCTCCCACTGCCGCAAAGGTtttg                                <  2:1116948/100‑1 (MQ=255)                                                                        tCACGCCGACTTTAATGTGGTTTGTATCTTTTTCATCCTGACCGCAGCCTACCAGTGCCAGTGATCCGATCAGGGCTCCCACTGCAGCAAAGGTtttg                                <  2:644735/98‑1 (MQ=255)                                                                        tCACGCCGACTTTAATGTGGTTTGGATCTTTTTCATCCTGACCGCAGCCTACCAGTGCCAGTGATCCGATCAGGGCTCCCACTGCCGCAAAGGTtttgaa                              >  2:1052762/1‑100 (MQ=255)                                                                          aCGCCGACTTTAATGTGGTTTGGATC‑TTTTCATCCTGACCGCAGCCTACCAGTGCCAGTGATCCGATCAGGGCTCCCACTGCCGCAAAGGTtttgaattt                           >  2:815529/1‑100 (MQ=255)                                                                           cGCCGACTTTAATGTGGTTTGGATCTTTTTCATCCTGACCGCAGCCTACCAGTGCCAGTGATCCGATCAGGGCTCCCACTGCCGCAAAGGTtttgaattt                           >  1:40008/1‑100 (MQ=255)                                                                           cGCCGACTTTAATGTGGTTTGGATC‑TTTTCATCCTGACCGCAGCCTACCAGTGCCAGTGATCCGATCAGGGCTCCCACTGCCGCAAAGGTtttgaatttg                          <  2:153980/100‑1 (MQ=255)                                                                           cGCCGACTTTAATGTGGTTTGGATC‑TTTTCATCCTGACCGCAGCCTACCAGTGCCAGTGATCCGATCAGGGCTCCCACTGCCGCAAAGGTtttgaatttg                          <  2:821031/100‑1 (MQ=255)                                                                                aCTTTAATGTGGTTTGGATCTTTTTCATCCTGACCGCAGCCTACCAGTGCCAGTGATCCGATCAGGGCTCCCACTGCCGCAAAGGTTTTGAATTTGAACg                      <  2:976839/100‑1 (MQ=255)                                                                                aCTTTAATGTGGTTTGGATC‑TTTTCATCCTGACCGCAGCCTACCAGTGCCAGTGATCCGATCAGGGCTCCCACTGCCGCAAAGGTTTTGAATTTGAACGc                     >  2:1191193/1‑100 (MQ=255)                                                                                 cTTTAATGTGGTTTGGATCTTTTTCATCCTGACCGCAGCCTACCAGTGCCAGTGATCCGATCAGGGCTCCCACTGCCGCAAAGGTTTTGAATTTGAACGc                     <  2:801564/100‑1 (MQ=255)                                                                                 cTTTAATGTGGTTTGGATC‑TTTTCATCCTGACCGCAGCCTACCAGTGCCAGTGATCCGATCAGGGCTCCCACTGCCGCAAAGGTTTTGAATTTGAACGcc                    >  2:3218/1‑100 (MQ=255)                                                                                  tttAATGTGTGTTGGATCTTTTTCATCCTTACCGCATCCTACCAGTGCCAGTGATCCGATCATGGCTCCCTCTGCCGCAAAGGTTTTGAATTTGAACGcc                    <  2:765396/100‑1 (MQ=255)                                                                                  tttAATGTGGTTTGGATCTTTTTCATCCTGACCGCAGCCTACCAGTGCCAGTGATCCGATCAGGGCTCCCACTGCCGCAAAGGTTTTGAATTTGAACGcc                    <  1:1566502/100‑1 (MQ=255)                                                                                       tgtgGTTTGGATCTTTTTCATCCTGACCGCAGCCTACCAGTGCCAGTGATCCGATCAGGGCTCCCACTGCCGCAAAGGTTTTGAATTTGAACGCCATAcc               >  2:1802645/1‑100 (MQ=255)                                                                                       tgtgGTTTGGATC‑TTTTCATCCTGACCGCAGCCTACCAGTGCCAGTGATCCGATCAGGGCTCCCACTGCCGCAAAGGTTTTGAATTTGAACGCCATACCt              >  1:1836065/1‑100 (MQ=255)                                                                                        atcttttGGATCTTTTTCATCCTGACCGCAGCCTACCAGTGCCAGTGATCCGATCAGGGCTCCCACTGCCGCAAAGGTTTTGAATTTGAACGCCATACCt              <  1:1950661/96‑1 (MQ=255)                                                                                         tgGTTTGGATCTTTTTCATCCTGACCGCAGCCTACCAGTGCCAGTGATCCGATCAGGGCTCCCACTGCCGCAAAGGTTTTGAATTTGAACGCCATACCtt             >  1:709050/1‑100 (MQ=255)                                                                                          ggTTTGGATC‑TTTTCATCCTGACCGCAGCCTACCAGTGCCAGTGATCCGATCAGGGCTCCCACTGCCGCAAAGGTTTTGAATTTGAACGCCATACCTTAt           >  1:999466/1‑100 (MQ=255)                                                                                          ggTTTGGATC‑TTTTCATCCTGACCGCAGCCTACCAGTGCCAGTGATCCGATCAGGGCTCCCACTGCCGCAAAGGTTTTGAATTTGAACGCCATACCTTAt           <  1:1438257/100‑1 (MQ=255)                                                                                           gTTTGGATC‑TTTTCATCCTGACCGCAGCCTACCAGTGCCAGTGATCCGATCAGGGCTCCCACTGCCGCAAAGGTTTTGAATTTGAACGCCATACCTTAtt          <  2:1485617/100‑1 (MQ=255)                                                                                           gTTTGGATC‑TTTTCATCCTGACCGCAGCCTACCAGTGCCAGTGATCCGATCAGGGCTCCCACTGCCGCAAAGGTTTTGAATTTGAACGCCATACCTTAtt          <  1:1480954/100‑1 (MQ=255)                                                                                            tttGGATCTTTTTCATCCTGACCGCAGCCTACCAGTGCCAGTGATCCGATCAGGGCTCCCACTGCCGCAAAGGTTTTGAATTTGAACGCCATACCtagat          >  2:1950648/1‑96 (MQ=255)                                                                                             ttGGATC‑TTTTCATCCTGACCGCAGCCTACGAGTGCCAGTGATCCGATCAGGGCTCCCACTGCCGCAAAGGTTTTGAATTTGAACGCCATACCTTATTcc        <  1:709094/100‑1 (MQ=255)                                                                                               ggATCTTTTTCATCCTGACCGCAGCCTACCAGTGCCAGTGATCCGATCAGGGCTCCCACTGCCGCAAAGGTTTTGAATTTGAACGCCATACCTTATTCct       >  1:1650830/1‑100 (MQ=255)                                                                                               ggATC‑TTTTCATCCTGACCGCAGCCTACCAGTGCCAGTGATCCGATCAGGGCTCCCACTGCCGCAAAGGTTTTGAATTTGAACGCCATACCTTATTCctt      >  1:235264/1‑100 (MQ=255)                                                                                               ggATC‑TTTTCATCCTGACCGCAGCCTACCAGTGCCAGTGATCCGATCAGGGCTCCCACTGCCGCAAAGGTTTTGAATTTGAACGCCATACCTTATTCctt      >  2:768759/1‑100 (MQ=255)                                                                                                gATCTTTTTCATCCTGACCGCAGCCTACCAGTGCCAGTGATCCGATCAGGGCTCCCACTGCCGCAAAGGTTTTGAATTTGAACGCCATACCTTATTCctt      <  1:2236282/100‑1 (MQ=255)                                                                                                gATCTTTTTCATCCTGACCGCAGCCTACCAGTGCCAGTGATCCGATCAGGGCTCCCACTGCCGCAAAGGTTTTGAATTTGAACGCCATACCTTATTCctt      >  2:1411487/1‑100 (MQ=255)                                                                                                  tCTTTTTCATCCTGACCGCAGCCTACCAGTGCCAGTGATCCGATCAGGGCTCCCACTGCCGCAAAGGTTTTGAATTTGAACGCCATACCTTATTCcttct    >  2:904949/1‑100 (MQ=255)                                                                                                    tttttCATCCTGACCGCAGCCTACCAGTGCCAGTGATCCGATCAGGGCTCCCACTGCCGCAAAGGTTTTGAATTTGAACGCCATACCTTATTCCttcttc  >  1:1875285/1‑100 (MQ=255)                                                                                                    |                                                                                                    AACGTCCAGGCCATATTTGTCTTTCGCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATCTTTTTCATCCTGACCGCAGCCTACCAGTGCCAGTGATCCGATCAGGGCTCCCACTGCCGCAAAGGTTTTGAATTTGAACGCCATACCTTATTCCTTCTTC  >  NC_000913/220746‑220944

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