Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,079,293	(T)6→5	100%	intergenic (‑411/‑12)	putA ← / → putP	fused DNA‑binding transcriptional regulator/proline dehydrogenase/pyrroline‑5‑carboxylate dehydrogenase/proline:sodium symporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,079,288	0	T	.	79.7%	195.6 / 65.2	79	intergenic (‑406/‑17)	putA/putP	fused DNA‑binding transcriptional regulator/proline dehydrogenase/pyrroline‑5‑carboxylate dehydrogenase/proline:sodium symporter
Reads supporting (aligned to +/- strand):  ref base T (7/9);  new base . (32/31);  total (39/40)
Fisher's exact test for biased strand distribution p-value = 7.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.54e-01

GCTTCCGGCAGGATACGGTCGCCCTGGTAAAACATAAACTCTGTTACCCCGTTCCGGTGGCAGATATAACGGCAAGTTTCGACATTGCCGATAATAATTTTTTGGAGACTTTAGATGGCTATTAGCACACCGATGTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGC  >  NC_000913/1079191‑1079383                                                                                                  |                                                                                                gcTTCCGGCAGGATACGGTCGCCCTGGTAAAACATAAACTCTGTTACCCCGTTCCGGTGGCAGATATAACGGCAAGTTTCGACATTGCCGCTaataat                                                                                                 >  2:1120521/1‑98 (MQ=255) gcTTCCGGCAGGATACGGTCGCCCTGGTAAAACATAAACTCTGTTACCCCGTTCCGGTGGCAGATATAACGGCAAGTTTCGACATTGCCGATAATAAttt                                                                                               <  1:537140/100‑1 (MQ=255) gcTTCCGGCAGGATACGGTCGCCCTGGTAAAACATAAACTCTGTTACCCCGTTCCGGTGGCAGATATAACGGCAAGTTTCGACATTGCCGATAATAAttt                                                                                               >  2:1783680/1‑100 (MQ=255)    tCCGGCAGGATACGGTCGCCCTGGTAAAACATAAACTCTGTTACCCCGTTCCGGTGGCAGATATAACGGCAAGTTTCGGCATTGCCGATAATAAtttttg                                                                                            <  1:1112755/100‑2 (MQ=255)    tCCGGCAGGATACGGTCGCCCTGGTAAAACATAAACTCTGTTACCCCGTTCCGGTGGCAGATATAACGGCAAGTTTCGACATTGCCGATAATAAtttttg                                                                                            <  1:2072253/100‑2 (MQ=255)      cGGCAGGATACGGTCGCCCTGGTAAAACATAAACTCTGTTACCCCGTTCCGGTGGCAGATATAACGGCAAGTTTCGACATTGCCGATAATAAtttttgga                                                                                          <  2:1229614/100‑4 (MQ=255)      cGGCAGGATACGGTCGCCCTGGTAAAACATAAACTCTGTTACCCCGTTCCGGTGGCAGATATAACGGCAAGTTTCGACATTGCCGATAATAAtttttgga                                                                                          <  1:203480/100‑4 (MQ=255)         cAGGATACGGTCGCCCTGGTAAAACATAAACTCTGTTACCCCGTTCCGGTGGCAGATATAACGGCAAGTTTCGACATTGCCGATAATAATTTTTTGgaga                                                                                       <  2:1610685/100‑1 (MQ=255)          aGGATACGGTCGCCCTGGTAAAACATAAACTCTGTTACCCCGTTCCGGTGGCAGATATAACGGCAAGTTTCGACATTGCCGATAATAA‑TTTTTGGAGACt                                                                                     <  2:2417839/100‑1 (MQ=255)           ggATACGGTCGCCCTGGTAAAACATAAACTCTGTTACCCCGTTCCGGTGGCAGATATAACGGCAAGTTTCGACATTGCCGATAATAA‑TTTTTGGAGACtt                                                                                    >  2:588575/1‑100 (MQ=255)           ggATACGGTCGCCCTGGTAAAACATAAACTCTGTTACCCCGTTCCGGTGGCAGATATAACGGCAAGTTTCGACATTGCCGATAATAA‑TTTTTGGAGACtt                                                                                    >  1:2500287/1‑100 (MQ=255)             aTACGGTCGCCCTGGTAAAACATAAACTCTGTTACCCCGTTCCGGTGGCAGATATAACGGCAAGTTTCGACATTGCCGATAATAA‑TTTTTGGAGACTTTa                                                                                  >  2:1394055/1‑100 (MQ=255)              tACGGTCGCCCTGGTAAAACATAAACTCTGTTACCCCGTTCCGGTGGCAGATATAACGGCAAGTTTCGACATTGCCGATAATAA‑TTTTTGGAGACTTTAg                                                                                 <  2:1393718/100‑1 (MQ=255)               aCGGTCGCCCTGGTAAAACATAAACTCTGTTACCCCGTTCCGGTGGCAGATATAACGGCAAGTTTCGACATTGCCGATAATAA‑TTTTTGGAGACTTTAGa                                                                                <  2:907778/100‑1 (MQ=255)                 ggTCGCCCTGGTAAAACATAAACTCTGTTACCCCGTTCCGGTGGCAGATATAACGGCAAGTTTCGACATTGCCGATAATAATTTTTTGGAGACTTTAGAt                                                                               <  2:2529336/100‑1 (MQ=255)                  gTCGCCCTGGTAAAACATAAACTCTGTTACCCCGTTCCGGTGGCAGATATAACGGCAAGTTTCGACATTGCCGATAATAATTTTTTGGAGACTTTAGATg                                                                              <  1:2513520/100‑1 (MQ=255)                      cccTGGTAAAACATAAACTCTGTTACCCCGTTCCGGTGGCAGATATAACGGCAAGTTTCGACATTGCCGATAATAATTTTTTGGAGACTTTAGATGGCTa                                                                          <  1:243570/100‑1 (MQ=255)                         tGGTAAAACATAAACTCTGTTACCCCGTTCCGGTGGCAGATATAACGGCAAGTTTCGACATTGCCGATAATAATTTTTTGGAGACTTTAGATGGCTATTa                                                                       >  2:534837/1‑100 (MQ=255)                         tGGTAAAACATAAACTCTGTTACCCCGTTCCGGTGGCAGATATAACGGCAAGTTTCGACATTGCCGATAATAA‑TTTTTGGAGACTTTAGATGGCTATTAg                                                                      >  2:1079853/1‑100 (MQ=255)                         tGGTAAAACATAAACTCTGTTACCCCGTTCCGGTGGCAGATATAACGGCAAGTTTCGACATTGCCGATAATAA‑TTTTTGGAGACTTTAGATGGCTATTAg                                                                      >  2:109952/1‑100 (MQ=255)                          ggTAAAACATAAACTCTGTTACCCCGTTCCGGTGGCAGATATAACGGCAAGTTTCGACATTGCCGATAATAA‑TTTTTGGAGACTTTAGATGGCTATTAGc                                                                     <  1:1141498/100‑1 (MQ=255)                          ggTAAAACATAAACTCTGTTACCCCGTTCCGGTGGCAGATATAACGGCAAGTTTCGACATTGCCGATAATAA‑TTTTTGGAGACTTTAGATGGCTATTAGc                                                                     <  1:239817/100‑1 (MQ=255)                             aaaaCATAAACTCTGTTACCCCGTTCCGGTGGCAGATATAACGGCAAGTTTCGACATTGCCGATAATAA‑TTTTTGGAGACTTTAGATGGCTATTAGcaca                                                                  >  1:2146147/1‑100 (MQ=255)                              aaaCATAAACTCTGTTACCCCGTTCCGGTGGCAGATATAACGGCAAGTTTCGACATTGCCGATAATAA‑TTTTTGGAGACTTTAGATGGCTATTAGCacac                                                                 <  1:1274727/100‑1 (MQ=255)                               aaCATAAACTCTGTTACCCCGTTCCGGTGGCAGATATAACGGCAAGTTTCGACATTGCCGATAATAA‑TTTTTGGAGACTTTAGATGGCTATTAGCACAcc                                                                <  1:776753/100‑1 (MQ=255)                                aCATAAACTCTGTTACCCCGTTCCGGTGGCAGATATAACGGCAAGTTTCGACATTGCCGATAATAA‑TTTTTGGAGACTTTAGATGGCTATTAGCACACCg                                                               <  2:864811/100‑1 (MQ=255)                                    aaaCTCTGTTACCCCGTTCCGGTGGCAGATATAACGGCAAGTTTCGACATTGCCGATAATAA‑TTTTTGGAGACTTTAGATGGCTATTAGCACACCGATg                                                            >  2:1576629/1‑99 (MQ=255)                                    aaaCTCTGTTACCCCGTTCCGGTGGCAGATATAACGGCAAGTTTCGACATTGCCGATAATAA‑TTTTTGGAGACTTTAGATGGCTATTAGCACACCGATGt                                                           >  1:1614099/1‑100 (MQ=255)                                    aaaCTCTGTTACCCCGTTCCGGTGGCAGATATAACGGCAAGTTTCGACATTGCCGATAATAA‑TTTTTGGAGACTTTAGATGGCTATTAGCACACCGATGt                                                           >  2:1795019/1‑100 (MQ=255)                                    aaaCTCTGTTACCCCGTTCCGGTGGCAGATATAACGGCAAGTTTCGACATTGCCGATAATAA‑TTTTTGGAGACTTTAGATGGCTATTAGCACACCGATGt                                                           >  1:793857/1‑100 (MQ=255)                                     aaCTCTGTTACCCCGTTCCGGTGGCAGATATAACGGCAAGTTTCGACATTGCCGATAATAA‑TTTTTGGAGACTTTAGATGGCTATTAGCACACCGATGtt                                                          <  1:2414042/100‑1 (MQ=255)                                       ctctGTTACCCCGTTCCGGTGGCAGATATAACGGCAAGTTTCGACATTGCCGATAATAATTTTTTGGAGACTTTAGATGGCTATTAGCACACCGATGTTg                                                         >  1:1861940/1‑100 (MQ=255)                                       ctctGTTACCCCGTTCCGGTGGCAGATATAACGGCAAGTTTCGACATTGCCGATAATAA‑TTTTTGGAGACTTTAGATGGCTATTAGCACACCGATGTTgg                                                        >  2:837049/1‑100 (MQ=255)                                       ctctGTTACCCCGTTCCGGTGGCAGATATAACGGCAAGTTTCGACATTGCCGATAATAA‑TTTTTGGAGACTTTAGATGGCTATTAGCACACCGATGTTgg                                                        <  1:2132220/100‑1 (MQ=255)                                         ctGTTACCCCGTTCCGGTGGCAGATATAACGGCAAGTTTCGACATTGCCGATAATAA‑TTTTTGGAGACTTTAGATGGCTATTAGCACACCGATGTTGGTg                                                      >  2:638365/1‑100 (MQ=255)                                          tGTTACCCCGTTCCGGTGGCAGATATAACGGCAAGTTTCGACATTGCCGATAATAA‑TTTTTGGAGACTTTAGATGGCTATTAGCACACCGATGTTGGTGa                                                     <  2:1921731/100‑1 (MQ=255)                                          tGTTACCCCGTTCCGGTGGCAGATATAACGGCAAGTTTCGACATTGCCGATAATAA‑TTTTTGGAGACTTTAGATGGCTATTAGCACACCGATGTTGGTGa                                                     >  1:1530582/1‑100 (MQ=255)                                          tGTTACCCCGTTCCGGTGGCAGATATAACGGCAAGTTTCGACATTGCCGATAATAA‑TTTTTGGAGACTTTAGATGGCTATTAGCACACCGATGTTGGTGa                                                     >  2:2046448/1‑100 (MQ=255)                                          tGTTACCCCGTTCCGGTGGCAGATATAACGGCAAGTTTCGACATTGCCGATAATAA‑TTTTTGGAGACTTTAGATGGCTATTAGCACACCGATGTTGGTGa                                                     >  2:165956/1‑100 (MQ=255)                                               ccccGTTCCGGTGGCAGATATAACGGCAAGTTTCGACATTGCCGATAATAA‑TTTTTGGAGACTTTAGATGGCTATTAGCACACCGATGTTGGTGACAttt                                                >  1:2400471/1‑100 (MQ=255)                                                 ccGTTCCGGTGGCAGATATAACGGCAAGTTTCGACATTGCCGATAATAA‑TTTTTGGAGACTTTAGATGGCTATTAGCACACCGATGTTGGTGACATTTtg                                              >  2:1247633/1‑100 (MQ=255)                                                 ccGTTCCGGTGGCAGATATAACGGCAAGTTTCGACATTGCCGATAATAA‑TTTTTGGAGACTTTAGATGGCTATTAGCACACCGATGTTGGTGACATTTtg                                              <  2:2131701/100‑1 (MQ=255)                                                   gTTCCGGTGGCAGATATAACGGCAAGTTTCGACATTGCCGATAATAA‑TTTTTGGAGACTTTAGATGGCTATTAGCACACCGATGTTGGTGACATTTtgtg                                            <  1:2293434/100‑1 (MQ=255)                                                    ttCAGGTGGCAGATATAACGGCAAGTTTCGTCATTGCCGATAATAA‑TTTTTGGAGACTTTATAGGGCTATTAGCACACCGATGTTGGTGACATTTtg                                              <  2:2366297/97‑1 (MQ=255)                                                     tCCGGTGGCAGATATAACGGCAAGTTTCGACATTGCCGATAATAATTTTTTGGAGACTTTAGATGGCTATTAGCACACCGATGTTGGTGACATTTTgtgt                                           >  1:2410425/1‑100 (MQ=255)                                                     tCCGGTGGCAGATATAACGGCAAGTTTCGACATTGCCGATAATAA‑TTTTTGGAGACTTTAGATGGCTATTAGCACACCGATGTTGGTGACATTTTGTGTc                                          <  1:415398/100‑1 (MQ=255)                                                       cGGTGGCAGATATAACGGCAATTTTCGACATTGCTGATAAT‑ATTTGTTGGAGACTTTAGATTGAGATTAGCACACAGATGTTGGTGACATTTTTTGTcga                                        <  1:10070/100‑3 (MQ=255)                                                       cGGTGGCAGATATAACAGTAAGAAACGACATTGCCGATAATAATTTTTTGGAGACTTTAGATGGCTATTAGCACACCGATGTTGGTGACATTTTGTGTCt                                         >  1:356184/1‑100 (MQ=255)                                                        ggTGGCAGATATAACGGCAAGTTTCGACATTGCCGATAATAA‑TTTTTGGAGACTTTAGATGGCTATTAGCACACCGATGTTGGTGACATTTTGTGTCtat                                       <  2:1990109/100‑1 (MQ=255)                                                          tGGCAGATATAACGGCAAGTTTCGACATTGCCGATAATAA‑TTTTTGGAGACTTTAGATGGCTATTAGCACACCGATGTTGGTGACATTTTGTGTCTatat                                     >  1:1019042/1‑100 (MQ=255)                                                           ggCAGATATAACGGCAAGTTTCGACATTGCCGATAATAA‑TTTTTGGAGACTTTAGATGGCTATTAGCACACCGATGTTGGTGACATTTTGTGTCTATATc                                    >  2:141479/1‑100 (MQ=255)                                                             cAGATATAACGGCAAGTTTCGACATTGCCGATAATAA‑TGTTTGGAGACTTTAGATGGGTATTAACACATCGATGTTGGTGACATTTTATGTTTATTTcct                                  >  2:1012854/1‑98 (MQ=255)                                                                atatAACGGCAAGTTTCGACATTGCCGATAATAATTTTTTGGAGACTTTAGATGGCTATTAGCACACCGATGTTGGTGACATTTTGTGTCTATATCTTTg                                <  1:534835/100‑1 (MQ=255)                                                                atatAACGGCAAGTTTCGACATTGCCGATAATAA‑TTTTTGGAGACTTTAGATGGCTATTAGCACACCGATGTTGGTGACATTTTGTGTCTATATCTTTgg                               >  1:839466/1‑100 (MQ=255)                                                                 tataACGGCAAGTTTCGACATTGCCGATAATAA‑TTTTTGGAGACTTTAGATGGCTATTAGCACACCGATGTTGGTGACATTTTGTGTCTATATCTTTGGc                              >  1:701453/1‑100 (MQ=255)                                                                  ataACGGCAAGTTTCGACATTGCCGATAATAA‑TTTTTGGAGACTTTAGATGGCTATTAGCACACTGATGTTGGTGACATTTTGTGTCTATATCTTTgtta                             <  1:1242603/100‑4 (MQ=255)                                                                  ataACGGCAAGTTTCGACATTGCCGATAATAA‑TTTTTGGAGACTTTAGATGGCTATTAGCACACCGATGTTGGTGACATTTTGTGTCTATATCTTTGGCa                             <  2:59675/100‑1 (MQ=255)                                                                   taACGGCAAGTTTCGACATTGCCGATAATAA‑TTTTTGGAGACTTTAGATGGCTATTAGCACACCGATGTTGGTGACATTTTGTGTCTATATCTTTGACAt                            <  1:2544798/100‑1 (MQ=255)                                                                    aaCGGCAAGTTTCGACATTGCCGATAATAA‑TTTTTGGAGACTTTAGATGGCTATTAGCACACCGATGTTGGTGACATTTTGTGTCTATATCTTTGGCATg                           <  2:1391527/100‑1 (MQ=255)                                                                    aaCGGCAAGTTTCGACATTGCCGATAATAA‑TTTTTGGAGACTTTAGATGGCTATTAGCACACCGATGTTGGTGACATTTTGTGTCTATATCTTTGGCATg                           <  2:608191/100‑1 (MQ=255)                                                                    aaCGGCAAGTTTCGACATTGCCGATAATAA‑TTTTTGGAGACTTTAGATGGCTATTAGCACACCGATGTTGGTGACATTTTGTGTCTATATCTTTGGCATg                           <  1:1732674/100‑1 (MQ=255)                                                                        gCAAGTTTCGACATTGCCGATAATAATTTTTTGGAGACTTTAGATGGCTATTAGCACACCGATGTTGGTGACATTTTGTGTCTATATCTTTGGCATGata                        <  2:128005/100‑1 (MQ=255)                                                                        gCAAGTTTCGACATTGCCGATAATAATTTTTTGGAGACTTTAGATGGCTATTAGCACACCGATGTTGGTGACATTATGTGTCTATATCTTTGGCATGata                        >  2:88830/1‑100 (MQ=255)                                                                        gCAAGTTTCGACATTGCCGATAATAA‑TTTTTGGAGACTTTAGATGGCTATTAGCACACCGATGTTGGTGACATTTTGTGTCTATATCTTTGGCATGatat                       >  1:1311995/1‑100 (MQ=255)                                                                        gCAAGTTTCGACATTGCCGATAATAA‑TTTTTGGAGACTTTAGATGGCTATTAGCACACCGATGTTGGTGACATTTTGTGTCTATATCTTTGGCATGatat                       >  2:373422/1‑100 (MQ=255)                                                                         cAAGTTTCGACATTGCCGATAATAA‑TTTTTGGAGACTTTAGATGGCTATTAGCACACCGATGTTGGTGACATTTTGTGTCTATATCTTTGGCATGATatt                      <  2:2220449/100‑1 (MQ=255)                                                                          aaGTTTCGACATTGCCGATAATTATTTTTTGGAGACTTTAGATGGCTATTAGCACACCGATGTTGGTGACATTTTGTGTCTATATCTTTGGCATGata                        <  1:361143/98‑1 (MQ=255)                                                                            gTTTCGACATTGCCGATAATAATTTTTTGGAGACTTTAGATGGCTATTAGCACACCGATGTTGGTGACATTTTGTGTCTATATCTTTGGCATGATattga                    >  2:527896/1‑100 (MQ=255)                                                                            gTTTCGACATTGCCGATAATAA‑TTTTTGGAGACTTTAGATGGCTATTAGCACACCGATGTTGGTGACATTTTGTGTCTATATCTTTGGCCTGATgttgat                   >  1:2000950/1‑100 (MQ=255)                                                                              ttCGACATTGCCGATAATAA‑TTTTTGGAGACTTTAGATGGCTATTAGCACACCGATGTTGGTGACATTTTGTGTCTATATCTTTGGCATGATattgattg                 <  2:1325150/100‑1 (MQ=255)                                                                                cGACATTGCCGATAATAA‑TTTTAGAAGACTTTAGATGGCTATTAGCACACCGATGTTGGTGACATATTGTGTCTATAGCTATGGCATGATattgaat                  >  1:1270399/1‑95 (MQ=255)                                                                                 gACATTGCCGATAATAA‑TTTTTGGAGACTTTAGATGGCTATTAGCACACCGATGTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGt              >  1:2536026/1‑100 (MQ=255)                                                                                 gACATTGCCGATAATAA‑TTTTTGGAGACTTTAGATGGCTATTAGCACACCGATGTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGt              >  2:826244/1‑100 (MQ=255)                                                                                   cATTGCCGATAATAA‑TTTTTGGAGACTTTAGATGGCTATTAGCACACCGATGTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGttt            >  1:1243390/1‑100 (MQ=255)                                                                                    aTTGCCGATAATAA‑TTTTTGGAGACTTTATATGGCTATTAGCACACCGATGTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTa           <  1:1886690/100‑1 (MQ=255)                                                                                    aTTGCCAATAACAA‑TTTTTGGAGACTTTAGATGGCTATTAGCACACCGATGTTGGTGAAATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGtt             <  1:426144/98‑1 (MQ=255)                                                                                      tGCCGATAATAA‑TTTTTGGAGACTTTAGATGGCTATTAGCACACCGATGTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATc         <  2:1663122/100‑1 (MQ=255)                                                                                       gCCGATAATAATTTTTTTTTTTTTTTATATGGCTATTAGCACACCGATGTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATc         >  1:80353/1‑100 (MQ=255)                                                                                       gCCGATAATAA‑TTTTTGGAGACTTTAGATGGCTATTAGCACACCGATGTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCg        >  2:831412/1‑100 (MQ=255)                                                                                       gCCGATAATAA‑TTTTTGGAGACTTTAGATGGCTATTAGCACACCGATGTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCg        <  2:538502/100‑1 (MQ=255)                                                                                        ccGATAATAA‑TTTTTGGAGACTTTAGATGGCTATTAGCACACCGATGTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGc       >  1:194996/1‑100 (MQ=255)                                                                                        ccGATAATAA‑TTTTTGGAGACTTTAGATGGCTATTAGCACACCGATGTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGc       <  2:2174357/100‑1 (MQ=255)                                                                                          gATAATAATTTTTTGGAGACTTTAGATGGCTATTAGCACACCGATGTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGcc      <  1:983875/100‑1 (MQ=255)                                                                                          gATAATAA‑TTTTTGGAGACTTTAGATGGCTATTAGCACACCGATGTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCt     >  1:1994076/1‑100 (MQ=255)                                                                                            taataa‑TTTTTGGAGACTTTAGATGGCTATTAGCACACCGATGTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTgg   <  1:831409/100‑1 (MQ=255)                                                                                            taataa‑TTTTTGGAGACTTTAGATGGCTATTAGCACACCGATGTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTgg   <  1:1095924/100‑1 (MQ=255)                                                                                             aataa‑tTTTTGGAGACTTTAGATGGCTATTAGCACACCGATGTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGc  <  2:2363627/100‑1 (MQ=255)                                                                                             aataa‑tTTTTGGAGACTTTAGATGGCTATTAGCACACCGATGTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGc  <  1:52004/100‑1 (MQ=255)                                                                                             aataa‑tTTTTGGAGACTTTAGATGGCTATTAGCACACCGATATTGGTGACATTTTGTGTCTATATCTTAGGCATGATATTGATTGGGTTTATCGCCTgg   >  1:1529253/1‑99 (MQ=255)                                                                                                  |                                                                                                GCTTCCGGCAGGATACGGTCGCCCTGGTAAAACATAAACTCTGTTACCCCGTTCCGGTGGCAGATATAACGGCAAGTTTCGACATTGCCGATAATAATTTTTTGGAGACTTTAGATGGCTATTAGCACACCGATGTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGC  >  NC_000913/1079191‑1079383

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