Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,953,071 A→C 14.4% intergenic (‑69/+126) PP_3486 ← / ← PP_3487 cytochrome c/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,953,0710AC14.4% 214.6 / 14.2 87intergenic (‑69/+126)PP_3486/PP_3487cytochrome c/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (30/43);  new base C (7/7);  total (37/50)
Fisher's exact test for biased strand distribution p-value = 5.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACGCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGT  >  NC_002947/3952922‑3953220
                                                                                                                                                     |                                                                                                                                                     
gaCGCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCgca                                                                                                                                                       <  1:257867/150‑1 (MQ=255)
  cGCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCgcagc                                                                                                                                                     <  2:1392760/150‑1 (MQ=255)
   gCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCgcagca                                                                                                                                                    <  1:1356758/150‑1 (MQ=255)
   gCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCgcagca                                                                                                                                                    <  1:1356785/150‑1 (MQ=255)
       tCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCTGCGCAGccc                                                                                                                                                <  2:2460393/150‑1 (MQ=255)
        cgcgGTAAAGGCGTTTGGCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCGAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGTAAGGGCTGCGCAGCAGcccc                                                                                                                                               <  2:141579/150‑1 (MQ=255)
        cgcgGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCTGCGCAGccct                                                                                                                                               <  1:2756023/150‑2 (MQ=255)
        cgcgGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCTGCGCAGccct                                                                                                                                               <  2:1217189/150‑2 (MQ=255)
         gcgGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCTGCGCAGCAGCCCCg                                                                                                                                              <  1:3348375/150‑1 (MQ=255)
          cgGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCTGCGCAGcccttt                                                                                                                                             <  2:2233236/150‑4 (MQ=255)
           ggTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCTGCGCAGccctttc                                                                                                                                            <  1:1448810/150‑5 (MQ=255)
           ggTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGc                                                                                                                                            <  2:402884/150‑1 (MQ=255)
            gTAACGGCGTGTGCCTGCGTGGTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCAGGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCa                                                                                                                                           <  2:2334037/150‑1 (MQ=39)
            gTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCa                                                                                                                                           <  2:1772804/150‑1 (MQ=255)
            gTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGCCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCa                                                                                                                                           <  1:750997/150‑1 (MQ=255)
                 ggCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTACGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGc                                                                                                                                      <  1:3497325/150‑1 (MQ=255)
                   cGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAg                                                                                                                                    <  2:1781469/150‑1 (MQ=255)
                    gTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGATGCGCAGCAGCCCCGGCAGTGGCAGc                                                                                                                                   <  2:62425/150‑1 (MQ=255)
                      ttGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTc                                                                                                                                 <  2:4485787/150‑1 (MQ=255)
                       tGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCTGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCa                                                                                                                                <  1:3124746/150‑1 (MQ=255)
                             cTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCa                                                                                                                          <  2:3662562/150‑1 (MQ=255)
                             cTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCa                                                                                                                          <  2:2780023/150‑1 (MQ=255)
                             cTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCa                                                                                                                          <  2:524062/150‑1 (MQ=255)
                                   cGTGGTCGGTCAGGTCCAGGGCCTGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCGAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAaga                                                                                                                    <  1:327093/150‑1 (MQ=255)
                                    gTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAagag                                                                                                                   <  2:2910974/150‑1 (MQ=255)
                                    gTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAagag                                                                                                                   <  1:276516/150‑1 (MQ=255)
                                     tGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCAGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCACTGGCAGCTCAAGCGCAGAAAGAgg                                                                                                                  <  2:2910609/150‑1 (MQ=255)
                                      ggtcggtcAGGTCCTGGGGCAGGGCTGAGGTCGAAGTTGCAATGCCAAGCACGCAGGTAAGCCATATCGATAGCCCTGTGGGTGCGGCCTTGCGGCGGGAAAGGGCTGAGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAggg                                                                                                                 <  2:3008868/150‑1 (MQ=255)
                                        tcggtcAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGa                                                                                                               <  2:2063646/150‑1 (MQ=255)
                                        tcggtcAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGa                                                                                                               <  2:2063513/150‑1 (MQ=255)
                                        tcggtcAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCGTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCCGCCGTGGCAGCACAAGCGCAGACAGAGGGGa                                                                                                               >  2:2610878/1‑150 (MQ=255)
                                        tcggtcAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCGTGCGGCGCGCAAGGGCGGGGCCGCAGCCCCCCCCCGGGCAGCGCAAGCGCACCAACAggggg                                                                                                               >  2:4317860/1‑149 (MQ=255)
                                            tcAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGGCGTGGGGGGCGGGAGGGCTGCGCCGGCGCCCCCCCCCGGGCAGGGCGCGCGCCGCAAGCGGGGAgcgc                                                                                                           >  1:4476218/1‑150 (MQ=255)
                                            tcAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCGTGCGGCGCGGAAGGGCTGCGCCGCAGCCCCCGCCCTGGCAGCTCAAGCGCCGCAAGAGGGgcgcgc                                                                                                           >  2:2346482/1‑150 (MQ=255)
                                               ggTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTgc                                                                                                        <  2:2521326/150‑1 (MQ=255)
                                               ggTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTgc                                                                                                        <  2:2855114/150‑1 (MQ=255)
                                               ggTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTgc                                                                                                        <  2:2521720/150‑1 (MQ=255)
                                               ggTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTgc                                                                                                        <  1:153355/150‑1 (MQ=255)
                                                gTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTgcg                                                                                                       <  1:4479520/150‑1 (MQ=255)
                                                 tCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTgcgc                                                                                                      <  1:1031822/150‑1 (MQ=255)
                                                 tCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTgcgc                                                                                                      >  2:3370787/1‑150 (MQ=255)
                                                   cagggccagggcTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTgcgcgc                                                                                                    <  2:2146029/150‑1 (MQ=255)
                                                   cagggccagggcTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTgcgcgc                                                                                                    <  1:2545240/150‑1 (MQ=255)
                                                   cagggccagggcTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCGTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCCGCCCTGGCAGCACAAGCGCAGCAACAGGGGAGCGCGgcgcgc                                                                                                    >  2:2954320/1‑150 (MQ=38)
                                                   cagggccagggcTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCGTGCGGCGCGCAACGGCGGCGCCGCCCCCCCGGCCGTGGCAGCGCCAGCGCCGACACAGGggcgcgcggcgcgc                                                                                                    >  2:2954619/1‑150 (MQ=255)
                                                    agggccagggcTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCt                                                                                                   <  1:4120449/150‑1 (MQ=255)
                                                     gggCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTc                                                                                                  >  2:2273400‑M2/1‑149 (MQ=255)
                                                         cAGGGCTGAGGTCGAAGTTGCCATGCCACGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGAGGAGCGCGCAGCGCTcccct                                                                                              >  2:186996‑M2/1‑145 (MQ=255)
                                                          aGGGCTGTGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttc                                                                                             >  1:463675‑M2/1‑144 (MQ=255)
                                                           gggCTGCGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttcg                                                                                            >  2:1895481‑M2/1‑143 (MQ=255)
                                                           gggCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCTTATCGCATGCCCTGTGGGAGCGACCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttcg                                                                                            >  1:3664265‑M2/1‑143 (MQ=255)
                                                           gggCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGACCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttcg                                                                                            >  1:3664289‑M2/1‑143 (MQ=255)
                                                                   gTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttcgcgacggaa                                                                                    <  2:3812200‑M2/150‑16 (MQ=255)
                                                                     cGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGGAACGGCTGCGCAGCAGCCCCGGCAGTGGCAGCGCAAGCGCAGAAAGGGGGGAGGGCTGCGCGCTccgtgcgcgacgccagg                                                                                  <  2:1796034‑M2/150‑18 (MQ=255)
                                                                      gAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcatttcgcgacgcaaggc                                                                                 >  2:4510741‑M2/1‑132 (MQ=255)
                                                                         gTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttcgcgacgcaaggccgc                                                                              >  1:3370787‑M2/1‑129 (MQ=255)
                                                                              catgccaAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCGTGCGTCGCGCAAGGGCGGCGCCGCCCCCCCCCCCCGGGGAGCCCCACCCCCCAAAGAGGGGGGCGCGGCCCGCGGCGGTCGCGGCGCCACGCCGCGcccc                                                                         >  2:81862/1‑149 (MQ=255)
                                                                              catgccaAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCGTGCGGCGCGCAAGGGCTGCGCACCAGCCCCCGCCGTGGCAGCGCAAGCGCCCAAAGAGGGGAGCGCAGCGCGCACCATTCGCGACGCCAGGCCCCAcccc                                                                         >  2:3767765/1‑149 (MQ=255)
                                                                              catgccaAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCGTGCGGCGCGCAAGGGCGGCGCCGCAGCCCCCGCCGTGGCCGCGCCAGCGCCCCACGAGGGGGGCGCTGCGCGCGCCATTCGCGGCGCCAGGCCGCGcccc                                                                         >  2:3675242/1‑149 (MQ=255)
                                                                                    aaGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCGTGCCTCGGGCACAGGCGGGGCACCAGCCCCGCCAGAGGCAGCGCCAGCCCAGAAAGAAGGGGGCGCGGCGCGCGACGTGCCCGACCCAACGCCCCGCCCACacaaa                                                                   >  2:2971204/1‑149 (MQ=255)
                                                                                                   ccATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGCAAGGGCTGCGCCGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGCGCGCGGCGCGCTCCGTTTGCGACGCAAGGCCGCTCCCCCCAAAGCGGATGCCCCAGAgg                                                    >  2:2532893/1‑150 (MQ=38)
                                                                                                              gCCTTGTGGGTGCGGCGTTGCGTGGGTAATGGGGTGCGCCGCACCCCGGGCCGGGGCCGCCAACGCGACCAACGTGGGGTGCGGTGGGCGCGCCTTGCGCGACGCATGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCa                                         <  2:1372338/150‑1 (MQ=255)
                                                                                                                             ccTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCACGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATAc                          >  2:3489797/1‑150 (MQ=255)
                                                                                                                              cTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCGACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACt                         >  2:2090176/1‑150 (MQ=255)
                                                                                                                              cTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCTCAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACt                         >  2:2767596/1‑150 (MQ=255)
                                                                                                                              cTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACt                         >  2:2767506/1‑150 (MQ=255)
                                                                                                                              cTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACt                         >  2:4358962/1‑150 (MQ=255)
                                                                                                                               ttGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTc                        >  1:4022526/1‑150 (MQ=255)
                                                                                                                                 gCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttcgcgacgcaaggccgctcccacagggcttgcgatatggcttacctgcgtgcttggcatggcaacttcgacct                      <  2:1901634‑M2/150‑78 (MQ=255)
                                                                                                                                   gTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttcgcgacgcaaggccgctcccacagggcttgcgatatggcttacctgcgtgcttggcatggcaacttcgacctca                    <  1:3683327‑M2/150‑80 (MQ=255)
                                                                                                                                      gcgAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttcgcgacgcaaggccgctcccacagggcttgcgatatggcttacctgcgtgcttggcatggcaacttcgacctcagcc                 <  1:1718328‑M2/150‑83 (MQ=255)
                                                                                                                                      gcgAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttcgcgacgcaaggccgctcccacagggcttgcgatatggcttacctgcgtgcttggcatggcaacttcgacctcagcc                 <  2:2864508‑M2/150‑83 (MQ=255)
                                                                                                                                      gcgAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCaccgacc                 >  1:3401010/1‑150 (MQ=255)
                                                                                                                                       cgAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttcgcgacgcaaggccgctcccacagggcttgcgatatggcttacctgcgtgcttggcatggcaacttcgacctcagccc                <  1:419868‑M2/150‑84 (MQ=255)
                                                                                                                                       cgAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCaccgaccg                >  1:3428247/1‑150 (MQ=255)
                                                                                                                                        gAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttcgcgacgcaaggccgctcccacagggcttgcgatatggcttacctgcgtgcttggcatggcaacttcgacctcagccct               <  2:2610221‑M2/150‑85 (MQ=255)
                                                                                                                                        gAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCgg               >  1:2815796/1‑150 (MQ=255)
                                                                                                                                        gAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCgg               >  2:2645361/1‑150 (MQ=255)
                                                                                                                                         aaaGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGc              >  2:207822/1‑150 (MQ=255)
                                                                                                                                          aaGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttcgcgacgcaaggccg                                                                               >  2:3987757‑M2/1‑64 (MQ=255)
                                                                                                                                          aaGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttcgcgacgcaaggccg                                                                               <  1:3987757‑M2/84‑21 (MQ=255)
                                                                                                                                           gggggggcgcCGCCGCCGGGGCTGGGGCACCACACGCGCAGACTGGGGGGTGGGCGGCGGGCGCTTTGGGCGCCGCAAGCCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCa            <  2:3501110/144‑1 (MQ=255)
                                                                                                                                            gggCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttcgcgacgcaaggccgctcccacagggcttgcgatatggcttacctgcgtgcttggcatggcaacttcgacctcagccctggcc           <  1:953627‑M2/150‑89 (MQ=255)
                                                                                                                                            gggCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttcgcgacgcaaggccgctcccacagggcttgcgatatggcttacctgcgtgcttggcatggcaacttcgacctcagccctggcc           <  1:1533431‑M2/150‑89 (MQ=255)
                                                                                                                                            gggCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttcgcgacgcaaggccgctcccacagggcttgcgatatggcttacctgcgtgcttggcatggcaacttcgacctcagccctggcc           <  1:953599‑M2/150‑89 (MQ=255)
                                                                                                                                            gggCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAg           >  1:1159058/1‑150 (MQ=255)
                                                                                                                                             ggCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttcgcgacgcaaggccgctcccacagggcttgcgatatggcttacctgcgtgcttggcatggcaacttcgacctcagccctggccc          <  1:3340846‑M2/150‑90 (MQ=255)
                                                                                                                                              gCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCg         >  1:674677/1‑150 (MQ=255)
                                                                                                                                               cTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCgg        >  1:353793/1‑150 (MQ=255)
                                                                                                                                               cTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCgg        >  1:2611732/1‑150 (MQ=255)
                                                                                                                                                tGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGc       >  1:1028402/1‑150 (MQ=255)
                                                                                                                                                     agcagcCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGt  >  2:2884954/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GACGCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGT  >  NC_002947/3952922‑3953220

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