Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,285,003 G→A 7.1% A302V (GCC→GTC)  PP_2887 ← catalase‑like peroxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,285,0030GA 7.1% 155.8 / 2.9 56A302V (GCC→GTC) PP_2887catalase‑like peroxidase
Reads supporting (aligned to +/- strand):  ref base G (28/24);  new base A (2/2);  total (30/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.17e-01

TCCGGTGCTGAGCGATTCGCGGCTACGGCGGTTGAAGGATTCCGAATAGGCGGCCGAGCGGGCGGCGAGGATCGGGTCGGCCGACGCCTCTATACCGCGC‑GGCAGGATCAGCGGGTCGAAGTTGAGGTCGCGGCAGGCCCCTTGCTCGGGGGCGTCTACCTGTTCCAGCACCAGGGTGCCGGCGTCCACGCTGCGCCGCTCGGCGGGCCACGGCCGGGCAGGGTCGTCCACGGCATCGCCTGGTTCGGCCAGTACCAGGCGCAGCGTCCAGCGCAATGGGCCTTGG  >  NC_002947/3284867‑3285152
                                                                                                                                         |                                                                                                                                                     
tCCGGTGCTGAGCGATTCGCGGCTACGGCGGTTGAAGGATTCCGAATAGGCGGCCGAGCGGGCGGCGAGGATCGGGTCGGCCGACGCCTCTATACCGCGC‑GGCAGGATCAGCGGGTCGAAGTTGAGGTCGCGGCAGGCCCCTTGCTCggg                                                                                                                                          <  2:222467/150‑1 (MQ=255)
tCCGGTGCTGAGCGATTCGCGGCTACGGCGGTTGAAGGATTCCGAATAGGCGGCCGAGCGGGCGGCGAGGATCGGGTCGGCCGACGCCTCTATACCGCGC‑GGCAGGATCAGCGGGTCGAAGTTGAGGTCGCGGCAGGCCCCTTGCTCggg                                                                                                                                          <  2:135914/150‑1 (MQ=255)
 ccGGTGCTGAGCGATTCGCGGCTACGGCGGTTGAAGGATTCCGAATAGGCGGCCGAGCGGGCGGCGAGGATCGGGTCGGCCGACGCCTCTATACCGCGC‑GGCAGGATCAGCGGGTCGAAGTTGAGGTCGCGGCAGGCCCCTTGCTCgggg                                                                                                                                         <  2:998707/150‑1 (MQ=255)
  cgctGCCGAGCGATTCGCGGCTACGGCGGATGAAGGATTCCGAATAGGCGGCCCAGCGGGCGGCGAGCATCGGGACGGCCGACGCCTCTATACCGCGC‑GGCAGGACCAGCGGGCCGCAGTTGAGGACGCGGCAGGCCCCTTGCTCggggg                                                                                                                                        >  2:135658/4‑150 (MQ=255)
        tGAGCGATTCGCGGCTACGGCGGTTGAAGGATTCCGAATAGGCGGCCGAGCGGGCGGCGAGGATCGGGTCGGCCGACGCCTCTATACCGCGC‑GGCAGGATCAGCGGGTCGAAGTTGAGGTCGCGGCAGGCCCCTTGCTCGGGGGCGTCTa                                                                                                                                  >  1:640237/1‑150 (MQ=255)
                 cgcgGCTACGGCGGTTGAAGGATTCCGAATAGGCGGCCGAGCGGGCGGCGAGGATCGGGTCGGCCGACGCCTCTATACCGCGC‑GGCAGGATCAGCGGGTCGAAGTTGAGGTCGCGGCAGGCCCCTTGCTCGGGGGCGTCTACCTGTTCCa                                                                                                                         <  1:64910/150‑1 (MQ=255)
                       tACGGCGGTTGAAGGATTCCGAATAGGCGGCCGAGCGGGCGGCGAGGATCGGGTCGGCCGACGCCTCTATACCGCGC‑GGCAGGATCAGCGGGTCGAAGTTGAGGTCGCGGCAGACCCCTTGCTCGGGGGCGTCTACCTGTTCCAGCACCa                                                                                                                   >  2:314952/1‑150 (MQ=255)
                       tACGGCGGTTGAAGGATTCCGAATAGGCGGCCGAGCGGGCGGCGAGGATCGGGTCGGCCGACGCCTCTATACCGCGC‑GGCAGGATCAGCGGGTCGAAGTTGAGGTCGCGGCAGACCCCTTGCTCGGGGGCGTCTACCTGTTCCAGCACCa                                                                                                                   >  2:314755/1‑150 (MQ=255)
                         cggcggTTGAAGGATTCCGAATAGGCGGCCGAGCGGGCGGCGAGGATCGGGTCGGCCGACGCCTCTATACCGCGC‑GGCAGGATCAGCGGGTCGAAGTTGAGGTCGCGGCAGGCCCCTTGCTCGGGGGCGTCTACCTGTTCCAGCACCAgg                                                                                                                 <  2:824753/150‑1 (MQ=255)
                         cggcggTTGAAGGATTCCGAATAGGCGGCCGAGCGGGCGGCGAGGATCGGGTCGGCCGACGCCTCTATACCGCGC‑GGCAGGATCAGCGGGTCGAAGTTGAGGTCGCGGCAGGCCCCTTGCTCGGGGGCGTCTACCTGTTCCAGCACCAgg                                                                                                                 >  1:1216136/1‑150 (MQ=255)
                          ggcggTTGAAGGATTCCGAATAGGCGGCCGAGCGGGCGGCGAGGATCGGGTCGGCCGACGCCTCTATACCGCGC‑GGCAGGATCAGCGGGTCGAAGTTGAGGTCGCGGCAGGCCCCTTGCTCGGGGGCGTCTACCTGTTCCAGCACCAggg                                                                                                                >  2:533566/1‑150 (MQ=255)
                          ggcggTTGAAGGATTCCGAATAGGCGGCCGAGCGGGCGGCGAGGATCGGGTCGGCCGACGCCTCTATACCGCGC‑GGCAGGATCAGCGGGTCGAAGTTGAGGTCGCGGCAGGCCCCTGCCTCGGGGGCGGCTACCTGTTCCAGCACCAggg                                                                                                                >  1:1161150/1‑150 (MQ=255)
                           gcggTTGAAGGATTCCGAATAGGCGGCCGAGCGGGCGGCGAGGATCGGGTCGGCCGACGCCTCTATACCGCGC‑GGCAGGATCAGCGGGTCGAAGTTGAGGTCGCGGCAGGCCCCTTGCTCGGGGGCGTCTACCTGTTCCAGCACCAGGGt                                                                                                               <  2:640237/150‑1 (MQ=255)
                               ttGAAGGATTCCGAATAGGCGGCCGAGCGGGCGGCGAGGATCGGGTCGGCCGACGCCTCTATACCGCGC‑GGCAGGATCAGCGGGTCGAAGTTGAGGTCGCGGCAGGCCCCTTGCTCGGGGGCGTCTACCTGTTCCAGCACCAGGGTGCCg                                                                                                           >  2:213627/1‑150 (MQ=255)
                                 gAAGGATTCCGAATAGGCGGCCGAGCGGGCGGCGAGGATCGGGTCGGCCGACGCCTCTATACCGCGC‑GGCAGGATCAGCGGGTCGAAGTTGAGGTCGCGGCAGGCCCCTTGCTCGGGGGCGTCTACCTGTTCCAGCACCAGGGTGCCGGc                                                                                                         >  2:245492/1‑150 (MQ=255)
                                 gAAGGATTCCGAATAGGCGGCCGAGCGGGCGGCGAGGATCGGGTCGGCCGACGCCTCTATACCGCGC‑GGCAGGATCAGCGGGTCGAAGTTGAGGTCGCGGCAGGCCCATCAACGCTGGGCGTCTACCTGTTCCAGCACCAGGGTGCCGGc                                                                                                         >  1:1142839/1‑150 (MQ=255)
                                   aGGATTCCGAATAGGCGGCCGAGCGGGCGGCGAGGATCGGGTCGGCCGACGCCTCTATACCGCGC‑GGCAGGATCAGCGGGTCGAAGTTGAGGTCGCGGCAGGCCCCTTGCTCgggg                                                                                                                                         <  1:1199166/116‑1 (MQ=255)
                                   aGGATTCCGAATAGGCGGCCGAGCGGGCGGCGAGGATCGGGTCGGCCGACGCCTCTATACCGCGC‑GGCAGGATCAGCGGGTCGAAGTTGAGGTCGCGGCAGGCCCCTTGCTCgggg                                                                                                                                         >  2:1199166/1‑116 (MQ=255)
                                   aGGATTCCGAATAGGCGGCCGAGCGGGCGGCGAGGATCGGGTCGGCCGACGCCTCTATACCGCGC‑GGCAGGATCAGCGGGTCGAAGTTGAGGTCGCGGCAGGCCCCTTGCTCGGGGGCGTCTACCTGTTCCAGCAcc                                                                                                                    >  2:1319913/1‑137 (MQ=255)
                                   aGGATTCCGAATAGGCGGCCGAGCGGGCGGCGAGGATCGGGTCGGCCGACGCCTCTATACCGCGC‑GGCAGGATCAGCGGGTCGAAGTTGAGGTCGCGGCAGGCCCCTTGCTCGGGGGCGTCTACCTGTTCCAGCAcc                                                                                                                    <  1:1319913/137‑1 (MQ=255)
                                    ggATTCCGAATAGGCGGCCGAGCGGGCGGCGAGGATCGGGTCGGCCGACGCCTCTATACCGCGC‑GGCAGGATCAGCGGGTCGAAGTTGAGGTCGCGGCAGGCCCCTTGCTCGGGGGCGTCTACCTGTTCCAGCACCAGGGTGCCGGCGTc                                                                                                      >  2:207482/1‑150 (MQ=255)
                                            aaTAGGCGGCCGAGCGGGCGGCGAGGTTCGGGTCGGGCGACGCCTCTATACCGCGC‑GGCAGGGTTAGCGGGTCGAAGTTGAGGTCGTGGCAGGCCCCTTGCTCGGGGGCGTCTACCTGTTCCAGCACCAGGGTGCCGGCGTCCACGCTgc                                                                                              <  1:135658/150‑1 (MQ=255)
                                            aaTAGGCGGCCGAGCGGGCGGCGAGGATCGGGTCGGCCGACGCCTCTATACCGCGC‑GGCAGGATCAGCGGGTCGAAGTTGAGGTCGCGGCAGACCCCTTGCTCGGGGGCGTCTACCTGTTCCAGCACCAGGGTGCCGGCGTCCACGCTgc                                                                                              <  1:314755/150‑1 (MQ=255)
                                            aaTAGGCGGCCGAGCGGGCGGCGAGGATCGGGTCGGCCGACGCCTCTATACCGCGC‑GGCAGGATCAGCGGGTCGAAGTTGAGGTCGCGGCAGACCCCTTGCTCGGGGGCGTCTACCTGTTCCAGCACCAGGGTGCCGGCGTCCACGCTgc                                                                                              <  1:314952/150‑1 (MQ=255)
                                                ggcggcCGAGCGGGCGGCGAGGATCGGGTCGGCCGACGCCTCTATACCGCGC‑GGCAGGATCAGCGGGTCGAAGTTGAGGTCGCGGCAGGCCCCTTGCTCGGGGGCGTCTACCTGTTCCAGCACCAGGGTGCCGGCGTCCACGCTGcgccg                                                                                          >  2:740964/1‑150 (MQ=255)
                                                 gcggcCGAGCGGGCGGCGAGGATCGGGTCGGCCGACGCCTCTATACCGCGC‑GGCAGGATCAGCGGGTCGAAGTTGAGGTCGCGGCAGGCCCCTTGCTCGGGGGCGTCTACCTGTTCCAGCACCAGCGTGCCGGCGTCCACGCTGcgccgc                                                                                         >  2:497079/1‑150 (MQ=255)
                                                  cggcCGAGCGGGCGGCGAGGATCGGGTCGGCCGACGCCTCTATACCGCGC‑GGCAGGATCAGCGGGTCGAAGTTGAGGTCGCGGCAGGCCCCTTGCTCGGGGGCGTCTACCTGTTCCAGCACCAGGGTGCCGGCGTCCACGCTGCGCCGCt                                                                                        <  1:159274/150‑1 (MQ=255)
                                                        aGCGGGCGGCGAGGATCGGGTCGGCCGACGCCTCTATACCGCGC‑GGCAGGATCAGCGGGTCGAAGTTGAGGTCGCGGCAGGCCCCTTGCTCGGGGGCGTCTACCTGTTCCAGCACCAGGGTGCCGGCGTCCACGCTGCGCCGCTcggcgg                                                                                  <  1:533566/150‑1 (MQ=255)
                                                               ggcgAGGATCGGGTCGGCCGACGCCTCTATACCGCGC‑GGCAGGATCAGCGGGTCGAAGTTGAGGTCGCGGCAGGCCCCTTGCTCGGGGGCGTCTACCTGTTCCAGCACCAGGGTGCCg                                                                                                           >  1:607406/1‑118 (MQ=255)
                                                               ggcgAGGATCGGGTCGGCCGACGCCTCTATACCGCGC‑GGCAGGATCAGCGGGTCGAAGTTGAGGTCGCGGCAGGCCCCTTGCTCGGGGGCGTCTACCTGTTCCAGCACCAGGGTGCCg                                                                                                           <  2:607406/118‑1 (MQ=255)
                                                               ggcgAGGATCGGGTCGGCCGACGCCTCTATACCGCGC‑GGCAGGATCAGCGGGTCGAAGTTGAGGTCGCGGCAGGCCCCTTGCTCGGGGGCGTCTACCTGTTCCAGCACCAGGGTGCCGGCGTCCACGCTGCGCCGCTCGGCGGGCCACgg                                                                           >  2:1376197/1‑150 (MQ=255)
                                                                  gAGGATCGGGTCGGCCGACGCCTCTATACCGCGC‑GGCAGGATCAGCGGGTCGAAGTTGAGGTCGCGGCAGGCCCCTTGCTCGGGGGCGTCTACCTGTTCCAGCACCAGGGTGCCGGCGTCCACGCTGCGCCGCTCGGCGGGCCACGGCCg                                                                        >  1:624624/1‑150 (MQ=255)
                                                                   aGGTTCGGGTCGGCCGACGCCTCTATACCGCGC‑GGCAGGATCAGCGGGTCGAAGTTGAGGTCGCGGCAGGCCCCTTGCTCGGGGGCGTCTACCTGTTCCAGCACCAGGGTGCCGGCGTCCACGCTGCGCCGCTCGGCGGGCCACGGCCgg                                                                       >  2:541100/1‑150 (MQ=255)
                                                                   aGGATCGGGTCGGCCGACGCCTCTATACCGCGC‑GGCAGGATCAGCGGGTCGAAGTTGAGGTCGCGGCAGGCCCCTTGCTCGGGGGCGTCTACCTGTTCCAGCACCAGGGTGCCGGCGTCCACGCTGCGCCGCTCGGCGGGCCACGGCCgg                                                                       <  1:1376197/150‑1 (MQ=255)
                                                                     gATCGGGTCGGCCGACGCCTCTATACCGCGC‑GGCAGGATCAGCGGGTCGAAGTTGAGGTCGCGGCAGGCCCCTTGCTCGGGGGCGTCTACCTGTTCCAGCACCAGGGTGCCGGCGTCCACGCTGCGCCGCTcggcg                                                                                   <  1:132543/136‑1 (MQ=255)
                                                                     gATCGGGTCGGCCGACGCCTCTATACCGCGC‑GGCAGGATCAGCGGGTCGAAGTTGAGGTCGCGGCAGGCCCCTTGCTCGGGGGCGTCTACCTGTTCCAGCACCAGGGTGCCGGCGTCCACGCTGCGCCGCTcggcg                                                                                   >  2:132543/1‑136 (MQ=255)
                                                                      aTCGGGTCGGCCGCCACCTCTCACCCGTGCAGGCCGG‑TCAGCGGGTCGAAGTTGAGGTCGCGGCAGGCCCCTTGCTCGGGGGCGTCTACCTGTTCCAGCACCAGGGTGCCGGCGTCCACGCTGCGCCGCTCGGCGGGCCACGGCCGGGCa                                                                    <  2:1142839/150‑1 (MQ=255)
                                                                        cgtggcGGCCGACGCCTCTATACCGCGC‑GGCAGGATCAGCGGGTCGAAGTTGAGGTCGCGGCAGGCCCCTTGCTCGGGGGCGTCTACCTGTTCCAGCACCAGGGTGCCGGCGTCCACGCTGCGCCGCTCGGCGGGCCACGGCCGGGCAgg                                                                  <  1:1059147/145‑1 (MQ=255)
                                                                        cGGGTCGGCCGACGCCTCTATACCGCGC‑GGCAGGATCAGCGGGTCGAAGTTGAGGTCGCGGCAGGCCCCTTGCTCGGGGGCGTCTACCTGTTCCAGCACCAGGGTGCCGGCGTCCACGCTGCGCCGCTCGGCGGGCCACGGCCGGGCAgg                                                                  <  2:1216136/150‑1 (MQ=255)
                                                                        cGGGTCGGCCGACGCCTCTATACCGCGC‑GGCAGGATCAGCGGGTCGAAGTTGAGGTCGCGGCAGGCCCCTTGCTCGGGGGCGTCTACCTGTTCCAGCACCAGGGTGCCGGCGTCCACGCTGCGCCGCTCGGCGGGCCACGGCCGGGCAgg                                                                  <  1:541100/150‑1 (MQ=255)
                                                                                 cGACGCCTCTATTCCGCGC‑GGCAGGATCAGCGGGTCGAAGTTGAGGTCGCGGCAGGCCCCTTGCTCGGGGGCGTCTACCTGTTCCAGCACCAGGGTGCCGGCGTCCACGCTGCGCCGCTCGGCGGGCCACGGCCGGGCAGGGTCGTCCAc                                                         <  1:574700/150‑1 (MQ=255)
                                                                                   aCGCCTCTATACCGCGC‑GGCAGGATCAGCGGGTCGAAGTTGAGGTCGCGGCAGGCCCCTTGCTCGGGGGCGTCTACCTGTTCCAGCACCAGGGTGCCGGCGTCCACGCTGCGCCGCTCGGCGGGCCACGGCCGGGCAGGGGCGTCCACgg                                                       >  1:879435/1‑150 (MQ=255)
                                                                                       ctctATACCGCGC‑GGCAGGATCAGCGGGTCGAAGTTGAGGTCGCGGCAGGCCCCTTGCTCGGGGGCGTCTACCTGTTCCAGCACCAGGGTGCCGGCGTCCACGCTGCGCCGCTCGGCGGGCCACGGCCgg                                                                       <  1:807387/130‑1 (MQ=255)
                                                                                       ctctATACCGCGC‑GGCAGGATCAGCGGGTCGAAGTTGAGGTCGCGGCAGGCCCCTTGCTCGGGGGCGTCTACCTGTTCCAGCACCAGGGTGCCGGCGTCCACGCTGCGCCGCTCGGCGGGCCACGGCCgg                                                                       >  2:807387/1‑130 (MQ=255)
                                                                                       ctctATACCGCGC‑GGCAGGATCAGCGGGTCGAAGTTGAGGTCGCGGCAGGCCCCTTGCTCGGGGGCGTCTACCTGTTCCAGCACCAGGGTGCCGGCGTCCACGCTGCGCCGCTCGGCGGGCCACGGCCGGGCAGGGTCGTCCACGGCATc                                                   >  2:645036/1‑150 (MQ=255)
                                                                                            taCCGCGC‑GGCAGGATCAGCGGGTCGAAGTTGAGGTCGCGGCAGGCCCCTTGCTCGGGGGCGTCTACCTGTTCCAGCACCAGGGTGCCGGCGTCCACGCTGCGCCGCTCGGCGGGCTACGGCCGGGCAGGGTCGTCCACGGCATCGCCTg                                              <  2:578278/150‑1 (MQ=255)
                                                                                                gcgc‑gGCAGGATCAGCGGGTCGAAGTTGAGGTCGCGGCAGGCCCCTTGCTCGGGGGCGTCTACCTGTTCCAGCACCAGGGTGCCGGCGTCCACGCTGCGCCGCTCGGCGGGCCACGGCCGGGCAGGGTCGTCCACGGCATCGCCTGGTTc                                          <  1:1276678/150‑1 (MQ=255)
                                                                                                         ggATCAGCGGGTCGAAGTTGAGGTCGCGGCAGGCCCCTTGCTCGGGGGCGTCTACCTGTTCCAGCACCAGGGTGCCGGCGTCCACGCTGCGCCGCTCGGCGGGCCACGGCCGGGCAGGGTCGTCCACGGCATCGCCTGGTTCGGCCAGTa                                  >  1:582435/1‑150 (MQ=255)
                                                                                                               gCGGGTCGAAGTTGAGGTCGCGGCAGGCCCCTTGCTCggg                                                                                                                                          >  1:970344/1‑40 (MQ=255)
                                                                                                               gCGGGTCGAAGTTGAGGTCGCGGCAGGCCCCTTGCTCggg                                                                                                                                          <  2:970344/40‑1 (MQ=255)
                                                                                                               gCGGGCCGACGTTGAGGTCGCGGCAGGCCCCTTGCTCGGGGGCGTCTACCTGTTCAAGCAACAGGGTGCCGGCGTCCACGCTGCGCCGCTCGGCGGGCCAAGGCCGGGCAGGGTCGACCACGGCAGCGACTGGTTCGGCCACTACGAGgc                            >  2:275967/1‑150 (MQ=255)
                                                                                                                     cGAAGTTGAGGTCGCGGCAGGCCCCTTGCTCGGGGGCGTCTACCTGTTCCAGCACCAGGGTGCCGGCGTCCACGCTGCGCCGCTCGGCGGGCCACGGCCGGGCAGGGTCGTCCACGGCATCGCCTGGTTCGGCCAGTACCAGGCGCAGCg                      >  2:564632/1‑150 (MQ=255)
                                                                                                                         gTTGAGGTCGCGGCAGGCCCCTTGCTCGGGGGCGTCTACCTGTTCCAGCACCAGGGTGCCGGCGTCCACGCTGCGCCGCTCGGCGGGCCACGGCCGGGCAGGGTCGTCCACGGCATCGCCTGGTTCGGCCAGTACCAGGCGCAGCGTCCa                  >  1:575582/1‑150 (MQ=255)
                                                                                                                           tGAGGTCGCGGCAGGCCCCTTGCTCGGGGGCGTCTACCTGTTCCAGCACCAGGGTGCCGGCGTCCACGCTGCGCCGCTCGGCGGGCCACGGCCGGGCAGGGTCGTCCACGGCATCGCCTGGTTCGGCCAGTACCAGGCGCAGCGTCCAgc                <  1:564632/150‑1 (MQ=255)
                                                                                                                             aGGTCGCGGCAGGCCCCTTGCTCGGGGGCGTCTACCTGTTCCAGCACCAGGGTGCCGGCGTCCACGCTGCGCCGCTCGGCGGGCCACGGCCGGGCAGGGTCGTCCACGGCATCGCCTGGTTCGGCCAGTACCAGGCGCAGCGTCCAgcgc              >  1:960637/1‑150 (MQ=255)
                                                                                                                                       aGGCCCCTTGCTCGGGGGCGTCTACCTGTTCCAGCACCAGGGTGCCGGCGTCCACGCTGCGCCGCTCGGCGGGCCACGGCCGGGCAGGGTCGTCCACGGCATCGCCTGGTTCGGCCAGTACCAGGCGCAGCGTCCAGCGCAATGGGCCtt    >  1:7697/1‑150 (MQ=255)
                                                                                                                                        ggCCCCTTGCTCGGGGGCGTCTACCTGTTCCAGCACCAGGGTGCCGGCGTCCACGCTGCGCCGCTCGGCGGGCCACGGCCGGGCAGGGTCGTCCACGGCATCGCCTGGTTCGGCCAGTACCAGGCGCAGCGTCCAGCGCAATGGGCCTTg   <  1:740964/150‑1 (MQ=255)
                                                                                                                                        ggCCCCTTGCTCGGGGGCGTCTACCTGTTCCAGCACCAGGGTGCCGGCGTCCACGCTGCGCCGCTCGGCGGGCCACGGCCGGGCAGGGTCGTCCACGGCATCGCCTGGTTCGGCCAGTACCAGGCGCAGCGTCCAGCGCAATGGGCCTTg   >  1:551852/1‑150 (MQ=255)
                                                                                                                                         gCCCCTTGCTCGGGGGCGTCTACCTGTTCCAGCACCAGGGTGCCGGCGTCCACGCTGCGCCGCTCGGCGGGCCACGGCCGGGCAGGGTCGTCCACGGCATCGCCTGGTTCGGCCAGt                                   >  1:1289929/1‑117 (MQ=255)
                                                                                                                                         gCCCCTTGCTCGGGGGCGTCTACCTGTTCCAGCACCAGGGTGCCGGCGTCCACGCTGCGCCGCTCGGCGGGCCACGGCCGGGCAGGGTCGTCCACGGCATCGCCTGGTTCGGCCAGt                                   <  2:1289929/117‑1 (MQ=255)
                                                                                                                                         gCCCCTTGCTCGGGGGCGTCTACCTGTTCCAGCACCAGGGTGCCGGCGTCCACGCTGCGCCGCTCGGCGGGCCACGGCCGGGCAGGGTCGTCCACGGCATCGCCTGGTTCGGCCAGTACCAGGCGCAGCGTCCAGCGCAATGGGCCTTgg  <  1:497079/150‑1 (MQ=255)
                                                                                                                                         gCCCCTTGCTCGGGGGCGTCTACCTGTTCCAGCACCAGGGTGCCGGCGTCCACGCTGCGCCGCTCGGCGGGCCACGGCCGGGCAGGGTCGTCCACGGCATCGCCTGGTTCGGCCAGTACCAGGCGCAGCGTCCAGCGCAATGGGCCTTgg  <  2:935006/150‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                                     
TCCGGTGCTGAGCGATTCGCGGCTACGGCGGTTGAAGGATTCCGAATAGGCGGCCGAGCGGGCGGCGAGGATCGGGTCGGCCGACGCCTCTATACCGCGC‑GGCAGGATCAGCGGGTCGAAGTTGAGGTCGCGGCAGGCCCCTTGCTCGGGGGCGTCTACCTGTTCCAGCACCAGGGTGCCGGCGTCCACGCTGCGCCGCTCGGCGGGCCACGGCCGGGCAGGGTCGTCCACGGCATCGCCTGGTTCGGCCAGTACCAGGCGCAGCGTCCAGCGCAATGGGCCTTGG  >  NC_002947/3284867‑3285152

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑