Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,730,937 C→G 7.9% R77R (CGC→CGG polB → DNA polymerase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,730,9370CG 7.9% 272.2 / 13.1 89R77R (CGC→CGGpolBDNA polymerase II
Reads supporting (aligned to +/- strand):  ref base C (39/43);  new base G (4/3);  total (43/46)
Fisher's exact test for biased strand distribution p-value = 7.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.47e-01

ACCGATGCCGGGCCGCGACTGTTGCGGCTGGCGGCGCAGGAGTCGGTGGCTTTCATTCCTCAGGCGCAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAGGCCGATATCCGGCCGCCGGAGCGTTACCTCATGGAGCGCTTCATCACGGCACCGGTGCAGTT  >  NC_002947/2730788‑2731086
                                                                                                                                                     |                                                                                                                                                     
aCCGATGCCGGGCCGCGACTGTTGCGGCTGGCGGCGCAGGAGTCGGTGGCTTTCATTCCTCAGGCGCAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGc                                                                                                                                                       <  2:1512587/150‑1 (MQ=255)
          ggCCGCGACTGTTGCGGCTGGCGGCGCAGGAGTCGGTGGCTTTCATTCCTCAGGCGCAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCg                                                                                                                                             <  1:390330/150‑1 (MQ=255)
                    gTTGCGGCTGGCGGCGCAGGAGTCGGTGGCTTTCATTCCTCAGGCGCAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATAc                                                                                                                                   <  2:963490/150‑1 (MQ=255)
                        cGGCTGGCGGCGCAGGAGTCGGTGGCTTTCATTCCTCAGGCGCAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGc                                                                                                                               <  2:1571227/150‑1 (MQ=255)
                          gCTGGCGGCGCAGGAGTCGGTGGCTTTCATTCCTCAGGCGCAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCa                                                                                                                             <  2:824210/150‑1 (MQ=255)
                                  cgcAGGAGTCGGTGGCTTTCATTCCTCAGGCGCAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGcc                                                                                                                     <  2:741127/150‑1 (MQ=255)
                                       gAGTCGGTGGCTTTCATTCCTCAGGCGCAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTg                                                                                                                >  2:546176/1‑150 (MQ=255)
                                          tCGGTGGCTTTCATTCCTCAGGCGCAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATg                                                                                                             <  2:1279068/150‑1 (MQ=255)
                                             gTGGCTTTCATTCCTCAGGCGCAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCaa                                                                                                          >  1:290933/1‑150 (MQ=255)
                                             gTGGCTTTCATTCCTCAGGCGCAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCaa                                                                                                          >  1:1443007/1‑150 (MQ=255)
                                              tGGCTTTCATTCCTCAGGCGCAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGAt                                                                                                              >  1:54192/1‑145 (MQ=255)
                                              tGGCTTTCATTCCTCAGGCGCAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGAt                                                                                                              <  2:54192/145‑1 (MQ=255)
                                                gCTTTCATTCCTCAGGCGCAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTg                                                                                                       <  2:904440/150‑1 (MQ=255)
                                                  tttCATTCCTCAGGCGCAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGc                                                                                                                               >  2:792702/1‑124 (MQ=255)
                                                  tttCATTCCTCAGGCGCAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGc                                                                                                                               <  1:792702/124‑1 (MQ=255)
                                                  tttCATTCCTCAGGCGCAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGa                                                                                                     <  1:1067857/150‑1 (MQ=255)
                                                  tttCATTCCTCAGGCGCAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGa                                                                                                     >  1:2276490/1‑150 (MQ=255)
                                                  tttCATTCCTCAGGCGCAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGa                                                                                                     >  1:1348437/1‑150 (MQ=255)
                                                  tttCATTCCTCAGGCGCAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTCCGAGCATCGGGCGATGCTCgaa                                                                                                                                           <  2:1756230/112‑3 (MQ=255)
                                                  tttCATTCCTCAGGCGCAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTCCGAGCATCGGGCGATGCTCgaa                                                                                                                                           >  1:1756230/1‑110 (MQ=255)
                                                  tttCATTCCTCAGGCGCAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTAAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGa                                                                                                     <  1:359476/150‑1 (MQ=255)
                                                  tttCATTCCTCAGGCGCAGGCAGAGCACGCCCGTCTGCTGCTGGCCAAAGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGa                                                                                                     >  2:155930/1‑150 (MQ=255)
                                                   ttCATTCCTCAGGCGCAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGag                                                                                                    <  2:1672847/150‑1 (MQ=255)
                                                       ttCCTCAGGCGCAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGagcagc                                                                                                <  2:2276490/150‑1 (MQ=255)
                                                       ttCCTCAGGCGCAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGagcagc                                                                                                >  1:57968/1‑150 (MQ=255)
                                                         ccTCAGGCGCAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATAACCGCCAGCATCCCCAAATAACGCACCTGAAGCAGCGt                                                                                              >  2:13784/1‑150 (MQ=255)
                                                          cTCAGGCGCAGGCCGAGCGCGCCCGTGTGCTGCTGGCCAACGAACACGGGACAAAGCTAAAGGCAGTCGACCTCAAGGACTTCGACAATAACCACGTGCTTCATCTCTATCTTCTTCACCATCGCCAATTTATGCAACTGGAGCAGCGTc                                                                                             <  2:1479263/150‑1 (MQ=255)
                                                           tCAGGCGCAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCtata                                                                                                                                    <  1:424135/110‑1 (MQ=255)
                                                           tCAGGCGCAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCtata                                                                                                                                    >  2:424135/1‑110 (MQ=255)
                                                           tCAGGCGCAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCt                                                                                            <  2:45343/150‑1 (MQ=255)
                                                                 gcAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGagca                                                                                                  >  1:2171510/1‑138 (MQ=255)
                                                                 gcAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGagca                                                                                                  <  2:2171510/138‑1 (MQ=255)
                                                                 gcAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCAc                                                                                      >  1:2061857/1‑150 (MQ=255)
                                                                    ggCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGAt                                                                                                                                                  >  1:1588661/1‑87 (MQ=255)
                                                                    ggCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGAt                                                                                                                                                  <  2:1588661/87‑1 (MQ=255)
                                                                    ggCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGc                                                                                   <  2:1443007/150‑1 (MQ=255)
                                                                      ccGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCg                                                                                 <  2:663558/150‑1 (MQ=255)
                                                                       cGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACTAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCCTCGCGTAGATCTGGGAGTGGAGCAGCGTCTGCGCACGGCCgg                                                                                >  2:194128/1‑150 (MQ=255)
                                                                         aGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTa                                                                              <  2:321780/150‑1 (MQ=255)
                                                                         aGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTa                                                                              <  2:57968/150‑1 (MQ=255)
                                                                          gCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGCGCAGCGTCTGCGCACGGCCGGTAt                                                                             >  2:1128315/1‑150 (MQ=255)
                                                                          gCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTAt                                                                             >  2:1860814/1‑150 (MQ=255)
                                                                           cACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATc                                                                            <  1:565820/150‑1 (MQ=255)
                                                                           cACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATc                                                                            >  1:625909/1‑150 (MQ=255)
                                                                             cGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGa                                                                          <  1:1860814/150‑1 (MQ=255)
                                                                             cGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGa                                                                          <  1:443972/150‑1 (MQ=255)
                                                                                ccGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACgt                                                                       <  2:625909/150‑1 (MQ=255)
                                                                                    gtgCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATTCCCGCCAGCATCGCCAATTGATCCAACTTGAGCGTCGCCTGCCCACGGCGCGTATCGCCGTGTTc                                                                   <  2:1061753/150‑1 (MQ=255)
                                                                                              ccAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGt                                                                                              <  1:1436620/113‑1 (MQ=255)
                                                                                              ccAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGt                                                                                              >  2:1436620/1‑113 (MQ=255)
                                                                                                aaCGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAGGCCGATATc                                                       >  2:1770237/1‑150 (MQ=255)
                                                                                                  cgaacgcGGGACGGAGCTTAAGCCAGTGCGCCTGAAGGACTTCAAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCACTTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAGGCCGATATCCg                                                     <  2:290933/144‑1 (MQ=255)
                                                                                                    aaaCCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAGGCCGATATCCgg                                                    <  2:1028357/149‑1 (MQ=255)
                                                                                                    aaaCCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAGGCCGATATCCgg                                                    >  1:1028357/1‑149 (MQ=255)
                                                                                                          ggACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAGGCCGATATCCGGCCGCCgg                                             >  1:2112434/1‑150 (MQ=255)
                                                                                                           gACGAAGCTCAAGCCACTGCGCCTGAAGGACTTGGTGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAGGCCGATATCCGGCCGCCGGa                                            >  1:1398677/1‑150 (MQ=255)
                                                                                                              gAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCt                                                                                            >  1:1653263/1‑99 (MQ=255)
                                                                                                              gAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCt                                                                                            <  2:1653263/99‑1 (MQ=255)
                                                                                                              gAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCg                                                                                 <  1:305000/110‑1 (MQ=255)
                                                                                                              gAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCg                                                                                 >  2:305000/1‑110 (MQ=255)
                                                                                                               aaGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAGGcc                                                             >  2:968844/1‑129 (MQ=255)
                                                                                                               aaGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAGGcc                                                             <  1:968844/129‑1 (MQ=255)
                                                                                                                  gtgctgCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAGGCCGATATCCGGCCGCCGGAGCGTTAc                                     >  2:2303491/6‑150 (MQ=255)
                                                                                                                  cTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGa                                                                 >  1:463483/1‑122 (MQ=255)
                                                                                                                  cTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGa                                                                 <  2:463483/122‑1 (MQ=255)
                                                                                                                   tCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAGGCCGATATCCGGCCGCCGGAGCGTTAcc                                    >  2:2095885/1‑150 (MQ=255)
                                                                                                                     aaGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCgg                                                                                <  1:1301595/104‑1 (MQ=255)
                                                                                                                     aaGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCgg                                                                                >  2:1301595/1‑104 (MQ=255)
                                                                                                                        ccACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAGGCCGATATCCGGCCGCCGGAGCGTTACCTCATg                               <  1:1770237/150‑1 (MQ=255)
                                                                                                                        ccACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAGGCCGATATCCGGCCGCCGGAGCGTTACCTCATg                               <  1:546176/150‑1 (MQ=255)
                                                                                                                           cTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACAAGTCATCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAGGCCGATATCCGGCCGCCGGAGCGTTACCTCATGGAg                            <  1:2329639/150‑1 (MQ=255)
                                                                                                                            tGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAGGCCGATATCCGGCCGCCGGAGCGTTACCTCATGGAgc                           <  1:155930/150‑1 (MQ=255)
                                                                                                                             gcgcCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAGGCCGATATCCGGCCGCCGGAGCGTTACCTCATGGAgcg                          >  2:850674/1‑150 (MQ=255)
                                                                                                                                  tGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAGGCCGATATCCGGCCGCCGGAGCGTTACCTCATGGAGCGCTtc                      <  1:559735/149‑1 (MQ=255)
                                                                                                                                  tGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAGGCCGATATCCGGCCGCCGGAGCGTTACCTCATGGAGCGCTtc                      >  2:559735/1‑149 (MQ=255)
                                                                                                                                       gaggccGAGCATCGGGCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAGGCCGATATCCGGCCGCCGGAGCGTTAc                                     >  1:570524/6‑129 (MQ=255)
                                                                                                                                       gaggccGAGCATCGGGCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAGGCCGATATCCGGCCGCCGGAGCGTTAc                                     <  2:570524/124‑1 (MQ=255)
                                                                                                                                         cTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAGGCCGATATCCGGCCGCCGGAGCGTTACCTCATg                               <  1:1952041/133‑1 (MQ=255)
                                                                                                                                         cTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAGGCCGATATCCGGCCGCCGGAGCGTTACCTCATg                               >  2:1952041/1‑133 (MQ=255)
                                                                                                                                          ttCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAgg                                                               <  1:630363/100‑1 (MQ=255)
                                                                                                                                          ttCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAGGCCGATATCCGGCCGCCGGAGCGTTACCTCATGGAGCGCTTCATCACGGCa             >  1:2014990/1‑150 (MQ=255)
                                                                                                                                          ttCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAGGCCGATATCCGGCCGCCGGAGCGTTACCTCATGGAGCGCTTCATCACGGCa             >  2:463718/1‑150 (MQ=255)
                                                                                                                                          ttCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAGGCCGATATCCGGCCGCCGGAGCGTTACCTCATGGAGCGCTTCATCACGGCa             >  2:2189827/1‑150 (MQ=255)
                                                                                                                                          ttCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGAACGGCCGGTATCGACGTGTTCGAgg                                                               >  2:630363/1‑100 (MQ=255)
                                                                                                                                           tCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAGGCCGATATCCGGCCGCCGGAGCGTTACCTCATGGAGCGCTTCATCACGGCAc            <  2:1348437/150‑1 (MQ=255)
                                                                                                                                           ccGAGCATCGGGCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAGGCCGATATcc                                                      >  2:1692501/2‑108 (MQ=255)
                                                                                                                                           ccGAGCATCGGGCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAGGCCGATATcc                                                      <  1:1692501/107‑1 (MQ=255)
                                                                                                                                            cGAGCATCGGGCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAGGCCGATATCCGGCCGCCGGAGCGTTACCTCATGGAGCGCTTCATCACGGCAcc           >  1:320385/1‑150 (MQ=255)
                                                                                                                                                  tCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAGGCCGATATCCGGCCGCCGGAGCGTTACCTCATGGAGCGCTTCATCACGGCACCGGTGCa     >  2:2085271/1‑150 (MQ=255)
                                                                                                                                                    gCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAGGCCGATATCCGGCCGCCGGAGCGTTACCTCATGGAGCGCTTCATCACGGCACCGGTGCAGt   >  2:168620/1‑150 (MQ=255)
                                                                                                                                                     cccGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAGGCCg                                                            <  2:1469544/92‑1 (MQ=255)
                                                                                                                                                     cccGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAGGCCg                                                            >  1:1469544/1‑92 (MQ=255)
                                                                                                                                                     cccGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAGGCCGATATCCGGCCGCCGGAGCGTTACCTCATGGAGCGCTTCATCACGGCACCGGTGCAGtt  >  1:1502825/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACCGATGCCGGGCCGCGACTGTTGCGGCTGGCGGCGCAGGAGTCGGTGGCTTTCATTCCTCAGGCGCAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAGGCCGATATCCGGCCGCCGGAGCGTTACCTCATGGAGCGCTTCATCACGGCACCGGTGCAGTT  >  NC_002947/2730788‑2731086

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