Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,881,128	C→A	100%	E228* (GAA→TAA)	dmlR ←	transcriptional activator of dmlA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,881,128	0	C	A	98.8%	311.7 / ‑3.8	95	E228* (GAA→TAA)	dmlR	transcriptional activator of dmlA
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base A (49/44);  total (50/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATCGCGCTCTTT  >  NC_000913/1881031‑1881224                                                                                                  |                                                                                                 aCTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTtatc                                                                                                <  2:1323800/100‑1 (MQ=255) aCTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTtatc                                                                                                <  1:226091/100‑1 (MQ=255) aCTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTtatc                                                                                                <  2:417700/100‑1 (MQ=255)  cTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATcc                                                                                               <  2:959199/100‑1 (MQ=255)  cTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATcc                                                                                               >  2:945649/1‑100 (MQ=255)    tGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCag                                                                                             >  1:959205/1‑99 (MQ=255)    tGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCAt                                                                                             <  1:1932259/100‑1 (MQ=255)      cACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGa                                                                                           <  2:507683/100‑1 (MQ=255)      cACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACAACACTGATGGACAATTTATCCATGa                                                                                           >  1:932846/1‑100 (MQ=255)       aCTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGAt                                                                                          >  1:512458/1‑100 (MQ=255)        cTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGAt                                                                                          <  1:370424/99‑1 (MQ=255)        cTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCAAGACACCACTGATGGACAATTTATCCATGAtt                                                                                         >  1:826356/1‑100 (MQ=255)         tAAATGACCACTGGCAATGTTTTCTCTGCCATCCCACGAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCTTTGATGGACAATTTATCCATGAt                                                                                          <  2:1763416/98‑1 (MQ=255)         tAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTa                                                                                        >  1:1512818/1‑100 (MQ=255)           aaTGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTa                                                                                        >  2:2956922/1‑98 (MQ=255)           aaTGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAga                                                                                      <  1:1332224/100‑1 (MQ=255)           aaTGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAga                                                                                      >  1:2154823/1‑100 (MQ=255)            aTGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAgac                                                                                     <  1:148380/100‑1 (MQ=255)            aTGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAgac                                                                                     <  1:945655/100‑1 (MQ=255)             tGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCATGCACATAGAAACACAACTGATGGACAATTTATCCATGATTaaa                                                                                      >  1:766833/1‑96 (MQ=255)                ccACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACAACTGATGGACAATTTATCCCTGATTAGACGACa                                                                                 >  2:2348012/1‑100 (MQ=255)                   cTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTg                                                                              <  1:2796902/100‑1 (MQ=255)                    tGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCACGCTTAGACGACAGTg                                                                              >  2:1673447/1‑99 (MQ=255)                       cAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGcc                                                                          <  2:49056/100‑1 (MQ=255)                         aTGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCgg                                                                        <  2:2561098/100‑1 (MQ=255)                           gTTTTCGCTAACATCCTACTAGGAACGCAGAGCAATGCCTTGTCCATCGCGACACCACTGTTGTACACTTTCTTCATGCTTGGCCGTGAGTggtcttgc                                                                       >  2:2330112/1‑93 (MQ=255)                            ttttCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGa                                                                     >  1:1203559/1‑100 (MQ=255)                             tttCGCTAACATCCCACCATGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGAc                                                                    <  2:888124/100‑1 (MQ=255)                              ttCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACt                                                                   >  1:1360868/1‑100 (MQ=255)                              ttCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACt                                                                   >  1:2717425/1‑100 (MQ=255)                                cGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACttt                                                                 <  1:2263950/100‑1 (MQ=255)                                cGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACttt                                                                 <  2:2976225/100‑1 (MQ=255)                                 gCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTa                                                                >  1:1036683/1‑100 (MQ=255)                                   tAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTaa                                                               <  2:2860505/99‑1 (MQ=255)                                    aaCATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATc                                                             <  2:78365/100‑1 (MQ=255)                                    aaCATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGAGACTTTAATc                                                             >  2:1520491/1‑100 (MQ=255)                                    aaCATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCAATGATTAGACGACAGTGAGCCGGTGACTTTAAt                                                              >  1:163235/1‑99 (MQ=255)                                     aCATCCCACCAGGAACGCAGAGCAATCCCCTGCGCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCg                                                            <  2:2655667/100‑1 (MQ=255)                                     aCATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCg                                                            <  2:1431279/100‑1 (MQ=255)                                     aCATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATACATGATTAGACGACAGTGGGCCGGTGACTTTAATCg                                                            >  1:780393/1‑100 (MQ=255)                                      cATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGc                                                           <  1:1021876/100‑1 (MQ=255)                                       aTCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCa                                                          >  2:1008353/1‑100 (MQ=255)                                        tCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCAt                                                         >  1:1837866/1‑100 (MQ=255)                                          ccaccaGGAACGCAGAGCAATCCCCTGCCCATCTAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATg                                                        >  1:1156262/1‑99 (MQ=255)                                          ccaccaGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATgg                                                       >  2:3039025/1‑100 (MQ=255)                                          ccaccaGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATAGCATgg                                                       >  1:1156730/1‑100 (MQ=255)                                           caccaGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTAGGCCGGTGACTTTAATCGCATggg                                                      >  1:2793058/1‑100 (MQ=255)                                            accaGGAACGCAGAGCAATCCCCTCCCCATCGAGACACCACTGATGGACAAATTATCCATGATTAGACGACAGTGGTCCGGTGACTTTAAGCGCATaggg                                                     >  2:1181563/1‑100 (MQ=255)                                             ccaGGAACGCAGAGCAATCCCCTGCCCATCGATACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATggggg                                                    <  1:2526275/100‑1 (MQ=255)                                             ccaGGAACGCAGAGCAATCCCCTGCCCATCGAGACTCCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATggg                                                      <  2:2293118/98‑1 (MQ=255)                                             ccaGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATggggg                                                    <  2:1582392/100‑1 (MQ=255)                                             ccaGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATggggg                                                    <  1:1008359/100‑1 (MQ=255)                                             ccaGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATggggg                                                    <  2:2717393/100‑1 (MQ=255)                                             ccaGGAACGCAGAGCAATCCCCTGCCCATCGAAACACCACTGATGGCCAATTTATACATAAAAAAACGAAAGTGGGCCGGTGACTTTAATCGCATgggg                                                     >  1:307401/1‑99 (MQ=255)                                              caGGAACGCAGAGCAATCCCCTGCCCATCGAGACACGACTGATGGACAATTTATCCGTGATTAGTTGTTAGTGGGTCGGTGGCTTTAATCGCATggggg                                                    >  2:2891764/1‑99 (MQ=255)                                              caGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATgggg                                                     >  2:1571581/1‑98 (MQ=255)                                               aGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGcc                                                  <  1:2346256/100‑1 (MQ=255)                                                ggAACGCAGAGCAATCCCCTGCCCATCGAGACACCTCTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCct                                                 <  1:695592/100‑1 (MQ=255)                                                ggAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCct                                                 <  2:255274/100‑1 (MQ=255)                                                ggAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCct                                                 >  1:687496/1‑100 (MQ=255)                                                   aCGCAGAGCAATCCCCTGCCCATCGAGACACCACTGAGGGACAATATATCCATGATTAGACGACAGTGGGCCGGTGAGTTTAATCGCATGAGGGCcatc                                               >  2:1173252/1‑99 (MQ=255)                                                   aCGCAGAGCAAGCCCCGGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCctt                                                <  2:1745464/98‑1 (MQ=255)                                                     gCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCttt                                            <  2:292411/100‑1 (MQ=255)                                                      cAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTa                                           >  1:1890608/1‑100 (MQ=255)                                                      cAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCACATGGGGGCCTTCTTTa                                           >  1:1594960/1‑100 (MQ=255)                                                        gagCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTAtt                                         >  1:1239843/1‑100 (MQ=255)                                                          gCAATCCCCTGCCCATCGAGACACCACTGATTGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTg                                        <  1:742343/99‑1 (MQ=255)                                                          gCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGc                                       <  1:1529831/100‑1 (MQ=255)                                                          gCAATCCCCTGCCCATCGAGACACCACTGATGGACAAGTTAGCCATGATTACACAACAGTGGGCCGGTGACTTTAATCGCATGGAGGCCTTCTTGATTGc                                       >  1:1648972/1‑100 (MQ=255)                                                             aTCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTa                                    <  1:2946691/100‑1 (MQ=255)                                                                 cttccgatctgagaCACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGCCTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTg                                <  2:1619184/90‑1 (MQ=255)                                                                   tGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACTACAGTGTGCCGGAGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGcc                              >  1:807649/1‑100 (MQ=255)                                                                    gCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCa                             <  1:2438618/100‑1 (MQ=255)                                                                     cccATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAg                            <  2:1685962/100‑1 (MQ=255)                                                                        aTCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGaca                         >  2:2368321/1‑100 (MQ=255)                                                                        aTCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGaca                         <  2:1225948/100‑1 (MQ=255)                                                                          cGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACAcc                       <  2:3037383/100‑1 (MQ=255)                                                                           gagaCACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTgagatcggaag                      >  1:1619203/1‑90 (MQ=255)                                                                              acacCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCaaaa                   <  1:2706582/100‑1 (MQ=255)                                                                               caccacTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAg                  <  2:2713000/100‑1 (MQ=255)                                                                                  cacTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGAGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGgat               >  1:2418442/1‑100 (MQ=255)                                                                                      gATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGgatgat            >  1:2674498/1‑99 (MQ=255)                                                                                      gATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATc           <  2:1203548/100‑1 (MQ=255)                                                                                       ctGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATcg          <  2:1551721/99‑1 (MQ=255)                                                                                       aTGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTCATTGCCAGACACCAAAAGGATGATcg          >  2:1339827/1‑100 (MQ=255)                                                                                       aTGGACAATTTAGCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATcg          <  2:164349/100‑1 (MQ=255)                                                                                        tGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATcgc         >  1:2486909/1‑100 (MQ=255)                                                                                         ggACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATcgcg        >  2:792631/1‑100 (MQ=255)                                                                                         ggACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATcgag        >  1:1551739/1‑98 (MQ=255)                                                                                         ggACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATcg          >  1:845387/1‑98 (MQ=255)                                                                                          gACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTACCAGACACCAAAAGGATGATCgcgc       >  1:636111/1‑100 (MQ=255)                                                                                           aCAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATCgcgc       >  1:2281856/1‑99 (MQ=255)                                                                                            cAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATCGCGctc     <  1:2747357/100‑1 (MQ=255)                                                                                             aaTTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGTGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATCGCGctct    <  2:2321579/100‑1 (MQ=255)                                                                                             aaTTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATCGCGctct    >  1:571114/1‑100 (MQ=255)                                                                                             aaTTTATCCATGATTAGACGACAGTGGGCCGGTGACTCTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATCGCGctct    >  1:1346235/1‑100 (MQ=255)                                                                                              aTTTATCCATGATTAGACGACAGTGGGCCGGTCACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATCGCGCTCtt   >  1:2884029/1‑100 (MQ=255)                                                                                               tttATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATCGCGCTCttt  <  1:2502153/100‑1 (MQ=255)                                                                                                  |                                                                                                 ACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATCGCGCTCTTT  >  NC_000913/1881031‑1881224

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