Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,255,625 T→C 5.4% L52L (CTA→CTG PP_2852 ← sulfatase domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,255,6250TC 5.4% 224.7 / 2.1 74L52L (CTA→CTGPP_2852sulfatase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (27/43);  new base C (2/2);  total (29/45)
Fisher's exact test for biased strand distribution p-value = 6.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.40e-01

GAATTCTCCATGAATTGCGGGTCGGTCAGGTAGTTAAGGTCGGCGCTGGTGGGCAGTCGGCCCACTGCGCTTACCAGTTCGGCCGAGGCCACCCACAGCGCCGCCCAGGCCAGCAGCACGGGCAGCGCCAGCCACCACGGCCTGCGGTGTAGCAGCACTATCAACAGGCTGCCGAGGGCGAG‑‑ATCGGACAGGTAACCGAACGGGTTGGACCAACCCAGCATGGCGCGCGTACCCAGTGGGATCACCAGCACCAGAGCCGCCAAAGTGGCAAGCCGGGTGCTTGGGTGGTCGGACAGG  >  NC_002947/3255476‑3255772
                                                                                                                                                     |                                                                                                                                                     
gAATTCTCCATGAATTGCGGGTCGGTCAGGTAGTTAAGGTGTGCGCTGGTGGGCAGTCGGCCCACTGCGCTTACCAGTTCGGCCGAGGCCACCCACAGCGCCGCCCAGGCCAGCAGCACGGGCAGTGCCAGCCACCACGGCCTGCGgtgt                                                                                                                                                       <  1:1261244/150‑1 (MQ=255)
gAATTCTCCATGAATTGCGGGTCGGTCAGGTAGTTAAGGTCGGCGCTGGTGGGCAGTCGGCCCACTGCGCTTACCAGTTCGGCCGAGGCCACCCACAGCGCCGCCCAGGCCAGCAGCACGGGCAGCGCCAGCCACCACGGCCTGCGgtgt                                                                                                                                                       <  1:829529/150‑1 (MQ=255)
  aTTCGCCATGAATTGCGGGTCGGTCAGGTAGTTAAGGTCGGCGCTGGTGGGCAGTCGGCCCACTGCGCTTACCAGTTCGGCCGAGGCCACCCATAGCGCCGCCCAGGCCAGCAGCACGGGCAGCGCCAGCCACCACGGCCTGCGGTGTag                                                                                                                                                     <  2:923075/150‑1 (MQ=255)
    tctcCATGAATTGCGGGTCGGTCAGGTAGTTAAGGTCGGCGCTGGTGGGCAGTCGGCCCACTGCGCTTACCAGTTCGGCCGAGGCCACCCACAGCGCCGCCCAGGCCAGCAGCACGGGCAGCGCCAGCCACCACGGCCTGCGGTGTagca                                                                                                                                                   <  2:766094/150‑1 (MQ=255)
    tctcCATGAATTCCCGGTCAACCAGGTAGTTAAGGTCGGCGCTGGTGGGCAGTCGGCCCACTGCGCTTACCAGTTCGGCCGAGGCCACCCACAGCGCCGCCCAGGCCAGCAGCACGGGCAGCGCCAGCCACCACGGCCTGCGGTGTagca                                                                                                                                                   >  2:464952/1‑150 (MQ=255)
          tGAATTGCGGGTCGGTCAGGTAGTTAAGGTCGGCGCTGGTGGGCAGTCGGCCCACTGCGCTTACCAGTTCGGCCGAGGCCACCCACAGCGCCGCCCAGGCCAGCAGCACGGGCAGCGCCAGCCTCCACGGCCTGCGGTGTAGCAGCACTa                                                                                                                                             >  1:1232840/1‑150 (MQ=255)
           gAATTGCGGGTCGGTCAGGTAGTTAAGGTCGGCGCTGGTGGGCAGTCGGCCCACTGCGCTTACCAGTTCGGCCGAGGCCACCCACAGCGCCGCCCAGGCCAGCAGCACGGGCAGCGCCAGCCACCACGGCCTGCGGTGTAGCAGCACTAt                                                                                                                                            <  2:1834408/150‑1 (MQ=255)
           gAATTGCGGGTCGGTCAGGTAGTTAAGGTCGGCGCTGGTGGGCAGTCGGCCCACTGCGCTTACCAGTTCGGCCGAGGCCACCCACAGCGCCGCCCAGGCCAGCAGCACGGGCAGCGCCAGCCACCACGGCCTGCGGTGTAGCAGCACTAt                                                                                                                                            <  1:145647/150‑1 (MQ=255)
              ttGCGGGTCGGTCAGGTAGTTAAGGTCGGCGCTGGTGGGCAGTCGGCCCACTGCGCTTACCAGTTCGGCCGAGGCCACCCACAGCGCCGCCCAGGCCAGCAGCACGGGCAGCGCCAGCCACCACGGCCTGCGGTGTAGCAGCACTATCaa                                                                                                                                         <  2:758721/150‑1 (MQ=255)
               tGCGGGTCGGTCAGGTAGTTAAGGTCGGCGCTGGTGGGCAGTCGGCCCACTGCGCTTACCAGTTCGGCCGAGGCCACCCACAGCGCCGCCCAGGCCAGCAGCACGGGCAGCGCCAGCCACCACGGCCTGCGGTGTAGCAGCACTATCAAc                                                                                                                                        <  1:1497316/150‑1 (MQ=255)
               tGCGGGTCGGTCAGGTAGTTAAGGTCGGCGCTGGTGGGCAGTCGGCCCACTGCGCTTACCAGTTCGGCCGAGGCCACCCACAGCGCCGCCCAGGCCAGCAGCACGGGCAGCGCCAGCCACCACGGCCTGCGGTGTAGCAGCACTATCAAc                                                                                                                                        <  2:395609/150‑1 (MQ=255)
                    gtcggtcAGGTAGTTAAGGTCGGCGCTGGTGGGCAGTCGGCCCACTGCGCTTACCAGTTCGGCCGAGGCCACCCACAGCGCCGCCCAGGCCAGCAGCACGGGCAGCGCCAGCCACCACGGCCTGCGGTGTAGCAGCACTATCAACAGGCt                                                                                                                                   >  2:328111/1‑150 (MQ=255)
                           aGGTAGTTAAGGTCGGCGCTGGTGGGCAGTCGGCCCACTGCGCTTCCCAGTTCGGCCGAGGCCACCCACAGCGCCGCCCAGGCCAGCAGCACGGGCAGCGCCAGCCACCACGGCCTGCGGTGTAGCAGCACTATCAACAGGCTGCCGAgg                                                                                                                            <  2:1789833/150‑1 (MQ=255)
                           aGGTAGTTAAGGTCGGCGCTGGTGGGCAGTCGGCCCACTGCGCTTACCAGTTCGGCCGAGGCCACCCACAGCGCCGCCCAGGCCAGCAGCACGGGCAGCGCCAGCCACCACGGCCTGCGGTGTAGCAGCACTATCAACAGGCTGCCGAgg                                                                                                                            <  2:853784/150‑1 (MQ=255)
                               aGTTAAGGTCGGCGCTGGTGGGCAGTCGGCCCACTGCGCTTACCAGTTCGGCCGAGGCCACCCACAGCGCCGCCCAGGCCAGCAGCACGGGCAGCGCCAGCCACCACGGCCTGCGGTGTAGCAGCACTATCAACAGGCTGCCGAGGGCga                                                                                                                        <  2:924175/150‑1 (MQ=255)
                                      gTCGGCGCTGGTGGGCAGTCGGCCCACTGCGCTTACCAGTTCGGCCGAGGCCACCCACAGCGCCGCCCAGGCCAGCAGCACGGGCAGCGCCAGCCACCACGGCCTGCGGTGTAGCAGCACTATCAACAGGCTGCCGAGGGCGAG‑‑ATCGGa                                                                                                               >  2:500722/1‑150 (MQ=255)
                                        cGGCGCTGGTGGGCAGTCGGCCCACTGCGCTTACCAGTTCGGCCGAGATTACCCACAGCGCCGCCCAGGCCAGCAGCACGGGCAGCGCCAGCCACCACGGCCTGCGGTGTAGCAGCACTATCAACAGGCTGCCGAGGGCGAG‑‑ATCGGACa                                                                                                             <  1:592362/150‑1 (MQ=255)
                                         ggCGCTGGTGGGCAGTCGGCCCACTGCGCTTACCAGTTCGGCCGAGGCCACCCACAGCGCCGCCCAGGCCAGCAGCACGGGCAGCGCCAGCCACCACGGCCTGCGGTGTAGCAGCACTATCAACAGGCTGCCGAGGGCGAG‑‑ATCGGACAg                                                                                                            >  1:817312/1‑150 (MQ=255)
                                              tggtggGCAGTCGGCCCACTGCGCTTACCAGTTCGGCCGAGGCCACCCACAGCGCCGCCCAGGCCAGCAGCACGGGCAGCGCCAGCCACCACGGCCTGCGGTGTAGCAGCACTATCAACAGGCTGCCGAGGGCGAG‑‑ATCGGACAGGTAAc                                                                                                       >  2:1297153/1‑150 (MQ=255)
                                              tggtggGCAGTCGGCCCACTGCGCTTACCAGTTCGGCCGAGGCCACCCACAGCGCCGCCCAGGCCAGCAGCACGGGCAGCGCCAGCCACCACGGCCTGCGGTGTAGCAGCACTATCAACAGGCTGCCGAGGGCGAG‑‑ATCGGACAGGTAAc                                                                                                       <  1:1137341/150‑1 (MQ=255)
                                               ggtggGCAGTCGGCCCACTGCGCTTACCAGTTCGGCCGAGGCCACCCACAGCGCCGCCCAGGCCAGCAGCACGGGCAGCGCCAGCCACCACGGCCTGCGGTGTAGCAGCACTATCAACAGGCTTCGTAGGGCGAG‑‑ATCGGACAGGTAAcc                                                                                                      <  1:258801/150‑1 (MQ=255)
                                                     cAGTCGGCCCACTGCGCTTACCAGTTCGGCCGAGGCCACCCACAGCGCCGCCCAGGCCAGCAGCACGGGCAGCGCCAGCCACCACGGCCTGCGGTGTAGCAGCACTATCAACAGGCTGCCGAGGGCGAG‑‑ATCGGACAGGTAACCGAACgg                                                                                                <  1:339193/150‑1 (MQ=255)
                                                         cgatcctcTGCGCGTACCAGTTCGGCCGAGGCCACCCACAGCGCCGCCCAGGCCAGCAGCACGGGCAGCGCCAGCCACCACGGCCTGCGGTGTAGCAGCACTATCAACAGGCTGCCGAGGGCGAG‑‑ATCGGACAGGTAACCGAACGGGTTg                                                                                            <  1:1461972/143‑1 (MQ=255)
                                                         cGGCCCACTGCGCTTACCAGTTCGGCCGAGGCCACCCACAGCGCCGCCCAGGCCAGCAGCACGGGCAGCGCCAGCCACCACGGCCTGCGGTGTAGCAGCACTATCAACAGGCTGCCGAGGGCGAG‑‑ATCGGACAGGTAACCGAACGGGTTg                                                                                            <  1:201237/150‑1 (MQ=255)
                                                         cGGCCCACTGCGCTTACCAGTTCGGCCGAGGCCACCCACAGCGCCGCCCAGGCCAGCAGCACAGGCAGCGCCAGCCACCACGGCCTGCGGTGTAGCAGCACTATCAACAGGCTGCCGAGGGCGAG‑‑ATCGGACAGGTAACCGAACGGGTTg                                                                                            <  2:824263/150‑1 (MQ=255)
                                                         cGGCCCACTGCGCTTACCAGTTCGGCCGAGGCCACCCACAGCGCCGCCCAGGCACGCCGCGAGCAGTTCGCCAGCCACCACGGCCTGGGGTGTAGCAGCACTATCAACAGGCTGCCGAGGGCGAG‑‑ATCGGACAGGTAACCGAACGGGTTg                                                                                            <  1:1875681/150‑1 (MQ=255)
                                                           gCCCACTGCGCTTACCAGTTCGGCCGAGGCCACCCACAGCGCCGCCCAGGCCAGCAGCACGGGCAGCGCCAGCCACCACGGCCTGCGGTGTAGCAGCACTGTCAACAGGCTGCCGAGGGCGAG‑‑ATCGGACAGGTAACCGAACGGGTTGGa                                                                                          <  1:189278/150‑1 (MQ=255)
                                                           gCCCACTGCGCTTACCAGTTCGGCCGAGGCCACCCACAGCGCCGCCCAGGCCAGCAGCACGGGCAGCGCCAGCCACCACGGCCTGCGGTGTAGCAGCACTATCAACAGGCTGCCGAGGGCGAG‑‑ATCGGACAGGTAACCGAACGGGTTGGa                                                                                          <  1:678310/150‑1 (MQ=255)
                                                           gCCCACTGCGCTTACCAGTTCGGCCGAGGCCACCCACAGCGCCGCCCAGGCCAGCAGCACGGGCAGCGCCAGCCACCACGGCCTGCGGTGTAGCAGCACTATCAACAGGCTGCCGAGGGCGAG‑‑ATCGGACAGGTAACCGAACGGGTTGGa                                                                                          <  2:623924/150‑1 (MQ=255)
                                                           gCCCACTGCGCTTACCAGTTCGGCCGAGGCCACCCACAGCGCCGCATCAATCAATAGCACGGGCAGCGCCAGCCACCACGGCCTGCGGTGTAGCAGCACTATCAACAGGCTGCCGAGGGCGAG‑‑ATCGGACAGGTAACCGAACGGGTTGGa                                                                                          <  1:1058304/150‑1 (MQ=255)
                                                                  gcgcTTACCAGTTCGGCCGAGGCCACCCACAGCGCCGCCCAGGCCAGCAGCACGGGCAGCGCCAGCCACCACGGCCTGCGGTGTAGCAGCACTATCAACAGGCTGCCGAGGGCGAG‑‑ATCGGACAGGTAACCGGACGGGTTGGACCAAccc                                                                                   <  2:221674/150‑1 (MQ=255)
                                                                   cgcTTACCAGTTCGGCCGAGGCCACCCACAGCGCCGCCCAGGCCAGCAGCACGGGCAGCGCCAGCCACCACGGCCTGCGGTGTAGCAGCACTATCAACAGGCTGCCGAGGGCGAG‑‑ATc                                                                                                                  <  2:1372067/118‑1 (MQ=255)
                                                                   cgcTTACCAGTTCGGCCGAGGCCACCCACAGCGCCGCCCAGGCCAGCAGCACGGGCAGCGCCAGCCACCACGGCCTGCGGTGTAGCAGCACTATCAACAGGCTGCCGAGGGCGAG‑‑ATc                                                                                                                  >  1:1372067/1‑118 (MQ=255)
                                                                   cgcTTACCAGTTCGGCCGAGGCCACCCACAGCGCCGCCCAGGCCAGCAGCACGGGCAGCGCCAGCCACCACGGCCTGCGGTGTAGCAGCACTATCAACAGGCTGCCGAGGGCGAG‑‑ATCGGACAGGTAACCGAACGGGTTGGACCAACCCa                                                                                  <  2:745876/150‑1 (MQ=255)
                                                                   cgcTTACCAGTTCGGCCGAGGCCACCCACAGCGCCGCCCAGGCCAGCAGCACGGGCAGCGCCAGCCACCACGGCCTGCGGTGTAGCAGCACTATCAACAGGCTGCCGAGGGCGAG‑‑ATCGGACAGGTAACCGAACGGGTTGGACCAACCCa                                                                                  <  2:93168/150‑1 (MQ=255)
                                                                    gcTTACCAGTTCGGCCGAGGCCACCCACAGCGCCGCCCAGGCCAGCAGCACGGGCAGCGCCAGCCACCACGGCCTGCGGTGTAGCAGCACTATCAACAGGCTGCCGAGGGCGAG‑‑ATCGGACAGGTAACCGAACGGGTTGGACCAACCCAg                                                                                 >  1:877597/1‑150 (MQ=255)
                                                                        aCCAGTTCGGCCGAGGCCACCCACAGCGCCGCCCAGGCCAGCAGCACGGGCAGCGCCAGCCACCACGGCCTGCGGTGTAGCAGCACTATCAACAGGCTGCCGAGGGCGAG‑‑ATCGGACAGGTAACCGAACGGGTTGGACCAACCCAGCATg                                                                             <  1:328111/150‑1 (MQ=255)
                                                                        aCCAGTTCGGCCGAGGCCACCCACAGCGCCGCCCAGGCCAGCAGCACGGGCAGCGCCAGCCACCACGGCCTGCGGTGTAGCAGCACTATCAACAGGCTGCCGAGGGCGAG‑‑ATCGGACAGGTAACCGAACGGGTTGGACCAACCCAGCATg                                                                             >  1:1478162/1‑150 (MQ=255)
                                                                            gTTCGGCCGAGGCCACCCACAGCGCCGCCCAGGCCAGCAGCACGGGCAGCGCCAGCCACCACGGCCTGCGGTGTAGCAGCACTATCAACAGGCTGCCGAGGGCGAG‑‑ATCGGACAGGTAACCGAACGGGTTGGACCAACCCAGCATGgcgc                                                                         >  1:842715/1‑150 (MQ=255)
                                                                             ttCGGCCGAGGCCACCCACAGCGCCGCCCAGGCCAGCAGCACGGGCAGCGCCAGCCACCACGGCCTGCGGTGTAGCAGCACTATCAACAGGCTGCCGAGGGCGAG‑‑ATCGGACAGGTAACCGAACGGGTTGGACCAACCCAGCATGgcgcg                                                                        >  1:1520481/1‑150 (MQ=255)
                                                                                     aGGCCACCCACCGCGCCGCCCATGCCAGCAGCACGGGCAGCGCCAGCCACCACGGCCTGCGGTGTAGCAGCACTATCAACAGGCTGCCGAGGGCGAG‑‑ATCGGACAGGTAACCGAACGGGTTGGACCAACCCAGCATGGCGCGCGTACCCa                                                                <  2:1295555/150‑1 (MQ=255)
                                                                                      ggCCACCCACAGCGCCGCCCAGGCCAGCAGCACGGGCAGCGCCAGCCACCACGGCCTGCGGTGTAGCAGCACTATCAACAGGCTGCCGAGGGCGAG‑‑ATCGGACAGGTAACCGAACGGGTTGGACCAACCCAGCATGGCGCGCGTACCCAg                                                               <  1:1297153/150‑1 (MQ=255)
                                                                                       gCCACCCACAGCGCCGCCCAGGCCAGCAGCACGGGCAGCGCCAGCCACCACGGCCTGCGGTGTAGCAGCACTATCAACAGGCTGCCGAGGGCGAG‑‑ATCGGACAGGTAACCGAACGGGTTGGACCAACCCAGCATGGCGCGCGTACCCAGt                                                              >  2:1490961/1‑150 (MQ=255)
                                                                                       gCCACCCACAGCGCCGCCCAGGCCAGCAGCACGGGCAGCGCCAGCCACCACGGCCTGCGGTGTAGCAGCACTATCAACAGGCTGCCGAGGGCGAG‑‑ATCGGACAGGTAACCGAACGCGTTGGACCAACCCAGCATGGCGCGCGTACCCAGt                                                              <  2:817312/150‑1 (MQ=255)
                                                                                        ccacccacAGCGCCGCCCAGGCCAGCAGCACGGGCAGCGCCAGCCACCACGGCCTGCGGTGTAGCAGCACTATCAACAGGCTGCCGAGGGCGAG‑‑ATCGGACAGGTAACCGAACGGGTTGGACCAACCCAGCATGGCGCGCGTACCCAGTg                                                             <  2:1478162/150‑1 (MQ=255)
                                                                                          acccacAGCGCCGCCCAGGCCAGCAGCACGGGCAGCGCCAGCCACCACGGCCTGCGGTGTAGCAGCACTATCAACAGGCTGCCGAGGGCGAG‑‑ATCGGACAGGTAACCGAAAGGGTTGGACCAACCCAGCATGGCGCGCGTACCCAGTggg                                                           >  1:1147935/1‑150 (MQ=255)
                                                                                            ccacAGCGCCGCCCAGGCCAGCAGCACGGGCAGCGCCAGCCACCACGGCCTGCGGTGTAGCAGCACTATCAACAGGCTGCCGAGGGCGAGCTAGCCGACAG‑‑CGCGTAAGCTGTTGGACCAACCCAGCATGGCGCGCGTACCCAGTGGGAt                                                         >  2:527624/1‑150 (MQ=255)
                                                                                            ccacAGCGCCGCCCAGGCCAGCAGCACGGGCAGCGCCAGCCACCACGGCCTGCGGTGTAGCAGCACTATCAACAGGCTGCCGAGGGCGAG‑‑ATCGGACAGGTAACCGAACGGGTTGGACCAACCCAGCATgg                                                                            >  2:340404/1‑131 (MQ=255)
                                                                                            ccacAGCGCCGCCCAGGCCAGCAGCACGGGCAGCGCCAGCCACCACGGCCTGCGGTGTAGCAGCACTATCAACAGGCTGCCGAGGGCGAG‑‑ATCGGACAGGTAACCGAACGGGTTGGACCAACCCAGCATgg                                                                            <  1:340404/131‑1 (MQ=255)
                                                                                            ccacAGCGCCGCCCAGGCCAGCAGCACGGGCAGCGCCAGCCACCACGGCCTGCGGTGTAGCAGCACTATCAACAGGCTGCCGAGGGCGAG‑‑ATCGGACAGGTAACCGAACGGGTTGGACCAACCCAGCATGGCGCGCGTACCCAGTGGGAt                                                         >  1:907765/1‑150 (MQ=255)
                                                                                            ccacAGCGCCGCCCAGGCCAGCAGCACGGGCAGCGCCAGCCACCACGGCCTGCGGTGTAGCAGCACTATCAACAGGCTGCCGAGGGCGAG‑‑ATCGGACAGGTAACCGAACGGGTTGGACCAACCCAGCATGGCGCGCGTACCCAGTGGGAt                                                         >  2:723919/1‑150 (MQ=255)
                                                                                            ccacAGCGCCGCCCAGGCCAGCAGCACGGGCAGCGCCAGCCACCACGGCCTGCGGTGTAGCAGCACTATCAACAGGCTGCCGAGGGCGAG‑‑ATCGGACAGGTAACCGAACGGGTTGGACCAACCCAGCATGGCGCGCGTACCCAGTGGGAt                                                         >  2:1226690/1‑150 (MQ=255)
                                                                                              acaGCGCCGCCCAGGCCAGCAGCACGGGCAGCGCCAGCCACCACGGCCTGCGGTGTAGCAGCACTATCAACAGGCTGCCGAGGGCGAG‑‑ATCGGACAGGTAACCGAACGGGTTGGACCAACCCAGCATGGCGCGCGTACCCAGTGGGATca                                                       >  2:363599/1‑150 (MQ=255)
                                                                                              acaGCGCCGCCCAGGCCAGCAGCACGGGCAGCGCCAGCCACCACGGCCTGCGGTGTAGCAGCACTATCAACAGGCTACCGAGGGCGAG‑‑ATCGGACAGGTAACCGAACGGGTTGGACCAACCCAGCATGGCGCGCGTACCCAGTGGGATca                                                       >  1:852836/1‑150 (MQ=255)
                                                                                              acaGCGCCGCCCAGGCCAGCAGCACGGGCAGCGCCAGCCACCACGGCCTGCGGTGCAGCAGCACTATc                                                                                                                                           >  1:341020/1‑68 (MQ=255)
                                                                                              acaGCGCCGCCCAGGCCAGCAGCACGGGCAGCGCCAGCCACCACGGCCTGCGGTGCAGCAGCACTATc                                                                                                                                           <  2:1284558/68‑1 (MQ=255)
                                                                                              acaGCGCCGCCCAGGCCAGCAGCACGGGCAGCGCCAGCCACCACGGCCTGCGGTGCAGCAGCACTATc                                                                                                                                           <  2:341020/68‑1 (MQ=255)
                                                                                              acaGCGCCGCCCAGGCCAGCAGCACGGGCAGCGCCAGCCACCACGGCCTGCGGTGCAGCAGCACTATc                                                                                                                                           >  1:1284558/1‑68 (MQ=255)
                                                                                                         cAGGCCAGCAGCACGGGCAGCGCCAGCCACCACGGCCTGCGGTGTAGCAGCACTATCAACAGGCTGCCGAGGGCGAG‑‑ATCGGACAGGTAACCGAACGGGTTGGACCAACCCAGCATGGCGCGCGTACCCAGTGGGATCACCAGCACCaga                                            <  2:877597/150‑1 (MQ=255)
                                                                                                         cAGGCCAGCAGCACGGGCAGCGCCAGCCACCACGGCCTGCATTTTAGCAGCACTATCAACAGGCTGCCGAGGGCGAG‑‑ATCGGACAGGTAACCGAACGGGTTGGACCAACCCAGCATGGCGCGCGTACCCAGTGGGATCACCAGCACCaga                                            <  1:497044/150‑1 (MQ=255)
                                                                                                          aGGCCAGCAGCACGGGCAGCGCCAGCCACCACGGCCTGCGGTGTAGCAGCACTATCAACAGGCTGCCGAGGGCGAG‑‑ATCGGACAGGTAACCGAACGGGTTGGACCAACCCAGCATGGCGCGCGTACCCAGTGGGATCACCAGCACCagag                                           >  2:644405/1‑150 (MQ=255)
                                                                                                              cagcagCACGGGCAGCGCCAGCCACCACGGCCTGCGGTGTAGCAGCACTATCAACAGGCTGCCGAGGGCGAG‑‑ATCGGACAGGTAACCGAACGGGTTGGACCAACCCAGCATGGCGCGCGTACCCAGTGGGATCACCAGCACCAGAgccgc                                       <  2:1520481/150‑1 (MQ=255)
                                                                                                                     aCGGGCAGCGCCAGCCACCACGGCCTGCGGTGTAGCAGCACTATCAACAGGCTGCCGAGGGCGAG‑‑ATCGGACAGGTAACCGAACGGGTTGGACCAACCCAGCATGGCGCGCGTACCCAGTGGGATCACCAGCACCAGAGCCGCCAAAGTg                                >  2:1593001/1‑150 (MQ=255)
                                                                                                                              gccagccaCCACGGCCTGCGGTGTAGCAGCACTATCAACAGGCTGCCGAGGGCGAG‑‑ATCGGACAGGTAACCGAACGGGTTGGACCAACCCAGCATGGCGCGCGTACGCAGTGGGATCACCAGCACCAGAGCCGCCAAAGTGGCAAGCCgg                       >  1:1023870/1‑150 (MQ=255)
                                                                                                                              gccagccaCCACGGCCTGCGGTGTAGCAGCACTATCAACAGGCTGCCGAGGGCGAG‑‑ATCGGACAGGTAACCGAACGGGTTGGACCAACCCAGCATGGCGCGCGTACCCAGTGGGATCACCAGCACCAGAGCCGCCAAAGTGGCAAGCCgg                       <  1:1226690/150‑1 (MQ=255)
                                                                                                                              gccagccaCCACGGCCTGCGGTGTAGCAGCACTATCAACAGGCTGCCGAGGGCGAG‑‑ATCGGACAGGTAACCGAACGGGTTGGACCAACCCAGCATGGCGCGCGTACCCAGTGGGATCACCAGCACCAGAGCCGCCAAAGTGGCAAGCCgg                       <  1:1850825/150‑1 (MQ=255)
                                                                                                                               ccagccaCCACGGCCTGCGGTGTAGCAGCACTATCAACAGGCTGCCGAGGGCGAG‑‑ATCGGACAGGTAACCGAACGGGTTGGACCAACCCAGCATGGCGCGCGTACCCAGTGGGATCACCAGCACCAGAGCCGCCAAAGTGGCAAGCCggg                      >  1:1591547/1‑150 (MQ=255)
                                                                                                                                cagccaCCACGGCCTGCGGTGTAGCAGCACTATCAACAGGCTGCCGAGGGCGAG‑‑ATCGGACAGGTAACCGAACGGGTTGGACCAACCCAGCATGGCGCGCGTACCCAGTGGGATCACCAGCACCAGAGCCGCCAAAg                                  <  2:1693434/137‑1 (MQ=255)
                                                                                                                                cagccaCCACGGCCTGCGGTGTAGCAGCACTATCAACAGGCTGCCGAGGGCGAG‑‑ATCGGACAGGTAACCGAACGGGTTGGACCAACCCAGCATGGCGCGCGTACCCAGTGGGATCACCAGCACCAGAGCCGCCAAAg                                  >  1:1693434/1‑137 (MQ=255)
                                                                                                                                cagccaCCACGGCCTGCGGTGTAGCAGCACTATCAACAGGCTGCCGAGGGCGAG‑‑ATCGGACAGGTAACCGAACGGGTTGGACCAACCCAGCATGGCGCGCGTACCCAGTGGGATCACCAGCACCAGAGCCGCCAAAGTGGCAAGCCGGGt                     >  1:1708815/1‑150 (MQ=255)
                                                                                                                                    caccacGGCCTGCGGTGTAGCAGCACTATCAACAGGCTGCCGAGGGCGAG‑‑ATCGGACAGGTAACCGAACGGGTTGGACCAACCCAGCATGGCGCGCGTACCCAGTGGGATCACCAGCACCAGAGCCGCCAAAGTGGCAAGCCGGGTGCtt                 <  2:842715/150‑1 (MQ=255)
                                                                                                                                      ccacGGCCTGCGGTGTAGCAGCACTATCAACAGGCTGCCGAGGGCGAG‑‑ATCGGACAGGTAACCGAACGGGTTGGACCAACCCAGCATGGCGCGCGTACCCAGTGGGATCACCAGCACCAGAGCCGCCAAAGTGGCAAGCCGGGTGCTTgg               <  1:1490961/150‑1 (MQ=255)
                                                                                                                                        acGGCCTGCGGTGTAGCAGCACTATCAACAGGCTGCCGAGGGCGAG‑‑ATCGGACAGGTAACCGAACGGGTTGGACCAACCCAGCATGGCGCGCGTACCCAGTGGGATCACCAGCACCAGAGCCGCCAAAGTGGCAAGCCGGGTGCTTGggt             <  2:1591547/150‑1 (MQ=255)
                                                                                                                                            ccTGCGGTGTAGCAGCACTATCAACAGGCTGCCGAGGGCGAG‑‑ATCGGACAGGTAACCGAACGGGTTGGACCAACCCAGCATGGCGCGCGTACCCAGTGGGATCACCAGCACCAGAGCCGCCAAAGTGGCAAGCCGGGTGCTTGGGTGGTc         <  1:500722/150‑1 (MQ=255)
                                                                                                                                                 ggTGTAGCAGCACTATCAACAGGCTGCCGAGGGCGAG‑‑ATCGGACAGGTAACCGAACGGGTTGGa                                                                                          >  1:977658/1‑64 (MQ=255)
                                                                                                                                                 ggTGTAGCAGCACTATCAACAGGCTGCCGAGGGCGAG‑‑ATCGGACAGGTAACCGAACGGGTTGGa                                                                                          <  2:977658/64‑1 (MQ=255)
                                                                                                                                                   tgtAGCAGCACTATCAACAGGCTGCCGAGGGCGAG‑‑ATCGGACAGGTAACCGAACGGGTTGGACCAACCCAGCATGGCGCGCGTACCCAGTGGGATCACCAGCACCAGAGCCGCCAAAACCTGTGGCCGGGTGCTTGGGTGGTCGGACAgg  <  1:527624/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GAATTCTCCATGAATTGCGGGTCGGTCAGGTAGTTAAGGTCGGCGCTGGTGGGCAGTCGGCCCACTGCGCTTACCAGTTCGGCCGAGGCCACCCACAGCGCCGCCCAGGCCAGCAGCACGGGCAGCGCCAGCCACCACGGCCTGCGGTGTAGCAGCACTATCAACAGGCTGCCGAGGGCGAG‑‑ATCGGACAGGTAACCGAACGGGTTGGACCAACCCAGCATGGCGCGCGTACCCAGTGGGATCACCAGCACCAGAGCCGCCAAAGTGGCAAGCCGGGTGCTTGGGTGGTCGGACAGG  >  NC_002947/3255476‑3255772

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