Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,786,066 C→G 5.5% A142G (GCC→GGC)  nikE → nickel ABC transporter ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,786,0660CG 5.5% 547.5 / 16.0 164A142G (GCC→GGC) nikEnickel ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (69/86);  new base G (7/2);  total (76/88)
Fisher's exact test for biased strand distribution p-value = 8.26e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

Only 100 of 172 total aligned reads displayed.

CTGGTGTTCCAGGACGCCCCTGGCGCCTTCAACCCGCAGCGCAGCATCGGCTGGAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCTCGACGAAGCGTTGTCGAACCTGGACCGGCTGTTGCAGCTGCAGATCCTGC  >  NC_002947/3785918‑3786215
                                                                                                                                                    |                                                                                                                                                     
cTGGTGTTCCAGGACGCCCCTGGCGCCTTCAACCCGCAGCGCAGCATCGGCTGGAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGcc                                                                                                                                                      >  2:3233562/1‑150 (MQ=255)
 tGGTGTTCCAGGACGCCCCTGGCGCCTTCAACCCGCAGCGCAGCATCGGCTGGAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCg                                                                                                                                                     >  2:1144434/1‑150 (MQ=255)
  ggTGTTCCAGGACGCCCCTGGCGCCTTCAACCCGCAGCGCAGCATCGGCTGGAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGa                                                                                                                                                    <  2:2067415/150‑1 (MQ=255)
  ggTGTTCCAGGACGCCCCTGGCGCCTTCAACCCGCAGCGCAGCATCGGCTGGAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGa                                                                                                                                                    <  1:3925262/150‑1 (MQ=255)
  ggTGTTCCAGGACGCCCCTGGCGCCTTCAACCCGCAGCGCAGCATCGGCTGGAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGa                                                                                                                                                    <  1:3855212/150‑1 (MQ=255)
  ggTGTTCCAGGACGCCCCTGGCGCCTTCAACCCGCAGCGCAGCATCGGCTGGAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGa                                                                                                                                                    <  1:2192840/150‑1 (MQ=255)
      ttCCAGGACGCCCCTGGCGCCTTCAACCCGCAGCGCAGCATCGGCTGGAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCAt                                                                                                                                                >  1:3097172/1‑150 (MQ=255)
      ttCCAGGACGCCCCTGGCGCCTTCAACCCGCAGCGCAGCATCGGCTGGAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCAt                                                                                                                                                >  1:2743523/1‑150 (MQ=255)
      ttCCAGGACGCCCATGGCGCCTTCAACCCGCAGCGCAGCATCGGCTGGAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCCCTGCTGGAGGAAATGGGCCTGCGCGCCGAGCAt                                                                                                                                                >  2:3941679/1‑150 (MQ=255)
        ccAGGACGCCCCTGGCGCCTTCAACCCGCAGCGCAGCATCGGCTGGAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGc                                                                                                                                              >  1:979681/1‑150 (MQ=255)
        ccAGGACGCCCCTGGCGCCTTCAACCCGCAGCGCAGCATCGGCTGGAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGc                                                                                                                                              >  1:1351185/1‑150 (MQ=255)
         cAGGACGCCCCTGGCGCCTTCAACCCGCAGCGCAGCATCGGCTGGAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGcc                                                                                                                                             >  2:66593/1‑150 (MQ=255)
         cAGGACGCCCCTGGCGCCTTCAACCCGCAGCGCAGCATCGGCTGGAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGcc                                                                                                                                             <  1:3593091/150‑1 (MQ=255)
         cAGGACGCCCCTGGCGCCTTCAACCCGCAGCGCAGCATCGGCTGGAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCATGCTCGGCGCGCAGGcc                                                                                                                                             <  2:2835823/150‑1 (MQ=255)
          aGGACGCCCCTGGCGCCTTCAACCCGCAGCGCAGCATCGGCTGGAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGccc                                                                                                                                            >  2:3380238/1‑150 (MQ=255)
            gACGCCCCTGGCGCCTTCAACCCGCAGCGCAGCATCGGCTGGAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCaa                                                                                                                                          >  2:3591466/1‑150 (MQ=255)
              cGCCCCTGGCGCCTTCAACCCGCAGCGCAGCATCGGCTGGAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCATGCTCGGCGCGCAGGCCCattt                                                                                                                                        <  2:549145/150‑4 (MQ=255)
                ccccTGGCGCCTTCAACCCGCAGCGCAGCATCGGCTGGAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGcc                                                                                                                                      <  2:1613271/150‑1 (MQ=255)
                    tGGCGCCTTCAACCCGCAGCGCAGCATCGGCTGGAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCTAGCATGCCCAAAGCCTGcc                                                                                                                                  <  2:5911/150‑1 (MQ=255)
                          cTTCAACCCGCAGCGCAGCATCGGCTGGAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCgcagca                                                                                                                            >  2:1576621/1‑150 (MQ=255)
                           ttCAACCCGCAGCGCAGCATCGGCTGGAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGcagcag                                                                                                                           >  2:303268/1‑150 (MQ=255)
                           ttCAACCCGCAGCGCAGCATCGGCTGGAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGcagcag                                                                                                                           <  1:3964497/150‑1 (MQ=255)
                             cAACCCGCAGCGCAGCATCGGCTGGAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGtt                                                                                                                         >  2:350529/1‑150 (MQ=255)
                              aaCCCGCAGCGCAGCATCGGCTGGAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTg                                                                                                                        >  2:456725/1‑150 (MQ=255)
                               aCCCGCAGCGCAGCATCGGCTGGAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGa                                                                                                                       >  1:830544/1‑150 (MQ=255)
                                 ccGCAGCGCAGCATCGGCTGGAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGt                                                                                                                     >  2:4507339/1‑150 (MQ=255)
                                 ccGCAGCGCAGCATCGGCTGGAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGt                                                                                                                     >  1:4320346/1‑150 (MQ=255)
                                   gcagcgcagcATCGGCTGGAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTgg                                                                                                                   >  2:2200516/1‑150 (MQ=255)
                                    cagcgcagcATCGGCTGGAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTggc                                                                                                                  <  2:2691845/150‑1 (MQ=255)
                                    cagcgcagcATCGGCTGGAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTggc                                                                                                                  <  1:1313478/150‑1 (MQ=255)
                                     agcgcagcATCGGCTGGAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTggcg                                                                                                                 <  1:2471037/150‑1 (MQ=255)
                                     agcgcagcATCGGCTGGAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTggcg                                                                                                                 <  1:260848/150‑1 (MQ=255)
                                                ggCTGGAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCCGTTGCAg                                                                                                      >  1:2184318/1‑150 (MQ=255)
                                                 gCTGGAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGc                                                                                                     <  1:1437531/150‑1 (MQ=255)
                                                  cTGGAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCg                                                                                                    >  2:1386948/1‑150 (MQ=255)
                                                      aGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGcc                                                                                                <  1:1735943/150‑1 (MQ=255)
                                                         aTCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAAc                                                                                             <  1:3632322/150‑1 (MQ=255)
                                                         aTCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAAc                                                                                             <  2:480773/150‑1 (MQ=255)
                                                          tCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACa                                                                                            >  2:4029579/1‑150 (MQ=255)
                                                            gCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACAtt                                                                                          >  2:3584417/1‑150 (MQ=255)
                                                              cGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGc                                                                                        <  1:2016657/150‑1 (MQ=255)
                                                               gAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGcc                                                                                       >  1:1183687/1‑150 (MQ=255)
                                                                 gCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCg                                                                                     >  2:1777461/1‑150 (MQ=255)
                                                                  ccGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCgg                                                                                    <  1:2341902/150‑1 (MQ=255)
                                                                  ccGCTGCGCCATCTGACCGACATGGATGAAGCAGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGACTGCGCGCCGAGCATGCGCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCgg                                                                                    >  2:2708349/1‑150 (MQ=255)
                                                                     cTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGgcg                                                                                 >  2:4140872/1‑150 (MQ=255)
                                                                      tGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGgcgc                                                                                >  1:2426507/1‑150 (MQ=255)
                                                                        cgcCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTg                                                                              >  1:1619483/1‑150 (MQ=255)
                                                                          ccATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGTAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGgc                                                                            <  1:3746565/150‑1 (MQ=255)
                                                                          ccATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGgc                                                                            <  2:750046/150‑1 (MQ=255)
                                                                           cATCTGAGCGACATGGAAGAAGGCGCCCGGCAGGCGCGCTCCCTTGCACTGCTGGTGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCTGTTGCTGCGTGCCAACATTGCCCGGGCGCTGgcc                                                                           <  2:3801383/150‑1 (MQ=255)
                                                                                gACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTc                                                                      <  1:3435096/150‑1 (MQ=255)
                                                                                gACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTc                                                                      <  2:4353430/150‑1 (MQ=255)
                                                                                     aCATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCcac                                                                 >  2:2619097/1‑150 (MQ=255)
                                                                                     aCATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCcac                                                                 >  1:1956319/1‑150 (MQ=255)
                                                                                      cATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCcaca                                                                >  1:4687878/1‑150 (MQ=255)
                                                                                       aTGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAg                                                               >  1:3918889/1‑150 (MQ=255)
                                                                                         ggATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCt                                                             >  1:1072944/1‑150 (MQ=255)
                                                                                          gATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTg                                                            <  2:1939817/150‑1 (MQ=255)
                                                                                          gATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTg                                                            >  2:2242729/1‑150 (MQ=255)
                                                                                          gATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTg                                                            <  1:4837155/150‑1 (MQ=255)
                                                                                           aTGTAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTgg                                                           >  2:4789434/1‑150 (MQ=255)
                                                                                           aTGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTgg                                                           <  2:3915823/150‑1 (MQ=255)
                                                                                             gAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTa                                                         <  1:1066804/150‑1 (MQ=255)
                                                                                                  cgccCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCt                                                    >  2:3281679/1‑150 (MQ=255)
                                                                                                  cgccCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCt                                                    <  1:192301/150‑1 (MQ=255)
                                                                                                      cGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCTcgac                                                <  1:1327591/150‑1 (MQ=255)
                                                                                                              gcgcACCCTTGCACTGCTGGAGGACATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCTCGACGAAGCGtt                                        <  2:1721383/150‑1 (MQ=255)
                                                                                                              gcgcACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCTCGACGAAGCGtt                                        >  2:959894/1‑150 (MQ=255)
                                                                                                              gcgcACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCTCGACGAAGCGtt                                        <  1:803435/150‑1 (MQ=255)
                                                                                                              gcgcACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCTCGACGAAGCGtt                                        >  1:878007/1‑150 (MQ=255)
                                                                                                              gcgcACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCTCGACGAAGCGtt                                        <  2:2209384/150‑1 (MQ=255)
                                                                                                                     cTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCTCGACGAAGCGTTGTCGAAc                                 <  1:4498643/150‑1 (MQ=255)
                                                                                                                     cTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCTCGACGAAGCGTTGTCGAAc                                 >  2:1851680/1‑150 (MQ=255)
                                                                                                                          aCTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCTCGACGAAGCGTTGTCGAACCTGGa                            <  1:4615783/150‑1 (MQ=255)
                                                                                                                          aCTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCTCGACGAAGCGTTGTCGAACCTGGa                            <  1:4615580/150‑1 (MQ=255)
                                                                                                                          aCTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCTCGACGAAGCGTTGTCGAACCTGGa                            <  1:2040176/150‑1 (MQ=255)
                                                                                                                            tgctgGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCTCGACGAAGCGTTGTCGAACCTGGAcc                          >  2:2519941/1‑150 (MQ=255)
                                                                                                                               tgGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCTCGACGAAGCGTTGTCGAACCTGGACCGGc                       >  2:2414950/1‑150 (MQ=255)
                                                                                                                               tgGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCTCGACGAAGCGTTGTCGAACCTGGACCGGc                       >  1:3370257/1‑150 (MQ=255)
                                                                                                                                 gaggaAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCTCGACGAAGCGTTGTCGAACCTGGACCGGCTg                     >  1:3484191/1‑150 (MQ=255)
                                                                                                                                 gaggaAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCTCGACGAAGCGTTGTCGAACCTGGACCGGCTg                     >  2:195539/1‑150 (MQ=255)
                                                                                                                                 gaggaAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCTCGACGAAGCGTTGTCGAACCTGGACCGGCTg                     >  1:911838/1‑150 (MQ=255)
                                                                                                                                  aggaAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCTCGACGAAGCGTTGTCGAACCTGGACCGGCTGt                    >  1:4202064/1‑150 (MQ=255)
                                                                                                                                   ggaAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCTCGACGAAGCGTTGTCGAACCTGGACCGGCTGtt                   <  1:1997858/150‑1 (MQ=255)
                                                                                                                                      aaTGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCTCGACGAAGCGTTGTCGAACCTGGACCGGCTGTTGCa                >  1:2521422/1‑150 (MQ=255)
                                                                                                                                       aTGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCTCGACGAAGCGTTGTCGAACCTGGACCGGCTGTTGCAg               <  1:4003243/150‑1 (MQ=255)
                                                                                                                                       aTGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCTCGACGAAGCGTTGTCGAACCTGGACCGGCTGTTGCAg               <  1:3233562/150‑1 (MQ=255)
                                                                                                                                          ggCATGCTCGGCGCGCAGGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCTCGACGAAGCGTTGTCGAACCTGGACCGGCTGTTGCAGCTg            >  2:3497342/1‑150 (MQ=255)
                                                                                                                                          ggCATGCTCGGCGCGCAGGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCTCGACGAAGCGTTGTCGAACCTGGACCGGCTGTTGCAGCTg            >  2:4082034/1‑150 (MQ=255)
                                                                                                                                           gcatGCTCGGCGCGCAGGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCTCGACGAAGCGTTGTCGAACCTGGACCGGCTGTTGCAGCTGc           >  1:2544230/4‑150 (MQ=255)
                                                                                                                                            ccTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCTGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCTCGACGAAGCGTTGTCGAACCTGGACCGGCTGTTGCAGCTGCa          <  2:2743523/150‑1 (MQ=255)
                                                                                                                                            ccTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCTCGACGAAGCGTTGTCGAACCTGGACCGGCTGTTGCAGCTGCa          <  2:3097172/150‑1 (MQ=255)
                                                                                                                                              tGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCTCGACGAAGCGTTGTCGAACCTGGACCGGCTGTTGCAGCTGCAGa        <  1:2565060/150‑1 (MQ=255)
                                                                                                                                              tGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCTCGACGAAGCGTTGTCGAACCTGGACCGGCTGTTGCAGCTGCAGa        <  2:4604432/150‑1 (MQ=255)
                                                                                                                                              tGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCTCGACGAAGCGTTGTCGAACCTGGACCGGCTGTTGCAGCTGCAGa        >  1:2940622/1‑150 (MQ=255)
                                                                                                                                                 gcgcCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCTCGACGAAGCGTTGTCGAACCTGGACCGGCTGTTGCAGCTGCAGATcc     >  2:4144281/1‑150 (MQ=255)
                                                                                                                                                 gcgcCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCTCGACGAAGCGTTGTCGAACCTGGACCGGCTGTTGCAGCTGCAGATcc     <  2:4222215/150‑1 (MQ=255)
                                                                                                                                                   gcCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCTCGACGAAGCGTTGTCGAACCTGGACCGGCTGTTGCAGCTGCAGATCCTg   >  1:2862012/1‑150 (MQ=255)
                                                                                                                                                   gcCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCTCGACGAAGCGTTGTCGAACCTGGACCGGCTGTTGCAGCTGCAGATCCTg   <  2:741522/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CTGGTGTTCCAGGACGCCCCTGGCGCCTTCAACCCGCAGCGCAGCATCGGCTGGAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCTCGACGAAGCGTTGTCGAACCTGGACCGGCTGTTGCAGCTGCAGATCCTGC  >  NC_002947/3785918‑3786215

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