Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,791,448 C→G 6.4% G283G (GGC→GGG)  PP_3353 → sulfatase‑modifying protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,791,4480CG 6.4% 291.2 / 9.7 94G283G (GGC→GGG) PP_3353sulfatase‑modifying protein
Reads supporting (aligned to +/- strand):  ref base C (39/49);  new base G (4/2);  total (43/51)
Fisher's exact test for biased strand distribution p-value = 4.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.04e-01

Only 100 of 103 total aligned reads displayed.

GATGGTTTTACCGGTACATCGCCCGTGGGTTGCTTCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGC  >  NC_002947/3791299‑3791596
                                                                                                                                                     |                                                                                                                                                    
  tGGTTTTACCGGTACATCGCCCGTGGGTTGCTTCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAggccc                                                                                                                                                    <  2:1995407/150‑1 (MQ=255)
    gTTTTACCGGTACATCGCCCGTGGGTTGCTTCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAggcccgg                                                                                                                                                  >  2:1799698/1‑150 (MQ=255)
     ttttACCGGTAGATCGCCCGTGGGTTGCTTCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAggcccggc                                                                                                                                                 >  2:1297397/1‑150 (MQ=255)
     ttttACCGGTACATCGCCCGTGGGTTGCTTCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAggcccggc                                                                                                                                                 >  1:370946/1‑150 (MQ=255)
       ttACCGGTACATCGCCCGTGGGTTGCTTCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAggcccggccc                                                                                                                                               <  2:1411261/150‑1 (MQ=255)
          ccGGTACATCGCCCGTGGGTTGCTTCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGa                                                                                                                                            <  2:210969/150‑1 (MQ=255)
          ccGGTACATCGCCCGTGGGTTGCTTCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGa                                                                                                                                            >  2:2151414/1‑150 (MQ=255)
          ccGGTACATCGCCCGTGGGTTGCTTCCCGGCCAACGGCTTTGGACTGTTCGATGCCGGCGGCAATATCGGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGAAGCCAAGCAAGAACCCTCAGGCGCGGCCCTGa                                                                                                                                            >  1:411310/1‑150 (MQ=255)
            ggTACATCGCCCGTGGGTTGCTTCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGt                                                                                                                                          <  2:2246331/150‑1 (MQ=255)
             gTACATCGCCCGTGGGTTGCTTCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGGCCGGGCCTGAggg                                                                                                                                         >  1:353163/1‑148 (MQ=255)
              tACATCGCCCGTGGGTTGCTGCCCGGTCAACGGCGTTGGCCTGTGCGATGCTGGCGGTAATGTCTGGGAGCTGCCTCGCCCGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGa                                                                                                                                        <  2:565244/150‑1 (MQ=255)
               aCATCGCCCGTGGGTTGCTTCCCGGCCAACGGCTTTGGCATGTTCGGTGACGGCGGCAGTGTCTGGGCGCTGACTAGCGCGGTTTAGCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGAc                                                                                                                                       <  2:273511/150‑1 (MQ=255)
                  tCGCCCGTGGGTTGCTTCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGc                                                                                                                                    <  1:1534154/150‑1 (MQ=255)
                  tCGCCCGTGGGTTGCTTCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGc                                                                                                                                    <  2:1254623/150‑1 (MQ=255)
                   cGCCCGTGGGTTGCTTCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCt                                                                                                                                   <  2:1235996/150‑1 (MQ=255)
                    gCCCGTGGGTTGCTTCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGGCGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCtt                                                                                                                                  <  2:125619/150‑1 (MQ=255)
                     cccGTGGGTTGCTTCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCCGCATCGCTGCCTTTATCGGCCAGTCCATGATGGGCAACGCGCCCTCTATCGCTACCCCCCAGGCCCGGCCCTGAGTGGCAGCTTc                                                                                                                                 >  1:1141988/1‑150 (MQ=255)
                      ccGTGGGTTGCTTCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCCTCg                                                                                                                                >  2:390731/1‑150 (MQ=255)
                        gTGGGTTGCTTCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGAc                                                                                                                              <  1:2151414/150‑1 (MQ=255)
                           ggTTGCTTCCCGTCCAACAGCTTTGGCCTGTTCGATGCGGGCTGCAATGTCTTGGAGCTGACTGGCCAGGGCTATCGGCCCGGCCATGACGCACAACGCGACGCCAAGCTCGACCCCCCAGGCCCTGCCCTGGGGAACAGCTGCGACCCg                                                                                                                           >  2:836697/1‑150 (MQ=255)
                                     cGGCCAACGGCTTTGGCCTGTTCGATGCCGGTGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACCGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCg                                                                                                                 <  2:936855/150‑1 (MQ=255)
                                     cGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCg                                                                                                                 >  1:2048930/1‑150 (MQ=255)
                                      ggCCAACGGCTTTGGCCTGTTCAATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCgg                                                                                                                <  2:1046819/150‑1 (MQ=255)
                                       gCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCACGACGCACAGCGCGACGCAAACATAGCGGACGAGCACCCGTCTCTGAGTGACAGCTTCGACCCGGCAGACCCCGGc                                                                                                               >  1:542056/1‑150 (MQ=255)
                                        ccAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGGCCGGGCCTGAggg                                                                                                                                         >  1:917385/1‑121 (MQ=255)
                                        ccAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGGCCGGGCCTGAggg                                                                                                                                         <  2:917385/123‑3 (MQ=255)
                                         cAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGt                                                                                                             <  2:1357829/150‑1 (MQ=255)
                                           aCGGCTTTGGCCTTTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGc                                                                                                           <  2:1053791/150‑1 (MQ=255)
                                                tttGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTg                                                                                                      >  2:1064760/1‑150 (MQ=255)
                                                 ttGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTgg                                                                                                     <  2:1187637/150‑1 (MQ=255)
                                                  tGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAggcccg                                                                                                                                                   <  1:2033449/103‑1 (MQ=255)
                                                  tGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAggcccg                                                                                                                                                   >  2:2033449/1‑103 (MQ=255)
                                                   ggCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCg                                                                                                   <  2:958541/150‑1 (MQ=255)
                                                    gCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCGCGGGCTATCGTCCAAGCCATGACGCACAGCTCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGGGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCgg                                                                                                  >  2:1898711/1‑150 (MQ=255)
                                                    gCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGAccc                                                                                                                            <  2:1508803/124‑1 (MQ=255)
                                                    gCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGAccc                                                                                                                            >  1:1508803/1‑124 (MQ=255)
                                                     ccTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGGCCGGGCCTGAggg                                                                                                                                         <  1:1790785/110‑3 (MQ=255)
                                                     ccTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGGCCGGGCCTGAggg                                                                                                                                         >  2:1790785/1‑108 (MQ=255)
                                                      cTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTa                                                                                                >  1:1848625/1‑150 (MQ=255)
                                                       tGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTaa                                                                                               >  2:2227559/1‑150 (MQ=255)
                                                       tGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGACAGACCCCGGCGTACCGGTGGCGGTaa                                                                                               <  1:1271566/150‑1 (MQ=255)
                                                         ttCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATc                                                                                             <  1:1860390/150‑1 (MQ=255)
                                                         ttCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATc                                                                                             >  2:1860390/1‑150 (MQ=255)
                                                              tGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAgg                                                                                        >  2:774525/1‑150 (MQ=255)
                                                              tGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAgg                                                                                        <  1:679840/150‑1 (MQ=255)
                                                               gCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAggc                                                                                       <  2:370946/150‑1 (MQ=255)
                                                                 cggcggCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAggcgg                                                                                     >  2:458130/1‑150 (MQ=255)
                                                                       cAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCa                                                                               <  2:1077790/150‑1 (MQ=255)
                                                                        aaTGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCAc                                                                              <  1:1922153/150‑1 (MQ=255)
                                                                           gTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCtg                                                                           >  1:1423673/1‑150 (MQ=255)
                                                                              tGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTgtgt                                                                        <  1:2236076/150‑1 (MQ=255)
                                                                               gggAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGTCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGtt                                                                       >  1:242645/1‑150 (MQ=255)
                                                                                   gCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTcggc                                                                   <  1:1297397/150‑1 (MQ=255)
                                                                                           ccACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCt                                                           <  1:1555366/149‑1 (MQ=255)
                                                                                             aCGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCa                                                                               >  2:1991445/1‑128 (MQ=255)
                                                                                             aCGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCa                                                                               <  1:1991445/128‑1 (MQ=255)
                                                                                               gggCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCAt                                                       <  1:777664/150‑1 (MQ=255)
                                                                                                 gCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATgc                                                     <  1:1631040/150‑1 (MQ=255)
                                                                                                   tATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCg                                                                                                   >  1:1036124/1‑102 (MQ=255)
                                                                                                   tATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCg                                                                                                   <  2:1036124/102‑1 (MQ=255)
                                                                                                      cGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTAc                                                <  1:989918/150‑1 (MQ=255)
                                                                                                         ccaggccaTGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGc                                             <  1:458130/150‑1 (MQ=255)
                                                                                                           aggccaTGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCgg                                                                                                  >  1:1046623/1‑95 (MQ=255)
                                                                                                           aggccaTGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCgg                                                                                                  <  2:1046623/95‑1 (MQ=255)
                                                                                                             gccaTGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCt                                         >  2:1987152/1‑150 (MQ=255)
                                                                                                               caTGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCTTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCg                                       <  2:921551/150‑1 (MQ=255)
                                                                                                               caTGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGCCCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGc                                             <  2:1112276/144‑1 (MQ=255)
                                                                                                               caTGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGc                                             >  1:1112276/1‑144 (MQ=255)
                                                                                                                 tGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGc                                     >  2:2016780/1‑150 (MQ=255)
                                                                                                                      cacaGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGGGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCg                                                                                                   <  2:523849/83‑1 (MQ=255)
                                                                                                                      cacaGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCg                                                                                                   >  1:523849/1‑83 (MQ=255)
                                                                                                                      cacaGCGCGACACCAAGTTCGCCCCCCCAGGCCAGGCCCTAAGGGGCAGCTTCGACCCAGCACACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGCCACGcc                                >  2:1235076/1‑150 (MQ=255)
                                                                                                                         aGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGAACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCcagc                             <  2:542056/150‑1 (MQ=255)
                                                                                                                            gcgACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCcagccgc                          >  2:411521/1‑150 (MQ=255)
                                                                                                                              gACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGGGGCGGTAATCAAAGGCGACTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCCCCCCTCGGCACGCcagccgccg                        >  1:981369/1‑148 (MQ=255)
                                                                                                                               aCGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCcagccgcagc                       <  1:390731/150‑1 (MQ=255)
                                                                                                                                 gCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCcagccgcagccg                     >  2:2261058/1‑150 (MQ=255)
                                                                                                                                 gCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCcagccgcagccg                     <  1:774525/150‑1 (MQ=255)
                                                                                                                                 accAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCCCCTGTGTTCGGCGGACCTCTGCATGCGCTACCGCCCCTCGGCACGCcagccgcagccg                     <  1:1502240/149‑1 (MQ=255)
                                                                                                                                  ccAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCgg                    >  2:1370936/1‑150 (MQ=255)
                                                                                                                                   cAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGTGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGgt                   >  1:1840720/1‑150 (MQ=255)
                                                                                                                                    aaGCTCGTCCCTTCAGGCCCGGCCCTAAGTGACATCTTCGTCTCTGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGgtg                  <  1:533757/150‑1 (MQ=255)
                                                                                                                                     aGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCg                                       >  2:1276877/1‑128 (MQ=255)
                                                                                                                                     aGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCg                                       <  1:1276877/128‑1 (MQ=255)
                                                                                                                                     aGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGgtgt                 <  1:1370936/150‑1 (MQ=255)
                                                                                                                                      gCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGtt                <  2:1455976/150‑1 (MQ=255)
                                                                                                                                       cTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTc               >  2:1212333/1‑150 (MQ=255)
                                                                                                                                        tCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTATGCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCa              >  1:2202826/1‑150 (MQ=255)
                                                                                                                                          tgtcacTCAGGGCCGGGCCTGAGGGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATg            >  1:2213697/6‑150 (MQ=255)
                                                                                                                                          gACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTgtgt                                                                        >  2:392532/1‑90 (MQ=255)
                                                                                                                                           aCCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGTATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATga           <  2:411310/150‑1 (MQ=255)
                                                                                                                                           aCCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATga           <  1:2227559/150‑1 (MQ=255)
                                                                                                                                            ccccTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATgac          <  1:1064760/150‑1 (MQ=255)
                                                                                                                                             cccTCAGGCCCGGCCCTGAGGGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGGGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGGTGCATGCGCTACCGCCCCTCGGCACCCCACCCGCAGCGGGTGTTCATgagg         >  1:1958798/1‑148 (MQ=255)
                                                                                                                                              ccTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCCGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATgacga        >  2:775151/1‑150 (MQ=255)
                                                                                                                                                tCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGc                                             <  1:561903/111‑1 (MQ=255)
                                                                                                                                                tCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGc                                             >  2:561903/1‑111 (MQ=255)
                                                                                                                                                 cAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCCGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATGACGACCt     <  1:398323/150‑1 (MQ=255)
                                                                                                                                                  aGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATGACGACCTc    >  2:186551/1‑150 (MQ=255)
                                                                                                                                                    gcccggcccTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGc  >  1:1647628/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GATGGTTTTACCGGTACATCGCCCGTGGGTTGCTTCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGC  >  NC_002947/3791299‑3791596

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