Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,983,761 C→T 100% E437K (GAA→AAA)  araG ← L‑arabinose ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,983,7610CT97.1% 331.6 / ‑0.3 103E437K (GAA→AAA) araGL‑arabinose ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base T (54/46);  total (55/48)
Fisher's exact test for biased strand distribution p-value = 5.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

Only 100 of 103 total aligned reads displayed.

CCGGTCGGCAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTCGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGCTTTTTGCT  >  NC_000913/1983666‑1983856
                                                                                               |                                                                                               
cccgTCTGCAACTCCGAGGACTTCAGGTTAGTCGCTGGATGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGc                                                                                             <  2:3637128/97‑1 (MQ=255)
ccGGTCGGCAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGc                                                                                             <  1:2082437/100‑1 (MQ=255)
ccGGTCGGCAACGCCGAGGACTTCAGGTAAGTCGCTCGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGc                                                                                             <  1:2904443/100‑1 (MQ=255)
    tCGGCAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAg                                                                                         <  2:3556999/100‑1 (MQ=255)
    tCGGCAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAg                                                                                         <  2:3557784/100‑1 (MQ=255)
     cGGCAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAgc                                                                                        >  1:983467/1‑100 (MQ=255)
      ggCAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAgcg                                                                                       <  2:3861413/100‑1 (MQ=255)
      ggCAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAgcg                                                                                       >  2:1517489/1‑100 (MQ=255)
       gCAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAgcgc                                                                                      <  1:2257883/100‑1 (MQ=255)
        cAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGcc                                                                                     >  1:2122891/1‑100 (MQ=255)
             ccGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACa                                                                                 >  1:3607653/1‑99 (MQ=255)
              cGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATc                                                                               <  2:1439146/100‑1 (MQ=255)
                   aCTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGc                                                                          <  1:630643/100‑1 (MQ=255)
                   aCTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGc                                                                          >  1:2813472/1‑100 (MQ=255)
                    cTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGcc                                                                         >  1:3264783/1‑100 (MQ=255)
                    cTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGcc                                                                         >  2:852693/1‑100 (MQ=255)
                    cTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGcc                                                                         <  1:3828591/100‑1 (MQ=255)
                     ttCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCcg                                                                        <  1:2222046/100‑1 (MQ=255)
                     ttCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTCGTGCTTAGCGCCAACATCAATGCcg                                                                        <  1:2284506/100‑1 (MQ=255)
                      tCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCcgc                                                                       >  2:3058894/1‑100 (MQ=255)
                       cAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGTATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCcgcg                                                                      <  1:3687835/100‑1 (MQ=255)
                       cAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCcgcg                                                                      <  2:3859270/100‑1 (MQ=255)
                         ggTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCcgcgcg                                                                    >  1:1864733/1‑100 (MQ=255)
                         ggTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCcgcg                                                                      >  2:2307120/1‑98 (MQ=255)
                          gTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGt                                                                   >  1:804840/1‑100 (MQ=255)
                          gTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGt                                                                   >  1:3066222/1‑100 (MQ=255)
                          gTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGt                                                                   >  2:2983905/1‑100 (MQ=255)
                           tAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGTGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTa                                                                  <  1:741764/100‑1 (MQ=255)
                            aaGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAg                                                                 >  1:73062/1‑100 (MQ=255)
                              gTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGt                                                               <  1:228568/100‑1 (MQ=255)
                              gTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGt                                                               >  1:2743457/1‑100 (MQ=255)
                              gTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGt                                                               <  2:1609719/100‑1 (MQ=255)
                                cggtGGAGGCAAACAGCACCGCCACGCCCTTCGCCGGCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGtt                                                              <  2:2783693/96‑1 (MQ=255)
                                cGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTc                                                             >  2:3960654/1‑100 (MQ=255)
                                 gCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCa                                                            <  1:159584/100‑1 (MQ=255)
                                 gCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCa                                                            >  1:2472065/1‑100 (MQ=255)
                                 gCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCa                                                            >  1:4103363/1‑100 (MQ=255)
                                   tGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGTGCATAAATTACGTTATATATTTCGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATc                                                          <  2:2782199/100‑1 (MQ=255)
                                    ggAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATcc                                                         >  1:3065453/1‑100 (MQ=255)
                                    ggAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATc                                                          <  1:2443645/99‑1 (MQ=255)
                                      aGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCATAGGTTCATCCAg                                                       <  2:2111499/100‑1 (MQ=255)
                                       ggCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGc                                                      <  2:2112675/100‑1 (MQ=255)
                                       ggCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGc                                                      >  1:1529657/1‑100 (MQ=255)
                                       ggCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTGATATATTTTGTGCTTAGCGCCAACATCAATGCCGCCCGTAGGTTCATCCAGc                                                      >  2:859349/1‑100 (MQ=255)
                                         cagaCAGCACCGTCACGCCCTGCGCCGCCAGCGACTAAAATACGGTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCa                                                     >  2:1987859/4‑99 (MQ=255)
                                         cAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCaa                                                    >  1:833528/1‑100 (MQ=255)
                                         cAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAAGTACGGTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTATGTTCATCCAGCaa                                                    >  2:1019116/1‑100 (MQ=255)
                                          aaaCAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCaa                                                    >  2:1136445/1‑99 (MQ=255)
                                           aaCAGCACCGCCACGCCCTGCGCCGCCAGCGCGTAAATTTCGATATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCaaaa                                                  <  2:2902641/100‑1 (MQ=255)
                                           aaCAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCaaaa                                                  >  1:178079/1‑100 (MQ=255)
                                            aCAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAAt                                                 >  1:1094482/1‑100 (MQ=255)
                                                cACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGAcc                                             >  2:1306084/1‑100 (MQ=255)
                                                 aCCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCt                                            <  1:109171/100‑1 (MQ=255)
                                                 aCCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCt                                            >  1:3497243/1‑100 (MQ=255)
                                                  ccGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCtt                                           >  1:1480918/1‑100 (MQ=255)
                                                  ccGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCtt                                           >  2:1258652/1‑100 (MQ=255)
                                                  ccGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGGAGGTTCATCCAGCAAAATGACCtt                                           >  2:3069659/1‑100 (MQ=255)
                                                    gCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCa                                         >  2:2419487/1‑100 (MQ=255)
                                                    gCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCa                                         >  2:3047064/1‑100 (MQ=255)
                                                    gCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCa                                         >  1:248295/1‑100 (MQ=255)
                                                     ccACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCAt                                        >  2:1044710/1‑100 (MQ=255)
                                                      cACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATACTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCAt                                        <  2:3826551/99‑1 (MQ=255)
                                                         gCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATctct                                    >  1:1420222/1‑100 (MQ=255)
                                                          cccTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCtt                                   >  1:1669237/1‑100 (MQ=255)
                                                           ccTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTc                                  <  1:3205300/100‑1 (MQ=255)
                                                           ccTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTc                                  >  2:4011163/1‑100 (MQ=255)
                                                            cTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTcc                                 <  2:2511718/100‑1 (MQ=255)
                                                             tGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCg                                >  1:1902061/1‑100 (MQ=255)
                                                              gcgcCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGa                               <  1:19575/100‑1 (MQ=255)
                                                               cgccgcCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGAt                              <  1:1343628/100‑1 (MQ=255)
                                                               cgccgcCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGAt                              <  2:1638087/100‑1 (MQ=255)
                                                                gccgccAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATa                             >  1:68048/1‑100 (MQ=255)
                                                                  cgccAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAAc                           <  1:2347257/100‑1 (MQ=255)
                                                                   gccAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAAcc                          >  2:1187169/1‑100 (MQ=255)
                                                                     cAGCGCATAAATTACGTTATATATTTTGTGCTTATCGCCAACATCAATGCCGCGCGTATGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAAc                           <  2:3753807/97‑1 (MQ=255)
                                                                      aGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGc                       <  2:3764252/100‑1 (MQ=255)
                                                                       gcgcATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCg                      >  2:3843145/1‑100 (MQ=255)
                                                                       gcacATGAATTACGTGGTAGATTTTTGGCTTATCTCCAGCATCAATGCGGCGCGTAGGTGCATCCAGCAAAATTAGCTTCAGCTCTTCCGATAACCAGCg                      >  2:3448861/4‑100 (MQ=255)
                                                                        cgcATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCgg                     >  1:950766/1‑100 (MQ=255)
                                                                        cgcATAAATTACGTTATATATTTCGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCgg                     >  2:376636/1‑100 (MQ=255)
                                                                                ttACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAAt             <  2:1529645/100‑1 (MQ=255)
                                                                                ttACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAAt             <  1:174611/100‑1 (MQ=255)
                                                                                ttACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAAt             <  1:1436174/100‑1 (MQ=255)
                                                                                 tACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATg            <  1:3825555/100‑1 (MQ=255)
                                                                                    gTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGCt         >  1:1498052/1‑100 (MQ=255)
                                                                                     ttATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGCtt        <  1:2271086/100‑1 (MQ=255)
                                                                                     ttATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGCtt        >  1:3408672/1‑100 (MQ=255)
                                                                                     ttATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCACAGAATGGCtt        >  2:1475974/1‑100 (MQ=255)
                                                                                     ttATATATTTTGTGCTCAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGCtt        >  2:3063287/1‑100 (MQ=255)
                                                                                      tatataTTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTCGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGAGGCCCAGACTGGCttt       >  1:3281765/1‑100 (MQ=255)
                                                                                       atatatTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGCtttt      <  2:2007744/100‑1 (MQ=255)
                                                                                        tatatTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGCttttt     <  2:2085709/100‑1 (MQ=255)
                                                                                         atatTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGCTTTTTg    <  1:2187347/100‑1 (MQ=255)
                                                                                         atatTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGCTTTTTg    <  2:804834/100‑1 (MQ=255)
                                                                                          tatTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGCTTTTTGc   <  1:2348856/100‑1 (MQ=255)
                                                                                          tatTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGCTTTTTGc   <  2:559742/100‑1 (MQ=255)
                                                                                           atTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGCTTTTTGCt  >  1:3220944/1‑100 (MQ=255)
                                                                                           atTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGCTTTTTGCt  >  1:1833784/1‑100 (MQ=255)
                                                                                           atTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGCTTTTTGCt  >  1:3958741/1‑100 (MQ=255)
                                                                                           atTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGCTTTTTGCt  >  1:4098402/1‑100 (MQ=255)
                                                                                               |                                                                                               
CCGGTCGGCAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTCGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGCTTTTTGCT  >  NC_000913/1983666‑1983856

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: