Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,956,096 G→T 5.9% intergenic (‑67/+128) 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,0960GT 5.9% 501.6 / 15.1 154intergenic (‑67/+128)PP_3486/PP_3487cytochrome c/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (84/60);  new base T (4/5);  total (89/65)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.52e-01

Only 100 of 166 total aligned reads displayed.

GACGCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGG  >  NC_002947/3955949‑3956245
                                                                                                                                                   |                                                                                                                                                     
gaCGCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCgca                                                                                                                                                     <  1:753244/150‑1 (MQ=255)
gaCGCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCgca                                                                                                                                                     <  2:5990258/150‑1 (MQ=255)
 aCGCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCgcag                                                                                                                                                    <  1:8766/150‑1 (MQ=255)
     ccTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCagcagc                                                                                                                                                <  2:6942999/150‑1 (MQ=255)
       tCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGccc                                                                                                                                              >  1:6108662/1‑150 (MQ=255)
       tCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGccc                                                                                                                                              >  2:3843535/1‑150 (MQ=255)
        cgcgGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGcccc                                                                                                                                             <  2:2573468/150‑1 (MQ=255)
          cgGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCTGCGCAGcccttt                                                                                                                                           <  1:5774808/150‑4 (MQ=255)
          cgGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCgg                                                                                                                                           <  2:2968763/150‑1 (MQ=255)
          cgGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGCCGCGAAAGGGCTGCGCAGCAGCCCCgg                                                                                                                                           <  2:2189033/150‑1 (MQ=255)
           ggTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAGGGGCTGCTGCGCAGccctttc                                                                                                                                          <  1:6066003/150‑5 (MQ=255)
           ggTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCTGCGCAGccctttc                                                                                                                                          <  2:3439678/150‑5 (MQ=255)
           ggTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCTGCGCAGccctttc                                                                                                                                          <  2:4568612/150‑5 (MQ=255)
             tAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAg                                                                                                                                        >  1:1415273/1‑150 (MQ=255)
              aaaGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGt                                                                                                                                       <  2:739356/150‑1 (MQ=255)
              aaaGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGt                                                                                                                                       <  2:4882416/150‑1 (MQ=255)
              aaaGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGt                                                                                                                                       <  2:4984345/150‑1 (MQ=255)
              aaaGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGt                                                                                                                                       <  2:6602196/150‑1 (MQ=255)
               aaGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGGGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTg                                                                                                                                      <  2:6150968/150‑1 (MQ=255)
                aGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGGCGTGGGTGGCGGAAGGGCAGCGCCGCAGCCCCGGccgggg                                                                                                                                     >  2:2296032/1‑145 (MQ=255)
                aGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTgg                                                                                                                                     >  2:6767626/1‑150 (MQ=255)
                  gCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGATGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGATCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCa                                                                                                                                   <  2:4304977/150‑1 (MQ=255)
                  gCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCa                                                                                                                                   <  1:2305830/150‑1 (MQ=255)
                      ttGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTc                                                                                                                               <  2:2393776/150‑1 (MQ=255)
                          cTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAgc                                                                                                                           >  1:5959392/1‑150 (MQ=255)
                          cTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAgc                                                                                                                           <  2:7408450/150‑1 (MQ=255)
                          cTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCGTTCGTGGCGAAAAGGCAGCGCCGCACCCCCCGCCAAGGCAGCGCAAgc                                                                                                                           >  2:1639657/1‑150 (MQ=255)
                          cTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCGTGGGGCGCGAAAGGGCAGCGCCGCAGCCCCCGCCCGGGCAGCTCAAgc                                                                                                                           >  2:1119504/1‑150 (MQ=255)
                           tGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAgcg                                                                                                                          <  2:1415273/150‑1 (MQ=255)
                           tGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCGTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCCGCCCTGGCAGCACCAgcg                                                                                                                          >  2:2422479/1‑150 (MQ=39)
                             cTTGTTCGTTGTCGTGCCGGTCCAGGGCCAGGGCTGAGGTCTACGTTGCCATGCCAAGCACGCAGTTAAGCAATATCTCAAGCCCTTTGTGAGCGGCCTTGCGTCGCGAAAGGTCTTCGCAGCAGACCCGGCAGTGGCAGCTCAAGCGCa                                                                                                                        <  1:3533323/150‑1 (MQ=255)
                             cTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCa                                                                                                                        <  2:5808044/150‑1 (MQ=255)
                                 ttcgacgtcgggcAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGaaa                                                                                                                    <  2:5932890/144‑1 (MQ=255)
                                 ttCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGaaa                                                                                                                    <  1:3843535/150‑1 (MQ=255)
                                  tCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCa                                                                                                                                         <  1:1027245/128‑1 (MQ=255)
                                     tGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAgg                                                                                                                <  2:6814092/150‑1 (MQ=255)
                                     tGGTCGGTCAGATCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAgg                                                                                                                <  2:4846672/150‑1 (MQ=255)
                                       gtgggtcAGGTCCTGGGCCAGGGCTGAGGTCGATGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCGTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAgggg                                                                                                              <  1:2477768/147‑1 (MQ=255)
                                         cggtcAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAg                                                                                                            <  2:5615847/150‑1 (MQ=255)
                                                 tCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTgcgc                                                                                                    <  2:5223731/150‑1 (MQ=255)
                                                  ccagggccagggCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTgcgcg                                                                                                   <  2:6711674/150‑1 (MQ=255)
                                                   cagggccagggcTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCACGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTgcgcgc                                                                                                  <  2:4353662/150‑1 (MQ=255)
                                                   cagggccagggcTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTgcgcgc                                                                                                  <  2:5620243/150‑1 (MQ=255)
                                                     gggCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTc                                                                                                >  2:5071439‑M2/1‑149 (MQ=255)
                                                         cAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTccttt                                                                                            >  2:1673881‑M2/1‑145 (MQ=255)
                                                         cAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTccttt                                                                                            >  2:2832543‑M2/1‑145 (MQ=255)
                                                          aGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttc                                                                                           >  1:1659308‑M2/1‑144 (MQ=255)
                                                          aGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttc                                                                                           >  2:7199693‑M2/1‑144 (MQ=255)
                                                          aGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttc                                                                                           >  2:7187474‑M2/1‑144 (MQ=255)
                                                          aGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttc                                                                                           >  2:5279769‑M2/1‑144 (MQ=255)
                                                          aGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttc                                                                                           >  2:4244595‑M2/1‑144 (MQ=255)
                                                          aGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCCAAAGGGCGGCGCACCAGCCCCCGCACTGGCAGCGCCACCGCAGCAAGAGGGGAGCGCTGCGCGCGcctata                                                                                           <  2:7335227‑M2/150‑7 (MQ=255)
                                                            ggCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttcgc                                                                                         >  1:3786862‑M2/1‑142 (MQ=255)
                                                               tGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttcgcgac                                                                                      >  1:6610836‑M2/1‑139 (MQ=255)
                                                               tGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttcgcgac                                                                                      >  1:6610750‑M2/1‑139 (MQ=255)
                                                                  ggTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCGTGGGGGGCGAAAGGGCAGCGCACCCGCCCCCCCACGGGCAGCACAAGCGCCCAAAGAGGGGAGCGCAGCGCGCGCCAGTAGCGACGCa                                                                                   >  2:3060704/1‑150 (MQ=255)
                                                                   gTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttcgcgacgcaa                                                                                  >  1:129279‑M2/1‑135 (MQ=255)
                                                                   gTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttcgcgacgcaa                                                                                  >  1:2635454‑M2/1‑135 (MQ=255)
                                                                   gTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCGTGCCGCGCGGAAGGGCGGCGCACCCCCCCCCGCACTGGCCGCCCAAGCGCAGCAAGACGGGGGCGCTGCCCGCGCCGGTcgcgccccca                                                                                  >  2:2316994/1‑144 (MQ=255)
                                                                     cGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttcgcgacgcaagg                                                                                >  1:4497079‑M2/1‑133 (MQ=255)
                                                                     cGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttcgcgacgcaagg                                                                                >  2:475807‑M2/1‑133 (MQ=255)
                                                                      gAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAACGGGCTGCGCAGCCGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttcgcgacgcaaggc                                                                               >  1:1551983‑M2/1‑132 (MQ=255)
                                                                      gAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttcgcgacgcaaggc                                                                               >  2:6448985‑M2/1‑132 (MQ=255)
                                                                      gAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttcgcgacgcaaggc                                                                               >  2:6107957‑M2/1‑132 (MQ=255)
                                                                      gAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttcgcgacgcaaggc                                                                               >  1:7396886‑M2/1‑132 (MQ=255)
                                                                      gAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCCTGCGTCGCGAAACGGCGGCGCCGCAGCCCCCGCCCGGGCCGCTCAAGCCCCCAAAGGGGGGCGCGCGGCGCCCGCCAGTCGCGACGCAAGGc                                                                               >  2:7107200/1‑150 (MQ=255)
                                                                        aGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttcgcgacgcaaggccg                                                                             >  2:4253599‑M2/1‑130 (MQ=255)
                                                                          ttGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGGCGTGCCGCGCGCCAGGGGGGCCGCGCAGCCCCCGCCGGGGCAGCGCCATCGGAGAAAGAGGGGGGCGCGGCGCGCGCCCGTCGCGCCGCAAGGccccg                                                                           >  2:6719296/1‑147 (MQ=255)
                                                                                tgccaAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGGCGTGGGGCGGGGGAGGGGGGCGCCGCAGCCCCCGCCCGGGGGGCTCAAACGCCGGCATAGGGGGGCGCGGGCCCCGCCCTTAGCGGCCCAAAGCCGCGcccccc                                                                     >  2:2155102/1‑147 (MQ=255)
                                                                                        aCGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCGGGCGGCGCGCAAGGGGTGCGCAGCCCCCCCCGCACGGGCAGCCCAAGCGCCCAAACACGGGAGCGCGGCGCGCGCCGTTGGCCACGCAAGGCCGCACCCACcaaacccga                                                             >  2:5428428/1‑145 (MQ=255)
                                                                                                                            gCCGTGCGTCGCTAACGGGCTGCGCCGCCGCCCGGCCTGGGGCCGCTCACGCGCCGAACGTGGGGCGGCCTGCGCGGTCCTGCCGCCCCCCACGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATa                         <  2:3554424/150‑1 (MQ=255)
                                                                                                                             ccTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATAc                        >  1:1284244/1‑150 (MQ=255)
                                                                                                                             ccTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATAc                        >  1:1565246/1‑150 (MQ=255)
                                                                                                                                   gTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttcgcgacgcaaggccgctcccacagggcttgcgatatggcttacctgcgtgcttggcatggcaacttcgacctca                  <  2:4354763‑M2/150‑80 (MQ=255)
                                                                                                                                   gTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttcgcgacgcaaggccgctcccacagggcttgcgatatggcttacctgcgtgcttggcatggcaacttcgacctca                  <  1:5550645‑M2/150‑80 (MQ=255)
                                                                                                                                   gTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttcgcgacgcaaggccgctcccacagggcttgcgatatggcttacctgcgtgattggcatggcaacttagacctca                  <  1:1490632‑M2/150‑80 (MQ=255)
                                                                                                                                   gTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCaccg                  >  1:1667712/1‑150 (MQ=255)
                                                                                                                                   gTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCaccg                  >  1:6723488/1‑150 (MQ=255)
                                                                                                                                    tCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCaccga                 >  1:6966263/1‑150 (MQ=255)
                                                                                                                                    tCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCaccga                 >  1:210503/1‑150 (MQ=255)
                                                                                                                                       cgAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctct                                                                                            <  1:6605313‑M2/72‑6 (MQ=255)
                                                                                                                                       cgAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctct                                                                                            >  2:6605313‑M2/1‑67 (MQ=255)
                                                                                                                                        gAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttcgcgacgcaaggccgctcccacagggcttgcgatatggcttacctgcgtgcttggcatggcaacttcgacctcagccct             <  2:5843503‑M2/150‑85 (MQ=255)
                                                                                                                                        gAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCgg             >  2:2490457/1‑150 (MQ=255)
                                                                                                                                         actggggTGCGCCCCAGCCCCGTCTGGGGCAGCTCACGCGAATAAGGTGGGGCGCGGGGCCCGCCCCTTCCGCGCCGCAGGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGc            <  2:2767051/147‑1 (MQ=255)
                                                                                                                                         aaaGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttcgcgacgcaaggccgctcccacagggcttgcgatatggcttacctgcgtgcttggcatggcaacttcgacctcagccctg            <  2:739798‑M2/150‑86 (MQ=255)
                                                                                                                                           tgggCTGCGCAGCAGCCCCGGCTGTGGCAGCTCACGCGCAGAAAGTGGGGCGCGCTGCGCGCACCTTTCGCGCCGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCa          <  2:4912396/149‑1 (MQ=255)
                                                                                                                                           ggggCTGCTGCGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCa          >  1:7115694/2‑150 (MQ=255)
                                                                                                                                            gggCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttcgcgacgcaaggccgctcccacagggcttgcgatatggcttacctgcgtgcttggcatggcaacttcgacctcagccctggcc         <  2:179235‑M2/150‑89 (MQ=255)
                                                                                                                                            gggCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTcctttcgcgacgcaaggccgctcccacagggcttgcgatatggcttacctgcgtgcttggcatggcaacttcgacctcagccctggcc         <  1:106168‑M2/150‑89 (MQ=255)
                                                                                                                                            gggCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAg         >  1:794364/1‑150 (MQ=255)
                                                                                                                                             ggCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGc        >  2:1167540/1‑150 (MQ=255)
                                                                                                                                              gCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCg       >  1:973552/1‑150 (MQ=255)
                                                                                                                                               cTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACt                       >  2:3607063/1‑133 (MQ=255)
                                                                                                                                                tGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGc     >  1:3579166/1‑150 (MQ=255)
                                                                                                                                                tGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGc     >  1:4353789/1‑150 (MQ=255)
                                                                                                                                                tGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGc     >  1:890546/1‑150 (MQ=255)
                                                                                                                                                tGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGc     >  1:6950488/1‑150 (MQ=255)
                                                                                                                                                  cgcAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCgg                                                              <  2:820178/91‑1 (MQ=255)
                                                                                                                                                  cgcAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTg   >  2:2474653/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GACGCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGG  >  NC_002947/3955949‑3956245

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: