Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,952,569 C→T 39.8% A29V (GCG→GTG)  ilvE → branched‑chain amino acid aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,952,5690CT39.8% 129.1 / 182.2 166A29V (GCG→GTG) ilvEbranched‑chain amino acid aminotransferase
Reads supporting (aligned to +/- strand):  ref base C (49/50);  new base T (31/35);  total (80/86)
Fisher's exact test for biased strand distribution p-value = 8.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 171 total aligned reads displayed.

AAGGAATATAAAAATGACCACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTCTGCA  >  NC_000913/3952471‑3952668
                                                                                                  |                                                                                                   
naGGAATATAAAAATGACCACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCAcgcg                                                                                                    >  1:2469446/2‑100 (MQ=255)
 aGGAATATAAAAATGACCACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACgtgc                                                                                                   <  2:1228081/100‑1 (MQ=255)
  ggAATATAAAAATGACCACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCt                                                                                                  <  1:2913234/100‑1 (MQ=255)
   gAATATAAAAATGACCACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTg                                                                                                 >  2:3319272/1‑100 (MQ=255)
    aaTATAAAAATGACCACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGc                                                                                                <  2:2657346/100‑1 (MQ=255)
    aaTATAAAAATGACCACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGc                                                                                                <  1:2545109/100‑1 (MQ=255)
     atatAAAAATGACCACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCa                                                                                               >  1:3265361/1‑100 (MQ=255)
      tataAAAATGACCACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCAc                                                                                              >  1:1530656/1‑100 (MQ=255)
      tataAAAATGACCAAGAAGAAAGCTGATTACATTAAGTACAATGGAGAGATAGCACGCTAGGAAGACGCAAAGCTGCACCTCATGTCGCACGACCTGccc                                                                                              <  1:2239693/100‑3 (MQ=255)
       ataAAAATGACCACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACt                                                                                             >  1:170079/1‑100 (MQ=255)
       ataAAAATGACCACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGATGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACt                                                                                             >  1:182853/1‑100 (MQ=255)
        taAAAATGACCACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTa                                                                                            <  1:547014/100‑1 (MQ=255)
        taAAAATGACCACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTa                                                                                            >  1:1784288/1‑100 (MQ=255)
        taAAAATGACCACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTa                                                                                            >  2:557309/1‑100 (MQ=255)
        taAAAATGACCACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTa                                                                                            >  1:743135/1‑100 (MQ=255)
        taAAAATGACCACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTa                                                                                            <  1:2193049/100‑1 (MQ=255)
          aaaaTGACCACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATg                                                                                          <  1:1509/100‑1 (MQ=255)
          aaaaTGACCACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATg                                                                                          >  2:1490760/1‑100 (MQ=255)
            aaTGACCACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGc                                                                                        <  1:1587085/100‑1 (MQ=255)
            aaTGACCACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGc                                                                                        <  2:1850649/100‑1 (MQ=255)
             aTGACCACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCa                                                                                       <  1:1618819/100‑1 (MQ=255)
             aTGACCACGAAGAAAGCGGAGTACATTTGGTGCAATGGGGAGATTGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATgg                                                                                         <  1:2729109/98‑1 (MQ=255)
              tGACCACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCAc                                                                                      <  2:3111764/100‑1 (MQ=255)
              tGACCACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCAc                                                                                      >  2:174757/1‑100 (MQ=255)
                 ccACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGCCGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTc                                                                                   <  1:835224/100‑1 (MQ=255)
                 ccACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTc                                                                                   >  1:2342911/1‑100 (MQ=255)
                 ccACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGAAGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTc                                                                                   >  2:397471/1‑100 (MQ=255)
                  cACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCg                                                                                  <  1:130745/100‑1 (MQ=255)
                   aCGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCgg                                                                                 >  2:2025318/1‑100 (MQ=255)
                   aCGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGAACGCGCTGCACTATGGCACTTCgg                                                                                 >  1:3444198/1‑100 (MQ=255)
                   aCGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTGCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCgg                                                                                 <  2:2952181/100‑1 (MQ=255)
                    cGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGt                                                                                >  1:197961/1‑100 (MQ=255)
                    cGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGt                                                                                >  2:609689/1‑100 (MQ=255)
                     gaagaaAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGtt                                                                               <  2:2771946/100‑1 (MQ=255)
                     gaagaaAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGtt                                                                               >  1:3429880/1‑100 (MQ=255)
                        gaaAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGttttt                                                                            <  1:3361281/100‑1 (MQ=255)
                        gaaAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTCTGGCACTTCGGttttt                                                                            <  2:2373956/100‑1 (MQ=255)
                          aaGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGa                                                                          >  2:3439748/1‑100 (MQ=255)
                          aaGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGa                                                                          >  1:1047872/1‑100 (MQ=255)
                            gCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAg                                                                        >  1:2050236/1‑100 (MQ=255)
                            gCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAg                                                                        >  1:1267628/1‑100 (MQ=255)
                             cTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAgg                                                                       <  2:3285419/100‑1 (MQ=255)
                             cTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAgg                                                                       <  2:331043/100‑1 (MQ=255)
                             cTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAgg                                                                       <  2:3226658/100‑1 (MQ=255)
                             cTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAgg                                                                       <  1:1490746/100‑1 (MQ=255)
                                aTTATATTTGGTTCAATGGTGAGATGGTTCGCTGGGAAGCCGCGAAGGTGCATGTGATGTCGCACGTGCTGCATTTTGGCACTTCGGTTTTTGAAGGCAt                                                                    <  1:952248/100‑1 (MQ=255)
                                aTTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCAt                                                                    >  1:944043/1‑100 (MQ=255)
                                     aTTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGtt                                                               >  1:2058506/1‑100 (MQ=255)
                                      tttGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTg                                                              <  1:1577920/100‑1 (MQ=255)
                                      tttGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTg                                                              <  2:1146843/100‑1 (MQ=255)
                                       ttGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGc                                                             >  2:2734887/1‑100 (MQ=255)
                                         ggTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTa                                                           >  1:2842761/1‑100 (MQ=255)
                                            tCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGa                                                        >  2:1072466/1‑100 (MQ=255)
                                               aTGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTc                                                     >  1:1100046/1‑100 (MQ=255)
                                               aTGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTc                                                     >  2:2120799/1‑100 (MQ=255)
                                                 ggggAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGc                                                   <  2:2527800/100‑1 (MQ=255)
                                                 ggggAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCg                                                    <  1:854593/99‑1 (MQ=255)
                                                  gggAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGca                                                  <  2:230041/100‑1 (MQ=255)
                                                   ggAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGcac                                                 >  1:2590642/1‑100 (MQ=255)
                                                   ggAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTGCGGGTTTTGAAGGCATCCGTTGCTACGACTCGcac                                                 >  1:2933008/1‑100 (MQ=255)
                                                      gaTGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACaaa                                              >  2:2566262/1‑100 (MQ=255)
                                                       aTGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAg                                             <  2:795221/100‑1 (MQ=255)
                                                       aTGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAg                                             <  2:1590308/100‑1 (MQ=255)
                                                       aTGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAg                                             <  2:3147610/100‑1 (MQ=255)
                                                        tGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAgg                                            <  2:2171156/100‑1 (MQ=255)
                                                         ggTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGa                                           >  2:431402/1‑100 (MQ=255)
                                                          gTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGAc                                          >  2:2802029/1‑100 (MQ=255)
                                                           ttCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGAcc                                         <  1:2137629/100‑1 (MQ=255)
                                                            tCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCg                                        <  1:3127656/100‑1 (MQ=255)
                                                            tCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCg                                        <  2:2050253/100‑1 (MQ=255)
                                                            tCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCg                                        <  1:91579/100‑1 (MQ=255)
                                                             cGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCgg                                       <  2:2809193/100‑1 (MQ=255)
                                                              gCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGt                                      <  1:1636779/100‑1 (MQ=255)
                                                                tGGGAAGACGCGAAGGTGCATGTGATGTCGCACTCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTg                                    <  1:2387115/100‑1 (MQ=255)
                                                                tGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTg                                    <  2:2058523/100‑1 (MQ=255)
                                                                tGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTg                                    <  1:2672030/100‑1 (MQ=255)
                                                                tGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTg                                    <  2:1788346/100‑1 (MQ=255)
                                                                 gggAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGt                                   >  1:2998830/1‑100 (MQ=255)
                                                                     aGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTc                               >  1:3233499/1‑100 (MQ=255)
                                                                     aGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTAtt                                >  2:3084149/1‑99 (MQ=255)
                                                                      gACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTcc                              <  1:721896/100‑1 (MQ=255)
                                                                      gACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTcc                              <  2:716498/100‑1 (MQ=255)
                                                                        cgcgAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGc                            >  2:2896618/1‑100 (MQ=255)
                                                                        cgcgAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGc                            >  2:2451545/1‑100 (MQ=255)
                                                                        cgcgAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGc                            >  1:2784043/1‑100 (MQ=255)
                                                                          cgAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCa                          >  2:1958996/1‑100 (MQ=255)
                                                                          cgAAGGTGCAGGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTTCGTCTCGCACAAAGGACCGGTTGTATTCCGCCa                          <  1:2279086/100‑1 (MQ=255)
                                                                           gAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCAt                         <  2:1917233/100‑1 (MQ=255)
                                                                             aGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCg                       >  1:1929138/1‑100 (MQ=255)
                                                                               gTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTg                     <  2:384511/100‑1 (MQ=255)
                                                                                 gCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAg                   <  1:309222/100‑1 (MQ=255)
                                                                                 gCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAg                   <  1:2826178/100‑1 (MQ=255)
                                                                                  cATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGc                  >  1:1229316/1‑100 (MQ=255)
                                                                                        aTGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGc            <  2:1267640/100‑1 (MQ=255)
                                                                                          gTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAg          >  2:547988/1‑100 (MQ=255)
                                                                                            cGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCg        >  2:1507252/1‑100 (MQ=255)
                                                                                            cGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGGGAGCGTATGCAGCg        >  2:1370550/1‑100 (MQ=255)
                                                                                            cGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCg        <  1:238946/100‑1 (MQ=255)
                                                                                            cGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCg        <  1:591414/100‑1 (MQ=255)
                                                                                                  cgcTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTCTGCa  >  2:180044/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
AAGGAATATAAAAATGACCACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTCTGCA  >  NC_000913/3952471‑3952668

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: