Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,953,071 A→C 6.9% 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,0710AC 6.9% 206.5 / 3.4 76intergenic (‑69/+126)PP_3486/PP_3487cytochrome c/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (38/31);  new base C (4/2);  total (42/34)
Fisher's exact test for biased strand distribution p-value = 6.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGT  >  NC_002947/3952924‑3953220
                                                                                                                                                   |                                                                                                                                                     
cGCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCgcagc                                                                                                                                                     <  2:2859847/150‑1 (MQ=255)
   ccTCGCGGTAAAGGCGTTTGCCTGCGTGTTCGTGGTCGGGCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCagcagc                                                                                                                                                  <  2:4504785/150‑1 (MQ=255)
     tCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGccc                                                                                                                                                <  1:1797253/150‑1 (MQ=255)
        cgGTAAAGGCGTTTGCCTGCTTGTTCGTTGTCTGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCTGCGCAGcccttt                                                                                                                                             <  2:1747324/150‑4 (MQ=255)
           tAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCGTGCGTCGCGCACGGGCTGCGCAGCACCCCCCGccc                                                                                                                                          >  2:4955744/1‑148 (MQ=255)
            aaaGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGt                                                                                                                                         <  2:2738323/150‑1 (MQ=255)
             aaGGCGTTTGCCTGCTTGTTCGTGCTCGGTCAGGACCAGGGCCAGGGCTGAGGTCGAATTTGCCATCCCAAGCACGCAGGAAACCCATATCCCAAGCCCTGTGGCAGCGGCCTTGCGTCGCCAAAGGGCTGCGCACCAGCCCCGGCAGTg                                                                                                                                        <  2:2987919/150‑1 (MQ=255)
                 cGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAg                                                                                                                                    <  2:1128009/150‑1 (MQ=255)
                  gTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGc                                                                                                                                   <  2:1268018/150‑1 (MQ=255)
                           cTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCTGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCa                                                                                                                          <  1:2988528/150‑1 (MQ=255)
                           cTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCa                                                                                                                          <  2:3139106/150‑1 (MQ=255)
                           cTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCa                                                                                                                          <  1:3756254/150‑1 (MQ=255)
                           cTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCa                                                                                                                          <  1:1977693/150‑1 (MQ=255)
                                 cGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGGAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAaga                                                                                                                    <  1:3216782/150‑1 (MQ=255)
                                 cGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAaga                                                                                                                    <  1:804506/150‑1 (MQ=255)
                                  gTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCTAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAagag                                                                                                                   <  1:861205/150‑1 (MQ=255)
                                  gTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGCGCGGCCTTGCGTCGCGAATGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAagag                                                                                                                   <  1:1466756/150‑1 (MQ=255)
                                   tGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAgg                                                                                                                  <  1:1645144/150‑1 (MQ=255)
                                    ggtcggtcAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAggg                                                                                                                 <  2:3108034/150‑1 (MQ=255)
                                     gtcgttcTGGTCATGGGCCTGGGCTGAGGTCGATGTTGCCAGGCCAAGCCCGCAGTTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTGGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAgggg                                                                                                                <  1:3006529/150‑1 (MQ=38)
                                          tcAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCGTTTGGTGGGGTATTGGGTTGGGATGAGGCCCGGCAGTTGGAGCTCAATCGCAGAAAGAtgggttggg                                                                                                           >  2:365579/1‑145 (MQ=255)
                                                ccagggccagggCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGGCGTGCGGCGCGGCAGGGCGGCGCCCCAGCCCCGGCCGGGGCCGCGCCATCGGAGCAAGAGGGGGGCgcggcggg                                                                                                     >  2:3903072/1‑148 (MQ=255)
                                                 cagggccagggcTGAGGTCTAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTgcgcgc                                                                                                    <  1:1524010/150‑1 (MQ=255)
                                                   gggCCAGGGCTGAGGTCGATGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGTGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTc                                                                                                  >  1:5113638‑M2/1‑149 (MQ=255)
                                                   gggCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTc                                                                                                  >  2:2434751‑M2/1‑149 (MQ=255)
                                                   gggCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTc                                                                                                  >  1:1678506‑M2/1‑149 (MQ=255)
                                                    ggccagggcTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTgcgccagtt                                                                                                 <  2:37112/150‑6 (MQ=255)
                                                    ggCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCGCAGCGCTcc                                                                                                 <  2:2437926‑M2/150‑3 (MQ=255)
                                                     gCCAGGGCTGAGGTCGAAGGTGCCAGGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCGCAGCGCTccc                                                                                                >  2:4781992‑M2/1‑147 (MQ=255)
                                                      gcaGGGCTGTGGTCGAAGTTTCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCGTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGATAGAGGGTAGCGCTCAGCGCTcccc                                                                                               >  2:2580493‑M2/2‑146 (MQ=255)
                                                       ccgggctGCGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCACGCCCGGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTccttt                                                                                              >  1:1923654‑M2/3‑145 (MQ=255)
                                                       cAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTccttt                                                                                              >  2:1693424‑M2/1‑145 (MQ=255)
                                                          ggCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttcgc                                                                                           >  2:2569580‑M2/1‑142 (MQ=255)
                                                             tGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttcgcgac                                                                                        >  2:4021611‑M2/1‑139 (MQ=255)
                                                             tGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttcgcgac                                                                                        >  1:4572923‑M2/1‑139 (MQ=255)
                                                                 gTCGAAGTTGCCATGCCACGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttcgcgacgcaa                                                                                    >  1:2215886‑M2/1‑135 (MQ=255)
                                                                 gTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGGCCTGGGGCGCGGAAGGGGGGGGCAGCCGCCCCCGCCGGGGCAGCCGCAACGCCGCCAAAGGGGTGGGCTGCGCGCGCCGGTAGCGCCGCaa                                                                                    >  2:2580834/1‑150 (MQ=255)
                                                                 gTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttcgcgacgcaa                                                                                    >  2:3243177‑M2/1‑135 (MQ=255)
                                                                    gAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCGTGCGTCGCGAAAGGGCTGCGCACCACCCCCCGCCGTGGCCCCGCAAGCCCCGAAACAGGGGAGCGCGGCGCGCACCAGTCGCGgcgccaagc                                                                                 >  2:3730502/1‑145 (MQ=255)
                                                                          gccatgccaAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCGTGCGTCGCGCCAGGGCGGCGCCGCCGCCCCCCCCCGGGCCGCGCCAGCGCCGCAAGCGGGGGGCGCGGCGCGCGCCGGTCGCGACGCAAGGCCGCGcc                                                                           >  2:5022734/1‑150 (MQ=255)
                                                                               gccaAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAggg                                                                                                                 >  1:1813749/1‑107 (MQ=255)
                                                                               gccaAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAggg                                                                                                                 <  2:1813749/107‑1 (MQ=255)
                                                                                                 ccATATCGCAAGCCCTGTGGGAGCGGCCGTTCGGCGCGAAAGGGCGGCGCACCAGCCCCCGCAGTGGCAGCGCAAGCGCCGAACGTGGGGAGCGCGGCGCGCTCCGGTCGCGACGCAAGGCCGCGCCCCCCAAAGCGGATGCCCCCGAgg                                                    >  2:4579868/1‑150 (MQ=255)
                                                                                                   atatCGCAAGCCCTGTGGGAGCGGCCATGCGTAGCGAAAGGGCTGCGCAGCAGCCCCCGCAGTGGCAGCTGAAGCCCAGAAAGAGGGGAGCGCGGCGCGCTCCGTTCGCGACCCAAGGCCGCGCCCCCAAAAGCGCAAGACACAGAGGcc                                                  >  2:2937874/1‑150 (MQ=255)
                                                                                                       cGCACGCCCTTTGGTAGCGGCGTTGCTTCGCGAATGGGCTGCGCTGCAGCCCCGGCTGGGGCAGCTCAAGCGCTGAATGAGGGGAGGGCTGGGCGCTCGTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAAtt                                              <  2:5045208/150‑1 (MQ=37)
                                                                                                                          gCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGACAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATa                           >  2:2948611/1‑150 (MQ=255)
                                                                                                                           ccTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAACAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGAtcc                          >  2:2809166/1‑148 (MQ=255)
                                                                                                                            cTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAg                                                                                                                              <  2:2428418/49‑1 (MQ=255)
                                                                                                                            cTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAg                                                                                                                              >  1:2428418/1‑49 (MQ=255)
                                                                                                                             ttGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTc                        >  1:2026237/1‑150 (MQ=255)
                                                                                                                             ttGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTc                        >  1:2026203/1‑150 (MQ=255)
                                                                                                                             ttGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGAGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTc                        >  1:2175992/1‑150 (MQ=255)
                                                                                                                                  tCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttcgcgacgcaaggccgctcccacagggcttgcgatatggcttacctgcgtgcttggcatggcaacttcgacctcag                   <  2:2759633‑M2/150‑81 (MQ=255)
                                                                                                                                    gcgAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttcgcgacgcaaggccgctcccacagggcttgcgatatggcttacctgcgtgcttggcatggcaacttcgacctcagcc                 <  2:4688281‑M2/150‑83 (MQ=255)
                                                                                                                                    gcgAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttcgcgacgcaaggccgctcccacagggcttgcgatatggcttacctgcgtgcttggcatggcaacttcgacctcagcc                 <  1:1138639‑M2/150‑83 (MQ=255)
                                                                                                                                    gcgAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCaccgacc                 >  2:1890718/1‑150 (MQ=255)
                                                                                                                                    gcgAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGAGACTCaccgacc                 >  2:4121106/1‑150 (MQ=255)
                                                                                                                                    gcgAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGATGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCaccgacc                 >  1:2372533/1‑150 (MQ=255)
                                                                                                                                      gAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttcgcgacgcaaggccgctcccacagggcttgcgatatggcttacctgcgtgcttggcatggcaacttcgacctcagccct               <  2:630313‑M2/150‑85 (MQ=255)
                                                                                                                                      gAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttcgcgacgcaaggccgctcccacagggcttgcgatatggcttacctgcgtgcttggcatggcaacttcgacctcagccct               <  2:2566652‑M2/150‑85 (MQ=255)
                                                                                                                                      gAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttcgcgacgcaaggccgctcccacagggcttgcgatatggcttacctgcgtgcttggcatggcaacttcgacctcagccct               <  2:2566555‑M2/150‑85 (MQ=255)
                                                                                                                                      gAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCgg               >  2:596617/1‑150 (MQ=255)
                                                                                                                                      gAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCgg               >  1:3828732/1‑150 (MQ=255)
                                                                                                                                        cggggCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGcc             >  1:3880205/3‑150 (MQ=255)
                                                                                                                                         ggggCTGCTGCGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCa            >  1:119398/2‑150 (MQ=255)
                                                                                                                                         ggggCTGCTGCGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATAATCACCGACCGGCCa            >  1:119429/2‑150 (MQ=255)
                                                                                                                                          gggCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttcgcgacgcaaggccgctcccacagggcttgcgatatggcttacctgcgtgcttggcatggcaacttcgacctcagccctggcc           <  1:3230794‑M2/150‑89 (MQ=255)
                                                                                                                                          gggCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttcgcgacgcaaggccgctcccacagggcttgcgatatggcttacctgcgtgcttggcatggcaacttcgacctcagccctggcc           <  1:3854810‑M2/150‑89 (MQ=255)
                                                                                                                                          gggCGGCGCCGCCGCCCGGGCCGTGGCTGCTCAAGCGCCGAATGTGGGGCGCGCGGCGCGCTCCTTTCGCGCCGCACGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAg           <  2:3332565/150‑1 (MQ=255)
                                                                                                                                           ggCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttcgcgacgcaaggccgctcccacagggcttgcgatatggcttacctgcgtgcttggcatggcaacttcgacctcagccctggccc          <  1:4959623‑M2/150‑90 (MQ=255)
                                                                                                                                           ggCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGc          >  2:2094271/1‑150 (MQ=255)
                                                                                                                                              tGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGTCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGc       >  1:727572/1‑150 (MQ=255)
                                                                                                                                              tGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGc       >  1:3633753/1‑150 (MQ=255)
                                                                                                                                              tGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGc       >  1:1947672/1‑150 (MQ=255)
                                                                                                                                              tGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGc       >  2:3371908/1‑150 (MQ=255)
                                                                                                                                               gcgcAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCt      >  1:2273957/1‑150 (MQ=255)
                                                                                                                                                cgcAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCcaca                                                        >  2:1881069/1‑99 (MQ=255)
                                                                                                                                                cgcAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCcaca                                                        <  1:1881069/99‑1 (MQ=255)
                                                                                                                                                 gcagcaGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTgg    >  2:1646047/1‑150 (MQ=255)
                                                                                                                                                  cagcagCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTggg   >  2:1734085/1‑150 (MQ=255)
                                                                                                                                                   agcagcCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttcgcgacgcaaggccgctcccacagggcttgcgatatggcttacctgcgtgcttggcatggcaacttcgacctcagccctggccctggacctg  <  1:3158371‑M2/150‑98 (MQ=255)
                                                                                                                                                   agcagcCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGt  >  2:2918559/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CGCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGT  >  NC_002947/3952924‑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: