Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,254,479 G→A 6.6% A420V (GCA→GTA)  PP_2852 ← sulfatase domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,254,4790GA 6.6% 178.7 / 2.2 60A420V (GCA→GTA) PP_2852sulfatase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base G (25/31);  new base A (2/2);  total (27/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.88e-01

GAAGCCCCAGGACGAGGCCAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAA‑GTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCACGGTGAGCAGCGTCAGCATCCATGGCTTGTCCTGCTTCTGCAACTGGCCCACATAGCCCAGGGCGCCTTCGAAGAAAGCGCG  >  NC_002947/3254334‑3254627
                                                                                                                                                  |                                                                                                                                                    
gAAGCCCCAGGACGAGGCCAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAA‑GTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCt                                                                                                                                                   >  2:27784/1‑149 (MQ=255)
 aaGCCCCAGTACGAGGCCAGGGGCACGCCATCGATGCCGTGAGACTCGTCGGAA‑GTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTTCGTCATCCAGGTAGGCGACGGCTGCCt                                                                                                                                                   <  6:311590/148‑1 (MQ=255)
   gCCCCAGGACGAGGCCAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAA‑GTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGtt                                                                                                                                                >  3:293009/1‑149 (MQ=255)
        aGGACGAGGCCAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAA‑GTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGgt                                                                                                                                           <  6:224275/149‑1 (MQ=255)
         ggACGAGGCCAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAA‑GTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGgtg                                                                                                                                          <  6:26388/149‑1 (MQ=255)
          gACGAGGCCAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAA‑GTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGCAGCCGTCGCCTACCTgga                                                                                                                                                <  5:336876/142‑3 (MQ=255)
          gACGAGGCCAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAA‑GTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGCAGCCGTCGCCTACCTgga                                                                                                                                                <  5:331295/142‑3 (MQ=255)
          gACGAGGCCAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAA‑GTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGCAGCCGTCGCCTACCTgga                                                                                                                                                >  6:331295/1‑140 (MQ=255)
          gACGAGGCCAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAA‑GTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGCAGCCGTCGCCTACCTgga                                                                                                                                                >  6:336876/1‑140 (MQ=255)
             gAGGCCAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAA‑GTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTc                                                                                                                                      <  1:344420/149‑1 (MQ=255)
             gAGGCCAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAA‑GTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTc                                                                                                                                      <  4:286238/149‑1 (MQ=255)
             gAGGCCAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAA‑GTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCACGTAGGCGACGGCTGCCTGTTTCGGTGTGTc                                                                                                                                      <  1:29819/149‑1 (MQ=255)
                gCCAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAA‑GTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCg                                                                                                                                     >  8:102851/1‑147 (MQ=255)
                  cAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAA‑GTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAgc                                                                                                                                 >  3:118325/1‑149 (MQ=255)
                  cAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAA‑GTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAgc                                                                                                                                 >  6:154612/1‑149 (MQ=255)
                         aCGCCATCGATGCCGTGAGACTCGTCGGAA‑GTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAg                                                                                                                          <  3:76025/149‑1 (MQ=255)
                          cGCCATCGATGCCGTGAGACTCGTCGGAA‑GTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAgg                                                                                                                         <  7:54585/149‑1 (MQ=255)
                            ccATCGATGCCGTGAGACTCGTCGGAAGGTCACCACCACCTGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAgg                                                                                                                         >  1:22570/1‑148 (MQ=255)
                                 gATGCCGTGAGACTCGTCGGAA‑GTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTcgg                                                                                                                  >  1:361201/1‑149 (MQ=255)
                                      cGTGAGACTCGTCGGAA‑GTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTcggccggc                                                                                                             <  7:77938/149‑1 (MQ=255)
                                      cGTGAGACTCGTCGGAA‑GTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTcggccggc                                                                                                             >  3:312004/1‑149 (MQ=255)
                                             cTCGTCGGAA‑GTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAAt                                                                                                      <  4:293009/149‑1 (MQ=255)
                                              tcgtcgGAA‑GTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGaa                                                                                                       >  6:281178/1‑147 (MQ=255)
                                                           accaccaccAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTggtgggt                                                                                         >  8:346012/1‑149 (MQ=255)
                                                           accaccaccAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTggtgggt                                                                                         <  4:118325/149‑1 (MQ=255)
                                                           accaccacaAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGACCGGCGCCGAATAGGGCTggtgggt                                                                                         >  2:80890/1‑149 (MQ=255)
                                                             caccaccaGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGc                                                                                       <  3:48078/149‑1 (MQ=255)
                                                                caccaGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCa                                                                                    >  7:317035/1‑149 (MQ=255)
                                                                 accaGTGTGTCGTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCAc                                                                                   <  2:361201/149‑1 (MQ=255)
                                                                   caGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGGCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCACgg                                                                                 <  1:27784/149‑1 (MQ=255)
                                                                     gtgtgtCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCACGGTg                                                                               <  1:313573/149‑1 (MQ=255)
                                                                               gAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGACCGGCGCCGAATAGGGCTGGTGGGTGCCCACGGTGAGCAGCGt                                                                       <  1:80890/147‑1 (MQ=255)
                                                                                 ggACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCACGGTGAGCAGCGTCAg                                                                    >  3:314429/1‑148 (MQ=255)
                                                                                            cgcTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCACGGTGAGCAGCGTCAGCATCCATGGCtt                                                        <  7:346012/149‑1 (MQ=255)
                                                                                               tCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCACGGTGAGCAGCGTCAGCATCCATGGCTTGTc                                                     >  7:242020/1‑149 (MQ=255)
                                                                                               tCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCACGGTGAGCAGCGTCAGCATCCATGGCTTGTc                                                     >  1:317728/1‑149 (MQ=255)
                                                                                                 gagaTTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCACGGTGAGCAGCGTCAGCATCCAt                                                             <  1:328264/139‑1 (MQ=255)
                                                                                                 gagaTTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCACGGTGAGCAGCGTCAGCATCCAt                                                             >  2:328264/1‑139 (MQ=255)
                                                                                                 gagaTTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCACGGTGAGCAGCGTCAGCATCCATGGCTTGTcc                                                    <  2:317728/148‑1 (MQ=255)
                                                                                                 gagaTTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCACGGTGAGCAGCGTCAGCATCCATGGCTTGTCCt                                                   >  1:208357/1‑149 (MQ=255)
                                                                                                 gagaTTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGGCGAAGCGCTCGAGGTACTCGGCCGTCGCCGAATAGGGCTGGTGGGTGCCCACGGTGAGCAGCGTCAGCATCCATGGCTTGTc                                                     <  2:22570/147‑1 (MQ=255)
                                                                                                       gTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCACGGTGAGCAGCGTCAGCATCCATGGCTTGTCCTGCTTc                                              <  4:314429/148‑1 (MQ=255)
                                                                                                        tCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCACGGTGAGCAGCGTCAGCATCCATGGCTTGTCCTGCTTCt                                             <  2:183648/148‑1 (MQ=255)
                                                                                                            gAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCACGGTGAGCAGCGTCAGCATCCATGGCTTGTCCTGCTTCTGCAAc                                        >  3:279005/1‑149 (MQ=255)
                                                                                                               aTGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCACGGTGAGCAGCGTCAGCATCCATGGCTTGTCCTGCTTCTGCAACTgg                                     <  5:154612/149‑1 (MQ=255)
                                                                                                                tGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCACGGTGAGCAGCGTCAGCATCCATGGCTTGTCCTGCTTCTGCAACTGGc                                    <  5:281178/149‑1 (MQ=255)
                                                                                                                tGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCACGGTGAGCAGCGTCAGCATCCATGGCTTGTCCTGCTTCTGCAACTGGc                                    <  8:317035/149‑1 (MQ=255)
                                                                                                                 gCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCACGGTGagca                                                                           <  2:1648/109‑1 (MQ=255)
                                                                                                                 gCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCACGGTGagca                                                                           >  1:1648/1‑109 (MQ=255)
                                                                                                                   cccGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCACGGTGAGCAGCGTCAGCATCCATGGCTTGTCCTGCTTCTGCAACTGGCCca                                 >  8:157411/1‑149 (MQ=255)
                                                                                                                    ccGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCACGGTGAGCAGCGTCAGCATCCATGGCTTGTCCTGCTTCTGCAACTGGCCca                                 >  6:375096/1‑148 (MQ=255)
                                                                                                                    ccGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCACGGTGAGCAGCGTCAGCATCCATGGCTTGTCCTGCTTCTGCAACTGGCCca                                 <  5:375096/148‑1 (MQ=255)
                                                                                                                                  cAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCACGGTGAGCAGCGTCAGCATCCATGGCTTGTCCTGCTTCTGCAACTGGCCCACATAGCCCAGGgcgc                  >  7:70327/1‑149 (MQ=255)
                                                                                                                                    ggTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCACGGTGAGCAGCGTCAGCATCCATGGCTTGTCCTGCTTCTGCAACTGGCCCACATAGCCCAGGGCGCCt                <  7:157411/149‑1 (MQ=255)
                                                                                                                                     gTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCACGGTGAGCAGCGTCAGCATCCATGGCTTGTCCTGCTTCTGCAACTGGCCCACATAGACCAGGGCGCCt                >  7:202810/1‑148 (MQ=255)
                                                                                                                                        ggCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCACGGTGAGCAGCGTCAGCATCCATGGCTTGTCCTGCTTCTGCAACTGGCCCACATAGCCCAGGGCGCCTTCg             >  6:364514/1‑148 (MQ=255)
                                                                                                                                         gCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCACGGTGAGCAGCGTCAGCATCCATGGCTTGTCCTGCTTCTGCAACTGGCCCACATAGCCCAGGGCGCCTTCgaa           <  2:208357/149‑1 (MQ=255)
                                                                                                                                         gCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCACGGTGAGCAGCGTCAGCATCCATGGCTTGTCCTGCTTCTGCAACTGGCCCACATAGCCCAGGGCGCCTTCgaa           <  4:351223/149‑1 (MQ=255)
                                                                                                                                         gCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCACGGTGAGCAGCGTCAGCATCCATGGCTTGTCCTGCTTCTGCAACTGGCCCACATAGCCCAGGGCGCCTTCgaa           <  8:70327/149‑1 (MQ=255)
                                                                                                                                         gCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCACGGTGAGCAGCGTCAGCATCCATGGCTTGTCCTGCTTCTGCAACTGGCCCACATAGCCCAGGGCGCCTTCg             >  3:351223/1‑147 (MQ=255)
                                                                                                                                                  gCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCACGGTGAGCAGCGTCAGCATCCATGGCTTGTCCTGCTTCTGCAACTGGCCCACATAGCCCAGGGCGCCTTCGAAGAAAGcgcg  >  8:282792/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GAAGCCCCAGGACGAGGCCAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAA‑GTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCACGGTGAGCAGCGTCAGCATCCATGGCTTGTCCTGCTTCTGCAACTGGCCCACATAGCCCAGGGCGCCTTCGAAGAAAGCGCG  >  NC_002947/3254334‑3254627

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