Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,098,512	C→A	100%	R259R (CGC→CGA)	rhaS →	transcriptional activator of rhaBAD and rhaT

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,098,512	0	C	A	98.9%	304.3 / ‑3.9	94	R259R (CGC→CGA)	rhaS	transcriptional activator of rhaBAD and rhaT
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (43/49);  total (43/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.01e-01

AAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCATATTGC  >  NC_000913/4098415‑4098608                                                                                                  |                                                                                                 aaGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACg                                                                                                <  2:2067715/100‑1 (MQ=255)   gCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACgag                                                                                              <  1:1209981/100‑1 (MQ=255)   gCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACgag                                                                                              <  2:1866388/100‑1 (MQ=255)   gCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACgag                                                                                              <  1:1396948/100‑1 (MQ=255)   gCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACgag                                                                                              <  2:900973/100‑1 (MQ=255)        aCATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGttt                                                                                         <  1:2483727/100‑1 (MQ=255)         cATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTa                                                                                        <  1:1337320/100‑1 (MQ=255)           tCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAAc                                                                                      >  1:312309/1‑100 (MQ=255)            cTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACt                                                                                     <  2:1774387/100‑1 (MQ=255)             tGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTg                                                                                    <  1:1840610/100‑1 (MQ=255)             tGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTg                                                                                    >  1:188810/1‑100 (MQ=255)              gCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTgg                                                                                   >  2:2256339/1‑100 (MQ=255)                tACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTTTGTATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTgg                                                                                   <  2:2547672/98‑1 (MQ=255)                 actccagaccGAGGCGAGCGGTACTGAGATTGCCTATGGTTGTGGATTCAGCGACAGTAACCACTTCGCGACGCTTTTTCGTCGAGAGTGTAACTGGTc                                                                                 <  2:2395939/90‑1 (MQ=255)                  cgccgcGGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCAc                                                                               <  2:760958/100‑1 (MQ=255)                  cGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCCTTTTCGACGAGAGTTTAACTGGTCAc                                                                               >  2:344563/1‑100 (MQ=255)                   gCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGCTTAACTGGTCAcc                                                                              >  2:2395426/1‑100 (MQ=255)                    ccACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACcg                                                                             >  1:214630/1‑100 (MQ=255)                    ccACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACcg                                                                             <  1:739271/100‑1 (MQ=255)                     cacaGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACcgc                                                                            >  2:790280/1‑100 (MQ=255)                       caGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGt                                                                          >  2:730493/1‑100 (MQ=255)                          cGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGat                                                                       <  1:1410955/100‑1 (MQ=255)                          cGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGATAAATAAACGAGAGTTTAACTGGTCACCGCGTGat                                                                       >  1:466866/1‑100 (MQ=255)                           gAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGata                                                                      <  1:803711/100‑1 (MQ=255)                           gAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGata                                                                      <  1:801661/100‑1 (MQ=255)                             ggCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCATCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATAtt                                                                    <  2:1391660/100‑1 (MQ=255)                             ggCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATAtt                                                                    <  2:2108409/100‑1 (MQ=255)                              gCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTc                                                                   <  1:284319/100‑1 (MQ=255)                               ccAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCg                                                                  >  1:1028247/1‑100 (MQ=255)                               ccAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCg                                                                  >  2:559993/1‑100 (MQ=255)                                cAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGc                                                                 >  2:795756/1‑100 (MQ=255)                                cAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGc                                                                 >  2:795651/1‑100 (MQ=255)                                  gCGTTACTGACATCGCCTATCGATGTGGGCTCAGCGACAGTAACCACTTTTCGACGCTTTTTGTACGAGAGGTTAACTGGTCACCGCGGGATATTCGCCa                                                               >  2:1616822/1‑100 (MQ=255)                                    gTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAgg                                                             >  2:716424/1‑100 (MQ=255)                                        cTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACg                                                         <  2:1096787/100‑1 (MQ=255)                                         tGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACgg                                                        <  2:214628/100‑1 (MQ=255)                                          naCATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACggg                                                       >  1:410890/2‑100 (MQ=255)                                          gACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACggg                                                       >  1:1910033/1‑100 (MQ=255)                                          gACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACggg                                                       <  2:925047/100‑1 (MQ=255)                                          gACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACggg                                                       >  2:2348969/1‑100 (MQ=255)                                            cATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGAt                                                     >  2:2356827/1‑100 (MQ=255)                                               cGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACTAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGc                                                  <  1:1524883/100‑1 (MQ=255)                                               cGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGc                                                  >  1:1021135/1‑100 (MQ=255)                                               cGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGc                                                  <  2:1677473/100‑1 (MQ=255)                                                gCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCt                                                 >  2:647568/1‑100 (MQ=255)                                                gCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCt                                                 >  2:2485747/1‑100 (MQ=255)                                                 ccTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCtt                                                <  2:1480370/100‑1 (MQ=255)                                                    aTCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCt                                             <  2:1930240/100‑1 (MQ=255)                                                      cGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGc                                           <  1:931870/100‑1 (MQ=255)                                                      cGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGc                                           <  2:417/100‑1 (MQ=255)                                                      cGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGc                                           <  1:805034/100‑1 (MQ=255)                                                      cGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGc                                           <  1:139125/100‑1 (MQ=255)                                                      cGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGc                                           <  2:684292/100‑1 (MQ=255)                                                        cTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCaa                                         <  2:126713/100‑1 (MQ=255)                                                        cTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCaa                                         >  1:2363677/1‑100 (MQ=255)                                                        cTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCaa                                         >  1:2443081/1‑100 (MQ=255)                                                         tgtgGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAAt                                        <  2:187761/100‑1 (MQ=255)                                                             gATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAAcg                                    >  1:1959191/1‑100 (MQ=255)                                                             gATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAAcg                                    <  1:2387944/100‑1 (MQ=255)                                                             gATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAAcg                                    <  1:800261/100‑1 (MQ=255)                                                              aTTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAAcgc                                   <  2:948587/100‑1 (MQ=255)                                                              aTTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAAcgc                                   >  1:1577230/1‑100 (MQ=255)                                                               ttCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAAcgcg                                  <  1:554988/100‑1 (MQ=255)                                                               ttCAGCGACAGTAACCACTTTTCGACACTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAAcgcg                                  >  1:654790/1‑100 (MQ=255)                                                                 cAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGACGGCTTTCTGCAATAAcgcg                                  >  2:2038267/1‑98 (MQ=255)                                                                  tccgatcttAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAAt                               <  1:1249107/92‑1 (MQ=255)                                                                  aGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAAt                               <  2:1788158/100‑1 (MQ=255)                                                                  aGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAAt                               <  1:2251816/100‑1 (MQ=255)                                                                  aGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATCCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAAt                               >  2:1898372/1‑100 (MQ=255)                                                                    cGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAAtct                             <  2:1466883/100‑1 (MQ=255)                                                                      aCAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAAtcttc                           <  1:2122687/100‑1 (MQ=255)                                                                       cAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAAtcttct                          <  1:409168/100‑1 (MQ=255)                                                                       cAGTAACCACTTTTCGACACTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAAtcttct                          <  2:654783/100‑1 (MQ=255)                                                                         gTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCa                        <  1:74721/100‑1 (MQ=255)                                                                         gTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCa                        >  2:1176070/1‑100 (MQ=255)                                                                         gTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCa                        >  2:1426410/1‑100 (MQ=255)                                                                         gTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTACTGCAATAACGCGAATCTTCTCa                        >  2:1929794/1‑100 (MQ=255)                                                                          tAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAAtagatcgg                        >  2:1249099/1‑92 (MQ=255)                                                                           aaCCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATTGCTTTTTTCAATAACGCTAATCTTTTCaa                       >  2:2147348/1‑99 (MQ=255)                                                                           aaCCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGCATCTTCTCAAc                      >  2:1953014/1‑100 (MQ=255)                                                                            aCCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACg                     <  1:1876433/100‑1 (MQ=255)                                                                            aCCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACg                     <  1:456883/100‑1 (MQ=255)                                                                             ccACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGt                    <  2:413034/100‑1 (MQ=255)                                                                                  tttCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTg               >  2:1253692/1‑100 (MQ=255)                                                                                  tccgatctcTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTg               <  2:1450997/92‑1 (MQ=255)                                                                                     cGACGCTTTTTCGACGCGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACCCGAATCTTCTCAACGTATTTGTAc            >  2:1779432/1‑100 (MQ=255)                                                                                     cGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTAc            >  1:196604/1‑100 (MQ=255)                                                                                       aCGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGc          >  1:359364/1‑100 (MQ=255)                                                                                       aCGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGc          <  1:380637/100‑1 (MQ=255)                                                                                         gCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCa        <  1:2046643/100‑1 (MQ=255)                                                                                          cTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTgagatcgga       >  1:1451008/1‑92 (MQ=255)                                                                                          cTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCat       >  2:1674013/1‑100 (MQ=255)                                                                                              atCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCATATTg   >  2:567328/2‑100 (MQ=255)                                                                                               tCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTCCGCACTATTGc  >  2:1059680/1‑100 (MQ=255)                                                                                                  |                                                                                                 AAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCATATTGC  >  NC_000913/4098415‑4098608

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