Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,098,512 C→A 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_0009134,098,5120CA96.8% 297.7 / ‑5.5 95R259R (CGC→CGA) rhaStranscriptional activator of rhaBAD and rhaT
Reads supporting (aligned to +/- strand):  ref base C (0/2);  major base A (50/40);  minor base T (2/1);  total (52/43)
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

AAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCATATTGC  >  NC_000913/4098415‑4098608
                                                                                                 |                                                                                                
aaGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACg                                                                                                >  1:1209290/1‑100 (MQ=255)
  gCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACgag                                                                                              >  2:535547/1‑100 (MQ=255)
   cccGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACgaga                                                                                             >  2:943596/1‑100 (MQ=255)
   cccGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACgaga                                                                                             <  1:2720985/100‑1 (MQ=255)
    ccGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGagag                                                                                            <  1:2872946/100‑1 (MQ=255)
    ccGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGagag                                                                                            <  1:745135/100‑1 (MQ=255)
      gACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGtt                                                                                          >  1:3004387/1‑100 (MQ=255)
       aCATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGttt                                                                                         >  1:2689510/1‑100 (MQ=255)
           cTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCATTTTTCGACGCTTTTTCGTCGAGAGTTTAACt                                                                                     <  2:2337983/100‑1 (MQ=255)
           cTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACt                                                                                     >  1:2935659/1‑100 (MQ=255)
             gCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTgg                                                                                   <  1:2854750/100‑1 (MQ=255)
             gCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCGATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTgg                                                                                   >  1:1856855/1‑100 (MQ=255)
              cTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGCTTCCGCGACAGCAACCACTTTTCGACGCTTTTTCGACGAGCGTTTAACTGGt                                                                                  >  2:1935996/1‑100 (MQ=255)
                gtgCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCa                                                                                >  1:2454159/3‑100 (MQ=255)
                 cGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCa                                                                                >  2:2916335/1‑99 (MQ=255)
                  gCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCAcc                                                                              >  1:792504/1‑100 (MQ=255)
                     acaGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACcgcg                                                                           >  2:217100/1‑100 (MQ=255)
                     acaGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACcgcg                                                                           <  1:770434/100‑1 (MQ=255)
                     acaGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACcgcg                                                                           <  2:499488/100‑1 (MQ=255)
                     acaGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACcgcg                                                                           <  1:33784/100‑1 (MQ=255)
                      caGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGt                                                                          >  2:71034/1‑100 (MQ=255)
                       aGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTg                                                                         <  1:2391682/100‑1 (MQ=255)
                       aGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTg                                                                         <  1:2621016/100‑1 (MQ=255)
                        gCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTg                                                                         >  2:3439730/1‑99 (MQ=255)
                        gCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGa                                                                        <  2:3029394/100‑1 (MQ=255)
                        gCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGa                                                                        <  2:1243008/100‑1 (MQ=255)
                        gCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGa                                                                        <  2:789070/100‑1 (MQ=255)
                         cGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGat                                                                       <  1:381197/100‑1 (MQ=255)
                         cGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGat                                                                       <  1:2424739/100‑1 (MQ=255)
                          gAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTg                                                                         <  1:957430/97‑1 (MQ=255)
                             gCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTg                                                                         <  1:956678/94‑1 (MQ=255)
                             gCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTc                                                                   >  2:1769451/1‑100 (MQ=255)
                              ccAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCg                                                                  >  2:199964/1‑100 (MQ=255)
                                aGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGcc                                                                <  2:1566905/100‑1 (MQ=255)
                                aGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGcc                                                                <  1:3059132/100‑1 (MQ=255)
                                aGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGcc                                                                <  2:2470203/100‑1 (MQ=255)
                                aGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGcc                                                                <  2:2982394/100‑1 (MQ=255)
                                     tACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGa                                                           >  2:697433/1‑100 (MQ=255)
                                       cTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACg                                                         >  2:990479/1‑100 (MQ=255)
                                        tGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACgg                                                        <  1:1587713/100‑1 (MQ=255)
                                         gACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTTGTCGATAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACggg                                                       >  1:1661428/1‑100 (MQ=255)
                                         gACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACggg                                                       <  1:3043378/100‑1 (MQ=255)
                                           cATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTTTTTTTTTTTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGAt                                                     >  1:795366/1‑100 (MQ=255)
                                           cATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGAt                                                     <  1:3299703/100‑1 (MQ=255)
                                            aTCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATg                                                    >  1:1690536/1‑100 (MQ=255)
                                             tCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATgg                                                   >  2:1447925/1‑100 (MQ=255)
                                             tCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCAACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATgn                                                   <  1:1897562/100‑2 (MQ=255)
                                              cGCCTATCGCTGTGGATTCAGCGACAGTAACCCCTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGc                                                  >  1:567252/1‑100 (MQ=255)
                                               gCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGTGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCt                                                 >  1:440803/1‑100 (MQ=255)
                                                ncTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCtt                                                >  1:1216902/2‑100 (MQ=255)
                                                 cTATCGCTTTTGATTCAGCGACATTAACCACTTTTCGACTCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGTGACGGGATGGCttt                                               <  1:478388/100‑1 (MQ=255)
                                                  tATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTc                                              <  1:1239978/100‑1 (MQ=255)
                                                  tATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTc                                              >  1:3112444/1‑100 (MQ=255)
                                                   aTCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCCCCGCGTGATATTCGCCAGGGACGGGATGGCTTTCt                                             >  1:699452/1‑100 (MQ=255)
                                                   aTCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCt                                             >  1:2911628/1‑100 (MQ=255)
                                                   aTCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCt                                             >  2:346532/1‑100 (MQ=255)
                                                    tCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTg                                            <  1:1447903/100‑1 (MQ=255)
                                                     cGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGc                                           <  2:2367859/100‑1 (MQ=255)
                                                     cGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGc                                           <  1:199961/100‑1 (MQ=255)
                                                      gCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCa                                          <  1:3337536/100‑1 (MQ=255)
                                                        tgtgGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAAt                                        <  2:1390305/100‑1 (MQ=255)
                                                        tgtgGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAAt                                        <  2:496223/100‑1 (MQ=255)
                                                        tgtgGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAAt                                        >  2:1088086/1‑100 (MQ=255)
                                                           ggATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAAc                                     >  1:2556308/1‑100 (MQ=255)
                                                             aTTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAAcgc                                   >  2:561942/1‑100 (MQ=255)
                                                              ttCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAAcgcg                                  >  1:142049/1‑100 (MQ=255)
                                                                cAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGaa                                >  2:621329/1‑100 (MQ=255)
                                                                cAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGaa                                >  1:1526711/1‑100 (MQ=255)
                                                                cAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGaa                                >  2:685030/1‑100 (MQ=255)
                                                                cAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGaa                                >  1:3283557/1‑100 (MQ=255)
                                                                 aGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAAt                               >  1:531191/1‑100 (MQ=255)
                                                                 aGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGGTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAAt                               >  1:3149098/1‑100 (MQ=255)
                                                                   cGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAAtct                             >  2:333427/1‑100 (MQ=255)
                                                                    gACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAAtctt                            >  2:2843044/1‑100 (MQ=255)
                                                                      cAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAAtcttct                          <  2:3476530/100‑1 (MQ=255)
                                                                       aGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCAAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTtctc                         >  2:3091785/1‑100 (MQ=255)
                                                                        gTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCa                        <  1:708028/100‑1 (MQ=255)
                                                                          aaCCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAAc                      <  1:2593687/100‑1 (MQ=255)
                                                                             cACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTa                   >  2:2729627/1‑100 (MQ=255)
                                                                              aCTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTAt                  >  2:904324/1‑100 (MQ=255)
                                                                               cTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTAtt                 >  2:2615792/1‑100 (MQ=255)
                                                                               cTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTAtn                 <  2:773721/100‑2 (MQ=255)
                                                                               cTTTTCGACGCTCTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTAtt                 >  2:2282113/1‑100 (MQ=255)
                                                                                 tttCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTg               <  2:933578/100‑1 (MQ=255)
                                                                                 tttCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTg               <  1:431822/100‑1 (MQ=255)
                                                                                    cGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTAc            <  2:101490/100‑1 (MQ=255)
                                                                                     gACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACg           >  1:238391/1‑100 (MQ=255)
                                                                                     gACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACg           <  2:139014/100‑1 (MQ=255)
                                                                                     gACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACg           <  1:832292/100‑1 (MQ=255)
                                                                                        gCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCa        >  2:693195/1‑100 (MQ=255)
                                                                                         cTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCat       >  1:722226/1‑100 (MQ=255)
                                                                                         cTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCat       >  1:3126826/1‑100 (MQ=255)
                                                                                           ttttCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCatat     >  1:1416319/1‑100 (MQ=255)
                                                                                           ttttCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCat       <  1:1524011/98‑1 (MQ=255)
                                                                                              tCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCATATTGc  >  1:1299591/1‑100 (MQ=255)
                                                                                                 |                                                                                                
AAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCATATTGC  >  NC_000913/4098415‑4098608

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