Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,791,453 C→G 6.1% A285G (GCC→GGC)  PP_3353 → sulfatase‑modifying protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,791,4530CG 6.1% 346.6 / 10.6 110A285G (GCC→GGC) PP_3353sulfatase‑modifying protein
Reads supporting (aligned to +/- strand):  ref base C (54/49);  new base G (3/4);  total (57/53)
Fisher's exact test for biased strand distribution p-value = 7.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

Only 100 of 124 total aligned reads displayed.

TTTACCGGTACATCGCCCGTGGGTTGCTTCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGCACGTGG  >  NC_002947/3791305‑3791602
                                                                                                                                                    |                                                                                                                                                     
tttACCGGTACATCGCCCGTGGGTTGCTTCCCGGCCAACGGCTTTGGTCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAggcccggcc                                                                                                                                                      <  1:2127583/150‑1 (MQ=255)
    ccGGTACATCGCCCGTGGGTTGCTTCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGa                                                                                                                                                  <  1:980128/150‑1 (MQ=255)
        tACATCGCCCGTGAGTTGCTACCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGa                                                                                                                                              <  2:1695442/150‑1 (MQ=255)
           aTCGCCCGTGGGTTGCTTCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAg                                                                                                                                           >  1:1458098/1‑150 (MQ=255)
            tCGCCCGTGGGTTGCTTCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGGCCGGGCCTGAgg                                                                                                                                                >  2:769391/1‑143 (MQ=255)
            tCGCCCGTGGGTTGCTTCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGGCCGGGCCTGAgg                                                                                                                                                <  1:769391/144‑2 (MQ=255)
            tCGCCCGTGGGTTGCTTCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCACTCAGGGCCGGGCCTGAggg                                                                                                                                               >  2:396863/1‑143 (MQ=255)
            tCGCCCGTGGGTTGCTTCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCACTCAGGGCCGGGCCTGAggg                                                                                                                                               <  1:396863/145‑3 (MQ=255)
                ccGTGGGTTGCTTCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGt                                                                                                                                                <  1:1231791/140‑1 (MQ=255)
                ccGTGGGTTGCTTCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGt                                                                                                                                                >  2:1231791/1‑140 (MQ=255)
                 cGTGGGTTGCTTCCCGGCAAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGCTCCGACCCTGAGTGTCAGCTTCGa                                                                                                                                     <  2:2577205/150‑1 (MQ=255)
                   tGGGTTGCTTCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGAcc                                                                                                                                   >  2:2545772/1‑150 (MQ=255)
                    gggTTGCTTCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGAccc                                                                                                                                  <  1:2979414/150‑1 (MQ=255)
                      gTTGCTTCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCgg                                                                                                                                >  2:2209135/1‑150 (MQ=255)
                      gTTGCTTCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCgg                                                                                                                                >  1:1730999/1‑150 (MQ=255)
                           ttCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGGCAGCTTCGACCCGGCAGAc                                                                                                                           >  2:984161/1‑150 (MQ=255)
                           ttCCCGGCCAACGGCTTTGGCCTGGTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCACGCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCCGAc                                                                                                                           >  1:941955/1‑150 (MQ=255)
                            tCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGAcc                                                                                                                          >  1:459381/1‑150 (MQ=255)
                              ccGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGAcccc                                                                                                                        >  1:274267/1‑150 (MQ=255)
                                  ccgaCGGCTTTGGCCGGTTCGACGCCGGCGGCGATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCGCTGAGTGACAGTTTCGGCCCGGCACACCCCggtg                                                                                                                    <  2:606438/147‑3 (MQ=255)
                                   cAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGt                                                                                                                   <  2:274267/150‑1 (MQ=255)
                                              gCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCgg                                                                                                        >  1:2873622/1‑150 (MQ=255)
                                               ccTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTTCCGGGGGCggg                                                                                                       >  2:2750803/1‑149 (MQ=255)
                                               ccTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGt                                                                                                       >  2:890047/1‑150 (MQ=255)
                                                cTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGa                                                                                                                                                  <  2:91739/106‑1 (MQ=255)
                                                cTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTa                                                                                                      >  2:1843391/1‑150 (MQ=255)
                                                 tGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTaa                                                                                                     <  2:194769/150‑1 (MQ=255)
                                                 tGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTaa                                                                                                     >  2:351481/1‑150 (MQ=255)
                                                 tGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTaa                                                                                                     <  1:2750803/150‑1 (MQ=255)
                                                   ttCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATc                                                                                                   >  1:1802012/1‑150 (MQ=255)
                                                   ttCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCAAAGCGCGACGCCAAGCTCGACCCCTCAGCCCCGGCCCAGAGTGACAGCTTCGACCCTGCAGACCCCCGCGTGCCGGTGGCGGTAATc                                                                                                   >  2:2490485/1‑150 (MQ=255)
                                                      gATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCaaa                                                                                                >  1:2008735/1‑150 (MQ=255)
                                                      gATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCaaa                                                                                                <  1:984161/150‑1 (MQ=255)
                                                       aTGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAg                                                                                               <  1:351481/150‑1 (MQ=255)
                                                        tactggcggcACTGTCTGGGAGCTGACTCGCACGCGCTATCGGCCAGGCCATGGCGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTCACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAgg                                                                                              <  2:1318047/146‑1 (MQ=255)
                                                        tGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAgg                                                                                              >  2:1234341/1‑150 (MQ=255)
                                                        tGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAgg                                                                                              >  1:264768/1‑150 (MQ=255)
                                                           cggcggCAATGTCTGGGAGCTTACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAggcgg                                                                                           <  2:419008/150‑1 (MQ=255)
                                                           cggcggCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAggcgg                                                                                           <  1:2209135/150‑1 (MQ=255)
                                                            ggcggcAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAggcggc                                                                                          >  1:2400759/1‑150 (MQ=255)
                                                              cggcAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTc                                                                                        <  1:1526004/150‑1 (MQ=255)
                                                               ggcAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCg                                                                                       >  2:1923748/1‑150 (MQ=255)
                                                                    tGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCt                                                                                  >  2:1378600/1‑150 (MQ=255)
                                                                      tCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGGGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCtgt                                                                                <  2:2783639/150‑1 (MQ=255)
                                                                         gggAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGtt                                                                             >  2:1843969/1‑150 (MQ=255)
                                                                             gCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGt                                                                                                             >  2:2902796/1‑114 (MQ=255)
                                                                             gCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTcggc                                                                         <  2:2873622/150‑1 (MQ=255)
                                                                                   tCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCg                                                                   >  2:1091604/1‑150 (MQ=255)
                                                                                     gCACGGGCTATCGGCCAGGTCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCt                                                                 >  1:1673152/1‑150 (MQ=255)
                                                                                     gCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCt                                                                 <  1:1843391/150‑1 (MQ=255)
                                                                                      cACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTg                                                                >  2:1826122/1‑150 (MQ=255)
                                                                                       aCGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACa                                                                                                                                            >  2:1393298/1‑73 (MQ=255)
                                                                                       aCGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACa                                                                                                                                            <  1:1393298/73‑1 (MQ=255)
                                                                                                cGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTAc                                                      <  1:2328225/150‑1 (MQ=255)
                                                                                                  gccaggccaTGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTaa                                                                                                     <  1:2865253/101‑1 (MQ=255)
                                                                                                  gccaggccaTGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCg                                                    >  1:1954833/1‑150 (MQ=255)
                                                                                                   ccaggccaTGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGc                                                   <  2:1912484/150‑1 (MQ=255)
                                                                                                    caggccaTGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGCTGGCGGCAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGcc                                                  >  2:1969407/1‑150 (MQ=255)
                                                                                                       gccaTGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCt                                               <  1:2183204/150‑1 (MQ=255)
                                                                                                       gccaTGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCt                                               >  2:1751488/1‑150 (MQ=255)
                                                                                                        ccaTGACGCAGAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGCAATGCCGGTGGCGGTAATCAAAGGCGGCCAGCACCTGTCTTCGGCGGACCGCTGCATGCGCTACCGCCCCTc                                              <  1:1751488/150‑1 (MQ=255)
                                                                                                          aTGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCCAAGGCCGCTCCCACCTGTGTTCGGCGGACCGCTACATGCGGCACCGCCCCTCgg                                            >  2:915742/1‑150 (MQ=255)
                                                                                                          aTGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCgg                                            <  1:448250/150‑1 (MQ=255)
                                                                                                            gACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCa                                          <  2:1458098/150‑1 (MQ=255)
                                                                                                             aCGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCAc                                         <  1:890047/150‑1 (MQ=255)
                                                                                                                 acaGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACa                                                                                                                                            <  2:1491291/47‑1 (MQ=255)
                                                                                                                 acaGCGCGACGCCAAGCTCGACCCCTCAGGCCAGGCCCTGAGTGACa                                                                                                                                            >  1:1491291/1‑47 (MQ=255)
                                                                                                                      gcgACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGc                                       <  1:853111/143‑1 (MQ=255)
                                                                                                                           gCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCcagccgcagccg                           >  1:2713816/1‑150 (MQ=255)
                                                                                                                             cAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCg                                                                                       <  2:570706/88‑1 (MQ=255)
                                                                                                                             cAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGgt                         >  2:2458538/1‑150 (MQ=255)
                                                                                                                             cAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGgt                         >  2:1535502/1‑150 (MQ=255)
                                                                                                                              aaGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTcg                                                                           >  1:2704389/1‑99 (MQ=255)
                                                                                                                              aaGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTcg                                                                           <  2:2704389/99‑1 (MQ=255)
                                                                                                                              aaGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATgcgc                                                         <  1:1494941/117‑1 (MQ=255)
                                                                                                                              aaGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATgcgc                                                         >  2:1494941/1‑117 (MQ=255)
                                                                                                                                 cTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCGTGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTc                     >  2:2098318/1‑150 (MQ=255)
                                                                                                                                 cTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTc                     <  1:2098318/150‑1 (MQ=255)
                                                                                                                                 cTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTc                     >  1:1171880/1‑150 (MQ=255)
                                                                                                                                 cTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTc                     <  2:1673152/150‑1 (MQ=255)
                                                                                                                                  tCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGCCCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTCCCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCa                    <  1:1151754/150‑1 (MQ=255)
                                                                                                                                  tCGACCCCTCAGGCCCGGCACTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCa                    >  1:52492/1‑150 (MQ=255)
                                                                                                                                   cGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCAt                   <  1:2581239/150‑1 (MQ=255)
                                                                                                                                    gACCACTCAGGCCCGGCCCTGAAGTTCAGCGGCGACCCTGCAGACCCCTTCGTGCCGGTGGCGGTAATCAAAGGAGGCTCGCACCTGTGTTTAGCCGAGCGCTCCAGGCCCGACCTCGCCCCTGCACGCCAGCCGCAGCCGGTGTTCATg                  >  2:2200817/1‑150 (MQ=255)
                                                                                                                                     cccccTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATga                 <  1:1535502/149‑1 (MQ=255)
                                                                                                                                     aCCCCTCAGGCCCGGCACTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATga                 <  2:52492/150‑1 (MQ=255)
                                                                                                                                        ccTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATgacga              >  2:2898770/1‑150 (MQ=255)
                                                                                                                                         ctcgggcccggcccTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATgacgac             <  1:2587285/146‑1 (MQ=255)
                                                                                                                                         cTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTc                                              >  1:2157989/1‑117 (MQ=255)
                                                                                                                                         cTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTc                                              <  2:2157989/117‑1 (MQ=255)
                                                                                                                                         cTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATgacgac             >  2:2951623/1‑150 (MQ=255)
                                                                                                                                         cTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATgacgac             >  1:1914501/1‑150 (MQ=255)
                                                                                                                                            aGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATGACGACCTc          <  2:1730999/150‑1 (MQ=255)
                                                                                                                                                ccggcccTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCa                                                              >  2:2439244/1‑94 (MQ=255)
                                                                                                                                                ccggcccTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATga                 >  2:1954793/1‑139 (MQ=255)
                                                                                                                                                ccggcccTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATga                 <  1:1954793/139‑1 (MQ=255)
                                                                                                                                                ccggcccTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGCAc      >  1:792560/1‑150 (MQ=255)
                                                                                                                                                    cccTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGCACGTgg  <  1:452500/150‑1 (MQ=255)
                                                                                                                                                    cccTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGCACGTgg  <  1:296021/150‑1 (MQ=255)
                                                                                                                                                    cccTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGCACGTgg  <  2:423651/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TTTACCGGTACATCGCCCGTGGGTTGCTTCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGCACGTGG  >  NC_002947/3791305‑3791602

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: ‑