Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,701,108 A→T 15.9% N51I (AAT→ATT)  malY → PLP‑dependent beta‑cystathionase and maltose regulon regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,701,1080AT15.9% 266.7 / 40.6 106N51I (AAT→ATT) malYPLP‑dependent beta‑cystathionase and maltose regulon regulator
Reads supporting (aligned to +/- strand):  ref base A (50/39);  new base T (11/6);  total (61/45)
Fisher's exact test for biased strand distribution p-value = 5.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.45e-01

Only 100 of 109 total aligned reads displayed.

GGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCG  >  NC_000913/1701009‑1701206
                                                                                                   |                                                                                                  
gggATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGaa                                                                                                    >  1:382191/1‑100 (MQ=255)
  gATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATc                                                                                                  >  1:1106706/1‑100 (MQ=255)
   aTTATGTCGCTGACCGTTTCGGCACTGCTGACCTTTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGATGCGCTGAATCa                                                                                                 <  1:1529876/100‑1 (MQ=255)
    ttATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAg                                                                                                <  1:529263/100‑1 (MQ=255)
      aTGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTACCGAGGCGCTGAATCAgcg                                                                                              >  2:136418/1‑100 (MQ=255)
       tGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAgcgc                                                                                             >  2:2909342/1‑100 (MQ=255)
        gTCGCTGACCGGTTCGGCACTGCTGACATGTTACCGTTCACGATTTCAGACGTGGAGTTTGCCACTGCCCCCTGCATTAGCTAGGCGCTGAATCAgcgc                                                                                             >  1:2431511/1‑99 (MQ=255)
         tCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCt                                                                                           >  2:1206991/1‑100 (MQ=255)
          cGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGATTCAGCGCCtc                                                                                          <  2:724576/100‑2 (MQ=255)
          cGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGATTCAGCGCCtc                                                                                          <  1:2492553/100‑2 (MQ=255)
          cGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGATTCAGCGCCtc                                                                                          <  2:880369/100‑2 (MQ=255)
           gCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGATTCAGCGCCtcg                                                                                         <  1:2598249/100‑3 (MQ=255)
            cTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGAt                                                                                        <  1:605276/100‑1 (MQ=255)
             tGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGATTCAGCGCCtcgat                                                                                       <  2:464078/100‑5 (MQ=255)
             tGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATg                                                                                       <  2:1545575/100‑1 (MQ=255)
              gACGGTTTCGGCACTGCTGACCTGTTACTGTTTACGATTTCAGACATGGATTTTGCCACTGCCCCGTGCATTATCGAGGCGCTGAATCAGCGCCTGATGc                                                                                      >  1:1280608/1‑100 (MQ=255)
              gACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGc                                                                                      >  1:836134/1‑100 (MQ=255)
              gACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTAAGACAAGGATTTTGCCACTGGCACCTGCATTATCGAGGCGCTGAATCAGCGCCTGagac                                                                                      >  1:1232101/1‑97 (MQ=255)
                ccGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCAc                                                                                    >  1:1062176/1‑100 (MQ=255)
                ccGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCAc                                                                                    >  2:1188979/1‑100 (MQ=255)
                  gTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGTCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACgg                                                                                  >  1:1335611/1‑100 (MQ=255)
                  gTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGATTCAGCGCCtcgataatgc                                                                                  <  2:1424316/100‑10 (MQ=38)
                  gTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACgg                                                                                  <  1:1954197/100‑1 (MQ=255)
                   tttCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGc                                                                                 <  2:1371321/100‑1 (MQ=255)
                   tttCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGc                                                                                 <  1:878628/100‑1 (MQ=255)
                   tttCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGc                                                                                 <  2:2300550/100‑1 (MQ=255)
                      cGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTa                                                                              >  1:1816294/1‑100 (MQ=255)
                         cACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTAttt                                                                           <  2:127152/100‑1 (MQ=255)
                            tgctgACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCTAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGc                                                                        <  1:340282/100‑1 (MQ=255)
                                gACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACa                                                                    >  1:1808533/1‑100 (MQ=255)
                                 aCCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCCGCGCCTGATGCACGGCGTATTTGGCTACAg                                                                   >  1:2736785/1‑100 (MQ=255)
                                 aCCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAg                                                                   >  2:2817907/1‑100 (MQ=255)
                                 aCCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAg                                                                   <  2:707497/100‑1 (MQ=255)
                                   cTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGcc                                                                 <  1:2555550/100‑1 (MQ=255)
                                    tGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCg                                                                >  1:247149/1‑100 (MQ=255)
                                       tACCGGTGACGATTTCAGACATGGATTTGGTCACTGCCCCCGGCAGTATCGAGGCGCTGGATCAGCGCCTGATGCACGGCGTATTTGGCTACAgctgct                                                              <  1:2759776/99‑1 (MQ=255)
                                        aCCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTgg                                                            >  2:93915/1‑100 (MQ=255)
                                        aCCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTgg                                                            >  1:2374166/1‑100 (MQ=255)
                                         ccGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGa                                                           <  1:338711/100‑1 (MQ=255)
                                          cGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGaa                                                          >  1:1317216/1‑100 (MQ=255)
                                           gTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGaaa                                                         >  1:474172/1‑100 (MQ=255)
                                           gTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGaaa                                                         >  2:2009979/1‑100 (MQ=255)
                                           gTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGaaa                                                         >  2:10128/1‑100 (MQ=255)
                                           gTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGaaa                                                         >  1:2470971/1‑100 (MQ=255)
                                            ttCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGaaaa                                                        <  2:474373/100‑1 (MQ=255)
                                             tCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGaaaaa                                                       <  1:1206986/100‑1 (MQ=255)
                                              cACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAAc                                                      <  2:1836486/100‑1 (MQ=255)
                                               aCGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACg                                                     >  2:2533863/1‑100 (MQ=255)
                                                ctaTTTCAGACATGGATTTTGCCACTGCCCCCTGCATGATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGa                                                    <  1:1734814/98‑1 (MQ=255)
                                                cGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGa                                                    <  2:309787/100‑1 (MQ=255)
                                                cGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGa                                                    >  2:2619604/1‑100 (MQ=255)
                                                cGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGa                                                    <  2:2401438/100‑1 (MQ=255)
                                                  aTTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATg                                                  <  2:549814/100‑1 (MQ=255)
                                                    ttCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAg                                                <  2:933925/100‑1 (MQ=255)
                                                    ttCAGACATGGATTTCGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAg                                                >  1:861526/1‑100 (MQ=255)
                                                       aGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGttt                                             <  2:1214717/100‑1 (MQ=255)
                                                          cATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTtctc                                          >  2:2833497/1‑100 (MQ=255)
                                                           aTGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTcg                                         >  1:2538033/1‑100 (MQ=255)
                                                           aTGGATTATGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTcg                                         >  2:1522045/1‑100 (MQ=255)
                                                            tGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTcgc                                        >  1:1393815/1‑100 (MQ=255)
                                                             ggATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTcgcg                                       <  2:1844844/100‑1 (MQ=255)
                                                                ttttGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCt                                    >  2:2214239/1‑100 (MQ=255)
                                                                  ttGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTAt                                  >  2:2095314/1‑100 (MQ=255)
                                                                     ccACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGc                               >  2:1427439/1‑100 (MQ=255)
                                                                     ccACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGc                               <  1:2909323/100‑1 (MQ=255)
                                                                      cACTGCCCCCTGCATTATCGAGGCGCTGAATCAGTGCCTGATGGACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGcc                              <  2:815615/100‑1 (MQ=255)
                                                                       aCTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGccc                             >  1:1866091/1‑100 (MQ=255)
                                                                         tGCCCCCTGCATTTTCGAGGCGCTGAAGCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCAc                           <  2:1697102/100‑1 (MQ=255)
                                                                          gCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACt                          >  2:1314101/1‑100 (MQ=255)
                                                                           cccccTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAACCGATGAGTTTCTCGCGGCTATTGCCCACTg                         >  1:2528052/1‑100 (MQ=255)
                                                                             cccTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTgg                        >  1:2643608/1‑99 (MQ=255)
                                                                             cccTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGt                       >  1:1847431/1‑100 (MQ=255)
                                                                               cTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGttt                     <  2:258605/100‑1 (MQ=255)
                                                                                gccgtgcatcaGGCGCTGATTCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGtttt                    >  1:2888394/11‑100 (MQ=38)
                                                                                gccgtgcatcaGGCGCTGATTCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGtttt                    >  2:591604/11‑100 (MQ=38)
                                                                                  cgtgcatcaGGCGCTGATTCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTcc                  >  2:891013/9‑100 (MQ=39)
                                                                                  cATTATCGAGGCGCTGAATCAAAGCCTGATGCACGGCGTATTTGGCTAAAAAAAAAAAAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTc                   <  1:2559713/99‑1 (MQ=255)
                                                                                   gtgcatcaGGCGCTGATTCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCa                 >  2:2795039/8‑100 (MQ=39)
                                                                                       tCGAGGCGCTGAATCAGCGCCTGATGCACGGTGTATTTGGCTACAGCCGCTGGAAAAACGACGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCa             <  1:1433352/100‑1 (MQ=255)
                                                                                       tCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCAccc              >  2:746833/1‑99 (MQ=255)
                                                                                       tCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCa             <  1:1936238/100‑1 (MQ=255)
                                                                                       tCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCa             <  2:918543/100‑1 (MQ=255)
                                                                                       atcaGGCGCTGATTCCGCGCCTGATGCACGACGAATTTGGCTACTGCCGCTGGAAAAACGATGAGTTTCCAGCGACTATCGCCCACTGCAGTTCCACCCa             >  1:1916185/4‑100 (MQ=255)
                                                                                       atcaGGCGCTGATTCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCa             >  1:1669799/4‑100 (MQ=255)
                                                                                       atcaGGCGCTGATTCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCa             >  1:1428218/4‑100 (MQ=255)
                                                                                        tcaGGCGCTGATTCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAg            >  2:1725879/3‑100 (MQ=255)
                                                                                        cGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAg            <  2:1542088/100‑1 (MQ=255)
                                                                                         gAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGc           >  1:124903/1‑100 (MQ=255)
                                                                                         gAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGc           >  1:818518/1‑100 (MQ=255)
                                                                                         gAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGc           <  1:466953/100‑1 (MQ=255)
                                                                                         gAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGc           <  1:181443/100‑1 (MQ=255)
                                                                                         caGGCGCTGATTCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGc           >  2:1867733/2‑100 (MQ=255)
                                                                                          aGGCGCTGATTCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCa          >  2:2024485/1‑100 (MQ=255)
                                                                                            gcgcTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCAtt        <  2:1816302/100‑1 (MQ=255)
                                                                                              gcTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTac      >  1:116934/1‑100 (MQ=255)
                                                                                                tGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTacac    >  1:587400/1‑100 (MQ=255)
                                                                                                tGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTacac    >  2:484068/1‑100 (MQ=255)
                                                                                                 gAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACAcc   <  2:474175/100‑1 (MQ=255)
                                                                                                 gAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACAcc   >  2:1273330/1‑100 (MQ=255)
                                                                                                  aaTCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCg  <  2:1539984/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
GGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCG  >  NC_000913/1701009‑1701206

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: