Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 87,956 G→C 5.2% R454R (CGC→CGG betC ← choline‑sulfatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294787,9560GC 5.2% 438.3 / 9.5 135R454R (CGC→CGGbetCcholine‑sulfatase
Reads supporting (aligned to +/- strand):  ref base G (69/59);  new base C (2/5);  total (71/64)
Fisher's exact test for biased strand distribution p-value = 2.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

Only 100 of 144 total aligned reads displayed.

AGGGGGGCGGGCTGTGGATAACGTGCCTTGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTCGCGGCTCAGGTCATA  >  NC_002947/87807‑88102
                                                                                                                                                     |                                                                                                                                                  
aGGGGGGCGGGCTGTGGATAACGTGCCTTGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCcggcgg                                                                                                                                                    >  1:323286/1‑150 (MQ=255)
      gcgggcTGTGGATAACGTGCCTTGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTCGCGCCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGgcgc                                                                                                                                              <  1:772504/150‑1 (MQ=255)
          gcTGTGGATAACGTGCCTTGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCAGCGCCGCCgccggc                                                                                                                                          <  1:322498/150‑1 (MQ=255)
            tgtgGATAACGTGCCTTGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCgctggctg                                                                                                                                        <  2:1491115/150‑1 (MQ=255)
               ggATAACGTGCCTTGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCAGCGCCGCCGccggctggtg                                                                                                                                     <  2:1939412/150‑3 (MQ=255)
               agATAACGTTCCTTGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGTCAGCGCTTCGGCCACCAGCCAGCGCCGCCGccggctggtg                                                                                                                                     <  1:1083737/149‑3 (MQ=255)
                gATAACGTGCCTTGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCa                                                                                                                                    >  1:1534059/1‑150 (MQ=255)
                gATAACGTGCCTTGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCa                                                                                                                                    <  2:363913/150‑1 (MQ=255)
                gATAACGTGCCTTGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCAGCGCCGCCGccggctggtgg                                                                                                                                    <  1:1484592/150‑4 (MQ=255)
                gATAACGTGCCTTGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGCGCAAGCAATTCCAAGCGTTTCAGCTTGCTGATGGACAGCGCCTCGGCAACCAGCCGGCGGCGGCGCTGGCTGGCCa                                                                                                                                    >  1:2223124/1‑150 (MQ=255)
                    aCGTGCCTTGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTAc                                                                                                                                <  2:2173160/150‑1 (MQ=255)
                    aCGTGCCTTGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTAc                                                                                                                                <  2:1002896/150‑1 (MQ=255)
                        gCCTTGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACctgc                                                                                                                            >  1:439481/1‑150 (MQ=255)
                           ttGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTgg                                                                                                                         >  1:1475047/1‑150 (MQ=255)
                                  cGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGGCCGGTTTGGCCAGTACCTGCTGGCGCAGGc                                                                                                                  <  1:675348/150‑1 (MQ=255)
                                   gAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTAGGCGGCTTGCCGATGGCAAGCGCTTAGGCTTTCTTCCGTCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCt                                                                                                                 >  1:2342541/1‑150 (MQ=255)
                                     ggTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTgg                                                                                                               >  1:365761/1‑150 (MQ=255)
                                       tCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTg                                                                                                             >  1:516872/1‑150 (MQ=255)
                                         aTCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGAt                                                                                                           >  1:1988057/1‑150 (MQ=255)
                                         aTCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGAt                                                                                                           >  2:629168/1‑150 (MQ=255)
                                             agaTCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTcc                                                                                                       >  1:602509/1‑150 (MQ=255)
                                              gaTCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGt                                                                                                              >  2:1230401/1‑142 (MQ=255)
                                                  gATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAgcg                                                                                                  >  1:1410309/1‑150 (MQ=255)
                                                      tgGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTcc                                                                                                       >  2:312681/1‑141 (MQ=255)
                                                       ggTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGtt                                                                                             <  1:794010/150‑1 (MQ=255)
                                                             gcATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCAGCGCCGCCGccggctggtggc                                                                                                                                   <  2:1918577/106‑5 (MQ=38)
                                                               aTGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGa                                                                                                            <  1:1945729/127‑1 (MQ=255)
                                                               aTGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGa                                                                                                            >  2:1945729/1‑127 (MQ=255)
                                                               aTGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCa                                                                                     >  1:769059/1‑150 (MQ=255)
                                                                tGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCAt                                                                                    <  1:2056557/150‑1 (MQ=255)
                                                                tGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCAt                                                                                    <  2:1475047/150‑1 (MQ=255)
                                                                     tGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACa                                                                               >  1:2075035/1‑150 (MQ=255)
                                                                     tGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACa                                                                               >  1:2396984/1‑150 (MQ=255)
                                                                     tGTTGGCTGGCGCCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGTTTGCCGGTGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCCCTCCCTGCTGGCGCAGGCTGGGGATGTCCCGGCGCTGTTTCGCTCCATCGAcc                                                                               >  2:319568/1‑149 (MQ=255)
                                                                       ttGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCgctg                                                                                                                                            <  2:2158675/87‑1 (MQ=255)
                                                                       ttGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCgctg                                                                                                                                            >  1:2158675/1‑87 (MQ=255)
                                                                                gTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGACACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTgcag                                                                    >  2:855499/1‑150 (MQ=255)
                                                                                  ccaccaGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTgccgca                                                                  >  2:92283/1‑150 (MQ=255)
                                                                                     ccaGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGtt                                                                                             >  1:276438/1‑120 (MQ=255)
                                                                                     ccaGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTgcagcagca                                                               >  2:1884886/1‑150 (MQ=255)
                                                                                      caGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCAc                                                              <  2:2075035/150‑1 (MQ=255)
                                                                                          ggttggtGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCtgg                                                          <  2:2278529/150‑1 (MQ=255)
                                                                                          ggttggtGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCtgg                                                          <  1:1924474/150‑1 (MQ=255)
                                                                                          ggttggtGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCtgg                                                          <  1:2362975/150‑1 (MQ=255)
                                                                                            ttggtGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGc                                                                          >  1:1615661/1‑132 (MQ=255)
                                                                                            ttggtGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCtggtg                                                        <  2:792933/150‑1 (MQ=255)
                                                                                            ttggtGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCtggtg                                                        <  1:58988/150‑1 (MQ=255)
                                                                                            ttggtGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCtggtg                                                        >  2:966035/1‑150 (MQ=255)
                                                                                            ttggtGGTCCCAGCTTTTCAGCTTGCCGATCGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGc                                                                          <  2:1615661/132‑1 (MQ=255)
                                                                                              ggtggtCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTcc                                                                                                       <  2:1712407/101‑1 (MQ=255)
                                                                                              ggtggtCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTcc                                                                                                       >  1:1712407/1‑101 (MQ=255)
                                                                                               gtggtCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCgctggctg                                                                                                                                        >  1:679525/1‑67 (MQ=255)
                                                                                               gtggtCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCgctggctg                                                                                                                                        <  2:679525/67‑1 (MQ=255)
                                                                                                 ggtCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCg                                                   <  2:454796/150‑1 (MQ=255)
                                                                                                    cccAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGc                                                >  2:2140773/1‑150 (MQ=255)
                                                                                                    cccAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGc                                                >  1:1299074/1‑150 (MQ=255)
                                                                                                     ccAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCt                                               >  2:493033/1‑150 (MQ=255)
                                                                                                     ccAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCt                                               >  1:2305894/1‑150 (MQ=255)
                                                                                                      cAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTg                                              >  2:360045/1‑150 (MQ=255)
                                                                                                      cAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGACCAGTACCTGCTGGCGCAGGCTGGTGAt                                                                                                           <  2:1795865/89‑1 (MQ=255)
                                                                                                          ttttCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCgg                                          <  2:1988057/150‑1 (MQ=255)
                                                                                                          ttttCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCgg                                          >  2:1857624/1‑150 (MQ=255)
                                                                                                           tttCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGt                                         <  2:769059/150‑1 (MQ=255)
                                                                                                             tCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGa                                       <  1:911796/150‑1 (MQ=255)
                                                                                                              cAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAg                                      <  1:455090/150‑1 (MQ=255)
                                                                                                              cAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAg                                      >  1:1702587/1‑150 (MQ=255)
                                                                                                                gCTTGGCGATGGCCAGCGCTTCGGCCACCCGCCGGCGGCGGCGTTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGt                                    <  1:1265751/150‑1 (MQ=255)
                                                                                                                 cTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCAt                                                                                    >  1:116217/1‑101 (MQ=255)
                                                                                                                  ttGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCaccgt                                                                                 >  1:1104538/1‑99 (MQ=255)
                                                                                                                  ttGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTtc                                  >  2:1921226/1‑150 (MQ=255)
                                                                                                                  ttGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTtc                                  >  2:718945/1‑150 (MQ=255)
                                                                                                                  ttGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCACCgt                                                                                 <  2:1104538/103‑2 (MQ=255)
                                                                                                                  ttGCCGATGGCCAGCGCTTCGCCCACCAGCGGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGGGAGGTCCGAGCGCTGTTGCGCTTAATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGCGGTtc                                  >  2:373717/1‑150 (MQ=255)
                                                                                                                    gCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTtctc                                <  2:741765/150‑1 (MQ=255)
                                                                                                                      cGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGAAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTgcagca                                                                  <  2:1639197/114‑1 (MQ=255)
                                                                                                                      cGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGAAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTgcagca                                                                  >  1:1639197/1‑114 (MQ=255)
                                                                                                                           gCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCAt                                                                                    <  2:2014497/91‑1 (MQ=255)
                                                                                                                           gCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCAt                                                                                    >  1:2014497/1‑91 (MQ=255)
                                                                                                                            ccAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGt                                                                                                              >  2:556605/1‑64 (MQ=255)
                                                                                                                            ccAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGt                                                                                                              <  1:556605/64‑1 (MQ=255)
                                                                                                                             cAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTg                       >  2:1312461/1‑150 (MQ=255)
                                                                                                                              aGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGc                      <  2:854647/150‑1 (MQ=255)
                                                                                                                               gcgcTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCg                     <  1:822772/150‑1 (MQ=255)
                                                                                                                               gcgcTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCg                     <  1:718945/150‑1 (MQ=255)
                                                                                                                                cgcTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCgg                    <  1:1312461/150‑1 (MQ=255)
                                                                                                                                  cTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCgggg                  >  1:1845307/1‑149 (MQ=255)
                                                                                                                                  cTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGt                  <  1:966035/150‑1 (MQ=255)
                                                                                                                                   ttCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCtggtg                                                        >  1:534864/1‑111 (MQ=255)
                                                                                                                                   ttCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCtggtg                                                        <  2:534864/111‑1 (MQ=255)
                                                                                                                                   ttCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGtt                 >  2:1331443/1‑149 (MQ=255)
                                                                                                                                     cGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCCGGCcgc               <  2:602509/150‑1 (MQ=255)
                                                                                                                                      ggCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCgggagaac              <  1:360045/150‑6 (MQ=255)
                                                                                                                                      ggCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTcgcg              >  1:1864199/1‑150 (MQ=255)
                                                                                                                                            caGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGa                                       <  2:563274/119‑1 (MQ=255)
                                                                                                                                                cggcggcggCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTCGCGGCTCAGGTCa    <  2:439481/150‑1 (MQ=255)
                                                                                                                                                  gcggcggcgCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTcc                                                    >  1:1947099/1‑100 (MQ=255)
                                                                                                                                                  gcggcggcgCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCTTGCGGGTCGCGGCTCAGGTCTTa  <  2:1299074/150‑1 (MQ=255)
                                                                                                                                                  gcggcggcgCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTCGCGGCTCAGGTCATa  <  1:866366/150‑1 (MQ=255)
                                                                                                                                                  gcggcggcgCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTCGCGGCTCAGGTCATa  <  2:1473262/150‑1 (MQ=255)
                                                                                                                                                  gcggcggcgCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTCGCGGCTCAGGTCATa  >  2:1706991/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
AGGGGGGCGGGCTGTGGATAACGTGCCTTGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTCGCGGCTCAGGTCATA  >  NC_002947/87807‑88102

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: