Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,880,549	T→C	26.4%	T71A (ACC→GCC)	leuE ←	leucine efflux protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,880,549	0	T	C	26.4%	140.7 / 58.7	96	T71A (ACC→GCC)	leuE	leucine efflux protein
Reads supporting (aligned to +/- strand):  ref base T (29/41);  new base C (10/15);  total (40/56)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.17e-01

ATTTTTACCCTTCAGGGTCGCGTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCATACC  >  NC_000913/1880451‑1880645                                                                                                   |                                                                                                 aTTTTTACACTTCAGTGACGCGTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGggtggt                                                                                                  <  1:1660248/99‑1 (MQ=255)  tttttACCCTTCAGGGTCGCGTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGTCt                                                                                                <  2:868254/100‑1 (MQ=255)      tACCCTTCAGGGTCGCGTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGTCttaat                                                                                            <  2:3384184/100‑1 (MQ=255)      tACCCTTCAGGGTCGCGTAAAGAATTTTACACCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATCATATCGGGGTGGTCttaat                                                                                            <  1:2192754/100‑1 (MQ=255)       aCCCTTCAGGGTCGCGTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGTCttaatt                                                                                           <  1:2467274/100‑1 (MQ=255)       aCCCTTCAGGGTCGCGTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCttaat                                                                                            <  1:2653049/99‑1 (MQ=255)        cccTTCAGGGTCGCGTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGTCttaatta                                                                                          >  2:2147729/1‑100 (MQ=255)        cccTTCAGGGTCGCGTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCttaatt                                                                                           <  1:1890773/99‑1 (MQ=255)         ccTTCAGGGTCGCGTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGTCttaattaa                                                                                         <  1:2562460/100‑1 (MQ=255)          cTTCAGGGTCGCGTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGTCTtaattaat                                                                                        >  2:1452173/1‑100 (MQ=255)          cTTCAGGGTCGCGTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTtaattaat                                                                                        >  2:3194674/1‑100 (MQ=255)          cTTCAGGGTCGCGTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGAGAGGGCattaatta                                                                                          >  1:2490049/1‑98 (MQ=255)           ttCAGGGTCGCGTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGTCTTAATTAATg                                                                                       <  2:1622563/100‑1 (MQ=255)              aGGGTCGCGTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGTCTTAATTAATGTCg                                                                                    <  1:2213210/100‑1 (MQ=255)              aGGGTCGCGTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGTCTTAATTAATGTCg                                                                                    <  2:1898552/100‑1 (MQ=255)              aGGGTCGCGTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGTCTTAATTAATGTCg                                                                                    <  2:415961/100‑1 (MQ=255)                ggTCGCGTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGTCTTAATTAATGTCGcc                                                                                  <  2:1148098/100‑1 (MQ=255)                ggTCGCGTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGTCTTAATTAATGTCGcc                                                                                  <  1:1096982/100‑1 (MQ=255)                 gTCGCGTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGTCTTAATTAATGTCGcc                                                                                  <  2:2961049/99‑1 (MQ=255)                    gcgTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTc                                                                              >  2:1264520/1‑100 (MQ=255)                    gcgTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTc                                                                              >  2:1199553/1‑100 (MQ=255)                    gcgTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTc                                                                              >  1:2294649/1‑100 (MQ=255)                    gcgTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTc                                                                              <  1:1179252/100‑1 (MQ=255)                    gcgTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTAGAATAAGATCGGGGTGGCCTtaattaat                                                                                        >  1:3014496/1‑90 (MQ=255)                     cgTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTcc                                                                             <  1:2612355/100‑1 (MQ=255)                     cgTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTcc                                                                             <  2:1666309/100‑1 (MQ=255)                      gTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTCCa                                                                            <  1:1899456/100‑1 (MQ=255)                       tAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTCCAg                                                                           <  2:1579027/100‑1 (MQ=255)                        aaaGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTCCAGc                                                                          <  2:2186783/100‑1 (MQ=255)                         aaGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTCCAGcc                                                                         >  1:1307305/1‑100 (MQ=255)                         aaGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGcc                                                                         >  2:3233492/1‑100 (MQ=255)                           gAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTCCAGCCCa                                                                       >  2:750050/1‑100 (MQ=255)                            aaTTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTCCAGCCCAt                                                                      <  1:3194203/100‑1 (MQ=255)                              ttttACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACCCCAGCCCACGc                                                                    >  1:3044090/1‑100 (MQ=255)                               tttACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTCCAGCCCATGcc                                                                   <  2:542835/100‑1 (MQ=255)                                 tACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAg                                                                 >  2:2485678/1‑100 (MQ=255)                                     ccccAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAAc                                                             <  2:3160654/100‑1 (MQ=255)                                     ccccAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAAc                                                             >  2:2257323/1‑100 (MQ=255)                                      cccAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACa                                                            >  2:3237593/1‑100 (MQ=255)                                        cAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATc                                                          <  2:1498705/100‑1 (MQ=255)                                        cAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATc                                                          <  1:3326901/100‑1 (MQ=255)                                        cAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATc                                                          <  1:1544280/100‑1 (MQ=255)                                        cAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATc                                                          <  2:357395/100‑1 (MQ=255)                                         aGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCa                                                         <  2:3240347/100‑1 (MQ=255)                                          gATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCaa                                                        >  2:707878/1‑100 (MQ=255)                                           aTAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAAt                                                       >  1:3002271/1‑100 (MQ=255)                                            tAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATa                                                      >  2:2800661/1‑100 (MQ=255)                                              gagCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATAcc                                                    >  2:184032/1‑100 (MQ=255)                                              gagCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATAcc                                                    <  2:1806089/100‑1 (MQ=255)                                              gagCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATAcc                                                    <  1:2480489/100‑1 (MQ=255)                                              gagCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATAcc                                                    >  1:2501501/1‑100 (MQ=255)                                               agCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTCAAGCCCATGCCAGAAACATCAATAcc                                                    >  1:1524308/1‑99 (MQ=255)                                               agCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTAGTCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATAcc                                                    <  1:2163720/99‑1 (MQ=255)                                                gCAAATAAAACGCACCAAGATAAGGTACAATGTTGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATAcc                                                    <  1:803630/98‑1 (MQ=255)                                                gCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGc                                                  >  2:1601846/1‑100 (MQ=255)                                                gCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGc                                                  >  2:1558762/1‑100 (MQ=255)                                                gCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGc                                                  <  1:858458/100‑1 (MQ=255)                                                 cAAATACATTGCACCAAGGTAACGTACAATGTTGAATTATATCGGGGTTTTCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCATTACCGc                                                  <  1:2016339/99‑1 (MQ=255)                                                 cAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCa                                                 <  1:681981/100‑1 (MQ=255)                                                 cAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCa                                                 >  2:71023/1‑100 (MQ=255)                                                 cAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCa                                                 <  1:702085/100‑1 (MQ=255)                                                   aaTAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCAt                                                >  2:1512198/1‑99 (MQ=255)                                                   aaTAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATc                                               >  2:426497/1‑100 (MQ=255)                                                    aTAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCg                                              <  1:2862447/100‑1 (MQ=255)                                                    aTAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCg                                              <  1:46484/100‑1 (MQ=255)                                                    aTAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCg                                              <  2:3115493/100‑1 (MQ=255)                                                     tAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGc                                             >  2:2633528/1‑100 (MQ=255)                                                      aaaaCGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGcc                                            <  2:1304938/100‑1 (MQ=255)                                                      aaaaCGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGcc                                            <  2:774922/100‑1 (MQ=255)                                                      aaaaCGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACAGCATCGcc                                            >  1:1676872/1‑100 (MQ=255)                                                      aaaaCGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGcc                                            >  2:3319202/1‑100 (MQ=255)                                                        aaCGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCaa                                          >  1:69760/1‑100 (MQ=255)                                                         aCGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCaat                                         <  1:2909997/100‑1 (MQ=255)                                                         aCGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCaat                                         <  2:184274/100‑1 (MQ=255)                                                           gCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCaataa                                       >  1:545535/1‑100 (MQ=255)                                                               cAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCACATGCCAGAAACATCAATACCGCATCGCCCATAAATac                                   >  1:87641/1‑100 (MQ=255)                                                                aaGATAACGTACAATGTTGAATAATATCGGGGTGGACTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATaca                                  >  2:1150418/1‑100 (MQ=255)                                                                  gATAACGTACAATGTTGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACAcc                                <  2:2112799/100‑1 (MQ=255)                                                                    tAACGTACAATGTTGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGc                              >  2:580323/1‑100 (MQ=255)                                                                    tAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGc                              >  2:3238100/1‑100 (MQ=255)                                                                      aCGTACAATGTTGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAg                            <  1:3086114/100‑1 (MQ=255)                                                                       cGTACAATGTTGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAgg                           >  2:3137386/1‑100 (MQ=255)                                                                        gTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGc                          <  1:2852445/100‑1 (MQ=255)                                                                        gTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGc                          >  1:2802198/1‑100 (MQ=255)                                                                          aCAATGTTGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCg                        <  1:2172428/100‑1 (MQ=255)                                                                             aTGTTGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCaa                     >  1:2135733/1‑100 (MQ=255)                                                                             aTGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCaa                     <  2:139706/100‑1 (MQ=255)                                                                              tGTTGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAg                    <  2:1875255/100‑1 (MQ=255)                                                                               gTTGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGa                   <  2:3002297/100‑1 (MQ=255)                                                                               gTTGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGa                   >  2:1900930/1‑100 (MQ=255)                                                                                ttGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCACTAAATACACCGCAGGCCGCAAGAt                  >  1:2545723/1‑100 (MQ=255)                                                                                 tGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATa                 <  1:1502125/100‑1 (MQ=255)                                                                                     taataTCGGGGTGGTCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCg             <  2:891037/100‑1 (MQ=255)                                                                                     taataTCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCg             <  1:3241442/100‑1 (MQ=255)                                                                                        tatCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCt          <  1:3306766/100‑1 (MQ=255)                                                                                          tCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCttt        <  1:1923042/100‑1 (MQ=255)                                                                                             gggTGGCCGTAATTAATGTCGCCACTCCAGCCCATGCCAGAACCCTCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCa      <  1:791308/99‑1 (MQ=255)                                                                                               gtggcCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCATAc   >  2:2314180/1‑100 (MQ=255)                                                                                                tggtCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCATAcc  <  2:1427891/100‑1 (MQ=255)                                                                                                   |                                                                                                 ATTTTTACCCTTCAGGGTCGCGTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCATACC  >  NC_000913/1880451‑1880645

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