Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,596,395 C→A 17.1% E3D (GAG→GAT livH ← branched‑chain amino acid ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,596,3950CA17.1% 484.5 / 76.8 181E3D (GAG→GATlivHbranched‑chain amino acid ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (81/69);  new base A (11/20);  total (92/89)
Fisher's exact test for biased strand distribution p-value = 7.59e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

Only 100 of 185 total aligned reads displayed.

GATAATGCCGTAAACCATGGTGTAGCCGATGGCTATCAGCGCGTAGGTACTGCCCAGCGTGACGCCGTTAAACATCTGCTGCAAGAAATACAAAAACTGCTCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCACTTGGCTGCCGTGGATGAACCGTCGGCGTGCCACTGGAAGA  >  NC_000913/3596296‑3596494
                                                                                                   |                                                                                                   
gATAATGCCGTAAACCATGGTGTAGCCGATGGCTATCAGCGCGTAGGTACTGCCCAGCGTGACGCCGTTAAACATCTGCTGCAAGAAATACAAAAACTGc                                                                                                     >  1:1716270/1‑100 (MQ=255)
  tAATGCCGTAAACCATGGTGTAGCCGATTGCTATCAGCGCGTAGGTACTGCCCAGCGTGACGCCGTTAAACATCTGCTGCAAGAAATACAAAAACTGCt                                                                                                    <  1:3752443/99‑1 (MQ=255)
  tAATGCCGTAAACCATGGTGTAGCCGATGGCTATCAGCGCGTAGGTACTGCCCAGCGTGACGCCGTTAAACATCTGCTGCAAGAAATACAAAAACTGCTc                                                                                                   <  1:377561/100‑1 (MQ=255)
  tAATGCCGTAAACCATGGTGTAGCCGATGGCTATCAGCGCGTAGGTACTGCCCAGCGTGACGCCGTTAAACATCTGCTGCAAGAAATACAAAAACTGCTc                                                                                                   <  2:1088535/100‑1 (MQ=255)
   aaTGCCGTAAACCATGGTGTCGCCGATGGCTATCAGCGCGTAGGTACTGCCCAGCGTGACGCCGTTAAACATCTGCTGCAAGAAATACAAAAACTGCTCa                                                                                                  >  2:2330784/1‑100 (MQ=255)
   aaTGCCGTAAACCATGGTGTAGCCGATGGCTATCAGCGCGTAGGTACTGCCCAGCGTGACGCCGTTAAACATCTGCTGCAAGAAATACAAAAACTGCTCa                                                                                                  >  2:1347471/1‑100 (MQ=255)
   aaTGCCGTAAACCATGGTGTAGCCGATGGCTATCAGCGCGTAGGTACTGCCCAGCGTGACGCCGTTAAACATCTGCTGCAAGAAATACAAAAACTGCTCa                                                                                                  <  2:1524425/100‑1 (MQ=255)
    aTGCCGTAAACCAGGGTGGAGCCGACAGCTATCAGCGCGTAGGTACTGCCCAGATTGACGCCGTTAAATATCTGCTGAAAAAAAAAAAAAAACTGATCAg                                                                                                 <  2:3735958/100‑1 (MQ=255)
     tGCCGTAAACCATGGTGTAGCCGATGGCTATCAGCGCGTAGGTACTGCCCAGCGTGACGCCGTTAAACATCTGCTGCAAGAAATACAAAAACTGATCAGa                                                                                                <  2:1116037/100‑1 (MQ=255)
      gCCGTAAACCATGGTGTAGCCGATGGCTATCAGCGCGTAGGTACTGCCCAGCGTGACGCCGTTAAACATCTGCTGCAAGAAATACAAAAACTGATCAGAc                                                                                               <  2:1288893/100‑1 (MQ=255)
      gCCGTAAACCATGGTGTAGCCGATGGCTATCAGCGCGTAGGTACTGCCCAGCGTGACGCCGTTAAACATCTGCTAAAAAAAAAAAAAAAACTGATCAGa                                                                                                <  2:615630/99‑1 (MQ=255)
          tAAACCATGGTGTAGCCGATGGCTATCAGCGCGTAGGTACTGCCCAGCGTGACGCCGTTAAACATCTGCTGCAAGAAATACAAAAACTGCTCAGACATaa                                                                                           <  2:1788957/100‑1 (MQ=255)
            aaCCATGGTGTAGCCGATGGCTATCAGCGCGTAGGTACTGCCCAGCGTGACGCCGTTAAACATCTGCTGCAAGAAATACAAAAACTGCTCAGACATAAgg                                                                                         >  1:2792553/1‑100 (MQ=255)
            aaCCATGGTGTAGCCGATGGCTATCAGCGCGTAGGTACTGCCCAGCGTGACGCCGTTAAACATCTGCTGCAAGAAATACAAAAACTGATCAGACATAAgg                                                                                         <  2:306999/100‑1 (MQ=255)
             aCCATGGTGTAGCCGATGGCTATCAGCGCGTAGGTACTGCCCAGCGTGACGCCGTTAAACATCTGCTGCAAGAAATACAAAAACTGCTCAGACATAAGGt                                                                                        <  2:3608631/100‑1 (MQ=255)
              ccATGGTGTAGCCGATGGCTATCAGCGCGTAGGTACTGCCCAGCGTGACGCCGTTAAACATCTGCTGCAAGAAATACAAAAACTGCTCAGACATAAGGTa                                                                                       >  2:252214/1‑100 (MQ=255)
               cATGGTGTAGCCGATGGCTATCAGCGCGTAGGTACTGCCCAGCGTGACGCCGTTAAACATCTGCTGCAAGAAATACAAAAACTGCTCAGACATAAGGTaa                                                                                      >  1:419849/1‑100 (MQ=255)
                aTGGTGTAGCCGATGGCTATCAGCGCGTAGGTACTGCCCAGCGTGACGCCGTTAAACATCTGCTAAAAAAAAAAAAAAAACTGCTCAGACATAAgg                                                                                         <  2:2269393/96‑1 (MQ=255)
                 tGGTGTAGCCGATGGCTATCAGCGCGTAGGTACTGCCCAGCGTGACGCCGTTAAACATCTGCTGCAAGAAATACAAAAACTGCTCAGACATAAGGTAAcc                                                                                    <  2:2744549/100‑1 (MQ=255)
                  ggTGTAGCCGATGGCTATCAGCGCGGAGGTACTGCCCAGCGTGACGCCGTTAAACATCTGCTGCAAGAAATACAAAAACTGATCAGACATAAGGTAAcc                                                                                    <  2:2830021/99‑1 (MQ=255)
                     gtAGCCGATGGCTATCAGCGCGTAGGTACTGCCCAGCGTGACGCCGTTAAACATCTGCTGCAAGAAATACAAAAACTGCTCAGACATAAGGTAACCTTTc                                                                                >  2:3303585/1‑100 (MQ=255)
                     gtAGCCGATGGCTATCAGCGCGTAGGTACTGCCCAGCGTGACGCCGTTAAACATCTGCTGCAAGAAATACAAAAACTGCTCAGACATAAGGTAACCTTTc                                                                                <  1:3744172/100‑1 (MQ=255)
                       aGCCGATGGCTATCAGCGCGTAGGTACTGCCCAGCGTGACGCCGTTAAACATCTGCTGCAAGAAATACAAAAACTGCTCAGACATAAGGTAACCTTTCTa                                                                              >  2:3558556/1‑100 (MQ=255)
                       aGCCGATGGCTATCAGCGCGTAGGTACTGCCCAGCGTGACGCCGTTAAACATCTGCTGCAAGAAATACAAAAACTGCTCAGACATAAGGTAACCTTTCTa                                                                              <  1:3138240/100‑1 (MQ=255)
                       aGCCGATGGCTATCAGCGCGTAGGTACTGCCCAGCGTGACGCCGTTAAACATCTGCTGCAAGAAATACAAAAACTGCTCAGACATAAGGTAACCTTTCTa                                                                              >  2:96628/1‑100 (MQ=255)
                        gCCGATGGCTATCAGCGCGTAGGTACTGCCCAGCGTGACGCCGTTAAACATCTGCTGCAAGAAATACAAAAACTGATCAGACATAAGGTAACCTTTCTaa                                                                             >  1:3399128/1‑100 (MQ=255)
                         ccGATGGCTATCAGCGCGTAGGTACTGCCCAGCGTGACGCCGTTAAACATCTGCTGCAAGAAATACAAAAACTGCTCAGACATAAGGTAACCTTTCTaaa                                                                            >  1:2472130/1‑100 (MQ=255)
                          cGATGGCTATCAGCGCGTAGGTACTGCCCAGCGTGACGCCGTTAAACATCTGCTGCAAGAAATACAAAAACTGATCAGACATAAGGTAACCTTTCTAAAc                                                                           <  2:261276/100‑1 (MQ=255)
                           gATGGCTATCAGCGCGTAGGTACTGCCCAGCGTGACGCCGTTAAACATCTGCTGCAAGAAATACAAAAACTGCTCAGACATAAGGTAACCTTTCTAAAcc                                                                          >  1:1752053/1‑100 (MQ=255)
                               gCTATCAGCGCGTAGGTACTGCCCAGCGTGACGCCGTTAAACATCTGCTGCAAGAAATACAAACACTGCTCAGACATAAGGTAACCTTTCTAAAcccgcc                                                                      <  2:692440/100‑1 (MQ=255)
                               gCTATCAGCGCGTAGGTACTGCCCAGCGTGACGCCGTTAAACATCTGCTGCAAGAAATACAAAAACTGCTCAGACATAAGGTAACCTTTCTAAAcccgcc                                                                      >  1:477283/1‑100 (MQ=255)
                               gCTATCAGCGCGTAGGTACTGCCCAGCGTGACGCCGTTAAACATCTGCTGCAAGAAATACAAAAACTGATCAGACATAAGGTAACCTTTCTAAAcccgcc                                                                      >  2:1827134/1‑100 (MQ=255)
                                 tATCAGCGCGTAGGTACTGCCCAGCGTGACGCCGTTAAACATCTGCTGCAAGAAATACAAAAACTGCTCAGACATAAGGTAACCTTTCTAAAcccgcccg                                                                    <  2:2489798/100‑1 (MQ=255)
                                     aGCGCGTAGGTACTGCCCAGCGTGACGCCGTTAAACATCTGCTGCAAGAAATACAAAAACTGCTCAGACATAAGGTAACCTTTCTAAACCCGCCCGCAtt                                                                >  1:3079640/1‑100 (MQ=255)
                                      gcgcGTAGGTACTGCCCAGCGTGACGCCGTTAAACATCTGCTGCAAGAAATACAAAAACTGATCAGACATAAGGTAACCTTTCTAAACCCGCCCGCAttt                                                               >  1:1575430/1‑100 (MQ=255)
                                            aGGTACTGCCCAGCGTGACGCCGTTAAACATCTGCTGCAAGAAATACAAAAACTGCTCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACggg                                                         <  2:2733930/100‑1 (MQ=255)
                                             ggTACTGCCCAGCGTGACGCCGTTAAACATCTGCTGCAAGAAATACAAAAACTGCTCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGc                                                        >  1:262396/1‑100 (MQ=255)
                                             ggTACTGCCCAGCGTGACGCCGTTAAACATCTGCTGCAAGAAATACAAAAACTGCTCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGc                                                        >  1:1519871/1‑100 (MQ=255)
                                              gTACTGCCCAGCGTGACGCCGTTAAACATCTGCTGCAAGAAATACAAAAACTGATCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCg                                                       <  2:3394327/100‑1 (MQ=255)
                                               tACTGCCCAGCGTGACGCCGTTAAACATCTGCTGCAAGAAATACAAAAACTGCTCAGGCACAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCgg                                                      >  2:2676310/1‑100 (MQ=255)
                                                 cTGCCCAGCGTGACGCCGTTAAACATCTGCTGCAAGAAATACAAAAACTGCTCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTg                                                    >  1:816615/1‑100 (MQ=255)
                                                 cTGCCCAGCGTGACGCCGTTAAACATCTGCTGCAAGAAATACAAAAACTGCTCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTg                                                    <  1:285590/100‑1 (MQ=255)
                                                  tGCCCAGCGTGACGCCGTTAAACATCTGCTGCAAGAAATACAAAAACTGCTCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTgg                                                   <  2:197543/100‑1 (MQ=255)
                                                     ccAGCGTGACGCCGTTAAACATCTGCTGCAAGAAATACAAAAACTGATCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGgat                                                >  1:2094481/1‑100 (MQ=255)
                                                      cAGCGTGACGCCGTTAAACATCTGCTGCAAGAAATACAAAAACTGCTCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGgatg                                               >  1:1666009/1‑100 (MQ=255)
                                                      cAGCGTGACGCCGTTAAACATCTGCTGCAAGAAATACAAAAACTGCTCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGgat                                                <  2:3218419/99‑1 (MQ=255)
                                                       aGCGTGACGCCGTTAAACATCTGCTGCAAGAAATACAAAAACTGCTCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGgatga                                              <  2:2498347/100‑1 (MQ=255)
                                                         cGTGACGCCGTTAAACATCTGCTGCAAGAAATACAAAAACTGCTCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATc                                            <  1:237951/100‑1 (MQ=255)
                                                         cGTGACGCCGTTAAACATCTGCTGCAAGAAATACAAAAACTGCTCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATc                                            <  1:3430156/100‑1 (MQ=255)
                                                           taaCGCCGTTAAAGATCTGCTGCAAGAAATACAAAAACTGCTCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCAc                                          >  1:2589913/3‑100 (MQ=255)
                                                           tGACGCCGTTAAACATCTGCTGCAAGAAATACAAAAACTGATCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCAc                                          <  2:1583329/100‑1 (MQ=255)
                                                              cGCCGTTAAACATCTGCTGCAAGAAATACAAAAACTGCTCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCACTTg                                       <  1:662958/100‑1 (MQ=255)
                                                              cGCCGTTAAACATCTGCTGCAAGAAATACAAAAACTGCTCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCACTTg                                       <  2:3470488/100‑1 (MQ=255)
                                                               gCCGTTAAACATCTGCTGCAAGAAATACAAAAACTGCTCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCACTTgg                                      >  2:467545/1‑100 (MQ=255)
                                                               ccgatctaaCATCTGCTGCAAGAAATACAAAAACTGCTCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCACTTgg                                      <  1:988267/93‑1 (MQ=255)
                                                                ccGTTAAACATCTGCTGCAAGAAATACAAAAACTGCTCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCACTTGGc                                     <  1:1049419/100‑1 (MQ=255)
                                                                ccGTTAAACATCTGCTGCAAGAAATACAAAAACTGATCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCACTTGGc                                     <  1:594710/100‑1 (MQ=255)
                                                                 cGTTAAACATCTGCTGCAAGAAATACAAAAACTGCTCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCACTTGGCt                                    <  2:3012569/100‑1 (MQ=255)
                                                                 cGTTAAACATCTGCTGCAAGAAATACAAAAACTGCTCAGACATAAGGTAACCTTTCTAAACCAGCCCGCATTTTACGGGCGGTGGGATGATCACTTGTCt                                    >  2:303606/1‑100 (MQ=255)
                                                                 cGTTAAACATCTGCTGCAAGAAATACAAAAACTGCTCAGACATAAGATAACCTTTCTAAAACCGCCCGCATTTTACGGGCGGTGGCATGATCACATg                                       >  2:1660488/1‑97 (MQ=255)
                                                                  gTTAAACATCTGCTGCAAGAAATACAAAAACTGCTCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCACTTGGCTg                                   <  1:2127246/100‑1 (MQ=255)
                                                                  gTTAAACATCTGCTGCAAGAAATACAAAAACTGCTCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCACTTGGCTg                                   >  1:907374/1‑100 (MQ=255)
                                                                     aaaCATCTGCTGCAAGAAATACAAAAACTGCTCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCACTTGGCTGcc                                 >  2:2533491/1‑99 (MQ=255)
                                                                     aaaCATCTGCTGCAAGAAATACAAAAACTGCTCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCACTTGGCTGCCg                                >  2:2974240/1‑100 (MQ=255)
                                                                      aaCATCTGCTGCAAGAAATACAAAAACTGATCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCACTTGGCTGCCGt                               <  2:340849/100‑1 (MQ=255)
                                                                       aCATCTGCTGCAAGAAATACAAAAACTGCTCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCACTTGGCTGCCGTg                              <  1:1030450/100‑1 (MQ=255)
                                                                       aCATCTGCTGCAAGAAATACAAAAACTGCTCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCACTTGGCTGCCGTg                              <  1:2743340/100‑1 (MQ=255)
                                                                       aCATCTGCTGCAAGAAATACAAAAACTGATCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCACTTGGCTGCCGTg                              >  1:1778624/1‑100 (MQ=255)
                                                                        cATCTGCTGCAAGAAATACAAAAACTGCTCAGACATAAGGTAACCTTTCTAAACACGCCCGCATTTATCGGGCGGTGGGATGATCACTTGGCTGCCGTgg                             >  2:1454903/1‑100 (MQ=255)
                                                                        cATCTGCTGCAAGAAATACAAAAACTGATCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCACTTGGCTGCCGTgg                             >  2:770192/1‑100 (MQ=255)
                                                                          tCTGCTGCAAGAAATACAAAAACTGCTCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCACTTGGCTGCCGTGGAt                           <  1:125900/100‑1 (MQ=255)
                                                                          tCTGCTGCAAGAAATACAAAAACTGCTCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCACTTGGCTGCCGTGGAt                           <  1:1638374/100‑1 (MQ=255)
                                                                          tCTGCTGCAAGAAATACAAAAACTGCTCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCACTTGGCTGCCGTGGAt                           >  1:142354/1‑100 (MQ=255)
                                                                             gctgcAAGAAATACAAAAACTGCTCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCACTTGGCTGCCGTGGATGaa                        >  2:2065188/1‑100 (MQ=255)
                                                                              ctgcAAGAAATACAAAAACTGCTCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCACTTGGCTGCCGTGGATGAAc                       >  1:1304227/1‑100 (MQ=255)
                                                                              ctgcAAGAAATACAAAAACTGCTCAGACATAAGGTAACCTTTCTAAAACCGCCCGCATTTTACGGGCGGTGGGATGATCACTTGGCTGCCGTGGATGaa                        >  1:3236514/1‑99 (MQ=255)
                                                                              ctgcAAGAAAGACAAAAACTGATCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCACTTGGCTGCCGTGGATGAAc                       <  2:2384548/100‑1 (MQ=255)
                                                                               tgcAAGAAATACAAAAACTGCTCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCACTTGGCTGCCGTGGATGAAcc                      <  1:3361768/100‑1 (MQ=255)
                                                                               tgcAAGAAATACAAAAACTGCTCAGACATAAGGTAACCTTTCGAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCACTTGGCTGCCGTGGATGAAcc                      <  1:3678631/100‑1 (MQ=255)
                                                                                gcAAGAAATACAAAAACTGCTCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCACTTGGCTGCCGTGGATGAACCg                     >  1:2315675/1‑100 (MQ=255)
                                                                                gcAAGAAATACAAAAACTGCTCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCACTTGGCTGCCGTGGATGAACCg                     >  1:2026527/1‑100 (MQ=255)
                                                                                gcAAGAAATACAAAAACTGATCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCACTTGGCTGCCGTGGATGAACCg                     >  1:2087269/1‑100 (MQ=255)
                                                                                  aaGAAATACAAAAACTGCTCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCACTTGGCTGCCGTGGATGAACCGTc                   >  2:533167/1‑100 (MQ=255)
                                                                                    gAAATACAAAAACTGATCAGACATAAGGTAACCTTGCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCACTTGGCTGCCGTGGATGAACCGTCgg                 <  2:2889735/100‑1 (MQ=255)
                                                                                       aTACAAAAACTGCTCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCACTTGGCTGCCGTGGATGAACCGTCGGCGt              <  1:3300807/100‑1 (MQ=255)
                                                                                        tACAAAAACTGCTCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCACTTGGCTGCCGTGGATGAACCGTCGGCGTg             <  2:1881870/100‑1 (MQ=255)
                                                                                        tACAAAAACTGCTCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCACTTGGCTGCCGTGGATGAACCGTCGGCGTg             >  1:1504117/1‑100 (MQ=255)
                                                                                        tACAAAAACTGCTCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCACTTGGCTGCCGTGGATGAACCGTCGGCGTg             >  1:2465819/1‑100 (MQ=255)
                                                                                         aCAATAACTGCTCAGACATAAGGTAACCTTTCTAAACCCACCCGCCTTCTACGGGCGGTGGGACGATCACTTCGCTGCCGTCGATGAAACGTCGGCGTg             >  2:2349217/1‑99 (MQ=255)
                                                                                         aCAAAAACTGCTCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCACTTGGCTGCCGTGGATGAACCGTCGGCGTGc            <  2:1376455/100‑1 (MQ=255)
                                                                                         aCAAAAACTGCTCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCACTTGGCTGCCGTGGATGAACCGTCGGCGTGc            >  2:434669/1‑100 (MQ=255)
                                                                                          cAAAAACTGCTCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCACTTGGCTGCCGTGGATGAACCGTCGGCGTGcc           >  2:2253016/1‑100 (MQ=255)
                                                                                          cAAAAACTGCTCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCACTTGGCTGCCGTGGATGAACCGTCGGCGTGcc           >  2:866293/1‑100 (MQ=255)
                                                                                          cAAAAACTGCTCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCACTTGGCTGCCGTGGATGAACCGTCGGCGTGcc           >  2:2429781/1‑100 (MQ=255)
                                                                                          cAAAAACTGCTCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCACTTGGCTGCCGTGGATGAACCGTCGGCGTGcc           >  1:1966219/1‑100 (MQ=255)
                                                                                           aaaaaCTGCTCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCACTTGGCTGCCGTGGATGAACCGTCGGCGTGCCa          >  1:3416037/1‑100 (MQ=255)
                                                                                             aaaCTGCTCAGACATAAGGTGACCTTTCTTAACCCGCCCGCATTTTACGGGCGGGGGGATGATCACTTGGCTGCCGTGGATGTACCGTCGGGGTGCCACt        <  2:3332929/100‑1 (MQ=255)
                                                                                             aaaCTGCTCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCACTTGGCTGCCGTGGATGAACCGTCGGCGTGCCACt        >  1:65519/1‑100 (MQ=255)
                                                                                              aaCTGCTCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCACTTGGCTGCCGTGGATGAACCGTCGGCGTGCCACTg       <  1:2735415/100‑1 (MQ=255)
                                                                                                cTGCTCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCACTTGGCTGCCGTGGATGAACCGTCGGCGTGCCACTGGa     <  2:2922266/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
GATAATGCCGTAAACCATGGTGTAGCCGATGGCTATCAGCGCGTAGGTACTGCCCAGCGTGACGCCGTTAAACATCTGCTGCAAGAAATACAAAAACTGCTCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCACTTGGCTGCCGTGGATGAACCGTCGGCGTGCCACTGGAAGA  >  NC_000913/3596296‑3596494

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