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	99.7%	303.3 / ‑6.7	94	E228* (GAA→TAA)	dmlR	transcriptional activator of dmlA
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base A (49/43);  total (51/43)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATCGCGCTC  >  NC_000913/1881031‑1881221                                                                                                  |                                                                                              aCTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTtatc                                                                                             <  1:2237137/100‑1 (MQ=255)  cTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATcc                                                                                            >  2:1280905/1‑100 (MQ=255)     gCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATg                                                                                         >  1:1090657/1‑100 (MQ=255)       aCTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGAt                                                                                       <  1:1828767/100‑1 (MQ=255)         tAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTa                                                                                     >  1:2469949/1‑100 (MQ=255)          aaaTGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAg                                                                                    >  1:1052341/1‑100 (MQ=255)           aaTTACCACTGGTAATGTTTTCGCTAATATCCCACCTGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAg                                                                                    <  1:1491686/99‑1 (MQ=255)            aTGACCACTGGCAGTGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAgac                                                                                  <  2:1283005/100‑1 (MQ=255)             tGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAgacg                                                                                 >  1:2223951/1‑100 (MQ=255)               aCCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAgacgac                                                                               <  1:77094/100‑1 (MQ=255)               aCCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAgacgac                                                                               <  1:1618981/100‑1 (MQ=255)               aCCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAgacgac                                                                               <  1:1618903/100‑1 (MQ=255)                  aCTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACCCCACTGACGGACAATTTCTCCATGATTAGACGACAg                                                                             >  2:774069/1‑99 (MQ=255)                  aCTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGt                                                                            <  1:751856/100‑1 (MQ=255)                   cTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTg                                                                           <  2:469142/100‑1 (MQ=255)                    tGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTgg                                                                          <  1:1362562/100‑1 (MQ=255)                       cAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGcc                                                                       >  1:1939639/1‑100 (MQ=255)                       cAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGcc                                                                       >  1:723091/1‑100 (MQ=255)                         aTGTTTTCGCTAACATCCCACCATGAACGCATAGCAATCTCCTGACCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCg                                                                      <  2:614665/99‑1 (MQ=255)                           gTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTg                                                                   <  2:1834948/100‑1 (MQ=255)                           gTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTg                                                                   >  1:1836550/1‑100 (MQ=255)                           gTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTg                                                                   <  1:608914/100‑1 (MQ=255)                           gTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATACCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAAACGACAGTg                                                                           >  2:1049489/1‑92 (MQ=255)                             tttCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGAc                                                                 <  1:196875/100‑1 (MQ=255)                              ttCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCTTGATTAGACGACAGTGGGCCGGt                                                                    >  2:2551708/1‑96 (MQ=255)                                 gCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTa                                                             >  2:1235367/1‑100 (MQ=255)                                  cTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGGCAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTaa                                                            >  2:2338187/1‑100 (MQ=255)                                  cTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACttt                                                              >  2:1706889/1‑98 (MQ=255)                                  cTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTaa                                                            >  1:1765023/1‑100 (MQ=255)                                    aaCATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATc                                                          <  2:205681/100‑1 (MQ=255)                                    aaCATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATc                                                          <  2:1052333/100‑1 (MQ=255)                                      cATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTATACGACAGTGGGCCGGTGACTTTAATCGc                                                        >  2:697048/1‑100 (MQ=255)                                      cATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGc                                                        <  2:910235/100‑1 (MQ=255)                                       aTCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCa                                                       >  1:14813/1‑100 (MQ=255)                                         cccACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGTTTCGACGACAGTGCGCGAGACACTATAATCGCATg                                                     >  2:2283528/1‑100 (MQ=255)                                          ccaccaGGAACTCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACGTTAATCGCa                                                       <  2:881648/97‑1 (MQ=255)                                          ccaccaGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATgg                                                    >  1:920219/1‑100 (MQ=255)                                          ccaccaGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATgg                                                    >  1:1469364/1‑100 (MQ=255)                                           caccaGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATggg                                                   >  2:1892696/1‑100 (MQ=255)                                             ccaGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGAGCCGGTGACTTTAATCGCATgg                                                    >  2:1952289/1‑97 (MQ=255)                                              caGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGc                                                >  1:2289027/1‑100 (MQ=255)                                              caGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGc                                                >  1:1221259/1‑100 (MQ=255)                                                 gAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGAGGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCctt                                             >  2:1422818/1‑100 (MQ=255)                                                  aaCGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATgggggn                                                <  2:1562921/96‑2 (MQ=255)                                                  aaCGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCcttc                                            <  1:1235375/100‑1 (MQ=255)                                                   aCGCAGAGCAATCCCCTGCCCATCGAGACACCACTCATGGACAATTTATCCATGACTAGACGCCACTTGGCCGGTGACCTTAATCGCATAAGGGCcttct                                           >  2:1115429/1‑100 (MQ=255)                                                    cGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCcttctt                                          <  1:1390700/100‑1 (MQ=255)                                                    cGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCcttctt                                          >  1:1612359/1‑100 (MQ=255)                                                    cGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCcttctt                                          <  1:1276395/100‑1 (MQ=255)                                                     gCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCttt                                         >  1:801499/1‑100 (MQ=255)                                                     gCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCttt                                         >  1:393744/1‑100 (MQ=255)                                                      cAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTa                                        >  2:1387795/1‑100 (MQ=255)                                                       agagCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTAt                                       <  1:1399898/100‑1 (MQ=255)                                                        gagCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTAtt                                      <  1:1969360/100‑1 (MQ=255)                                                         agCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTg                                     <  2:172399/100‑1 (MQ=255)                                                         agCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTg                                     <  2:2081026/100‑1 (MQ=255)                                                         agCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTg                                     <  1:1589622/100‑1 (MQ=255)                                                         agCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTg                                     <  1:771325/100‑1 (MQ=255)                                                          gCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTAtt                                      >  1:699581/1‑98 (MQ=255)                                                          gCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGc                                    <  2:1899729/100‑1 (MQ=255)                                                            aaTCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGCCAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGt                                  <  2:1275048/100‑1 (MQ=255)                                                            aaTCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGt                                  <  2:1384785/100‑1 (MQ=255)                                                            aaTCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGt                                  <  2:1765010/100‑1 (MQ=255)                                                               ccccTGCCCATCGAGACACCACTGATGGACAATTTCTCCATCATTAGACGCCAGTGGGCCGGTGACTTTAATCGCTTGGCGGCCTTCTTTATTGCg                                   >  2:2018125/1‑96 (MQ=255)                                                               ccccTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAAt                               >  1:685785/1‑100 (MQ=255)                                                               ccccTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGCCCGGTGACTGTAATCGCATGGGGGCCTTCTTTATTGCGTAAt                               >  2:1799648/1‑100 (MQ=255)                                                                cccTGCCCATCGAGACACCACTGATGGCCAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAAtt                              >  2:903998/1‑100 (MQ=255)                                                                cccTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAAtt                              >  1:1487157/1‑100 (MQ=255)                                                                 ccTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTg                             <  1:1978564/100‑1 (MQ=255)                                                                 ccTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTg                             >  2:1341285/1‑100 (MQ=255)                                                                  cTGCCCATCGAGACACCACTGATGGAGAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGc                            <  2:2130536/100‑1 (MQ=255)                                                                   tGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGcc                           <  2:2140791/100‑1 (MQ=255)                                                                   tGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGcc                           <  1:1103819/100‑1 (MQ=255)                                                                    gCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCa                          >  1:1122554/1‑100 (MQ=255)                                                                    gCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCa                          <  2:2469935/100‑1 (MQ=255)                                                                       cATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGac                       <  2:1969570/100‑1 (MQ=255)                                                                        aTCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGaca                      >  1:2157655/1‑100 (MQ=255)                                                                          cGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACAcc                    <  1:2040670/100‑1 (MQ=255)                                                                          cGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACAcc                    >  2:1869370/1‑100 (MQ=255)                                                                           gagaCACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCa                   >  1:1299564/1‑100 (MQ=255)                                                                               caccacTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAg               <  2:1939625/100‑1 (MQ=255)                                                                               caccacTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAg               <  1:1554611/100‑1 (MQ=255)                                                                                accacTGATGGACAATTTATCCATGATTAGACGACAGCGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAgg              >  2:814916/1‑100 (MQ=255)                                                                                 ccacTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGga             >  1:2366213/1‑100 (MQ=255)                                                                                  cacTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGgat            >  2:262924/1‑100 (MQ=255)                                                                                     tGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGgatgat         <  1:346886/100‑1 (MQ=255)                                                                                      gATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTGTTTATTGCGTAATTGCCAGACACCAAAAGga             >  2:1097760/1‑95 (MQ=255)                                                                                       aTGGACAGTTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATcg       >  1:1341667/1‑100 (MQ=255)                                                                                        tGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATcgc      >  2:632056/1‑100 (MQ=255)                                                                                         ggACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATcgcg     >  1:1756221/1‑100 (MQ=255)                                                                                         ggACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATcgcg     >  1:1952779/1‑100 (MQ=255)                                                                                         ggACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATcgcg     >  1:2133103/1‑100 (MQ=255)                                                                                          gACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATCgcgc    <  1:68571/100‑1 (MQ=255)                                                                                          gACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATCgcgc    >  1:816107/1‑100 (MQ=255)                                                                                            cAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATCGCGctc  <  2:393740/100‑1 (MQ=255)                                                                                                  |                                                                                              ACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATCGCGCTC  >  NC_000913/1881031‑1881221

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