Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,036,649 C→G 5.6% L332V (CTG→GTG)  dfp → bifunctional 4'‑phosphopantothenoylcysteine decarboxylase/phosphopantothenoylcysteine synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,036,6490CG 5.6% 537.0 / 16.1 162L332V (CTG→GTG) dfpbifunctional 4'‑phosphopantothenoylcysteine decarboxylase/phosphopantothenoylcysteine synthetase
Reads supporting (aligned to +/- strand):  ref base C (76/77);  new base G (3/6);  total (79/83)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

Only 100 of 167 total aligned reads displayed.

AAGTGGTTGCCACACAAAAACTGAAGAAAGATCCTACGACCGGCGACGGCATGCTGCTGCAGATGGTGCGTAATCCCGATATCCTTGCCACCATTGCTGGCCGCGCCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAGGACAAGAACCTCGACCTGATCGTCGCCAATGATGTGGCCAACCCCAGCATCGGCTTCAACAGCGAAGAAAACGCCTTGACCGTGATCGACCGCCAGCAGCACCAGACCCTCTTCGCGCAG  >  NC_002947/6036500‑6036798
                                                                                                                                                     |                                                                                                                                                     
aaGTGGTTGCCACACAAAAACTGAAGAAAGATCCTACGACCGGCGACGGCATGCTGCTGCAGATGGTGCGTAATCCCGATATCCTTGCCACCATTGCTGGCCGCGCCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCAcc                                                                                                                                                       <  1:4559716/150‑1 (MQ=255)
   tGGTTGCCACACAAAAACTGAAGAAAGATCCTACGACCGGCGACGGCATGCTGCTGCAGATGGTGCGTAATCCCGATATCCTTGCCACCATTGCTGGCCGCGCCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGc                                                                                                                                                    <  1:4655106/150‑1 (MQ=255)
   tGGGTGCCACACAAAAACTGAAGCACGATCCTACGACCGGCGACGGCATGCTGCTGCAGATGGTGCGTAATCCCGATATCCTTGCCACCATTGCTGGCCGCGCCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGc                                                                                                                                                    <  2:1728989/150‑1 (MQ=255)
    ggTTGCCACACAAAAACTGAAGAAAGATCCTACGACCGGCGACGGCATGCTGCTGCAGATGGTGCGTAATCCCGATATCCTTGCCACCATTGCTGGCCGCGCCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCt                                                                                                                                                   >  1:3822184/1‑150 (MQ=255)
    ggTTGCCACACAAAAACTGAAGAAAGATCCTACGACCGGCGACGGCATGCTGCTGCAGATGGTGCGTAATCCCGATATCCTTGCCACCATTGCTGGCCGCGCCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCt                                                                                                                                                   <  1:3793126/150‑1 (MQ=255)
     gTTGCCACACAAAAACTGAAGAAAGATCCTACGACCGGCGACGGCATGCTGCTGCAGATGGTGCGTAATCCCGATATCCTTGCCACCATTGCTGGCCGCGCCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTc                                                                                                                                                  <  2:1374135/150‑1 (MQ=255)
           acacaAAAACTGAAGAAAGATCCTACGACCGGCGACGGCATGCTGCTGCAGATGGTGCGTAATCCCGATATCCTTGCCACCATTGCTGGCCGCGCCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAATCGAGCAGGTGCTCGGtttc                                                                                                                                            <  2:3028790/150‑3 (MQ=255)
           acacaAAAACTGAAGAAAGATCCTACGACCGGCGACGGCATGCTGCTGCAGATGGTGCGTAATCCCGATATCCTTGCCACCATTGCTGGCCGCGCCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTAc                                                                                                                                            >  2:1488857/1‑150 (MQ=255)
            cacaAAAACTGAAGAAAGATCCTACGACCGGCGACGGCATGCTGCTGCAGATGGTGCGTAATCCCGATATCCTTGCCACCATTGCTGGCCGCGCCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACg                                                                                                                                           >  1:3223071/1‑150 (MQ=255)
               aaaaaCTGAAGAAAGATCCTACGACCGGCGACGGCATGCTGCTGCAGATGGTGCGTAATCCCGATATCCTTGCCACCATTGCTGGCCGCGCCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGTAATCGAGCAGGTGCTCGGTTTCggcg                                                                                                                                        <  1:784740/150‑4 (MQ=255)
                  aaCTGAAGAAAGATCCTACGACCGGCGACGGCATGCTGCTGCAGATGGTGCGTAATCCCGATATCCTTGCCACCATTGCTGGCCGCGCCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCAcgc                                                                                                                                     <  1:4987029/150‑1 (MQ=255)
                  aaCTGAAGAAAGATCCTACGACCGGCGACGGCATGCTGCTGCAGATGGTGCGTAATCCCGATATCCTTGCCACCATTGCTGGCCGCGCCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCAcgc                                                                                                                                     <  1:4986931/150‑1 (MQ=255)
                       aagaaagaTCCTACGACCGGCGACGGCATGCTGCTGCAGATGGTGCGTAATCCCGATATCCTTGCCACCATTGCTGGCCGCGCCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAg                                                                                                                                <  2:3822184/150‑1 (MQ=255)
                          aaagaTCCTACGACCGGCGACGGCATGCTGCTGCAGATGGTGCGTAATCCCGATATCCTTGCCACCATTGCTGGCCGCGCCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTc                                                                                                                             >  2:2407980/1‑150 (MQ=255)
                          aaagaTCCTACGACCGGCGACGGCATGCTGCTGCAGATGGTGCGTAATCCCGATATCCTTGCCACCATTGCTGGCCGCGCCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTc                                                                                                                             >  2:2295267/1‑150 (MQ=255)
                           aagaTCCTACGACCGGCGACGGCATGCTGCTGCAGATGGTGCGTAATCCCGATATCCTTGCCACCATTGCTGGCCGCGCCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCa                                                                                                                            >  1:4052925/1‑150 (MQ=255)
                             gaTCCTACGACCGGCGACGGCATGCTGCTGCAGATGGTGCGTAATCCCGATATCCTTGCCACCATTGCTGGCCGCGCCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAg                                                                                                                          <  1:4443967/150‑1 (MQ=255)
                             gaTCCTACGACCGGCGACGGCATGCTGCTGCAGATGGTGCGTAATCCCGATATCCTTGCCACCATTGCTGGCCGCGCCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAg                                                                                                                          <  1:3333796/150‑1 (MQ=255)
                             gaTCCTACGACCGGCGACGGCATGCTGCTGCAGATGGTGCGTAATCCCGATATCCTTGCCACCATTGCTGGCCGCGCCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAg                                                                                                                          <  2:2854316/150‑1 (MQ=255)
                             gaTCCTACGACCGGCGACGGCATGCTGCTGCAGATGGTGCGTAATCCCGATATCCTTGCCACCATTGCTGGCCGCGCCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAg                                                                                                                          <  1:4928108/150‑1 (MQ=255)
                                ccTACGACCGGCGACGGCATGCTGCTGCAGATGGTGCGTAATCCCGATATCCTTGCCACCATTGCTGGCCGCGCCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAGGAc                                                                                                                       >  1:161760/1‑150 (MQ=255)
                                ccTACGACCGGCGACGGCATGCTGCTGCAGATGGTGCGTAATCCCGATATCCTTGCCACCATTGCTGGCCGCGCCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAGGAc                                                                                                                       >  1:161651/1‑150 (MQ=255)
                                 cTACGACCGGCGACGGCATGCTGCTGCAGATGGTGCGTAATCCCGATATCCTTGCCACCATTGCTGGCCGCGCCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCAcgc                                                                                                                                     <  1:961110/135‑1 (MQ=255)
                                  tACGACCGGCGACGGCATGCTGCTGCAGATGGTGCGTAATCCCGATATCCTTGCCACCATTGCTGGCCGCGCCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAGGACaa                                                                                                                     >  2:4037400/1‑150 (MQ=255)
                                  tACGACCGGCGACGGCATGCTGCTGCAGATGGTGCGTAATCCCGATATCCTTGCCACCATTGCTGGCCGCGCCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAGGACaa                                                                                                                     >  2:3153704/1‑150 (MQ=255)
                                      aCCGGCGACGGCATGCTGCTGCAGATGGTGCGTAATCCCGATATCCTTGCCACCATTGCTGGCCGCGCCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAGGACAAGAAc                                                                                                                 >  1:239845/1‑150 (MQ=255)
                                         ggCGACGGCATGCTGCTGCAGATGGTGCGTAATCCCGATATCCTTGCCACCATTGCTGGCCGCGCCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAGGACAAGAACCTc                                                                                                              <  1:1008387/150‑1 (MQ=255)
                                           cGACGGCATGCTGCTGCAGATGGTGCGTAATCCCGATATCCTTGCCACCATTGCTGGCCGCGCCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAGGACAAGAACCTCGa                                                                                                            <  2:239845/150‑1 (MQ=255)
                                                gCATGCTGCTGCAGATGGTGCGTAATCCCGATATCCTTGCCACCATTGCTGGCCGCGCCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAGGACAAGAACCTCGACCTGa                                                                                                       <  2:4794524/150‑1 (MQ=255)
                                                gCATGCTGCTGCAGATGGTGCGTAATCCCGATATCCTTGCCACCATTGCTGGCCGCGCCGACCGCCCGTTCAGCGTCCGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAGGACAAGAACCTCGACCTGa                                                                                                       >  2:2312260/1‑150 (MQ=255)
                                                  aTGCTGCTGCAGATGGTGCGTAATCCCGATATCCTTGCCACCATTGCTGGCCGCGCCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAGGACAAGAACCTCGACCTGAtc                                                                                                     <  1:2312260/150‑1 (MQ=255)
                                                  aTGCCGCTGCAGATGGTGCGTAATCCCGATATCCTTGCCACCATTGCTGGCCGCGCCGACCGCCCGTTCAGCGTCGGCTTTGCAGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAGGACAAGAACCTCGACCTGAtc                                                                                                     <  2:581134/150‑1 (MQ=255)
                                                      tgctgcAGATGGTGCGTAATCCCGATATCCTTGCCACCATTGCTGGCCGCGCCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAATCGAGCAGGTGCTCGGTTTCgg                                                                                                                                          <  2:3268960/109‑2 (MQ=255)
                                                      tgctgcAGATGGTGCGTAATCCCGATATCCTTGCCACCATTGCTGGCCGCGCCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAGGACAAGAACCTCGACCTGAtcgtcg                                                                                                 >  2:4379728/1‑150 (MQ=255)
                                                      tgctgcAGATGGTGCGTAATCCCGATATCCTTGCCACCATTGCTGGCCGCGCCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAGGACAAGAACCTCGACCTGAtcgtcg                                                                                                 >  1:1276837/1‑150 (MQ=255)
                                                          gcAGATGGTGCGTAATCCCGATATCCTTGCCACCATTGCTGGCCGCGCCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAGGACAAGAACCTCGACCTGATCGTCGCCaa                                                                                             >  1:4460743/1‑150 (MQ=255)
                                                          gcAGATGGTGCGTAATCCCGATATCCTTGCCACCATTGCTGGCCGCGCCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAGGACAAGAACCTCGACCTGATCGTCGCCaa                                                                                             <  1:2023196/150‑1 (MQ=255)
                                                          gcAGATGGTGCGTAATCCCGATATCCTTGCCACCATTGCTGGCCGCGCCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAGGACAAGAACCTCGACCTGATCGTCGCCaa                                                                                             <  1:2022929/150‑1 (MQ=255)
                                                            aGATGGTGCGTAATCCCGATATCCTTGCCACCATTGCTGGCCGCGCCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAGGACAAGAACCTCGACCTGATCGTCGCCAatg                                                                                           >  1:4161952/1‑150 (MQ=255)
                                                              aTGGTGCGTAATCCCGATATCCTTGCCACCATTGCTGGCCGCGCCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAGGACAAGAACCTCGACCTGATCGTCGCCAatgat                                                                                         >  2:534696/1‑150 (MQ=255)
                                                                ggTGCGTAATCCCGATATCCTTGCCACCATTGCTGGCCGCGCCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAGGACAAGAACCTCGACCTGATCGTCGCCAATGAtgt                                                                                       <  2:2690872/150‑1 (MQ=255)
                                                                 gTGCGTAATCCCGATATCCTTGCCACCATTGCTGGCCGCGCCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAGGACAAGAACCTCGACCTGATCGTCGCCAATGAtgtg                                                                                      <  1:1453334/150‑1 (MQ=255)
                                                                 gTGCGTAATCCCGATATCCTTGCCACCATTGCTGGCCGCGCCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAGGACAAGAACCTCGACCTGATCGTCGCCAATGAtgtg                                                                                      <  1:27690/150‑1 (MQ=255)
                                                                   gCGTAATCCCGATATCCTTGCCACCATTGCTGGCCGCGCCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAGGACAAGAACCTCGACCTGATCGTCGCCAATGATGTGGc                                                                                    <  1:687346/150‑1 (MQ=255)
                                                                    cGTAATCCCGATATCCTTGCCACCATTGCTGGCCGCGCCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAGGACAAGAACCTCGACCTGATCGTCGCCAATGATGTGGcc                                                                                   <  1:1222184/150‑1 (MQ=255)
                                                                      tAATCCCGATATCCTTGCCACCATTGCTGGCCGCGCCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAGGACAAGAACCTCGACCTGATCGTCGCCAATGATGTGGCCaa                                                                                 <  2:5277054/150‑1 (MQ=255)
                                                                      tAATCCCGATATCCTTGCCACCATTGCTGGCCGCGCCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAGGACAAGAACCTCGACCTGATCGTCGCCAATGATGTGGCCaa                                                                                 <  2:2551001/150‑1 (MQ=255)
                                                                       aaTCCCGATATCCTTGCCACCATTGCTGGCCGCGCCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAGGACAAGAACCTCGACCTGATCGTCGCCAATGATGTGGCCAAc                                                                                >  2:2064651/1‑150 (MQ=255)
                                                                         tCCCGATATCCTTGCCACCATTGCTGGCCGCGCCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAGGACAAGAACCTCGACCTGATCGTCGCCAATGATGTGGCCAAccc                                                                              >  1:3574823/1‑150 (MQ=255)
                                                                          cccGATATCCTTGCCACCATTGCTGGCCGCGCCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAGGACAAGAACCTCGACCTGATCGTCGCCAATGATGTGGCCAAcccc                                                                             <  1:363417/150‑1 (MQ=255)
                                                                           ccGATATCCTTGCCACCATTGCTGGCCGCGCCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAGGACAAGAACCTCGACCTGATCGTCGCCAATGATGTGGCCAACCCCa                                                                            >  1:5380925/1‑150 (MQ=255)
                                                                           ccGATATCCTTGCCACCATTGCTGGCCGCGCCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAGGACAAGAACCTCGACCTGATCGTCGCCAATGATGTGGCCAACCCCa                                                                            >  2:1595411/1‑150 (MQ=255)
                                                                           ccGATATCCTTGCCACCATTGCTGGCCGCGCCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAGGACAAGAACCTCGACCTGATCGTCGCCAATGATGTGGCCAACCCCa                                                                            >  2:1595409/1‑150 (MQ=255)
                                                                            cGATATCCTTGCCACCATTGCTGGCCGCGCCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAGGACAAGAACCTCGACCTGATCGTCGCCAATGATGTGGCCAACCCCAg                                                                           >  2:2594759/1‑150 (MQ=255)
                                                                                atCCTTGCCACCATTGCTGGCCGCGCCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAGGACAAGAACCTCGACCTGATCGTCGCCAATGATGTGGCCAACCCCAGCATc                                                                       >  1:4962091/1‑150 (MQ=255)
                                                                                   cTTGCCACCATTGCTGGCCGCGCCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAGGACAAGAACCTCGACCTGATCGTCGCCAATGATGTGGCCAACCCCAGCATCGGc                                                                    <  2:1276837/150‑1 (MQ=255)
                                                                                    ttGCCACCATTGCTGGCCGCGCCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAGGACAAGAACCTCGACCTGATCGTCGCCAATGATGTGGCCAACCCCAGCATCGGCt                                                                   >  1:4984815/1‑150 (MQ=255)
                                                                                        caccaTTGCTGGCCGCGCCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAGGACAAGAACCTCGACCTGATCGTCGCCAATGATGTGGCCAACCCCAGCATCGGCTTCaa                                                               >  2:1070915/1‑150 (MQ=255)
                                                                                         accaTTGCTGGCCGCGCCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAATCGAGCAGGTGCTCGGTTTCggcg                                                                                                                                        <  2:66378/76‑4 (MQ=25)
                                                                                         accaTTGCTGGCCGCGCCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCAcgc                                                                                                                                     <  1:1743402/79‑1 (MQ=255)
                                                                                          ccaTTGCTGGCCGCGCCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAGGACAAGAACCTCGACCTGATCGTCGCCAATGATGTGGCCAACCCCAGCATCGGCTTCAACa                                                             <  1:1595409/150‑1 (MQ=255)
                                                                                           caTTGCTGGCCGCGCCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAGGACAAGAACCTCGACCTGATCGTCGCCAATGATGTGGCCAACCCCAGCATCGGCTTCAACAg                                                            <  2:4878752/150‑1 (MQ=255)
                                                                                            aTTGCTGGCCGCGCCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAGGACAAGAACCTCGACCTGATCGTCGCCAATGATGTGGCCAACCCCAGCATCGGCTTCAACAGc                                                           >  2:1765702/1‑150 (MQ=255)
                                                                                             ttGCTGGCCGCGCCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAGGACAAGAACCTCGACCTGATCGTCGCCAATGATGTGGCCAACCCCAGCATCGGCTTCAACAGCg                                                          >  1:2584756/1‑150 (MQ=255)
                                                                                               gCTGGCCGCGCCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAGGACAAGAACCTCGGCCTGATCGTCGCCAATGATGTGGCCAACCCCAGCATCGGCTTCAACAGCgaa                                                        <  2:2027318/150‑1 (MQ=255)
                                                                                               gCTGGCCGCGCCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAGGACAAGAACCTCGACCTGATCGTCGCCAATGATGTGGCCAACCCCAGCATCGGCTTCAACAGCgaa                                                        <  2:3886370/150‑1 (MQ=255)
                                                                                               gCTGGCCGCGCCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAGGACAAGAACCTCGACCTGATCGTCGCCAATGATGTGGCCAACCCCAGCATCGGCTTCAACAGCgaa                                                        <  1:2275157/150‑1 (MQ=255)
                                                                                                 tGGCCGCGCCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAGGACAAGAACCTCGACCTGATCGTCGCCAATGATGTGGCCAACCCCAGCATCGGCTTCAACAGCgaaga                                                      >  2:1783062/1‑150 (MQ=255)
                                                                                                 tGGCCGCGCCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAGGACAAGAACCTCGACCTGATCGTCGCCAATGATGTGGCCAACCCCAGCATCGGCTTCAACAGCgaaga                                                      >  2:501433/1‑150 (MQ=255)
                                                                                                 tGGCCGCGCCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAGGACAAGAACCTCGACCTGATCGTCGCCAATGATGTGGCCAACCCCAGCATCGGCTTCAACAGCgaa                                                        <  1:1969181/148‑1 (MQ=255)
                                                                                                 tGGCCGCGCCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAGGACAAGAACCTCGACCTGATCGTCGCCAATGATGTGGCCAACCCCAGCATCGGCTTCAACAGCgaa                                                        >  2:1969181/1‑148 (MQ=255)
                                                                                                       cgcCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGACACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAGGACAAGAACCTCGACCTGATCGTCGCCAATGATGTGGCCAACCCCAGCATCGGCTTCAACAGCGAAGAAAACGc                                                <  2:4695350/150‑1 (MQ=255)
                                                                                                       cgcCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAGGACAAGAACCTCGACCTGATCGTCGCCAACGATGTGGCCAACCCCAGCATCGGCTTCAACAGCGAAGAAAACGc                                                <  1:4853771/150‑1 (MQ=255)
                                                                                                        gcCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAGGACAAGAACCTCGACCTGATCGTCGCCAATGATGTGGCCAACCCCAGCATCGGCTTCAACAGCGAAGAAAACGcc                                               <  2:3389570/150‑1 (MQ=255)
                                                                                                        gcCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAGGACAAGAACCTCGACCTGATCGTCGCCAATGATGTGGCCAACCCCAGCATCGGCTTCAACAGCGAAGAAAACGcc                                               >  2:4097145/1‑150 (MQ=255)
                                                                                                          cGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAGGACAAGAACCTCGACCTGATCGTCGCCAATGATGTGGCCAACCCCAGCATCGGCTTCAACAGCGAAGAAAACGCCtt                                             <  1:2872013/150‑1 (MQ=255)
                                                                                                              cGCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAGGACAAGAACCTCGACCTGATCGTCGCCAATGATGTGGCCAACCCCAGCATCGGCTTCAACAGCGAAGAAAACGCCTTGAcc                                         <  2:1390828/150‑1 (MQ=255)
                                                                                                              cGCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAGGACAAGAACCTCGACCTGATCGTCGCCAATGATGTGGCCAACCCCAGCATCGGCTTCAACAGCGAAGAAAACGCCTTGAcc                                         <  2:3268184/150‑1 (MQ=255)
                                                                                                               gCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAGGACAAGAACCTCGACCTGATCGTCGCCAATGATGTGGCCAACCCCAGCATCGGCTTCAACAGCGAAGAAAACGCCTTGACCg                                        >  2:2784628/1‑150 (MQ=255)
                                                                                                               gCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAGGACAAGAACCTCGACCTGATCGTCGCCAATGATGTGGCCAACCCCAGCATCGGCTTCAACAGCGAAGAAAACGCCTTGACCg                                        >  2:1238837/1‑150 (MQ=255)
                                                                                                                cccGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAGGACAAGAACCTCGACCTGATCGTCGCCAATGATGTGGCCAACCCCAGCATCGGCTTCAACAGCGAAGAAAACGCCTTGACCGt                                       >  2:4100251/1‑150 (MQ=255)
                                                                                                                  cGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAGGACAAGAACCTCGACCTGATCGTCGCCAATGATGTGGCCAACCCCAGCATCGGCTTCAACAGCGAAGAAAACGCCTTGACCGTGa                                     <  1:501433/150‑1 (MQ=255)
                                                                                                                     tCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAGGACAAGAACCTCGACCTGATCGTCGCCAATGATGTGGCCAACCCCAGCATCGGCTTCAACAGCGAAGAAAACGCCTTGACCGTGATCg                                  <  2:833194/150‑1 (MQ=255)
                                                                                                                      cAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAGGACAAGAACCTCGACCTGATCGTCGCCAATGATGTGGCCAACCCCAGCATCGGCTTCAACAGCGAAGAAAACGCCTTGACCGTGATCGa                                 >  1:1487403/1‑150 (MQ=255)
                                                                                                                       aGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAGGACAAGAACCTCGACCTGATCGTCGCCAATGATGTGGCCAACCCCAGCa                                                                         <  1:165498/109‑1 (MQ=255)
                                                                                                                        gCGTCGGCTTTGCCGCCGAATCCGAGCACCTTCTCGATTACGCCACGCGCAAGCTCAAGGACAAGAACCTCGACCTGATCGTCGCCAATGATGTGGCCAACCCCAGCATCGGCTTCAACAGCGAAGAAAACGCCTTGACCGTGATCGAcc                               <  2:4588126/150‑1 (MQ=255)
                                                                                                                           tCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAGGACAAGAACCTCGACCTGATCGTCGCCAATGATGTGGCCAACCCCAGCATCGGCTTCAACAGCGAAGAAAACGCCTTGACCGTGATCGACCGcc                            >  1:831618/1‑150 (MQ=255)
                                                                                                                             ggCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAGGACAAGAACCTCGACCTGATCGTCGCCAATGATGTGGCCAACCCCAGCATCGGCTTCAACAGCGAAGAAAACGCCTTGACCGTGATCGACCGCcag                          <  2:4662587/150‑1 (MQ=255)
                                                                                                                               cTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAGGACAAGAACCTCGACCTGATCGTCGCCAATGATGTGGCCAACCCCAGCATCGGCTTCAACAGCGAAGAAAACGCCTTGACCGTGATCGACCGCcagca                        <  2:3128469/150‑1 (MQ=255)
                                                                                                                                tttGCCGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAGGACAAGAACCTCGACCTGATCGTCGCCAATGATGTGGCCAACCCCAGCATCGGCTTCAACAGCGAAGAAAACGCCTTGACCGTGATCGAc                                >  1:2087576/1‑141 (MQ=255)
                                                                                                                                     cgccgAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAGGACAAGAACCTCGACCTGATCGTCGCCAATGATGTGGCCAACCCCAGCATCGGCTTCAACAGCGAAGAAAACGCCTTGACCGTGATCGACCGCCAGCAGCACCa                  >  1:1123762/1‑150 (MQ=255)
                                                                                                                                     cgccgAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAGGACAAGAACCTCGACCTGATCGTCGCCAATGATGTGGCCAACCCCAGCATCGGCTTCAACAGCGAAGAAAACGCCTTGACCGTGATCGACCGCCAGCAGCACCa                  >  2:5203922/1‑150 (MQ=255)
                                                                                                                                     cgccgAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAGGACAAGAACCTCGACCTGATCGTCGCCAATGATGTGGCCAACCCCAGCATCGGCTTCAACAGCGAAGAAAACGCCTTGACCGTGATCGACCGCCAGCAGCACCa                  >  2:5228894/1‑150 (MQ=255)
                                                                                                                                      gccgAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAGGACAAGAACCTCGACCTGATCGTCGCCAATGATGTGGCCAACCCCAGCATCGGCTTCAACAGCGAAGAAAACGCCTTGACCGTGATCGACCGCCAGCAGCACCAg                 >  1:541312/1‑150 (MQ=255)
                                                                                                                                       gcgtaatcGAGCAGGTGCTCGGTTTCGGCACGCGCAAGCTCAAGGACAAGAACCTCGACCTGATCGTCGCCAATGATGTGGCCAACCCCAGCATCGGCTTCAACAGCGAAGAAAACGCCTTGACCGTGATCGACCGCCAGCAGCACCAGa                >  1:120940/8‑150 (MQ=255)
                                                                                                                                         gAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAGGACAAGAACCTCGACCTGATCGTCGCCAATGATGTGGCCAACCCCAGCATCGGCTTCAACAGCGAAGAAAACGCCTTGACCGTGATCGACCGCCAGCAGCACCAGAcc              <  1:1783062/150‑1 (MQ=255)
                                                                                                                                               gAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAGGACAAGAACCTCGACCTGATCGTCGCCAATGATGTGGCCAACCCCAGCATCGGCTTCAACAGCGAAGAAAACGCCTTGACCGTGATCGACCGCCAGCAGCACCAGACCCTCTTc        >  1:1021051/1‑150 (MQ=255)
                                                                                                                                                    ccTGCTCGATTACGCCACGCGCAAGCTCAAGGACAAGAACCTCGACCTGATCGTCGCCAATGATGTGGCCAACCCCAGCATCGGCTTCAACAGCGAAGAAAACGCCTTGACCGTGATCGACCGCCAGCAGCACCAGACCCTCTTCGCGCa   <  1:2594759/150‑1 (MQ=255)
                                                                                                                                                    ccTGCTCGATTACGCCACGCGCAAGCTCAAGGACAAGAACCTCGACCTGATCGTCGCCAATGATGTGGCCAACCCCAGCATCGGCTTCAACAGCGAAGAAAACGCCTTGACCGTGATCGACCGCCAGCAGCACCAGACCCTCTTCGCGCa   >  1:3443432/1‑150 (MQ=255)
                                                                                                                                                     cTGCTCGATTACGCCACGCGCAAGCTCAAGGACAAGAACCTCGACCTGATCGTCGCCAATGATGTGGCCAACCCCAGCATCGGCGTCAACAGCGAAGAAAACGCCTTGACCGTGATCGACCGCCAGCAGCACCAGACCCTCTTCGCGCAg  >  2:1898604/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AAGTGGTTGCCACACAAAAACTGAAGAAAGATCCTACGACCGGCGACGGCATGCTGCTGCAGATGGTGCGTAATCCCGATATCCTTGCCACCATTGCTGGCCGCGCCGACCGCCCGTTCAGCGTCGGCTTTGCCGCCGAAACCGAGCACCTGCTCGATTACGCCACGCGCAAGCTCAAGGACAAGAACCTCGACCTGATCGTCGCCAATGATGTGGCCAACCCCAGCATCGGCTTCAACAGCGAAGAAAACGCCTTGACCGTGATCGACCGCCAGCAGCACCAGACCCTCTTCGCGCAG  >  NC_002947/6036500‑6036798

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: