Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,665,704 A→T 5.2% L69Q (CTG→CAG)  PP_1459 ← CBS domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,665,7040AT 5.2% 582.7 / 14.0 173L69Q (CTG→CAG) PP_1459CBS domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (78/86);  new base T (6/3);  total (84/89)
Fisher's exact test for biased strand distribution p-value = 3.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

Only 100 of 195 total aligned reads displayed.

TCAGCAGGCTGGCACCAAGGCTCACGAACGCCTGCGGGTTACGCCGGGCCAGGCGCCGCGGCAGATACTCGGCAAACACCAGCAAGGTGAGGGTCGCGGCCAGGCCTGCCAGCCAGAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAGCCCCAGCACCAGCAGCTTGCCCAGGCTGGCGCACAGCACCAAGGCCTGGGCTGGCAGCACGGGATGGCCGTCTTCGCCGGCGCGCAGTGTGCCATTGAGCTGCAGGCGGGCAGCCTCCACGGCAGTGAACAGTGCTGAC  >  NC_002947/1665555‑1665853
                                                                                                                                                     |                                                                                                                                                     
tCAGCAGGCTGGCACCAAGGCTCACGAACGCCTGCGGGTTACGCCGGGCCAGGCGCCGCGGCAGATACTCGGCAAACACCAGCAAGGTGAGGGTCGCGGCCAGGCCTGCCAGCCAGAAGCCGTGCTCGCCACTCTGGCGCTGGCccacca                                                                                                                                                       <  1:1977591/150‑1 (MQ=255)
  agcagGCTGGCACCAAGGCTCACGAACGCCTGCGGGTTACGCCGGGCCAGGCGCCGCGGCAGATACTCGGCAAACACCAGCAAGGTGAGGGTCGCGGCCAGGCCTGCCAGCCAGAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACcagg                                                                                                                                                     <  1:3988027/150‑1 (MQ=255)
      ggCTGGCACCAAGGCTCACGAACGCCTGCGGGTTACGCCGGGCCAGGCGCCGCGGCAGATACTCGGCAAACACCAGCAAGGTGAGGGTCGCGGCCAGGCCTGCCAGCCAGAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACcaggcagg                                                                                                                                                 <  1:1036963/150‑1 (MQ=255)
        cTGGCACCAAGGCTCACGAACGCCTGCGGGTTACGCCGGGCCAGGCGCCGCGGCAGATACTCGGCAAACACCAGCAAGGTGAGGGTCGCGGCCAGGCCTGCCAGCCAGAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGcc                                                                                                                                               >  1:3423354/1‑150 (MQ=255)
         tGGCACCAAGGCTCACGAACGCCTGCGGGTTACGCCGGGCCAGGCGCCGCGGCAGATACTCGGCAAACACCAGCAAGGTGAGGGTCGCGGCCAGGCCTGCCAGCCAGAAGCCGTACTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCa                                                                                                                                              <  2:1434426/150‑1 (MQ=255)
         tGGCACCAAGGCTCACGAACGCCTGCGGGTTACGCCGGGCCAGGCGCCGCGGCAGATACTCGGCAAACACCAGCAAGGTGAGGGTCGCGGCCAGGCCTGCCAGCCAGAAGCCGTACTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCa                                                                                                                                              <  2:1434677/150‑1 (MQ=255)
          ggCACCAAGGCTCACGAACGCCTGCGGGTTACGCCGGGCCAGGCGCCGCGGCAGATACTCGGCAAACACCAGCAAGGTGAGGGTCGCGGCCAGGCCTGCCAGCCAGAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAg                                                                                                                                             <  2:2219068/150‑1 (MQ=255)
          ggCACCAAGGCTCACGAACGCCTGCGGGTTACGCCGGGCCAGGCGCCGCGGCAGATACTCGGCAAACACCAGCAAGGTGAGGGTCGCGGCCAGGCCTGCCAGCCAGAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAg                                                                                                                                             <  2:3761582/150‑1 (MQ=255)
          ggCACCAAGGCTCACGAACGCCTGCGGGTTACGCCGGGCCAGGCGCCGCGGCAGATACTCGGCAAACACCAGCAAGGTGAGGGTCGCGGCCAGGCCTGCCAGCCAGAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAg                                                                                                                                             <  2:4614877/150‑1 (MQ=255)
           gCACCAAGGCTCACGAACGCCTGCGGGTTACGCCGGGCCAGGCGCCGCGGCAGATACTCGGCAAACACCAGCAAGGTGAGGGTCGCGGCCAGGCCTGCCAGCCAGAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAGc                                                                                                                                            <  1:4258211/150‑1 (MQ=255)
            cACCAAGGCTCACGAACGCCTGCGGGTTACGCCGGGCCAGGCGCCGCGGCAGATACTCGGCAAACACCAGCAAGGTGAGGGTCGCGGCCAGGCCTGCCAGCCAGAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAGcc                                                                                                                                           <  2:3423354/150‑1 (MQ=255)
               cAAGGCTCACGAACGCCTGCGGGTTACGCCGGGCCAGGCGCCGCGGCAGATACTCGGCAAACACCAGCAAGGTGAGGGTCGCGGCCAGGCCTGCCAGCCAGAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAGCCcca                                                                                                                                        <  1:841812/150‑1 (MQ=255)
                 aGGCTCACGAACGCCTGCGGGTTACGCCGGGCCAGGCGCCGCGGCAGATACTCGGCAAACACCAGCAAGGTGAGGGTCGCGGCCAGGCCTGCCAGCCAGAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAGCCccagc                                                                                                                                      >  1:2026473/1‑150 (MQ=255)
                     tCACGAACGCCTGCGGGTTACGCCGGGCCAGGCGCCGCGGCAGATACTCGGCAAACACCAGCAAGGTGAGGGTCGCGGCCAGGCCTGCCAGCCAGAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCTGGCAGGCCAGCCccagcacca                                                                                                                                  >  2:3603401/1‑150 (MQ=255)
                     tCACGAACGCCTGCGGGTTACGCCGGGCCAGGCGCCGCGGCAGATACTCGGCAAACACCAGCAAGGTGAGGGTCGCGGCCAGGCCTGCCAGCCAGAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAGCCccagcacca                                                                                                                                  >  1:94113/1‑150 (MQ=255)
                      cACGAACGCCTGCGGGTTACGCCGGGCCAGGCGCCGCGGCAGATACTCGTCAAACACCAGCAAGGTGAGGGTCGCGGCCAGGCCTGCCAGCCAGAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAGCCccagcaccag                                                                                                                                 <  2:1993447/150‑1 (MQ=255)
                         gAACGCCTGCGGGTTACGCCGGGCCAGGCGCCGCGGCAGATACTCGGCAAACACCAGCACGGTGAGGGTCGCGGCCAGGCCTGCCAGCCAGAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAGCCCCAGCACcagcag                                                                                                                              >  2:3678282/1‑150 (MQ=255)
                         gAACGCCTGCGGGTTACGCCGGGCCAGGCGCCGCGGCAGATACTCGGCAAACACCAGCAAGGTGAGGGTCGCGGCCAGGCCTGCCAGCCAGAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAGCCCCAGCACcagcag                                                                                                                              >  1:353350/1‑150 (MQ=255)
                         gAACGCCTGCGGGTTACGCCGGGCCAGGCGCCGCGGCAGATACTCGGCAAACACCAGCAAGGTGAGGGTCGCGGCCAGGCCTGCCAGCCAGAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAGCCCCAGCACcagcag                                                                                                                              >  2:2142814/1‑150 (MQ=255)
                          aaCGCCTGCGGGTTACGCCGGGCCAGGCGCCGCGGCAGATACTCGGCAAACACCAGCAAGGTGAGGGTCGCGGCCAGGCCTGCCAGCCAGAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAGCCCCAGCACCagcagc                                                                                                                             <  2:1146458/150‑1 (MQ=255)
                          aaCGCCTGCGGGTTACGCCGGGCCAGGCGCCGCGGCAGATACTCGGCAAACACCAGCAAGGTGAGGGTCGCGGCCAGGCCTGCCAGCCAGAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAGCCCCAGCACCagcagc                                                                                                                             <  2:4853993/150‑1 (MQ=255)
                               cTGCGGGTTACGCCGGGCCAGGCGCCGCGGCAGATACTCGGCAAACACCAGCAAGGTGAGGGTCGCGGCCAGGCCTGCCAGCCAGAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAGCCccagcaccag                                                                                                                                 >  2:331711/1‑141 (MQ=255)
                                tGCGGGTTACGCCGGGCCAGGCGCCGCGGCAGATACTCGGCAAACACCAGCAAGGTGAGGGTCGCGGCCAGGCCTGCCAGCCAGAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAGCCCCAGCACCAGCAGCTTGccc                                                                                                                       <  1:4180028/150‑1 (MQ=255)
                                tGCGGGTTACGCCGGGCCAGGCGCCGCGGCAGATACTCGGCAAACACCAGCAAGGTGAGGGTCGCGGCCAGGCCTGCCAGCCAGAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAGCCCCAGCACCAGCAGCTTGccc                                                                                                                       >  1:3069369/1‑150 (MQ=255)
                                 gCGGGTTACGCCGGGCCAGGCGCCGCGGCAGATACTCGGCAAACACCAGCAAGGTGAGGGTCGCGGCCAGGCCTGCCAGCCAGAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAGCCCCAGCACCAGCAGCTTGCCCa                                                                                                                      <  2:1292792/150‑1 (MQ=255)
                                  cGGGTTACGCCGGGCCAGGCGCCGCGGCAGATACTCGGCAAACACCAGCAAGGTGAGGGTCGCGGCCAGGCCTGCCAGCCAGAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAGCCCCAGCACCAGCAGCTTGCCCAg                                                                                                                     <  1:2351032/150‑1 (MQ=255)
                                    ggTTACGCCGGGCCAGGCGCCGCGGCAGATACTCGGCAAACACCAGCAAGGTGAGGGTCGCGGCCAGGCCTGCCAGCCAGAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAGCCCCAGCACCAGCAGCTTGCCCAGGc                                                                                                                   <  1:3349474/150‑1 (MQ=255)
                                     gTTACGCCGGGCCAGGCGCCGCGGCAGATACTCGGCAAACACCAGCAAGGTGAGGGTCGCGGCCAGGCCTGCCAGCCAGAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAGCCCCAGCACCAGCAGCTTGCCCAGGCt                                                                                                                  >  1:2817038/1‑150 (MQ=255)
                                      ttACGCCGGGCCAGGCGCCGCGGCAGATACTCGGCAAACACCAGCAAGGTGAGGGTCGCGGCCAGGCCTGCCAGCCAGAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAGCCCCAGCACCAGCAGCTTGCCCAGGCTg                                                                                                                 >  1:629389/1‑150 (MQ=255)
                                       tACGCCGGGCCAGGCGCCGCGGCAGATACTCGGCAAACACCAGCAAGGTGAGGGTCGCGGCCAGGCCTGCCAGCCAGAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAGCCCCAGCACCAGCAGCTTGCCCAGGCTgg                                                                                                                <  1:1518896/150‑1 (MQ=255)
                                          gCCGGGCCAGGCGCCGCGGCAGATACTCGGCAAACACCAGCAAGGTGAGGGTCGCGGCCAGGCCTGCCAGCCAGAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAGCCCCAGCACCAGCAGCTTGCCCAGGCTGgcgc                                                                                                             >  1:2980383/1‑150 (MQ=255)
                                          gCCGGGCCAGGCGCCGCGGCAGATACTCGGCAAACACCAGCAAGGTGAGGGTCGCGGCCAGGCCTGCCAGCCAGAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAGCCCCAGCACCAGCAGCTTGCCCAGGCTGgcgc                                                                                                             >  1:721841/1‑150 (MQ=255)
                                             gggCCAGGCGCCGCGGCAGATACTCGGCAAACACCAGCAAGGTGAGGGTCGCGGCCAGGCCTGCCAGCCAGAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAGCCCCAGCACCAGCAGCTTGCCCAGGCTGGCGcaca                                                                                                          <  1:4530553/150‑1 (MQ=255)
                                             gggCCAGGCGCCGCGGCAGATACTCGGCAAACACCAGCAAGGTGAGGGTCGCGGCCAGGCCTGCCAGCCAGAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAGCCCCAGCACCAGCAGCTTGCCCAGGCTGGCGcaca                                                                                                          >  2:75552/1‑150 (MQ=255)
                                             gggCCAGGCGCCGCGGCAGATACTCGGCAAACACCAGCAAGGTGAGGGTCGCGGCCAGGCCTGCCAGCCAGAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAGCCCCAGCACCAGCAGCTTGCCCAGGCTGGCGcaca                                                                                                          <  1:2268133/150‑1 (MQ=255)
                                              ggCCAGGCGCCGCGGCAGATACTCGGCAAACACCAGCAAGGTGAGGGTCGCGGCCAGGCCTGCCAGCCAGAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAGCCCCAGCACCAGCAGCTTGCCCAGGCTGGCGCACAg                                                                                                         <  1:2335737/150‑1 (MQ=255)
                                                ccAGGCGCCGCGGCAGATACTCGGCAAACACCAGCAAGGTGAGGGTCGCGGCCAGGCCTGCCAGCCAGAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAGCCCCAGCACCAGCAGCTTGCCCAGGCTGGCGCACAGCa                                                                                                       >  2:448991/1‑150 (MQ=255)
                                                 cAGGCGCCGCGGCAGATACTCGGCAAACACCAGCAAGGTGAGGGTCGCGGCCAGGCCTGCCAGCCAGAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAGCCCCAGCACCAGCAGCTTGCCCAGGCTGGCGCACAGCAc                                                                                                      >  1:1784163/1‑150 (MQ=255)
                                                  aGGCGCCGCGGCAGATACTCGGCAAACACCAGCAAGGTGAGGGTCGCGGCCAGGCCTGCCAGCCAGAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAGCCCCAGCACCAGCAGCTTGCCCAGGCTGGCGCACAGCAcc                                                                                                     >  2:177086/1‑150 (MQ=255)
                                                   ggCGCCGCGGCAGATACTCGGCAAACACCAGCAAGGTGAGGGTCGCGGCCAGGCCTGCCAGCCAGAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAGCCCCAGCACCAGCAGCTTGCCCAGGCTGGCGCACAGCACCa                                                                                                    >  2:249144/1‑150 (MQ=255)
                                                    gcgcCGCGGCAGATACTCGGCAAACACCAGCAAGGTGAGGGTCGCGGCCAGGCCTGCCAGCCAGAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAGCCCCAGCACCAGCAGCTTGCCCAGGCTGGCGCACAGCACCaa                                                                                                   >  2:4473275/1‑150 (MQ=255)
                                                     cgccgcGGCAGATACTCGGCAAACACCAGCAAGGTGAGGGTCGCGGCCAGGCCTGCCAGCCAGAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAGCCCCAGCACCAGCAGCTTGCCCAGGCTGGCGCACAGCACCAAg                                                                                                  <  1:4044077/150‑1 (MQ=255)
                                                      gccgcGGCAGATACTCGGCAAACACCAGCAAGGTGAGGGTCGCGGCCAGGCCTGCCAGCCAGAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAGCCCCAGCACCAGCAGCTTGCCCAGGCTGGCGCACAGCACCAAgg                                                                                                 >  1:756340/1‑150 (MQ=255)
                                                          cgGCAGATACTCGGCAAACACCAGCAAGGTGAGGGTCGCGGCCAGGCCTGCCAGCCAGAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAGCCCCAGCACCAGCAGCTTGCCCAGGCTGGCGCACAGCACCAAGGCCTg                                                                                             <  1:2304098/150‑1 (MQ=255)
                                                           ggCAGATACTCGGCAAACACCAGCAAGGTGAGGGTCGCGGCCAGGCCTGCCAGCCGGAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAGCCCCAGCACCAGCAGCTTGCCCAGGCTGGCGCACAGCACCAAGGCCTgg                                                                                            <  1:2767541/150‑1 (MQ=255)
                                                           ggCAGATACTCGGCAAACACCAGCAAGGTGAGGGTCGCGGCCAGGCCTGCCAGCCAGAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAGCCCCAGCACCAGCAGCGTGCCCAGGCTGGCGCACAGCACCAAGGCCTgg                                                                                            <  2:596387/150‑1 (MQ=255)
                                                              aGATACTCGGCAAACACCAGCAAGGTGAGGGTCGCGGCCAGGCCTGCCAGCCAGAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACcaggcagg                                                                                                                                                 >  1:5203321/1‑94 (MQ=255)
                                                                 tACTCGGCAAACACCAGCAAGGTGAGGGTCGCGGCCAGGCCTGCCAGCCAGAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAGCCCCAGCACCAGCAGCTTGCCCAGGCTGGCGCACAGCACCAAGGCCTGGGCTGgc                                                                                      >  2:4176162/1‑150 (MQ=255)
                                                                 tACTCGGCAAACACCAGCAAGGTGAGGGTCGCGGCCAGGCCTGCCAGCCAGAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAGCCCCAGCACCAGCAGCTTGCCCAGGCTGGCGCACAGCACCAAGGCCTGGGCTGgc                                                                                      >  2:1264711/1‑150 (MQ=255)
                                                                    tCGGCAAACACCAGCAAGGTGAGGGTCGCGGCCAGGCCTGCCAGCCAGAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAGCCCCAGCACCAGCAGCTTGCCCAGGCTGGCGCACAGCACCAAGGCCTGGGCTGgcagc                                                                                   >  1:4985914/1‑150 (MQ=255)
                                                                          aaCACCAGCAAGGTGAGGGTCGCGGCCAGGCCTGCCAGCCAGAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAGCCCCAGCACCAGCAGCTTGCCCAGGCTGGCGCACAGCACCAAGGCCTGGGCTGGCAGCACGGGa                                                                             >  1:2445272/1‑150 (MQ=255)
                                                                             acCAGCAAGGTGAGGGTCGCGGCCAGGCCTGCCAGCCAGAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAGCCccagcacca                                                                                                                                  <  1:1969919/94‑1 (MQ=255)
                                                                             acCAGCAAGGTGAGGGTCGCGGCCAGGCCTGCCAGCCAGAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAGCCCCAGCACCAGCAGCTTGCCCAGGCTGGCGCACAGCACCAAGGCCTGGGCTGGCAGCACGGGATgg                                                                          >  1:1237382/1‑150 (MQ=255)
                                                                             acCAGCAAGGTGAGGGTCGCGGCCAGGCCTGCCAGCCAGAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAGCCCCAGCACCAGCAGCTTGCCCAGGCCGGCGCACAGCACCAAGGCCTGGGCTGGCAGCACGGGATgg                                                                          <  2:3543610/150‑1 (MQ=255)
                                                                                 gCAAGGTGAGGGTCGCGGCCAGGCCTGCCAGCCAGAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAGCCCCTGCACCAGCAGCTTGCCCAGGCTGGCGCACAGCACCAAGGCCTGGGCTGGCAGCAGGGGATGGCCgg                                                                      >  2:213850/1‑149 (MQ=255)
                                                                                 gCAAGGTGAGGGTCGCGGCCAGGCCTGCCAGCCAGAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAGCCCCAGCACCAGCAGCTTGCCCAGGCTGGCGCACAGCACCAAGGCCTGGGCTGGCAGCACGGGATGGCCGt                                                                      >  1:2280620/1‑150 (MQ=255)
                                                                                 gCAAGGTGAGGGTCGCGGCCAGGCCTGCCAGCCAGAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAGCCCCAGCACCAGCAGCTTGCCCAGGCTGGCGCACAGCACCAAGGCCTGGGCTGGCAGCACGGGATGGCCGt                                                                      >  1:2537034/1‑150 (MQ=255)
                                                                                     ggTGAGGGTCGCGGCCAGGCCTGCCAGCCAGAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAGCCCCAGCACCAGCAGCTTGCCCAGGCTGGCGCACAGCACCAAGGCCTGGGCTGGCAGCACGGGTTGGCCGTCTTc                                                                  <  2:5252145/150‑1 (MQ=255)
                                                                                     ggTGAGGGTCGCGGCCAGGCCTGCCAGCCAGAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAGCCCCAGCACCAGCAGCTTGCCCAGGCTGGCGCACAGCACCAAGGCCTGGGCTGGCAGCACGGGATGGCCGTCTTc                                                                  <  2:770913/150‑1 (MQ=255)
                                                                                       tGAGGGTCGCGGCCAGGCCTGCCAGCCAGAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAGCCCCAGCACCAGCAGCTTGCCCAGGCTGGCGCACAGCACCAAGGCCTGGGCTGGCAGCACGGGATGGCCGTCTTCGc                                                                <  2:1968619/150‑1 (MQ=255)
                                                                                       tGAGGGTCGCGGCCAGGCCTGCCAGCCAGAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAGCCCCAGCACCAGCAGCTTGCCCAGGCTGGCGCACAGCACCAAGGCCTGGGCTGGCAGCACGGGATGGCCGTCTTCGc                                                                <  2:3754002/150‑1 (MQ=255)
                                                                                         aGGGTCGCGGCCAGGCCTGCCAGCCAGAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAGCCCCAGCACCAGCAGCTTGCCCAGGCTGGCGCACAGCACCAAGGCCTGGGCTGGCAGCACGGGATGGCCGTCTTCGCCg                                                              <  2:5295447/150‑1 (MQ=255)
                                                                                            gTCGCGGCCAGGCCTGCCAGCCAGAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAGCCCCAGCACCAGCAGCTTGCCCAGGCTGGCGCACAGCACCAAGGCCTGGGCTGGCAGCACGGGATg                                                                           >  2:3447720/1‑134 (MQ=255)
                                                                                            gTCGCGGCCAGGCCTGCCAGCCAGAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAGCCCCAGCACCAGCAGCTTGCCCAGGCTGGCGCACAGCACCAAGGCCTGGGCTGGCAGCACGGGATg                                                                           <  1:3447720/134‑1 (MQ=255)
                                                                                             tCGCGGCCAGGCCTGCCAGCCAGAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAGCCCCAGCACCAGCAGCTTGCCCAGGCTGGCGCACAGCACCAAGGCCTGGGCTGGCAGCACGGGATGGCCGTCTTCGCCGgcgc                                                          <  2:2658950/150‑1 (MQ=255)
                                                                                                  gCCAGGCCTGCCAGCCAGAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAGCCCCAGCACCAGCAGCTTGCCCAGGCTGGCGCACAGCACCAAGGCCGGGGCTGGCAGCACGGGATGGCCGTCTTCGCCGGCGCGCAgt                                                     >  1:3005470/1‑150 (MQ=255)
                                                                                                      ggCCTGCCAGCCAGAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAGCCCCAGCACCAGCAGCTTGCCCAGGCTGGCGCACAGCACCAAGGCCTGGGCTGGCAGCACGGGATGGCCGTCTTCGCCGGCGCGCAGTGTGc                                                 <  1:3314551/150‑1 (MQ=255)
                                                                                                      ggCCTGCCAGCCAGAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAGCCCCAGCACCAGCAGCTTGCCCAGGCTGGCGCACAGCACCAAGGCCTGGGCTGGCAGCACGGGATGGCCGTCTTCGCCGGCGCGCAGTGTGc                                                 <  1:4238813/150‑1 (MQ=255)
                                                                                                          tGCCAGCCAGAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAGCCCCAGCACCAGCAGCTTGCCCAGGCTGGCGCACAGCACCAAGGCCTGGGCTGGCAGCACGGGATGGCCGTCTTCGCCGGCGCGCAGTGTGCCAtt                                             <  1:4354910/150‑1 (MQ=255)
                                                                                                          tGCCAGCCAGAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAGCCCCAGCACCAGCAGCTTGCCCAGGCTGGCGCACAGCACCAAGGCCTGGGCTGGCAGCACGGGATGGCCGTCTTCGCCGGCGCGCAGTGTGCCAtt                                             <  2:3246869/150‑1 (MQ=255)
                                                                                                          tGCCAGCCAGAAGCCGTGCTCGCCACTCGGGCGCTGGCCCACCAGGCAGGCCAGCCCCAGCACCAGCAGCTTGCCCAGGCTGGCGCACAGCACCAAGGCCTGGGCTGGCAGCACGGGATGGCCGTCTTCGCCGGCGCGCAGTGTGCCAtt                                             <  2:3866061/150‑1 (MQ=255)
                                                                                                            ccagccagAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAGCCCCAGCACCAGCAGCTTGCCCAGGCTGGCGCACAGCACCAAGGCCTGGGCTGGCAGCACGGGATGGCCGTCTTCGCCGGCGCGCAGTGTGCCATTGa                                           >  1:3595146/1‑150 (MQ=255)
                                                                                                            ccagccagAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAGCCCCAGCACCAGCAGCTTGCCCAGGCTGGCGCACAGCACCAAGGCCTGGGCTGGCAGCACGGGATGGCCGTCTTCGCCGGCGCGCAGTGTGCCATTGa                                           <  1:4861689/150‑1 (MQ=255)
                                                                                                            ccagccagAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAGCCCCAGCACCAGCAGCTTGCCCAGGCTGGCGCACAGCACCAAGGCCTGGGCTGGCAGCACGGGATGGCCGTCTTCGCCGGCGCGCAGTGTGCCATTGa                                           <  2:4372256/150‑1 (MQ=255)
                                                                                                            ccagccagAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAGCCCCAGCACCAGCAGCTTGCCCAGGCTGGCGCACAGCACCAAGGCCTGGGCTGGCAGCACGGGATGGCCGTCTTCGCCGGCGCGCAGTGTGCCATTGa                                           <  1:2967150/150‑1 (MQ=255)
                                                                                                             cagccagAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAGCCCCAGCACCAGCAGCTTGCCCAGGCTGGCGCACAGCACCAAGGCCTGGGCTGGCAGCACGGGATGGCCGTCTTCGCCGGCGCGCAgtgt                                                   <  2:2984591/141‑1 (MQ=255)
                                                                                                               gccagAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAGCCCCAGCACCAGCAGCTTGCCCAGGCTGGCGCACAGCACCAAGGCCTGGGCTGGCAGCACGGGATGGCCGTCTTCGCCGGCGCGCAGTGTGCCa                                               >  2:516496/1‑143 (MQ=255)
                                                                                                                 cagAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAGCCCCAGCACCAGCAGCTTGCCCAGGCTGGCGCACAGCACCAAGGCCTGGGCTGGCAGCACGGGATGGCCGTCTTCGCCGGCGCGCAGTGTGCCATTGAGCTGc                                      >  1:1439679/1‑150 (MQ=255)
                                                                                                                 cagAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAGCCCCAGCACCAGCAGCTTGCCCAGGCTGGCGCACAGCACCAAGGCCTGGGCTGGCAGCACGGGATGGCCGTCTTCGCCGGCGCGCAGTGTGCCATTGAGCTGc                                      >  1:5055240/1‑150 (MQ=255)
                                                                                                                       ccGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAGCCCCAGCACCAGCAGCTTGCCCAGGCTGGCGCACAGCACCAAGGCCTGGGCTGGCAGCACGGGATGGCCGTCTTCGCCGGCGCGCAGTGTGCCATTGAGCTGCAGGCgg                                <  2:1854455/150‑1 (MQ=255)
                                                                                                                           gCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAGCCCCAGCACCAGCAGCTTGCCCAGGCTGGCGCACAGCACCAAGGCCTGGGCTGGCAGCACGGGATGGCCGTCTTCGCCGGCGCGCAGTGTGCCATTGAGCTGCAGGCGGGCAg                            >  1:4590890/1‑150 (MQ=255)
                                                                                                                            cTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAGCCCCAGCACCAGCAGCTTGCCCAGGCTGGCGCACAGCACCAAGGCCTGGGCTGGCAGCACGGGATGGCCGTCTTCGCCGGCGCGCAGTGTGCCATTGAGCTGCAGGCGGGCAGc                           >  1:4946315/1‑150 (MQ=255)
                                                                                                                             tCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAGCCCCAGCACCAGCAGCTTGCCCAGGCTGGCGCACAGCACCCAGGCCTGGGCTGGCAGCACGGGag                                                                            >  1:1281336/1‑99 (MQ=255)
                                                                                                                                ccACTCTGGCGCTGGCCCACCAGGCAGGCCAGCCCCAGCACCAGCAGCTTGCCCAGGCTGGCGCACAGCACCAAGGCCTGGGCTGGCAGCACGGGATGGCCGTCTTCGCCGGCGCGCAGTGTGCCATTGAGCTGCAGGCGGGCAGCCTcc                       >  2:3786801/1‑150 (MQ=255)
                                                                                                                                 cACTCTGGCGCTGGCCCACCAGGCAGGCCAGCCCCAGCACCAGCAGCTTGCCCAGGCTGGCGCACAGCACCAAGGCCTGGGCTGGCAGCACGGGATGGCCGTCTTCGCCGGCGCGCAGTGTGCCATTGAGCTGCAGGCGGGCAGCCTCCa                      >  2:1898200/1‑150 (MQ=255)
                                                                                                                                  aCTCTGGCGCTGGCCCACCAGGCAGGCCAGCCCCAGCACCAGCAGCTTGCCCAGGCTGGCGCACAGCACCAAGGCCTGGGCTGGCAGCACGTGATGGCCGTCTTCGCCGGCGCGCAGTGTGCCATTGAGCTGCAGGCGGGCAGCCTCCAc                     <  2:4773519/150‑1 (MQ=255)
                                                                                                                                   ctctGGCGCTGGCCCACCAGGCAGGCCAGCCCCAGCACCAGCAGCTTGCCCAGGCTGGCGCACAGCACCAAGGCCTGGGCTGGCAGCACGGGATGGCCGTCTTCGCCGGCGCGCAGTGTGCCATTGAGCTGCAGGCGGGCAGCCTCCACg                    <  2:2375375/150‑1 (MQ=255)
                                                                                                                                     ctGGCGCTGGCCCACCAGGCAGGCCAGCCCCAGCACCAGCAGCTTGCCCAGGCTGGCGCACAGCACCAAGGCCTGGGCTGGCAGCACGGGATGGCCGTCTTCGCCGGCGCGCAGTGTGCCATTGAGCTGCAGGCGGGCAGCCTCCACGGc                  <  1:4473275/150‑1 (MQ=255)
                                                                                                                                      tGGCGCTGGCCCACCAGGCAGGCCAGCCCCAGCACCAGCAGCTTGCCCAGGCTGGCGCACAGCACCAAGGCCTGGGCTGGCAGCACGGGATGGCCGTCTTCGCCGGCGCGCAGTGTGCCATTGAGCTGCAGGCGGGCAGCCTCCACGGCa                 >  1:1915249/1‑150 (MQ=255)
                                                                                                                                        gggcTGGCCTGCCTGGTGGGCCAGCGCCAGCACCAGCAGCTTGCCCAGGCTGGCGCACAGCACCAAGGCCTGGGCTGGCAGCACGGGATGGCCGTCTTCGCCGGCGCGCAGTGTGCCATTGAGCTGCAGGCGGGCAGCCTCCACGGCAGt               >  1:2078870/3‑150 (MQ=255)
                                                                                                                                        gggcTGGCCTGCCTGGTGGGCCAGCGCCAGCACCAGCAGCTTGCCCAGGCTGGCGCACAGCACCAAGGCCTGGGCTGGCAGCACGGGATGGCCGTCTTCGCCGGCGCGCAGTGTGCCATTGAGCTGCAGGCGGGCAGCCTCCACGGCAGt               >  1:1816148/3‑150 (MQ=255)
                                                                                                                                        gcgcTGGCCCACCAGGCAGGCCAGCCCCAGCACCAGCAGCTTGCCCAGGCTGGCGCACAGCACCAAGGCCTGGGCTGGCAGCACGGGATGGCCGTCTTCGCCGGCGCGCAGTGTGCCATTGAGCTGCAGGCGGGCAGCCTCCACGGCAGt               >  2:3923991/1‑150 (MQ=255)
                                                                                                                                            tGGCCCACCAGGCAGGCCAGCCCCAGCACCAGCAGCTTGCCCAGGCTGGCGCACAGCACCAAGGCCTGGGCTGGCAGCACGGGATGGCCGTCTTCGCCGGCGCGCAGTGTGCCATTGAGCTGCAGGCGGGCAGCCTCCACGGCAGTGAAc           >  2:1322064/1‑150 (MQ=255)
                                                                                                                                              gCCCACCAGGCAGGCCAGCCCCAGCACCAGCAGCTTGCCCAGGCTGGCGCACAGCACCAAGGCCTGGGCTGGCAGCACGGGATGGCCGTCTTCGCCGGCGCGCAGTGTGCCATTGAGCTGCAGGCGGGCAGCCTCCACGGCAGTGAACAg         >  1:4400201/1‑150 (MQ=255)
                                                                                                                                                  accaGGCAGGCCAGCCCCAGCACCAGCAGCTTGCCCAGGCTGGCGCACAGCACCAAGGCCTGGGCTGGCAGCACGGGATGGCCGTCTTCGCCGGCGCGCAGTGTGCCATTGAGCTGCAGGCGGGCAGCCTCCACGGCAGTGAACAGTGCt     >  1:1165443/1‑150 (MQ=255)
                                                                                                                                                  accaGGCAGGCCAGCCCCAGCACCAGCAGCTTGCCCAGGCTGGCGCACAGCACCAAGGCCTGGGCTGGCAGCACGGGATGGCCGTCTTCGCCGGCGCGCAGTGTGCCATTGAGCTGCAGGCGGGCAGCCTCCACGGCAGTGAACAGTGCt     <  2:3580118/150‑1 (MQ=255)
                                                                                                                                                  accaGGCAGGCCAGCCCCAGCACCAGCAGCTTGCCCAGGCTGGCGCACAGCACCAAGGCCTGGGCTGGCAGCACGGGATGGCCGTCTTCGCCGGCGCGCAGTGTGCCATTGAGCTGCAGGCGGGCAGCCTCCACGGCAGTGAACAGTGCt     <  1:3879209/150‑1 (MQ=255)
                                                                                                                                                   ccaGGCAGGCCAGCCCCAGCACCAGCAGCTTGCCCAGGCTGGCGCACAGCACCAAGGCCTGGGCTGGCAGCACGGGATGGCCGTCTTCGCCGGCGCGCAGTGTGCCATTGAGCTGCAGGCGGGCAGCCTCCACGGCAGTGAACAGTGCTg    >  2:102256/1‑150 (MQ=255)
                                                                                                                                                    caggcaggCCAGCCCCAGCACCAGCAGCTTGCCCAGGCTGGCGCACAGCACCAAGGCCTGGGCTGGCAGCACGGGATGGCCGTCTTCGCCGGCGCGCAGTGTGCCATTGAGCTGCAGGCGGGCAGCCTCCACGGCAGTGAACAGTGcctc   <  2:1612672/150‑4 (MQ=255)
                                                                                                                                                    caggcaggCCAGCCCCAGCACCAGCAGCTTGCCCAGGCTGGCGCACAGCACCAAGGCCTGGGCTGGCAGCACGGGATGGCCGTCTTCGCCGGCGCGCAGTGTGCCATTGAGCTGCAGGCGGGCAGCCTCCACGGCAGTGAACAGTGCTGa   <  2:3069369/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCAGCAGGCTGGCACCAAGGCTCACGAACGCCTGCGGGTTACGCCGGGCCAGGCGCCGCGGCAGATACTCGGCAAACACCAGCAAGGTGAGGGTCGCGGCCAGGCCTGCCAGCCAGAAGCCGTGCTCGCCACTCTGGCGCTGGCCCACCAGGCAGGCCAGCCCCAGCACCAGCAGCTTGCCCAGGCTGGCGCACAGCACCAAGGCCTGGGCTGGCAGCACGGGATGGCCGTCTTCGCCGGCGCGCAGTGTGCCATTGAGCTGCAGGCGGGCAGCCTCCACGGCAGTGAACAGTGCTGAC  >  NC_002947/1665555‑1665853

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