Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,825,133 G→C 5.4% R55G (CGC→GGC)  mutS ← DNA mismatch repair protein MutS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,825,1330GC 5.4% 467.6 / 11.4 143R55G (CGC→GGC) mutSDNA mismatch repair protein MutS
Reads supporting (aligned to +/- strand):  ref base G (56/79);  new base C (5/3);  total (61/82)
Fisher's exact test for biased strand distribution p-value = 2.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

Only 100 of 156 total aligned reads displayed.

CTCCACCGGGCCCTTGCTGGTGGCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGC‑‑‑‑‑GAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTACTGCTGCATCATCGGGGTGTGTGCGG  >  NC_002947/1824984‑1825282
                                                                                                                                                          |                                                                                                                                                     
 tCCACCGGGCCCTTGCTGGTGGCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGC‑‑‑‑‑GAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCgcgg                                                                                                                                                      >  2:778663/1‑150 (MQ=255)
  ccACCGGGCCCTTGCTGGTGGCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGC‑‑‑‑‑GAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCgcggg                                                                                                                                                     <  1:1548073/150‑1 (MQ=255)
   cACCGGGCCCTTGCTGGTGGCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGC‑‑‑‑‑GAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCgcgggc                                                                                                                                                    >  2:1088896/1‑150 (MQ=255)
    aCCGGGCCCTTGCTGGTGGCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGC‑‑‑‑‑GAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCgcgggcg                                                                                                                                                   <  2:704093/150‑1 (MQ=255)
    aCCGGGCCCTTGCTGGTGGCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGC‑‑‑‑‑GAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCgcgggcg                                                                                                                                                   <  1:1782176/150‑1 (MQ=255)
        ggCCCGTGCTGGTGGCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGC‑‑‑‑‑GAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCa                                                                                                                                               <  1:3719825/150‑1 (MQ=255)
          cccTTGCTGGTGGCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGC‑‑‑‑‑GAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAgg                                                                                                                                             <  1:2164982/150‑1 (MQ=255)
           ccTTGCTGGTGGCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGC‑‑‑‑‑GAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCCggg                                                                                                                                            >  1:799790/1‑150 (MQ=255)
           ccTTGCTGGTGGCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGC‑‑‑‑‑GAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCCCGCGGTCAgtc                                                                                                                                            <  1:2664279/150‑3 (MQ=255)
            cgtGCTGGTGGCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGC‑‑‑‑‑GAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGt                                                                                                                                           <  2:2412991/148‑1 (MQ=255)
            cTTGCTGGTGGCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGC‑‑‑‑‑GAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGt                                                                                                                                           <  2:4284212/150‑1 (MQ=255)
                  ggtggCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGC‑‑‑‑‑GAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATc                                                                                                                                     <  1:3608276/150‑1 (MQ=255)
                   gtggCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGC‑‑‑‑‑GAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATcc                                                                                                                                    <  2:5120095/150‑1 (MQ=255)
                    tggCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGC‑‑‑‑‑GAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCca                                                                                                                                   <  1:3168426/150‑1 (MQ=255)
                       cGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGC‑‑‑‑‑GAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCcagca                                                                                                                                <  2:4890711/150‑1 (MQ=255)
                        ggggTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGC‑‑‑‑‑GAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCcagcag                                                                                                                               >  1:3045405/1‑150 (MQ=255)
                                      gCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGC‑‑‑‑‑GAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCtttt                                                                                                                 <  2:5035483/150‑1 (MQ=255)
                                        tCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGC‑‑‑‑‑GAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTc                                                                                                               >  2:565063/1‑150 (MQ=255)
                                        tCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGC‑‑‑‑‑GAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTc                                                                                                               >  2:132669/1‑150 (MQ=255)
                                        tCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGC‑‑‑‑‑GAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTc                                                                                                               <  2:3391158/150‑1 (MQ=255)
                                        tCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGC‑‑‑‑‑GAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTc                                                                                                               >  1:4939347/1‑150 (MQ=255)
                                        tCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGC‑‑‑‑‑GAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTc                                                                                                               >  1:348311/1‑150 (MQ=255)
                                          acaGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGC‑‑‑‑‑GAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCCGGGTGATATACAGCAGTTTCGCGGCCGTTTTCGc                                                                                                             >  1:2285230/1‑150 (MQ=255)
                                           caGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGC‑‑‑‑‑GAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCa                                                                                                            <  2:2212618/150‑1 (MQ=255)
                                            aGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGC‑‑‑‑‑GAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCAt                                                                                                           >  2:2535196/1‑150 (MQ=255)
                                                caccacCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGC‑‑‑‑‑GAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGGCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTc                                                                                                       >  1:1361800/1‑150 (MQ=255)
                                                  ccaccGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGC‑‑‑‑‑GAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGt                                                                                                                                                 <  1:4756767/106‑1 (MQ=255)
                                                  ccaccGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGC‑‑‑‑‑GAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGt                                                                                                                                                 >  2:4756767/1‑106 (MQ=255)
                                                   caccGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGC‑‑‑‑‑GAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTa                                                                                                    <  1:2664989/150‑1 (MQ=255)
                                                   caccGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGC‑‑‑‑‑GAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTa                                                                                                    <  1:1722933/150‑1 (MQ=255)
                                                   caccGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGC‑‑‑‑‑GAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTa                                                                                                    >  2:2574168/1‑150 (MQ=255)
                                                   caccGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGC‑‑‑‑‑GAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTa                                                                                                    >  1:4897751/1‑150 (MQ=255)
                                                      cGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGC‑‑‑‑‑GAATGGAACGGAATCCCGCACCTGGGGATGGACTGGCCGGCCGACTGACCCCGGGCGGTCCGGGTGATATCCAGCCGTTTCGCGGCCTTTTTCGCATCTTCGTagaa                                                                                                 >  1:4008148/1‑150 (MQ=255)
                                                      cGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGC‑‑‑‑‑GAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTagaa                                                                                                 <  2:2683824/150‑1 (MQ=255)
                                                       gACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGC‑‑‑‑‑GAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTagaag                                                                                                >  1:2914501/1‑150 (MQ=255)
                                                          tCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGC‑‑‑‑‑GAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGAtc                                                                                             <  1:2950537/150‑1 (MQ=255)
                                                           cGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCAAATGGAACGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTagaa                                                                                                 <  2:1203643/150‑1 (MQ=255)
                                                            gCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGC‑‑‑‑‑GAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGAtctc                                                                                           <  2:4170736/150‑1 (MQ=255)
                                                            gCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGC‑‑‑‑‑GAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTAGCAACTTCGTAGAAGAtatc                                                                                           >  1:2720485/1‑150 (MQ=255)
                                                            gCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGC‑‑‑‑‑GAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACGGACCGCCGGCGGTCCGGGGGGGATCCAGCAGTTTCGCGGCCGTTTTCGCATCTTCGTAGAAGAtctc                                                                                           >  1:1467942/1‑150 (MQ=255)
                                                              cGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGC‑‑‑‑‑GAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGt                                                                                         >  1:1686713/1‑150 (MQ=255)
                                                               gAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGC‑‑‑‑‑GAATTGAACGGAATCCCGAACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGTGTGATATCCAGCAGTTTCTCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTa                                                                                        <  1:2651367/150‑1 (MQ=255)
                                                               gAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGC‑‑‑‑‑GAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTa                                                                                        <  1:2784054/150‑1 (MQ=255)
                                                               gAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGC‑‑‑‑‑GAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTa                                                                                        <  1:2784089/150‑1 (MQ=255)
                                                               gAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGC‑‑‑‑‑GAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTa                                                                                        <  2:1748054/150‑1 (MQ=255)
                                                               gAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGC‑‑‑‑‑GAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTa                                                                                        <  2:1696055/150‑1 (MQ=255)
                                                               gAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGC‑‑‑‑‑GAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTa                                                                                        <  1:4392345/150‑1 (MQ=255)
                                                               gAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGC‑‑‑‑‑GAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTa                                                                                        <  2:4193574/150‑1 (MQ=255)
                                                                  cTTGACCAGCTTGGCCAGGTAGCCTTCCAGC‑‑‑‑‑GAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGaa                                                                                     <  1:67641/150‑1 (MQ=255)
                                                                  cTTGACCAGCTTGGCCAGGTAGCCTTCCAGC‑‑‑‑‑GAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGaa                                                                                     <  1:610978/150‑1 (MQ=255)
                                                                   ttGACCAGCTTGGCCAGGTAGCCTTCCAGC‑‑‑‑‑GAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGGGGTGAGGGTGTTATACAGCAGTTTCGCGGCCTTTTTCGCATATTCGTAGAAGATATAGTATAAg                                                                                    >  1:870774/1‑150 (MQ=255)
                                                                     gACCAGCTTGGCCAGGTAGCCTTCCAGC‑‑‑‑‑GAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTc                                                                                  >  1:341773/1‑150 (MQ=255)
                                                                     gACCAGCTTGGCCAGGTAGCCTTCCAGC‑‑‑‑‑GAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTAGTAGAAGATCTCGTAGAAGTc                                                                                  >  1:4303688/1‑150 (MQ=255)
                                                                      aCCAGCTTGGCCAGGTAGCCTTCCAGC‑‑‑‑‑GAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCg                                                                                 <  2:5019069/150‑1 (MQ=255)
                                                                       ccAGCTTGGCCAGGTAGCCTTCCAGC‑‑‑‑‑GAATGGAACGGAATCCCGCACATGGGGAAGGACTGGCCGGCCGACTGACCGCGGGCGGTAAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGc                                                                                >  2:3215371/1‑150 (MQ=255)
                                                                        cAGCTTGGCCAGGTAGCCTTCCAGC‑‑‑‑‑GAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGcc                                                                               <  2:2402044/150‑1 (MQ=255)
                                                                                 cAGGTAGCCTTCCAGC‑‑‑‑‑GAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTa                                                                      <  1:878198/150‑1 (MQ=255)
                                                                                 cAGGTAGCCTTCCAGC‑‑‑‑‑GAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTAGCATCTTCGTAGAAGATCTAGTAGAAGTAGCCCATGCGGTa                                                                      >  1:1025368/1‑150 (MQ=255)
                                                                                    gTAGCCTTCCAGC‑‑‑‑‑GAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTTGCCCATGCGGTAGaa                                                                   <  2:4329407/150‑1 (MQ=255)
                                                                                    gTAGCCTTCCAGC‑‑‑‑‑GAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCCTGCGGTAGaa                                                                   >  2:17927/1‑150 (MQ=255)
                                                                                    gTAGCCTTCCAGC‑‑‑‑‑GAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGaa                                                                   <  1:1344034/150‑1 (MQ=255)
                                                                                          ttCCAGC‑‑‑‑‑GAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTTACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAg                                                             <  1:1518335/150‑1 (MQ=255)
                                                                                          ttCCAGC‑‑‑‑‑GAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAg                                                             >  2:3113043/1‑150 (MQ=255)
                                                                                           tCCAGC‑‑‑‑‑GAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGc                                                            <  1:2058790/150‑1 (MQ=255)
                                                                                             cAGC‑‑‑‑‑GAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTg                                                          <  1:2452413/150‑1 (MQ=255)
                                                                                             cAGC‑‑‑‑‑GAATGGAACGGAATACCGCACATGGGGATGGACTGGCCGGCCGACTGAACCCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTAGCATCATCGTGGAAGATCTCCGAGAAGTCGCCCCGGCGGGAGAacaacaccgg                                                          >  1:967769/1‑146 (MQ=255)
                                                                                              aGC‑‑‑‑‑GAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTgg                                                         <  1:5001992/150‑1 (MQ=255)
                                                                                              aGC‑‑‑‑‑GAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATACAGCAGTTTCGCGGCCTTTTTAGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTgg                                                         >  1:1427384/1‑150 (MQ=255)
                                                                                               gC‑‑‑‑‑GAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGGCAGGGTGATATACAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCCTGCGGTAGAACATCAGCTggg                                                        >  1:1286019/1‑149 (MQ=255)
                                                                                                       aaTGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCg                                                     <  1:778663/150‑1 (MQ=255)
                                                                                                        aTGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCgg                                                    <  2:733855/150‑1 (MQ=255)
                                                                                                        aTGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCgg                                                    <  2:3962000/150‑1 (MQ=255)
                                                                                                         tGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCggg                                                   >  2:2440665/1‑150 (MQ=255)
                                                                                                         tGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCggg                                                   <  1:1411168/150‑1 (MQ=255)
                                                                                                         tGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCggg                                                   >  2:378563/1‑150 (MQ=255)
                                                                                                         tGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGt                                                  >  1:4771349/1‑150 (MQ=255)
                                                                                                          ggAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATACAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGt                                                  >  1:3890500/1‑150 (MQ=255)
                                                                                                             aCGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCt                                               <  2:903611/150‑1 (MQ=255)
                                                                                                                  aTCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGttt                                          <  2:187849/150‑1 (MQ=255)
                                                                                                                  aTCCCGCACATGGGGACGGACTGGCCGGCCGACTGACCGCGGGCGGTAAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGttt                                          <  1:3215371/150‑1 (MQ=255)
                                                                                                                    cccGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGttttt                                        <  2:1943378/150‑1 (MQ=255)
                                                                                                                     ccGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTc                                       <  2:3499072/150‑1 (MQ=255)
                                                                                                                     ccGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTc                                       <  2:3499094/150‑1 (MQ=255)
                                                                                                                     ccGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCATTTTAGCATCATCGTAGAAGATCTCGTAGAAGTAGCCCAAGCGGTAGAACATCAGCAGGTCCGGGTGCTGGttttta                                       >  1:4444355/1‑149 (MQ=255)
                                                                                                                      cGCACTTGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCCTGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCa                                      <  2:204144/150‑1 (MQ=255)
                                                                                                                                gATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTa                            >  2:826217/1‑150 (MQ=255)
                                                                                                                                   ggACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTActg                         <  2:3717561/150‑1 (MQ=255)
                                                                                                                                     aCTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTActgct                       <  2:1720300/150‑1 (MQ=255)
                                                                                                                                     aCTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTActgct                       >  2:2468287/1‑150 (MQ=255)
                                                                                                                                     aCTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATACAACAGTTTAGCGGCCATTTTAGCATCTTCGTAGAAGATAAAGTATAAGTAGCCCAAGCCGTAGAAAAAAAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTAcagct                       >  1:2216576/1‑150 (MQ=255)
                                                                                                                                         gccggccgACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTACTGCTGcat                   <  1:1764037/150‑1 (MQ=255)
                                                                                                                                              cacccTGACCGCCCGCGGTCAGTGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTACTGCTGCATCATCg              >  2:2310505/5‑150 (MQ=255)
                                                                                                                                                cccTGACCGCCCGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTACTGCTGCATCATCggg            >  2:1203759/3‑150 (MQ=255)
                                                                                                                                                   tGACCGCGGGCTGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTACTGCTGCATCATCGGGgtg         >  1:2840674/1‑150 (MQ=255)
                                                                                                                                                   tGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTACTGCTGCATCATCgggggg         >  2:2127038/1‑148 (MQ=255)
                                                                                                                                                   tGACCGCCCGCGGTCAGTGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTACTGCTGCATCATCGGGgtg         >  2:1349016/1‑150 (MQ=255)
                                                                                                                                                     aCCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTACTGCTGCATCATCGGGgtgtg       >  1:3273611/1‑150 (MQ=255)
                                                                                                                                                        gcgggcggTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTACTGCTGCATCATCGGGGTGTGTGc    >  1:4951064/1‑150 (MQ=255)
                                                                                                                                                         cgggcggTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTACTGCTGCATCATCGGGGTGTGTGCg   <  2:4281654/150‑1 (MQ=255)
                                                                                                                                                         cgggcggTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTACTGCTGCATCATCGGGGTGTGTGCg   <  2:3247206/150‑1 (MQ=255)
                                                                                                                                                          |                                                                                                                                                     
CTCCACCGGGCCCTTGCTGGTGGCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGC‑‑‑‑‑GAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTACTGCTGCATCATCGGGGTGTGTGCGG  >  NC_002947/1824984‑1825282

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