Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,786,063 G→T 5.4% R141L (CGC→CTC)  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,0630GT 5.4% 556.8 / 14.3 165R141L (CGC→CTC) nikEnickel ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (70/86);  new base T (7/2);  total (77/88)
Fisher's exact test for biased strand distribution p-value = 8.37e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 168 total aligned reads displayed.

CTGGTGTTCCAGGACGCCCCTGGCGCCTTCAACCCGCAGCGCAGCATCGGCTGGAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCTCGACGAAGCGTTGTCGAACCTGGACCGGCTGTTGCAGCTGCAGATCC  >  NC_002947/3785918‑3786212
                                                                                                                                                 |                                                                                                                                                     
 tGGTGTTCCAGGACGCCCCTGGCGCCTTCAACCCGCAGCGCAGCATCGGCTGGAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCg                                                                                                                                                  >  2:1144434/1‑150 (MQ=255)
 tGGTGTTCCAGGACGCCCCTGGCGCCTTCAACCCGCAGCGCAGCATCGGCTGGAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCg                                                                                                                                                  >  2:3500329/1‑150 (MQ=255)
  ggTGTTCCAGGACGCCCCTGGCGCCTTCAACCCGCAGCGCAGCATCGGCTGGAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGa                                                                                                                                                 <  1:3925262/150‑1 (MQ=255)
  ggTGTTCCAGGACGCCCCTGGCGCCTTCAACCCGCAGCGCAGCATCGGCTGGAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGa                                                                                                                                                 <  1:3855212/150‑1 (MQ=255)
  ggTGTTCCAGGACGCCCCTGGCGCCTTCAACCCGCAGCGCAGCATCGGCTGGAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGa                                                                                                                                                 <  1:2192840/150‑1 (MQ=255)
   gtgtTCCAGGACGCCCCTGGCGCCTTCAACCCGCAGCGCAGCATCGGCTGGAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAg                                                                                                                                                >  2:3584350/1‑150 (MQ=255)
      ttCCAGGACGCCCCTGGCGCCTTCAACCCGCAGCGCAGCATCGGCTGGAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCAt                                                                                                                                             >  1:3097172/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                                                                                                                                          <  1:3593091/150‑1 (MQ=255)
         cAGGACGCCCCTGGCGCCTTCAACCCGCAGCGCAGCATCGGCTGGAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGcc                                                                                                                                          >  2:2734176/1‑150 (MQ=255)
         cAGGACGCCCCTGGCGCCTTCAACCCGCAGCGCAGCATCGGCTGGAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCGGGAGGAACTGGGCCTGCGCGCCGAGCCTGcc                                                                                                                                          >  2:2674894/1‑150 (MQ=255)
          aGGACGCCCCTGGCGCCTTCAACCCGCAGCGCAGCATCGGCTGGAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGccc                                                                                                                                         >  1:1954671/1‑150 (MQ=255)
          aGGACGCCCCTGGCGCCTTCAACCCGCAGCGCAGCATCGGCTGGAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGccc                                                                                                                                         >  2:3380238/1‑150 (MQ=255)
              cGCCCCTGGCGCCTTCAACCCGCAGCGCAGCATCGGCTGGAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCATGCTCGGCGCGCAGGCCCattt                                                                                                                                     <  2:549145/150‑4 (MQ=255)
                ccccTGGCGCCTTCAACCCGCAGCGCAGCATCGGCTGGAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGcc                                                                                                                                   <  2:1613271/150‑1 (MQ=255)
                 cccTGGCGCCTTCAACCCGCAGCGCAGCATCGGCTGGAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCt                                                                                                                                  >  2:3130213/1‑150 (MQ=255)
                 cccTGGCGCCTTCAACCCGCAGCGCAGCATCGGCTGGAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCt                                                                                                                                  >  1:2386801/1‑150 (MQ=255)
                    tGGCGCCTTCAACCCGCAGCGCAGCATCGGCTGGAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCTAGCATGCCCAAAGCCTGcc                                                                                                                               <  2:5911/150‑1 (MQ=255)
                          cTTCAACCCGCAGCGCAGCATCGGCTGGAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCgcagca                                                                                                                         >  2:1576621/1‑150 (MQ=255)
                           ttCAACCCGCAGCGCAGCATCGGCTGGAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGcagcag                                                                                                                        >  1:4604432/1‑150 (MQ=255)
                           ttCAACCCGCAGCGCAGCATCGGCTGGAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGcagcag                                                                                                                        <  1:3964497/150‑1 (MQ=255)
                           ttCAACCCGCAGCGCAGCATCGGCTGGAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGcagcag                                                                                                                        >  2:303268/1‑150 (MQ=255)
                             cAACCCGCAGCGCAGCATCGGCTGGAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGtt                                                                                                                      <  1:3598433/150‑1 (MQ=255)
                             cAACCCGCAGCGCAGCATCGGCTGGAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGtt                                                                                                                      <  2:2559960/150‑1 (MQ=255)
                              aaCCCGCAGCGCAGCATCGGCTGGAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTg                                                                                                                     >  2:456725/1‑150 (MQ=255)
                               aCCCGCAGCGCAGCATCGGCTGGAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGa                                                                                                                    >  2:2818378/1‑150 (MQ=255)
                               aCCCGCAGCGCAGCATCGGCTGGAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGa                                                                                                                    <  2:2274844/150‑1 (MQ=255)
                                cccGCAGCGCAGCATCGGCTGGAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAg                                                                                                                   >  2:3301646/1‑150 (MQ=255)
                                 ccGCAGCGCAGCATCGGCTGGAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGt                                                                                                                  >  1:4320346/1‑150 (MQ=255)
                                 ccGCAGCGCAGCATCGGCTGGAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGt                                                                                                                  >  2:4507339/1‑150 (MQ=255)
                                   gcagcgcagcATCGGCTGGAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTgg                                                                                                                >  2:2200516/1‑150 (MQ=255)
                                   gcagcgcagcATCGGCTGGAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCATGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTgg                                                                                                                >  2:62177/1‑150 (MQ=255)
                                    cagcgcagcATCGGCTGGAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTggc                                                                                                               <  1:1313478/150‑1 (MQ=255)
                                     agcgcagcATCGGCTGGAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTggcg                                                                                                              <  1:2471037/150‑1 (MQ=255)
                                                 gCTGGAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGc                                                                                                  <  1:1437531/150‑1 (MQ=255)
                                                  cTGGAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCg                                                                                                 >  2:1386948/1‑150 (MQ=255)
                                                    ggAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTg                                                                                               <  2:326310/150‑1 (MQ=255)
                                                      aGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGcc                                                                                             <  1:1735943/150‑1 (MQ=255)
                                                         aTCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAAc                                                                                          <  2:480773/150‑1 (MQ=255)
                                                            gCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACAtt                                                                                       >  2:3584417/1‑150 (MQ=255)
                                                              cGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGc                                                                                     <  1:2016657/150‑1 (MQ=255)
                                                               gAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGcc                                                                                    >  1:1183687/1‑150 (MQ=255)
                                                               gAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGcc                                                                                    >  2:1997795/1‑150 (MQ=255)
                                                                 gCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCg                                                                                  >  2:1777461/1‑150 (MQ=255)
                                                                  ccGCTGCGCCATCTGACCGACATGGATGAAGCAGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGACTGCGCGCCGAGCATGCGCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCgg                                                                                 >  2:2708349/1‑150 (MQ=255)
                                                                     cTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGgcg                                                                              >  2:4140872/1‑150 (MQ=255)
                                                                        cgcCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTg                                                                           <  1:2036035/150‑1 (MQ=255)
                                                                          ccATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGTAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGgc                                                                         <  1:3746565/150‑1 (MQ=255)
                                                                          ccATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGgc                                                                         <  2:750046/150‑1 (MQ=255)
                                                                          ccATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGgc                                                                         <  1:3869582/150‑1 (MQ=255)
                                                                            aTCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTgaaag                                                                       <  2:2235513/150‑5 (MQ=255)
                                                                                gACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTc                                                                   <  1:3435096/150‑1 (MQ=255)
                                                                                gACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTc                                                                   <  2:4353430/150‑1 (MQ=255)
                                                                                  ccGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGc                                                                 >  2:3728396/1‑150 (MQ=255)
                                                                                      cATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCcaca                                                             >  1:4687855/1‑150 (MQ=255)
                                                                                         ggATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCt                                                          >  1:1072944/1‑150 (MQ=255)
                                                                                          gATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTg                                                         <  2:1939817/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)
                                                                                                gccgccCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTg                                                   >  2:4171405/1‑150 (MQ=255)
                                                                                                  cgccCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCt                                                 <  1:192301/150‑1 (MQ=255)
                                                                                                  cgccCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCt                                                 <  1:4813516/150‑1 (MQ=255)
                                                                                                     ccGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCTcga                                              >  2:1815605/1‑150 (MQ=255)
                                                                                                      cGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCTcgac                                             <  1:1327591/150‑1 (MQ=255)
                                                                                                         cAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCTCGACGaa                                          <  2:3884330/150‑1 (MQ=255)
                                                                                                         cAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCTCGACGaa                                          <  1:403101/150‑1 (MQ=255)
                                                                                                              gcgcACCCTTGCACTGCTGGAGGACATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCTCGACGAAGCGtt                                     <  2:1721383/150‑1 (MQ=255)
                                                                                                              gcgcACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCTCGACGAAGCGtt                                     <  1:803435/150‑1 (MQ=255)
                                                                                                              gcgcACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCTCGACGAAGCGtt                                     >  1:878007/1‑150 (MQ=255)
                                                                                                              gcgcACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCTCGACGAAGCGtt                                     <  2:2209384/150‑1 (MQ=255)
                                                                                                              gcgcACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCTCGACGAAGCGtt                                     >  2:959894/1‑150 (MQ=255)
                                                                                                                 cACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCTCGACGAAGCGTTGTc                                  <  1:2734176/150‑1 (MQ=255)
                                                                                                                        gCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCTCGACGAAGCGTTGTCGAACCTg                           <  1:4463230/150‑1 (MQ=255)
                                                                                                                          aCTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCTCGACGAAGCGTTGTCGAACCTGGa                         <  1:4615783/150‑1 (MQ=255)
                                                                                                                          aCTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCTCGACGAAGCGTTGTCGAACCTGGa                         <  2:2853443/150‑1 (MQ=255)
                                                                                                                           ctgctgGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCTCGACGAAGCGTTGTCGAACCTGGAc                        <  1:3986649/150‑1 (MQ=255)
                                                                                                                            tgctgGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCTCGACGAAGCGTTGTCGAACCTGGAcc                       >  2:2519941/1‑150 (MQ=255)
                                                                                                                               tgGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCTCGACGAAGCGTTGTCGAACCTGGACCGGc                    >  1:3370257/1‑150 (MQ=255)
                                                                                                                                 gaggaAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCTCGACGAAGCGTTGTCGAACCTGGACCGGCTg                  >  2:195539/1‑150 (MQ=255)
                                                                                                                                 gaggaAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCTCGACGAAGCGTTGTCGAACCTGGACCGGCTg                  >  1:911838/1‑150 (MQ=255)
                                                                                                                                 gaggaAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCTCGACGAAGCGTTGTCGAACCTGGACCGGCTg                  >  1:3484191/1‑150 (MQ=255)
                                                                                                                                   ggaAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCTCGACGAAGCGTTGTCGAACCTGGACCGGCTGtt                <  1:1997858/150‑1 (MQ=255)
                                                                                                                                      aaTGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCTCGACGAAGCGTTGTCGAACCTGGACCGGCTGTTGCa             <  1:2794124/150‑1 (MQ=255)
                                                                                                                                      aaTGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCTCGACGAAGCGTTGTCGAACCTGGACCGGCTGTTGCa             <  1:5047349/150‑1 (MQ=255)
                                                                                                                                      aaTGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCTCGACGAAGCGTTGTCGAACCTGGACCGGCTGTTGCa             <  1:2200516/150‑1 (MQ=255)
                                                                                                                                       aTGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCTCGACGAAGCGTTGTCGAACCTGGACCGGCTGTTGCAg            <  1:3233562/150‑1 (MQ=255)
                                                                                                                                         gggCATGCTCGGCGCGCAGGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCTCGACGAAGCGTTGTCGAACCTGGACCGGCTGTTGCAGCt          >  1:2495528/1‑150 (MQ=255)
                                                                                                                                          ggCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCTCGACGAAGCGTTGTCGAACCTGGACCGGCTGTTGCAGCTg         <  1:1537915/150‑1 (MQ=255)
                                                                                                                                          ggCATGCTCGGCGCGCAGGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCTCGACGAAGCGTTGTCGAACCTGGACCGGCTGTTGCAGCTg         >  2:4081987/1‑150 (MQ=255)
                                                                                                                                          ggCATGCTCGGCGCGCAGGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCTCGACGAAGCGTTGTCGAACCTGGACCGGCTGTTGCAGCTg         >  2:4082034/1‑150 (MQ=255)
                                                                                                                                           gcatGCTCGGCGCGCAGGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCTCGACGAAGCGTTGTCGAACCTGGACCGGCTGTTGCAGCTGc        >  1:2544166/4‑150 (MQ=255)
                                                                                                                                           gcatGCTCGGCGCGCAGGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCTCGACGAAGCGTTGTCGAACCTGGACCGGCTGTTGCAGCTGc        >  1:2544255/4‑150 (MQ=255)
                                                                                                                                              tGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCTCGACGAAGCGTTGTCGAACCTGGACCGGCTGTTGCAGCTGCAGa     <  2:4604432/150‑1 (MQ=255)
                                                                                                                                              tGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCTCGACGAAGCGTTGTCGAACCTGGACCGGCTGTTGCAGCTGCAGa     <  1:2565060/150‑1 (MQ=255)
                                                                                                                                              tGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCTCGACGAAGCGTTGTCGAACCTGGACCGGCTGTTGCAGCTGCAGa     >  1:2940622/1‑150 (MQ=255)
                                                                                                                                                 gcgcCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCTCGACGAAGCGTTGTCGAACCTGGACCGGCTGTTGCAGCTGCAGATcc  <  1:2102618/150‑1 (MQ=255)
                                                                                                                                                 gcgcCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCTCGACGAAGCGTTGTCGAACCTGGACCGGCTGTTGCAGCTGCAGATcc  <  2:4222215/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
CTGGTGTTCCAGGACGCCCCTGGCGCCTTCAACCCGCAGCGCAGCATCGGCTGGAGCATCGCCGAGCCGCTGCGCCATCTGACCGACATGGATGAAGCCGCCCGGCAGGCGCGCACCCTTGCACTGCTGGAGGAAATGGGCCTGCGCGCCGAGCATGCCCAACGCCTGCCGCAGCAGTTGAGTGGCGGCCAGTTGCAGCGTGCCAACATTGCCCGGGCGCTGGCCGCCTCGCCACAGCTGGTAGTGCTCGACGAAGCGTTGTCGAACCTGGACCGGCTGTTGCAGCTGCAGATCC  >  NC_002947/3785918‑3786212

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: