Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,902,432 G→A P265L (CCG→CTG)  bglF ← PTS system beta‑glucosides‑specific enzyme IIABC component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,902,4320GA95.5% 589.5 / 7.0 176P265L (CCG→CTG) bglFPTS system beta‑glucosides‑specific enzyme IIABC component
Reads supporting (aligned to +/- strand):  ref base G (2/6);  new base A (111/57);  total (113/63)
Fisher's exact test for biased strand distribution p-value = 2.55e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.21e-01

Only 100 of 181 total aligned reads displayed.

CGCCCATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCA  >  NC_000913/3902343‑3902521
                                                                                         |                                                                                         
cgccCATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGcc                                                                                        >  1:2279045/1‑93 (MQ=255)
 gccCATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGccc                                                                                       >  2:687250/1‑93 (MQ=255)
  cccATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGcccc                                                                                      >  2:4162470/1‑93 (MQ=255)
     aTTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCccacc                                                                                   >  1:3408103/1‑93 (MQ=255)
     aTTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCccacc                                                                                   >  1:4308674/1‑93 (MQ=255)
     aTTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCccacc                                                                                   >  2:829094/1‑93 (MQ=255)
     aTTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCccacc                                                                                   >  2:4136365/1‑93 (MQ=255)
      ttACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCccacca                                                                                  >  2:3818117/1‑93 (MQ=255)
      ttACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCccacca                                                                                  >  1:167604/1‑93 (MQ=255)
       tACCGCGCCCGCAAATGCAGGAACCACCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCACCCCCACcag                                                                                 >  1:998230/1‑93 (MQ=255)
         ccGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACcagca                                                                               >  1:3930016/1‑93 (MQ=255)
          cgcgCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACcagcag                                                                              <  1:2348276/93‑1 (MQ=255)
           gcgcCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGa                                                                             >  1:1018798/1‑93 (MQ=255)
           gcgcCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGa                                                                             >  1:1018797/1‑93 (MQ=255)
            cgcCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGaa                                                                            <  1:4134122/93‑1 (MQ=255)
            cgcCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGaa                                                                            >  2:311618/1‑93 (MQ=255)
             gcCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGaaa                                                                           >  1:1673390/1‑93 (MQ=255)
             gcCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGaaa                                                                           <  2:2500565/93‑1 (MQ=255)
               ccGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACGAGCAGAAAgg                                                                         >  1:1500804/1‑93 (MQ=255)
               ccGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAgg                                                                         >  1:2627168/1‑93 (MQ=255)
               ccGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAgg                                                                         >  2:4101794/1‑93 (MQ=255)
                cGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAggg                                                                        >  1:550470/1‑92 (MQ=255)
                cGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGt                                                                        <  1:1546064/93‑1 (MQ=255)
                 gCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATGGGAGCCCCACCAGCAGAAGGGTg                                                                       >  1:4060223/1‑93 (MQ=255)
                 gCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTg                                                                       >  2:4241977/1‑93 (MQ=255)
                  cAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGa                                                                      <  2:3397154/93‑1 (MQ=255)
                   aaaTGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGAc                                                                     >  2:2158762/1‑93 (MQ=255)
                   aaaTGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGAc                                                                     >  1:126627/1‑93 (MQ=255)
                     aTGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACgg                                                                   <  2:932229/93‑1 (MQ=255)
                      tGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACggg                                                                  >  2:4384635/1‑93 (MQ=255)
                      agCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACggg                                                                  <  2:2950825/92‑1 (MQ=255)
                          ggAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGtgtg                                                              >  2:2570148/1‑93 (MQ=255)
                          ggAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGGGACGGGtgtg                                                              >  1:1256468/1‑93 (MQ=255)
                            aaCCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGAt                                                            >  1:3183446/1‑93 (MQ=255)
                             aCCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATa                                                           >  2:3129889/1‑93 (MQ=255)
                             aCCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATa                                                           <  1:1540178/93‑1 (MQ=255)
                              ccGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATaa                                                          >  2:1608778/1‑93 (MQ=255)
                              ccGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATaa                                                          >  1:2151867/1‑93 (MQ=255)
                               cGCCTGATAAAGCCAGAGATAACCGGCGGCAATGAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAAc                                                         <  2:2084391/93‑1 (MQ=255)
                                gCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATaaac                                                        >  1:2335314/1‑91 (MQ=255)
                                gCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAAcc                                                        >  2:1326258/1‑93 (MQ=255)
                                 ccTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCTACCAGCAGAAAGGTGACGGGTGGGATAACCa                                                       >  2:3795955/1‑93 (MQ=255)
                                  cTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCAt                                                      >  1:4531233/1‑93 (MQ=255)
                                   tGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATc                                                     <  1:4019458/93‑1 (MQ=255)
                                   tGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATc                                                     <  2:1593004/93‑1 (MQ=255)
                                   tGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATc                                                     <  2:1123049/93‑1 (MQ=255)
                                   tGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATc                                                     <  1:949597/93‑1 (MQ=255)
                                      tAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCGAGGTTGATAGCAGCCCCACCAGCAGAAAGTTGACGGGTGTGATAACCATCAGa                                                  <  2:2889969/93‑1 (MQ=255)
                                       aaaGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCACCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGAc                                                 <  2:4732054/93‑1 (MQ=255)
                                        aaGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACa                                                <  2:2877756/93‑1 (MQ=255)
                                         aGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACAt                                               >  2:2711431/1‑93 (MQ=255)
                                         aGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACAt                                               >  2:2311659/1‑93 (MQ=255)
                                         aGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACAt                                               <  2:1154660/93‑1 (MQ=255)
                                         aGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACAt                                               <  1:364427/93‑1 (MQ=255)
                                           ccAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAg                                             >  1:1237544/1‑93 (MQ=255)
                                             agagATAACCGGCGGCAATCAGTTCGCTTATCCCGGTTGATAGCAGCCCCACCAGCAGAAAAGTGACGGGTGTGATAACCATCAGACATACCa                                           >  1:3461317/1‑93 (MQ=255)
                                              gagaTAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCaa                                          >  1:3326314/1‑93 (MQ=255)
                                              gagaTAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCACCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCaa                                          >  1:3275769/1‑93 (MQ=255)
                                               agaTAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAAt                                         >  2:2508219/1‑93 (MQ=255)
                                               agaTAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAAt                                         >  1:4174598/1‑93 (MQ=255)
                                               agaTAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAAt                                         <  1:3747525/93‑1 (MQ=255)
                                                 aTAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATgg                                       <  2:2761178/93‑1 (MQ=255)
                                                 aTAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATgg                                       >  2:2668228/1‑93 (MQ=255)
                                                 aTAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATgg                                       >  1:2831508/1‑93 (MQ=255)
                                                 aTAACCGGCGGCAATCAGTTCGCTTATCCAGGGTGATAGCAGCCCCACCAGCAGGAAGGTGACGGGTGTGATAACCATCAGCCATAGCACTgg                                       >  2:2344290/1‑93 (MQ=255)
                                                  tAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGgt                                      >  1:778072/1‑93 (MQ=255)
                                                  tAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGgt                                      >  1:3010797/1‑93 (MQ=255)
                                                  tAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGgt                                      <  1:1752770/93‑1 (MQ=255)
                                                  tAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGgt                                      <  1:3963475/93‑1 (MQ=255)
                                                   aaCCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGgtg                                     >  1:2232559/1‑93 (MQ=255)
                                                   aaCCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGgtg                                     <  1:1135153/93‑1 (MQ=255)
                                                   aaCCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGgtg                                     <  1:4749401/93‑1 (MQ=255)
                                                      cggcggCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTga                                  <  1:669322/93‑1 (MQ=255)
                                                      cggcggCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTga                                  >  1:719128/1‑93 (MQ=255)
                                                       ggcggcAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTgaa                                 >  1:3254416/1‑93 (MQ=255)
                                                       ggcggcAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTgaa                                 <  2:3965284/93‑1 (MQ=255)
                                                          ggcAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTgaagaa                              >  1:3342566/1‑93 (MQ=255)
                                                             aaTCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAAtt                           >  2:1680404/1‑93 (MQ=255)
                                                             aaTCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAAtt                           <  1:4614661/93‑1 (MQ=255)
                                                              aTCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAAttt                          >  1:1207465/1‑93 (MQ=255)
                                                              aTCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAAttt                          >  1:3153599/1‑93 (MQ=255)
                                                              aTCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAAttt                          <  2:208831/93‑1 (MQ=255)
                                                                  gTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGa                      >  2:1885226/1‑93 (MQ=255)
                                                                   ttCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGAt                     >  2:3879731/1‑93 (MQ=255)
                                                                    tCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGAtt                    <  2:4415186/93‑1 (MQ=255)
                                                                       cTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGcc                 >  2:3278512/1‑93 (MQ=255)
                                                                        ttATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCg                >  1:2892575/1‑93 (MQ=255)
                                                                          aTCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGAc              >  1:4650178/1‑93 (MQ=255)
                                                                           tCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACg             >  2:3218095/1‑93 (MQ=255)
                                                                             cAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGt           >  2:1483723/1‑93 (MQ=255)
                                                                             cAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGt           <  2:3914418/93‑1 (MQ=255)
                                                                              aGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTa          >  2:3599400/1‑93 (MQ=255)
                                                                                gTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAAc        >  2:2258212/1‑93 (MQ=255)
                                                                                 ttGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAAcc       >  2:140117/1‑93 (MQ=255)
                                                                                   gATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCAc     >  1:4557300/1‑93 (MQ=255)
                                                                                   gATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCAc     >  1:2916559/1‑93 (MQ=255)
                                                                                    aTAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACg    >  2:4102819/1‑93 (MQ=255)
                                                                                    aTAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACg    >  1:1690938/1‑93 (MQ=255)
                                                                                    aTAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACg    <  2:1541864/93‑1 (MQ=255)
                                                                                      agcagcCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCa  >  1:741866/1‑93 (MQ=255)
                                                                                         |                                                                                         
CGCCCATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCA  >  NC_000913/3902343‑3902521

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