Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 606,299 A→T 15.8% V362E (GTG→GAG)  ybdK ← weak gamma‑glutamyl:cysteine ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913606,2990AT15.8% 419.5 / 60.3 172V362E (GTG→GAG) ybdKweak gamma‑glutamyl:cysteine ligase
Reads supporting (aligned to +/- strand):  ref base A (63/80);  new base T (9/18);  total (72/100)
Fisher's exact test for biased strand distribution p-value = 3.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-01

Only 100 of 174 total aligned reads displayed.

AACATTGTTAAATGTAAGTTAAAACACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCAATCGCG  >  NC_000913/606200‑606397
                                                                                                   |                                                                                                  
aaCATTGTTAAATGTAAGTTAAAACACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCa                                                                                                    <  2:4151743/100‑1 (MQ=255)
 aCATTGTTAAATGTAAGTTAAAACACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCAc                                                                                                   >  1:1074535/1‑100 (MQ=255)
 aCATTGTTAAATGTAAGTTAAAACACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCAc                                                                                                   >  1:516029/1‑100 (MQ=255)
  cATTGTTAAATGTAAGTTAAAACACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCAcc                                                                                                  >  2:2813593/1‑100 (MQ=255)
    ttGTTAAATGTAAGTTAAAACACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAg                                                                                                >  2:1784033/1‑100 (MQ=255)
     tGTTAAATGTAAGTTAACACACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGc                                                                                               <  1:1365314/100‑1 (MQ=255)
       ttAAATGTAAGTTAAAACACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGccc                                                                                             <  2:1389993/100‑1 (MQ=255)
       ttAAATGTAAGTTAAAACACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGccc                                                                                             <  1:728194/100‑1 (MQ=255)
       ttAAATGTAAGTTAAAACACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGccc                                                                                             <  1:4024768/100‑1 (MQ=255)
         aagtGTAAGTTAAAACACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCGTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCaa                                                                                           >  2:3538174/4‑100 (MQ=255)
          aaTGTAAGTTAAAACACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAAt                                                                                          >  2:2253508/1‑100 (MQ=255)
           aTGTAAGTTAAAACACCATTGTCATGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATc                                                                                         <  2:846884/100‑1 (MQ=255)
           aTGTAAGTTAAAACACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATc                                                                                         <  1:872971/100‑1 (MQ=255)
              tAAGTTAAAACACCATTGTCAGGGATATTCTTTTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGc                                                                                      <  2:232221/100‑1 (MQ=255)
              tAAGTTAAAACACCATTGTCAGGGATATTCTTCTGTACGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGc                                                                                      >  2:2226378/1‑100 (MQ=255)
               aaGTTAAAACACCATTGTCAGGGATATTCTGCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCg                                                                                     <  1:285351/100‑1 (MQ=255)
                     aaaCACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAgccgc                                                                               <  1:1096768/100‑1 (MQ=255)
                     aaaCACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAgccgc                                                                               >  1:2252510/1‑100 (MQ=255)
                       acacCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAgccgcc                                                                              <  1:1800520/99‑1 (MQ=255)
                       acacCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCCTTTTCACCAGCCCAATCAGCGAGCCGCCa                                                                             >  2:1064536/1‑100 (MQ=255)
                        cacCATTGTCAGGGCTACTCTTCTGTAAGGGTCTTCCCGTCTTAGCCACCGGCCCAGCTCTCACAATGCTTTTTGTCCAGCCCAATCAGCGAGCCGCCa                                                                             <  2:3891226/99‑1 (MQ=255)
                          ccATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCg                                                                          >  1:3877432/1‑100 (MQ=255)
                           cATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCgg                                                                         <  1:48760/100‑1 (MQ=255)
                           cATTGTCAGGGATATGCTTCTGTAAGGCAATTCCAGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCa                                                                             <  2:3286294/96‑1 (MQ=255)
                             ttGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCggag                                                                       <  1:2818774/100‑3 (MQ=255)
                               gTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTTATTTGTCTAATCAGCGAGCCGCCATTGGcgac                                                                     >  2:1246618/1‑100 (MQ=255)
                               gTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGcgac                                                                     <  1:2156797/100‑1 (MQ=255)
                                tCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTTACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGcgacg                                                                    <  2:2311451/100‑1 (MQ=255)
                                tCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGcgacg                                                                    <  1:2058884/100‑1 (MQ=255)
                                tCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGcgacg                                                                    <  1:914706/100‑1 (MQ=255)
                                tCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGcgacg                                                                    <  1:353313/100‑1 (MQ=255)
                                 cAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGcgacga                                                                   <  2:541058/100‑1 (MQ=255)
                                   gggATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGaaa                                                                 >  2:1275462/1‑100 (MQ=255)
                                   gggATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGaaa                                                                 <  1:998486/100‑1 (MQ=255)
                                     gATATTCTTCTGAAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATc                                                               >  1:2739076/1‑100 (MQ=255)
                                      atatTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATcg                                                              <  1:1908857/100‑1 (MQ=255)
                                      atatTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATcg                                                              <  1:4109169/100‑1 (MQ=255)
                                        atTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATcgcg                                                            <  2:1394046/100‑1 (MQ=255)
                                        atTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATcgcg                                                            <  1:3442901/100‑1 (MQ=255)
                                         ttcttcTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCgcgc                                                           >  1:4007959/1‑100 (MQ=255)
                                         ttcttcTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCgcgc                                                           <  1:3134888/100‑1 (MQ=255)
                                          tcttctGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCa                                                          >  1:2463355/1‑100 (MQ=255)
                                              ctGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAg                                                      >  2:4045953/1‑100 (MQ=255)
                                                gTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCt                                                    >  1:424350/1‑100 (MQ=255)
                                                 tAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTg                                                   >  1:1727986/1‑100 (MQ=255)
                                                  aaGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTgc                                                  >  1:2517922/1‑100 (MQ=255)
                                                  aaGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTgc                                                  <  1:2590758/100‑1 (MQ=255)
                                                   aGGAAATTCACGGCGTGGTCACCGGCCCAGATCCCCCAACGCCTTGTCATGAGCCCAATCAGCGAGGCGCCATAGGCGACAACATCgcgc                                                           <  2:3297212/90‑1 (MQ=255)
                                                    ggCAATTCCCGGCTTAGTCACCGGCCCAGGTCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCTGCCATCGGCGACGAAATCGCGCATTAGCTgcgc                                                <  2:2499094/100‑1 (MQ=255)
                                                    ggCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTgc                                                  >  2:2622132/1‑98 (MQ=255)
                                                     gCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCt                                               <  1:2382413/100‑1 (MQ=255)
                                                     gCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCt                                               <  2:3509542/100‑1 (MQ=255)
                                                       aaTTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCttn                                             <  1:1518248/100‑2 (MQ=255)
                                                       aaTTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCttc                                             <  2:666213/100‑1 (MQ=255)
                                                       aaTTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCttc                                             <  1:1683554/100‑1 (MQ=255)
                                                        aTTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCttca                                            >  1:2720377/1‑100 (MQ=255)
                                                        aTTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCttca                                            >  1:1988863/1‑100 (MQ=255)
                                                        aTTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCttca                                            >  1:90539/1‑100 (MQ=255)
                                                         ttCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCttcat                                           <  2:1710126/100‑1 (MQ=255)
                                                         ttCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCGTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCttcat                                           <  2:1655889/100‑1 (MQ=255)
                                                           cccGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCttcattc                                         <  1:3507168/100‑1 (MQ=255)
                                                            ccGGCTTCGTCACCGGCCAAGATCTGACAATGCTCTTTCACCCGCGCAGCCAGCGAGCCACCATAGACGACGAAATCGCGCATTAGCTGCCCttcattca                                        >  2:3718280/1‑100 (MQ=255)
                                                            ccGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCttcattca                                        >  1:614328/1‑100 (MQ=255)
                                                            ccGGCTTAGTCACCGGCCCAGATCTCACAATGCGTTTTCACCAGCCCAATTAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCttcattca                                        <  1:1256640/100‑1 (MQ=255)
                                                              ggCTTAGTCACCGGCCCAGATCTCACAATCCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAg                                       <  2:3833064/99‑1 (MQ=255)
                                                               gCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGAc                                     <  2:4173237/100‑1 (MQ=255)
                                                                cTTAGTCACCGGCCCAGATCTCACCATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAACTCGCGCATTAGCTGCGCTTCATTCAGAcc                                    <  1:1230307/100‑1 (MQ=255)
                                                                cTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGAcc                                    >  2:13274/1‑100 (MQ=255)
                                                                cTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGAc                                     <  2:4111111/99‑1 (MQ=255)
                                                                  tAGTCACCGGCCCAGATCTCACAATGCTTTTTTTTTTTTTTTATCATCGAGCCGCCATCGGCAAAGTAATTACGCATTAGCTGCGCTTTATTCAGAcc                                    >  2:1892552/1‑98 (MQ=255)
                                                                    gTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTg                                >  2:348380/1‑100 (MQ=255)
                                                                    gTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTg                                >  2:3584706/1‑100 (MQ=255)
                                                                    gTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTg                                >  1:1424084/1‑100 (MQ=255)
                                                                    gTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTg                                <  2:182737/100‑1 (MQ=255)
                                                                      cACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTgac                              <  1:3918069/100‑1 (MQ=255)
                                                                      cACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTgac                              >  1:3808486/1‑100 (MQ=255)
                                                                      cACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTgac                              >  2:2552039/1‑100 (MQ=255)
                                                                      cACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGAGCTTCATTCAGACCGCTgac                              >  2:1525253/1‑100 (MQ=255)
                                                                       acgggCCCAGATCTCACAATGCTTTTTCACCAGTCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTgacg                             <  2:3503461/97‑1 (MQ=255)
                                                                       aCCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTgacg                             <  1:32334/100‑1 (MQ=255)
                                                                       aCCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTgacg                             <  1:3873133/100‑1 (MQ=255)
                                                                         cGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTgacg                             <  2:784825/98‑1 (MQ=255)
                                                                          ggCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGAcc                          <  1:2600787/100‑1 (MQ=255)
                                                                             ccAGCTCTCACAATGCTTTTTCTCCAGCCCACTCCGCGAGCCGCCATCGGCGACGAATTCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTgg                       <  2:3682676/100‑1 (MQ=255)
                                                                              cAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGc                      <  2:2053997/100‑1 (MQ=255)
                                                                              cAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGc                      <  2:1565984/100‑1 (MQ=255)
                                                                                gATCTCTCAATGCTTTTTCACCAGCCCAATCAGCTAGCCGCCATCGGCGACGAAATCGCGCATTAGCTTCGCTGCATTCAGACCGCTGACGACCTgg                       <  2:1219752/97‑1 (MQ=255)
                                                                                       caATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAggg             <  2:1396710/100‑1 (MQ=255)
                                                                                       caATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAggg             >  2:253756/1‑100 (MQ=255)
                                                                                       caATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAggg             >  2:3366006/1‑100 (MQ=255)
                                                                                        aaTGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGc            <  2:126413/100‑1 (MQ=255)
                                                                                        aaTGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGc            >  1:2263470/1‑100 (MQ=255)
                                                                                         atgttttttCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGcc           <  1:1647053/96‑1 (MQ=255)
                                                                                           gCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTc         <  1:504627/100‑1 (MQ=255)
                                                                                              ttttCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCAAt      >  1:3609375/1‑100 (MQ=255)
                                                                                               tttCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCAATc     <  1:1754310/100‑1 (MQ=255)
                                                                                                ttCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCAATcg    <  1:2200368/100‑1 (MQ=255)
                                                                                                  cACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCAATcgcg  >  1:116989/1‑100 (MQ=255)
                                                                                                  cACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCAATcgcg  >  1:3225735/1‑100 (MQ=255)
                                                                                                  cACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCAATcgcg  >  1:2077941/1‑100 (MQ=255)
                                                                                                   |                                                                                                  
AACATTGTTAAATGTAAGTTAAAACACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCAATCGCG  >  NC_000913/606200‑606397

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