Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 966,368 A→G 25.5% A683A (GCA→GCG ycaI → ComEC family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913966,3680AG25.5% 184.0 / 72.6 117A683A (GCA→GCGycaIComEC family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (42/45);  new base G (17/13);  total (59/58)
Fisher's exact test for biased strand distribution p-value = 5.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.97e-01

Only 100 of 119 total aligned reads displayed.

TGGCAGCAAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACAT  >  NC_000913/966270‑966467
                                                                                                  |                                                                                                   
tGGCAGCAAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCAc                                                                                                    >  2:3160046/1‑100 (MQ=255)
tGGCAGCAAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCAc                                                                                                    >  2:1783626/1‑100 (MQ=255)
 ggCAGCAAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACt                                                                                                   >  1:2554958/1‑100 (MQ=255)
 ggCAGCAAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACt                                                                                                   <  1:616817/100‑1 (MQ=255)
   cagcaAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCCTTAATTCAGCGAGTGAATGGAAAAGTGGCAc                                                                                                    >  2:1594750/1‑97 (MQ=255)
   cagcaAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCg                                                                                                 >  1:1937258/1‑100 (MQ=255)
   cagcaAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCg                                                                                                 <  1:2271083/100‑1 (MQ=255)
   cagcaAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCg                                                                                                 <  1:1276534/100‑1 (MQ=255)
    agcaAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCGTCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCgc                                                                                                <  1:3176668/100‑1 (MQ=255)
    agcaAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCgc                                                                                                >  1:170791/1‑100 (MQ=255)
     gcaAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCgca                                                                                               >  1:2709598/1‑100 (MQ=255)
        aGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCgcatcg                                                                                            <  1:1929012/100‑1 (MQ=255)
        aGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCgcatcg                                                                                            >  2:1019545/1‑100 (MQ=255)
          tGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCgcatcggc                                                                                          <  2:2799529/100‑1 (MQ=255)
            cAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCgcatcggcan                                                                                        <  1:1812218/100‑2 (MQ=255)
            cAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCgcatcggcat                                                                                        >  1:3579718/1‑100 (MQ=255)
             aGGCAACATTGCTTCAGGTACCTCGCGATGGCAGTGATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAAGGGAAAAGTGGCACTCgcatcggcatc                                                                                       <  2:173487/100‑1 (MQ=255)
              ggCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCgcatcggcatcg                                                                                      >  1:3805821/1‑100 (MQ=255)
                cAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATcgcg                                                                                    >  1:2879947/1‑100 (MQ=255)
                cAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATcgcg                                                                                    <  2:2867536/100‑1 (MQ=255)
                   cATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCta                                                                                 >  1:3597577/1‑100 (MQ=255)
                     ttGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCtata                                                                               >  1:2366348/1‑100 (MQ=255)
                      tGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATaa                                                                              >  1:2937888/1‑100 (MQ=255)
                       gCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAAc                                                                             >  2:3208078/1‑100 (MQ=255)
                        cTTCAGGTACCTCACCATGTCAGTAATACCTCATCATCGTTGCCATTAATTCAGCAAATGGATGGAAAAGTGGCACTCGCATCGACATCGCGCTATGACg                                                                            >  2:657768/1‑100 (MQ=255)
                        cTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACg                                                                            >  1:3856829/1‑100 (MQ=255)
                        cTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGGGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAAc                                                                             >  1:2897360/1‑99 (MQ=255)
                          tCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCa                                                                          >  1:3989712/1‑100 (MQ=255)
                          tCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCa                                                                          <  2:1505165/100‑1 (MQ=255)
                              gTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGc                                                                      <  2:291096/100‑1 (MQ=255)
                              gTAACTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGGGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGc                                                                      <  2:3729721/100‑1 (MQ=255)
                                aCCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGa                                                                    >  2:1646542/1‑100 (MQ=255)
                                 ccTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGAc                                                                   >  2:3402113/1‑100 (MQ=255)
                                  cTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGATCGCATCGGCATCGCGCTATAACGCATGGCGACt                                                                  >  2:1216991/1‑100 (MQ=255)
                                  cTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACt                                                                  <  1:1804624/100‑1 (MQ=255)
                                    cACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGc                                                                <  2:2210747/100‑1 (MQ=255)
                                    cACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGc                                                                >  1:2781838/1‑100 (MQ=255)
                                      ccATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGccc                                                              <  1:2325117/100‑1 (MQ=255)
                                         tGGCAGTAATACCTTATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCctct                                                           <  2:2813202/100‑1 (MQ=255)
                                          ggCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTa                                                          <  2:2464827/100‑1 (MQ=255)
                                           gCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTaa                                                         >  2:3771466/1‑100 (MQ=255)
                                           gCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTaa                                                         >  1:1790140/1‑100 (MQ=255)
                                           gCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTaa                                                         <  1:1099274/100‑1 (MQ=255)
                                           gCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTaa                                                         >  2:3668060/1‑100 (MQ=255)
                                            cAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAAc                                                        <  2:1458090/100‑1 (MQ=255)
                                            cAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAAc                                                        <  2:1655378/100‑1 (MQ=255)
                                               tAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACaaa                                                     <  1:179734/100‑1 (MQ=255)
                                                aaTACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTTCCCTCTAActttg                                                    <  2:2559617/100‑5 (MQ=255)
                                                   aCCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAACAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTa                                                 >  1:3889071/1‑100 (MQ=255)
                                                   aCCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTa                                                 <  1:3661746/100‑1 (MQ=255)
                                                   aCCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACcaa                                                     >  2:3963118/1‑96 (MQ=255)
                                                   aCCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTa                                                 >  1:1681912/1‑100 (MQ=255)
                                                    ccTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTaa                                                <  1:2123271/100‑1 (MQ=255)
                                                      tcatcaTCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAg                                               <  2:1557300/99‑1 (MQ=255)
                                                       catcatCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCa                                             >  2:2937381/1‑100 (MQ=255)
                                                        atcatcGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCAt                                            <  2:1105162/100‑1 (MQ=255)
                                                        atcatcGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCAt                                            >  1:3258101/1‑100 (MQ=255)
                                                         tcatcGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATc                                           <  2:2128541/100‑1 (MQ=255)
                                                          catcGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATc                                           <  1:38390/99‑1 (MQ=255)
                                                               ttGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATc                                     <  1:466217/100‑1 (MQ=255)
                                                                tGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGTGACTGCCATCTAACAAAGTTAAGCATCGCTATc                                     <  2:949983/99‑1 (MQ=255)
                                                                 gCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCaa                                   >  2:358615/1‑100 (MQ=255)
                                                                  ccATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAAc                                  <  1:1898741/100‑1 (MQ=255)
                                                                   cATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACt                                 >  2:3638186/1‑100 (MQ=255)
                                                                   cATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACt                                 >  2:2751561/1‑100 (MQ=255)
                                                                    aTTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTg                                <  1:2813942/100‑1 (MQ=255)
                                                                     ttAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACGGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGc                               <  2:1714354/100‑1 (MQ=255)
                                                                       aaTTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCaa                             >  2:335092/1‑100 (MQ=255)
                                                                       aaTTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCaa                             <  1:2384081/100‑1 (MQ=255)
                                                                         ttCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAgg                           <  1:1797919/100‑1 (MQ=255)
                                                                           cAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGat                         <  1:1019551/100‑1 (MQ=255)
                                                                               gAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCa                     <  2:3872663/100‑1 (MQ=255)
                                                                                  tGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATg                  <  1:358618/100‑1 (MQ=255)
                                                                                   gAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATgg                 >  1:1508094/1‑100 (MQ=255)
                                                                                   gAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATgg                 >  1:3427999/1‑100 (MQ=255)
                                                                                   gAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGTCTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATgg                 >  2:1235415/1‑100 (MQ=255)
                                                                                   gAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATgg                 >  2:3884864/1‑100 (MQ=255)
                                                                                   gAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATgg                 >  1:1513568/1‑100 (MQ=255)
                                                                                    aaTGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGa                >  1:1248125/1‑100 (MQ=255)
                                                                                    aaTGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGa                >  1:362504/1‑100 (MQ=255)
                                                                                     aTGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGAt               >  1:1592857/1‑100 (MQ=255)
                                                                                      tGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGAtt              >  1:1320048/1‑100 (MQ=255)
                                                                                      tGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGAtt              <  2:1681900/100‑1 (MQ=255)
                                                                                      tGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGAtt              <  1:2654470/100‑1 (MQ=255)
                                                                                      tGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGAtt              <  1:3788237/100‑1 (MQ=255)
                                                                                      tGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGAtt              <  1:4055237/100‑1 (MQ=255)
                                                                                      tGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGAtt              <  1:488784/100‑1 (MQ=255)
                                                                                        gAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGa            >  2:722487/1‑100 (MQ=255)
                                                                                        gAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGa            >  2:1859952/1‑100 (MQ=255)
                                                                                         aaaaGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGAt           <  1:1081116/100‑1 (MQ=255)
                                                                                          aaaGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATa          >  1:2241479/1‑100 (MQ=255)
                                                                                           aaGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATAc         <  2:989363/100‑1 (MQ=255)
                                                                                            aGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACt        >  2:3428675/1‑100 (MQ=255)
                                                                                            aGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACt        >  1:1490521/1‑100 (MQ=255)
                                                                                              tgacACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTcc      <  1:1504688/97‑1 (MQ=255)
                                                                                              tGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGGCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTcc      <  2:4054576/100‑1 (MQ=255)
                                                                                               ggCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACACAGTTAAGCATCGCTATCCACAGCAAGGATATCAATGGATTGATACTc       >  2:599393/1‑98 (MQ=255)
                                                                                                gCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCcac    <  1:4086171/100‑1 (MQ=255)
                                                                                                gCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATCCATTGATatccttg    <  2:3207025/100‑7 (MQ=255)
                                                                                                  aCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACAt  <  2:362501/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
TGGCAGCAAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACAT  >  NC_000913/966270‑966467

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