Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,971,690 A→C 100% Q144P (CAG→CCG)  wecC → UDP‑N‑acetyl‑D‑mannosaminuronic acid dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,971,6900AC97.9% 319.0 / ‑2.3 95Q144P (CAG→CCG) wecCUDP‑N‑acetyl‑D‑mannosaminuronic acid dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base C (47/45);  total (47/48)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

GCTGGTGATCCTTGAATCCACCTCGCCGGTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCAGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGT  >  NC_000913/3971598‑3971786
                                                                                            |                                                                                                
gcTGGTGATCCTTGAATCCACCTCGCCGGTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTg                                                                                           <  1:119668/100‑1 (MQ=255)
    gTGATCCTTGAATCCACCTCGCCGGTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCg                                                                                       <  1:394545/100‑1 (MQ=255)
    gTGATCCTTGAATCCACCTCGCCGGTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCg                                                                                       <  1:539562/100‑1 (MQ=255)
    gTGATCCTTGAATCCACCTCGCCGGTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCg                                                                                       <  1:1712718/100‑1 (MQ=255)
    gTGATCCTTGAATCCACCTCGCCGGTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCg                                                                                       <  2:2339224/100‑1 (MQ=255)
     tGATCCTTGAATCCACCTCGCCGGTGTGGTCAACCGCGAAGACGGCAGCACGGGTAGCAGAGAGTCGTCCGCCCCTCACTTCCCCGacgcag                                                                                              >  2:2218210/1‑92 (MQ=255)
     tGATCCTTGAATCCACCTCGCCGGTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGa                                                                                      >  1:968891/1‑100 (MQ=255)
      gATCGTTGAGTCCACCTCGCCGGTGGGGTCAACCGAGAGGATGGCAGCATGGGTAGCAGGGATGGGTCCGGTTCTCACTTATCCGACGCAggggtgc                                                                                        >  2:59442/1‑92 (MQ=255)
      gATCCTTGAATCCACCTCGCCGGTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAg                                                                                     <  1:812822/100‑1 (MQ=255)
      gATCCTTGAATCCACCTCGCCGGTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCACAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAg                                                                                     <  2:1400883/100‑1 (MQ=255)
      gATCCTTGAATCCACCTCGCCGGTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCACAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAg                                                                                     <  2:1400919/100‑1 (MQ=255)
       aTCCTTGAATCCACCTCGCCGGTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCg                                                                                       >  2:1872903/1‑97 (MQ=255)
        tCCTTGAATCCACCTCGCCGGTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCa                                                                                   <  1:868759/100‑1 (MQ=255)
        tCCTTGAATCCACCTCGCCGGTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCa                                                                                   <  2:1030059/100‑1 (MQ=255)
        tCCTTGAATCCACCTCGCCGGTGGGGTCAACCGAGAAGATGGCACAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCTCAGGTGGGCGAGc                                                                                    <  2:1315338/99‑1 (MQ=255)
         ccTTGAATCCACCTCGCCGGTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAg                                                                                  >  2:508238/1‑100 (MQ=255)
          cTTGAATCCACCTCGCCGGTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAgg                                                                                 >  1:1524217/1‑100 (MQ=255)
           ttGATTCCACCTCGCCGGTGGGGTCAGCCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTAGGCGAGCAGGc                                                                                <  2:1167579/100‑1 (MQ=255)
           ttGAATCCACCTCGCCGGTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGc                                                                                <  1:2337129/100‑1 (MQ=255)
           ttGAATCCACCTCGCCGGTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGc                                                                                <  1:1290079/100‑1 (MQ=255)
           ttGAATCCACCTCGCCGGTGGGGTCAACCGAGAAGAGGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGc                                                                                <  1:784668/100‑1 (MQ=255)
             gAATCCACCTCGCCGGTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCgg                                                                              >  2:370586/1‑100 (MQ=255)
               aTCCACCTCGCCGGTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGa                                                                             <  1:362025/99‑1 (MQ=255)
                tCCACCTCGCCTGTTGGGTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACg                                                                           <  2:426302/100‑1 (MQ=255)
                nccACCTCGCCGGTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACg                                                                           >  1:571076/2‑100 (MQ=255)
                  cACCTCGCCGGTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATGGCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTc                                                                         <  2:2121532/100‑1 (MQ=255)
                    ccTCGCCGGTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCaa                                                                       >  1:1481566/1‑100 (MQ=255)
                       cGCCGGTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACAt                                                                    >  1:2199482/1‑100 (MQ=255)
                       cGCCGGTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCCCACTTTCCCGCCGCAGGTGGGCGCGCAGGCGGACGTCAACAt                                                                    >  2:1794592/1‑100 (MQ=255)
                            gTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGGTCTCCCTGTCCCGCCGCAGGTGGGCGACCAGGCGGACGTCGACATTGCt                                                                >  2:1635527/1‑99 (MQ=255)
                               gggTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACt                                                            >  1:1047554/1‑100 (MQ=255)
                                ggTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTg                                                           >  1:1967257/1‑100 (MQ=255)
                                ggTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTg                                                           >  1:424870/1‑100 (MQ=255)
                                 gTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGc                                                          >  1:737654/1‑100 (MQ=255)
                                 gTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGc                                                          <  2:968883/100‑1 (MQ=255)
                                   cAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGccc                                                        >  1:996629/1‑100 (MQ=255)
                                   cAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGccc                                                        <  2:1462584/100‑1 (MQ=255)
                                   cAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGccc                                                        <  1:2444756/100‑1 (MQ=255)
                                   cAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGcc                                                         >  2:1838207/1‑99 (MQ=255)
                                   cAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGGGGGCGAGCAGGCGGACGTCAACATGGCTTACTGcc                                                         >  2:1140103/1‑99 (MQ=255)
                                     aCCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTGACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCgg                                                      >  2:1523225/1‑100 (MQ=255)
                                     aCCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCgg                                                      >  1:1110737/1‑100 (MQ=255)
                                     aCCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCgg                                                      >  1:34869/1‑100 (MQ=255)
                                     aCCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCgg                                                      <  2:1482355/100‑1 (MQ=255)
                                     aCCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCg                                                       <  2:2045706/99‑1 (MQ=255)
                                      ccGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGAAGTGAACATCGCTTACGGCCCGGa                                                     >  2:1845109/1‑100 (MQ=255)
                                       cggtaagaTGGCAGATTGGTGAGCAGAGATGCCGCCCGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCgg                                                      <  2:2295498/94‑1 (MQ=255)
                                       cGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGaa                                                    <  1:289722/100‑1 (MQ=255)
                                         agaagaTGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAAcg                                                  >  1:2246493/1‑100 (MQ=255)
                                           aagaTGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAAcgcg                                                <  1:276513/100‑1 (MQ=255)
                                              aTGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCgtgt                                             <  1:996648/100‑1 (MQ=255)
                                               tGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGtt                                            >  1:790484/1‑100 (MQ=255)
                                                 gCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTAc                                          >  1:1124677/1‑100 (MQ=255)
                                                 gCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTAc                                          >  1:1076918/1‑100 (MQ=255)
                                                  cAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTAcc                                         >  1:689800/1‑100 (MQ=255)
                                                   aGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCa                                        <  2:29600/100‑1 (MQ=255)
                                                    gAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAg                                       <  2:1939847/100‑1 (MQ=255)
                                                    gAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCAGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAg                                       <  1:375093/100‑1 (MQ=255)
                                                     aaTGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAgg                                      <  1:370590/100‑1 (MQ=255)
                                                      aTGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGa                                     >  1:1675198/1‑100 (MQ=255)
                                                        ggTTAGCAGGGATGCGCCCGGATCTCACTTTCCAGCGGCAGGTGGGCGAGCAGGCGGGCGTCAACATAGCTTACTGCCCGGAACGCGTGTTACCAGGACa                                   <  2:72427/100‑1 (MQ=255)
                                                         gTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGAcn                                   <  2:72230/99‑2 (MQ=255)
                                                         gTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAg                                  >  2:1223411/1‑100 (MQ=255)
                                                           tAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGt                                >  1:1626803/1‑100 (MQ=255)
                                                            aGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGCGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTa                               >  2:846817/1‑100 (MQ=255)
                                                             gCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGCGTTACCAGGACAGGTaa                              >  2:360852/1‑100 (MQ=255)
                                                              cAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTaa                              >  2:1711225/1‑99 (MQ=255)
                                                               agagATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATg                            <  2:2478671/100‑1 (MQ=255)
                                                                gagaTGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATgg                           <  1:2192064/100‑1 (MQ=255)
                                                                 agaTGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGt                          <  2:1808713/100‑1 (MQ=255)
                                                                  gaTGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTc                         >  1:1030180/1‑100 (MQ=255)
                                                                    tGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCgn                       <  1:920519/100‑2 (MQ=255)
                                                                     gCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAg                      <  1:2001198/100‑1 (MQ=255)
                                                                      tccgatcttctcACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGc                     <  2:321136/92‑1 (MQ=255)
                                                                      cGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGc                     >  1:1310543/1‑100 (MQ=255)
                                                                       gTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCt                    >  1:389678/1‑100 (MQ=255)
                                                                          cGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGAt                 <  1:1523236/100‑1 (MQ=255)
                                                                          cGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGAt                 <  1:457456/100‑1 (MQ=255)
                                                                          cGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGAt                 <  1:777307/100‑1 (MQ=255)
                                                                          cGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGAt                 <  1:845472/100‑1 (MQ=255)
                                                                           ggATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGAtt                >  2:1894577/1‑100 (MQ=255)
                                                                             aTCTCACTTTCCCGCAGCAGGTGGGCTAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGAtt                <  1:594410/98‑1 (MQ=255)
                                                                              tctcACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTaaa             <  2:1953027/100‑1 (MQ=255)
                                                                              tctcACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTaaa             >  1:2177854/1‑100 (MQ=255)
                                                                              tctcACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCAGAtcgga             >  1:321140/1‑96 (MQ=255)
                                                                               ctcACTTTCCCGCCGCAGGAGGGCGAGCAGGCGGACGTCAACCTTGCTTACTGCACGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTaaaa            >  1:243556/1‑100 (MQ=255)
                                                                                tcACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTACTGGTCGAGCTGATTaaaa            >  2:1238571/1‑99 (MQ=255)
                                                                                tcACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTaaaca           >  1:1802506/1‑98 (MQ=255)
                                                                                tcACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTaa              >  2:883895/1‑97 (MQ=255)
                                                                                 cACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAGCATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAAc          >  1:2101176/1‑100 (MQ=255)
                                                                                 cACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAAc          >  2:1167255/1‑100 (MQ=255)
                                                                                 cACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAAc          <  1:684800/100‑1 (MQ=255)
                                                                                 cACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAAc          >  1:697830/1‑100 (MQ=255)
                                                                                     ttCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAGCATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATc      <  2:2101162/100‑1 (MQ=255)
                                                                                     ttCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATc      >  1:848536/1‑100 (MQ=255)
                                                                                       cccGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATcgc    <  2:1512400/100‑1 (MQ=255)
                                                                                        ccgccgCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATcgcg   >  1:1264957/1‑100 (MQ=255)
                                                                                         cgccgcAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGt  >  1:1333907/1‑100 (MQ=255)
                                                                                            |                                                                                                
GCTGGTGATCCTTGAATCCACCTCGCCGGTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCAGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGT  >  NC_000913/3971598‑3971786

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