Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,956,098 A→C 9.5% 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,956,0980AC 9.5% 419.7 / 14.8 147intergenic (‑69/+126)PP_3486/PP_3487cytochrome c/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (72/56);  new base C (7/10);  total (79/68)
Fisher's exact test for biased strand distribution p-value = 3.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 164 total aligned reads displayed.

GACGCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGG  >  NC_002947/3955949‑3956246
                                                                                                                                                     |                                                                                                                                                    
gaCGCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCgca                                                                                                                                                      <  2:6402426/150‑1 (MQ=255)
 aCGCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCgcag                                                                                                                                                     <  1:6805057/150‑1 (MQ=255)
   gCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCgcagca                                                                                                                                                   <  2:5863105/150‑1 (MQ=255)
   gCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCgcagca                                                                                                                                                   <  2:2314029/150‑1 (MQ=255)
   gCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCgcagca                                                                                                                                                   <  1:3130607/150‑1 (MQ=255)
   gCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCgcagca                                                                                                                                                   <  1:3129901/150‑1 (MQ=255)
   gCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCGTGCGTCGGGACAGGGGTGAgcagca                                                                                                                                                   >  2:5093158/1‑150 (MQ=255)
       tCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCGTGCGGCGCGGAACGGCGGCGCAGCAcccc                                                                                                                                               >  2:2131180/1‑150 (MQ=255)
        cgcgGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGcccc                                                                                                                                              <  2:883945/150‑1 (MQ=255)
         gcgGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCg                                                                                                                                             <  2:5860241/150‑1 (MQ=255)
          cgGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCgg                                                                                                                                            <  1:3662505/150‑1 (MQ=255)
          cgGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCgg                                                                                                                                            <  1:4905892/150‑1 (MQ=255)
              acaGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCTATGCCACGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGCGCGGCCTTGCGTCGCGACAGGGCTGCGCAGCAGCCCCGGCAGt                                                                                                                                        <  2:4897742/148‑1 (MQ=255)
              aaaGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGt                                                                                                                                        <  2:4742013/150‑1 (MQ=255)
               aaGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCGTTGGGCGCGCAAGGGCTGGGCAGCCCCCCCCGCCGTg                                                                                                                                       >  1:4331541/1‑150 (MQ=255)
                aGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGGCCTGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGTTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTgg                                                                                                                                      <  1:439207/150‑1 (MQ=255)
                 ggCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCGTGGGGCGCGGAAGGGCTGCGCCCCAGCCCCCGCacgggg                                                                                                                                     >  2:5258087/1‑145 (MQ=255)
                  gCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCa                                                                                                                                    <  2:6015435/150‑1 (MQ=255)
                  gCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCa                                                                                                                                    <  2:6015288/150‑1 (MQ=255)
                   cGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCTAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAg                                                                                                                                   <  1:5471959/150‑1 (MQ=255)
                      ttGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTc                                                                                                                                <  2:3585712/150‑1 (MQ=255)
                         ccTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAg                                                                                                                             <  2:3264270/150‑1 (MQ=255)
                          cTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAgc                                                                                                                            <  2:6862863/150‑1 (MQ=255)
                          cTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAgc                                                                                                                            <  1:3344918/150‑1 (MQ=255)
                           tGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAgcg                                                                                                                           <  2:1860552/150‑1 (MQ=255)
                              gtGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAg                                                                                                                        <  1:4713762/149‑1 (MQ=255)
                               tGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCGGTGCGTCGCGTCAGGGCTGCGCCCCACCCCCCCCCCGGGCCGCTCAACCGCacc                                                                                                                       >  2:3709541/1‑148 (MQ=255)
                                 ttCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGaaa                                                                                                                     <  2:5538093/150‑1 (MQ=255)
                                 ttCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGaaa                                                                                                                     <  2:3101298/150‑1 (MQ=255)
                                    gTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAagag                                                                                                                  <  2:1809943/150‑1 (MQ=255)
                                    gTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAagag                                                                                                                  <  2:2733715/150‑1 (MQ=255)
                                    gTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAagag                                                                                                                  <  1:5662627/150‑1 (MQ=255)
                                      ggtcgttcAGGTCCTGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCAATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAggg                                                                                                                <  1:6501304/150‑1 (MQ=255)
                                      ggtcggtcAGGTCCAGGGCCTGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAggg                                                                                                                <  1:1671343/150‑1 (MQ=255)
                                      ggtcggtcAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAggg                                                                                                                <  2:727504/150‑1 (MQ=255)
                                      ggtcggtcAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAggg                                                                                                                <  2:6266867/150‑1 (MQ=255)
                                        tcggtcAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCGTTCGTCGCGGCAGGGCTGCGCCGCACCCCCCCCCGTGGCAGCGCCAGCGCCCGAAGAggggg                                                                                                              >  2:1517986/1‑149 (MQ=255)
                                                   cagggccagggcTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTgcgcgc                                                                                                   <  2:4569654/150‑1 (MQ=255)
                                                     gggCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTc                                                                                                 >  2:6418592‑M2/1‑149 (MQ=255)
                                                       gccagggcTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCGTTCGGCGCGGAAGGGGGGCGCAGCAGCCCCCGCCGGGGCCACGCAAGCCCAGAAAGAGGGGGGCGCTGCGCGCGccc                                                                                               >  2:4362034/1‑149 (MQ=255)
                                                        ccagggcTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCGTGCGGCGCGCAACGGCGGCGCCGCCCCCCCCCGAGTGGCCGCGCAAGCCCCCCAAGAGGGGggcgcggcgcccgccgg                                                                                              >  2:1577267/1‑143 (MQ=255)
                                                          agggcTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCCTGGGGCCCGAGAGGGCGGCGCCCCCCCCCCCCCCGAGGCACCCCAACCGCCCACAGAGGGGGGCGCGGCGCGCGcccgtg                                                                                            >  2:6351231/1‑146 (MQ=255)
                                                          aGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttc                                                                                            >  1:4428882‑M2/1‑144 (MQ=255)
                                                           gggCTGAGGTCGAAGTTGCCAGGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttcg                                                                                           >  2:2591039‑M2/1‑143 (MQ=255)
                                                            ggCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttcgc                                                                                          >  1:1861039‑M2/1‑142 (MQ=255)
                                                            ggCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttcgc                                                                                          >  1:4940895‑M2/1‑142 (MQ=255)
                                                            ggCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGACGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttcgc                                                                                          >  1:3508964‑M2/1‑142 (MQ=255)
                                                               tGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttcgcgac                                                                                       >  2:3246078‑M2/1‑139 (MQ=255)
                                                                    tCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttcgcgacgcaag                                                                                  >  1:3808996‑M2/1‑134 (MQ=255)
                                                                    tCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCACAAGCGCCGCAAGAGGGGAGCGCAGCGCGCAccggtcgccacgcaag                                                                                  <  2:3951517‑M2/150‑17 (MQ=255)
                                                                    tCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCATTGCGTCGCTAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttcgcgacgcaag                                                                                  >  1:3090891‑M2/1‑134 (MQ=255)
                                                                     cGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttcgcgacgcaagg                                                                                 >  2:6260379‑M2/1‑133 (MQ=255)
                                                                     cGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCGTGCGTCGCGGAAGGGCTGCGCCGCAGCCCCCCCCGTGGCCGCGCAAGCGCAGAAAGAGGGGAGCGCAGCGCGCGccggtagcgacgcaaag                                                                                 <  2:5074573‑M2/150‑18 (MQ=255)
                                                                      gAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttcgcgacgc                                                                                     <  2:3444957‑M2/145‑14 (MQ=255)
                                                                           tgccttgccaAGCACGCTGGTAAGCAATATCGAAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCGATGCAGCCCCGGCAGTGGCAGCTCAAGCGCATAAAGAGGGGAGCGCTGCGGGCTCCTTTCGCGACGCAAGGCCGCTc                                                                           <  2:4554117/150‑1 (MQ=37)
                                                                                    aaGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCGTGCGTCGCGGAAGGGCGGCGCAGCACCCCCCCCAGGGGCCGCACAACCGCCGCAAGAAGGGGGCGCGGCGCCCGCCGGGCGCGACCCAAGGCCGCGCCCCCcaaac                                                                  >  2:4121957/1‑149 (MQ=255)
                                                                                     aGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAgtg                                                                                                                           <  2:6751134/92‑3 (MQ=255)
                                                                                     aGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAgtg                                                                                                                           >  1:6751134/1‑90 (MQ=255)
                                                                                                          ccacgCCTTGTGGTCGGGGCCTTGGGTCGCGATCGGGCTGCGCCGCACCCCCGGCTGTGGCAGCTAAAGCGCACAATGCGGGTCGCGCGGCGCCCCCCTTTCGCCCCGCACGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTc                                            <  2:2137432/146‑1 (MQ=255)
                                                                                                             cgCCCTGGGGGAGCGGCCTTGCGCCGCCACCGGGCGGCGCCGCAGCCCCGGCCGTGGCAGCCCAAGCGCCCAACGAGGGGAGCGCGGCGCGCACCTTCCGCGCCGCACGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTcccc                                         <  2:2809130/149‑1 (MQ=255)
                                                                                                                         gCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATg                             >  2:287572/1‑150 (MQ=255)
                                                                                                                             ccTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTTCCCAGGCGGCCCATGATAc                         >  2:2650004/1‑150 (MQ=255)
                                                                                                                             ccTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATAc                         >  1:2690255/1‑150 (MQ=255)
                                                                                                                              cTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACt                        >  2:1290819/1‑150 (MQ=255)
                                                                                                                              cTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACt                        >  2:4248831/1‑150 (MQ=255)
                                                                                                                               ttGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTc                       >  2:3845711/1‑150 (MQ=255)
                                                                                                                                   gTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttcgcgacgcaaggccgctcccacagggcttgcgatatggcttacctgcgtgcttggcatggcaacttcgacctca                   <  1:1409035‑M2/150‑80 (MQ=255)
                                                                                                                                   gTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttcgcgacgcaaggccgctcccacagggcttgcgatatggcttacctgcgtgcttggcatggcaacttcgacctca                   <  1:3465031‑M2/150‑80 (MQ=255)
                                                                                                                                   gTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttcgcgacgcaaggccgctcccacagggcttgcgatatggcttacctgcgtgcttggcatggcaacttcgacctca                   <  1:3434716‑M2/150‑80 (MQ=255)
                                                                                                                                   gTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttcgcgacgcaaggccgctcccacagggcttgcgatatggcttacctgcgtgcttggcatggcaacttcgacctca                   <  1:22006‑M2/150‑80 (MQ=255)
                                                                                                                                   gTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCaccg                   >  1:871168/1‑150 (MQ=255)
                                                                                                                                    tCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCaccga                  >  2:1664517/1‑150 (MQ=255)
                                                                                                                                     cgcgAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttcgcgacgcaaggccgctcccacagggcttgcgatatggcttacctgcgtgcttggcatggcaacttcgacctcagc                 <  2:2575454‑M2/150‑82 (MQ=255)
                                                                                                                                      gcgAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCaccgacc                >  1:3873025/1‑150 (MQ=255)
                                                                                                                                      gcgAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCaccgacc                >  1:6750934/1‑150 (MQ=255)
                                                                                                                                        gAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGTGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCgg              >  1:5575399/1‑150 (MQ=255)
                                                                                                                                        gAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCgg              >  2:5182597/1‑150 (MQ=255)
                                                                                                                                        gAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCgg              >  2:4119568/1‑150 (MQ=255)
                                                                                                                                        gAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCgg              >  1:3559655/1‑150 (MQ=255)
                                                                                                                                        gAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAAGTCCCCAGGCGGCCCATGATACTCACCGACCgg              >  2:6436989/1‑150 (MQ=255)
                                                                                                                                         aaaGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGACCGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACCAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGAAACTCACCGACCGGc             >  2:440597/1‑150 (MQ=255)
                                                                                                                                         aaaGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGc             >  1:4848792/1‑150 (MQ=255)
                                                                                                                                           ggggcTGCTGCGCAGCCCTGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttcgcgacgcaaggccgctcccacagggcttgcgatatggcttacctgcgtgcttggc                                     <  1:1029533‑M2/123‑62 (MQ=255)
                                                                                                                                           ggggcTGCTGCGCAGCCCTGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttcgcgacgcaaggccgctcccacagggcttgcgatatggcttacctgcgtgcttggc                                     >  2:1029533‑M2/2‑63 (MQ=255)
                                                                                                                                           ggggCTGCTGCGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCCCCGACCGGccc           >  2:1934176/2‑149 (MQ=255)
                                                                                                                                            gggCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttcgcgacgcaaggccgctcccacagggcttgcgatatggcttacctgcgtgcttggcatggcaacttcgacctcagccctggcc          <  2:306291‑M2/150‑89 (MQ=255)
                                                                                                                                            gggCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttcgcgacgcaaggccgctcccacagggcttgcgatatggcttacctgcgtgcttggcatggcaacttcgacctcagccctggcc          <  1:2652451‑M2/150‑89 (MQ=255)
                                                                                                                                             ggCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGc         >  2:3777718/1‑150 (MQ=255)
                                                                                                                                             ggCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGc         >  1:6183276/1‑150 (MQ=255)
                                                                                                                                               cTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttcgcgacgcaaggccgctcccacagggcttgcgatatggcttacctgcgtgcttggcatggcaacttcgacctcagccctggccctg       >  1:1006896‑M2/1‑59 (MQ=255)
                                                                                                                                               cTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCgg       >  2:2376954/1‑150 (MQ=255)
                                                                                                                                                tGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGc      >  1:1247125/1‑150 (MQ=255)
                                                                                                                                                tGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGc      >  1:2362946/1‑150 (MQ=255)
                                                                                                                                                tGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGc      >  2:2289444/1‑150 (MQ=255)
                                                                                                                                                tGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGc      >  1:6626240/1‑150 (MQ=255)
                                                                                                                                                tGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGc      >  1:2362942/1‑150 (MQ=255)
                                                                                                                                                tGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGc      >  1:249929/1‑150 (MQ=255)
                                                                                                                                                   gcagcaGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTgg   <  2:2690255/150‑1 (MQ=255)
                                                                                                                                                   gcagcaGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTgg   >  1:1826201/1‑150 (MQ=255)
                                                                                                                                                    cagcagCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTggg  >  1:4978781/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GACGCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGG  >  NC_002947/3955949‑3956246

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: