Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,830,082 G→C 5.6% 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,830,0820GC 5.6% 951.3 / 31.0 290R55G (CGC→GGC) mutSDNA mismatch repair protein MutS
Reads supporting (aligned to +/- strand):  ref base G (149/123);  new base C (12/5);  total (161/129)
Fisher's exact test for biased strand distribution p-value = 2.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 311 total aligned reads displayed.

TCCACCGGGCCCTTGCTGGTGGCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTACTGCTGCATCATCGGGGTGTGTGCG  >  NC_002947/1829934‑1830230
                                                                                                                                                    |                                                                                                                                                    
 ccACCGGGCCCTTGCTGGTGGCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCgcggg                                                                                                                                                    >  2:4145089/1‑150 (MQ=255)
 ccACCGGGCCCTTGCTGGTGGCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCgcggg                                                                                                                                                    <  1:859862/150‑1 (MQ=255)
  cACCGGGCCCTTGCTGGTGGCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCgcgggc                                                                                                                                                   >  2:4489181/1‑150 (MQ=255)
   aCCGGGCCCTTGCTGGTGGCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCATCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCgcgggcg                                                                                                                                                  <  2:6849716/150‑1 (MQ=255)
   aCCGGGCCCTTGCTGGTGGCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCgcgggcg                                                                                                                                                  >  1:4308649/1‑150 (MQ=255)
   aCCGGGCCCTTGCTGGTGGCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCgcgggcg                                                                                                                                                  <  2:5025853/150‑1 (MQ=255)
     cGGGCCCTTGCTGGTGGCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGt                                                                                                                                                >  2:5402921/1‑150 (MQ=255)
        gCCCTTGCTGGTGGCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCCCGCGGTCAg                                                                                                                                             <  2:4935099/150‑1 (MQ=255)
         cccTTGCTGGTGGCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAgg                                                                                                                                            <  1:652901/150‑1 (MQ=255)
           cTTGCTGGTGGCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACCTGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCCGGGt                                                                                                                                          >  1:6475074/1‑150 (MQ=255)
            ttGCAGGTGGCGGGGTCGCCGATCTGCTCACAGAACACCACCGACTCGCCGAGCTAGACCAGCTTGGCCAGGTAGCCTTCGAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTg                                                                                                                                         >  2:1029343/1‑150 (MQ=255)
            gtGCTGGTGGCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTg                                                                                                                                         <  1:2477788/149‑1 (MQ=255)
                tggtggCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGatat                                                                                                                                     >  1:4216098/1‑150 (MQ=255)
                tggtggCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGatat                                                                                                                                     >  1:2578229/1‑150 (MQ=255)
                  gtggCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATcc                                                                                                                                   >  2:3152679/1‑150 (MQ=255)
                   tggCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCca                                                                                                                                  >  1:1716142/1‑150 (MQ=255)
                      cGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCcagca                                                                                                                               <  1:4129082/150‑1 (MQ=255)
                           tCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTc                                                                                                                          <  2:5540402/150‑1 (MQ=255)
                             gcagATCTGCTCACAGATCACCACCGACTCGCCGAGATTGTCCAGCTTGGCCAGGTAGCCTTCCAGAGAATGGAACGGACTCCCGCACATGTGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCGGGGTGATATCCAGCAGTGTcgc                                                                                                                        >  2:6714509/4‑150 (MQ=255)
                             gCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGACTGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTcgc                                                                                                                        <  2:5689929/150‑1 (MQ=255)
                             gCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTcgc                                                                                                                        <  2:3940957/150‑1 (MQ=255)
                             gCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTcgc                                                                                                                        >  1:5015131/1‑150 (MQ=255)
                             gCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACAGGGGGATGGACTGGCCGGCCGACGGCCCGCGGGCGGGCCGGGGGGTATCCAGCAGGGTcgc                                                                                                                        >  1:1153108/1‑150 (MQ=255)
                                gATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGc                                                                                                                     >  1:1238449/1‑150 (MQ=255)
                                gATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGc                                                                                                                     >  2:2752344/1‑150 (MQ=255)
                                 aTCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGcc                                                                                                                    >  2:1917826/1‑150 (MQ=255)
                                      cTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCttttt                                                                                                               >  1:1838736/1‑150 (MQ=255)
                                       tCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTc                                                                                                              >  2:3645298/1‑150 (MQ=255)
                                        cacaGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCg                                                                                                             >  2:6634674/1‑150 (MQ=255)
                                          caGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCTGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCa                                                                                                           >  1:1478730/1‑150 (MQ=255)
                                           aGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCAt                                                                                                          >  2:2013513/1‑150 (MQ=255)
                                           aGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCAt                                                                                                          >  1:167814/1‑150 (MQ=255)
                                            gATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATc                                                                                                         >  1:2388262/1‑150 (MQ=255)
                                               caccacCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTc                                                                                                      <  1:2650885/150‑1 (MQ=255)
                                                  caccGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGACCGGAATCCCGCACCTGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTGTCGCATCTTCGtc                                                                                                   >  2:6691226/1‑149 (MQ=255)
                                                  caccGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTa                                                                                                   >  2:4082475/1‑150 (MQ=255)
                                                  caccGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTa                                                                                                   >  2:3661961/1‑150 (MQ=255)
                                                  caccGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTa                                                                                                   <  1:2180319/150‑1 (MQ=255)
                                                  caccGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTa                                                                                                   >  2:1128496/1‑150 (MQ=255)
                                                  caccGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTa                                                                                                   <  1:5996904/150‑1 (MQ=255)
                                                  caccGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTa                                                                                                   >  1:4542329/1‑150 (MQ=255)
                                                     cGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTagaa                                                                                                <  2:4661683/150‑1 (MQ=255)
                                                      gACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTagaag                                                                                               >  2:681522/1‑150 (MQ=255)
                                                            ccGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCg                                                                                         >  1:1613148/1‑150 (MQ=255)
                                                              gAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTa                                                                                       <  2:3624074/150‑1 (MQ=255)
                                                              gAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTa                                                                                       <  1:5259432/150‑1 (MQ=255)
                                                                gCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGa                                                                                     <  1:4726604/150‑1 (MQ=255)
                                                                 cTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGaa                                                                                    <  2:4279284/150‑1 (MQ=255)
                                                                 cTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGaa                                                                                    <  2:6191006/150‑1 (MQ=255)
                                                                    gACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTc                                                                                 >  1:3884963/1‑150 (MQ=255)
                                                                    gACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTc                                                                                 >  2:6126068/1‑150 (MQ=255)
                                                                    gACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATACAGCAGTTTCGCGGCCATTTTAGCATATTAGGAGAAGATCTAGTAGAAGTc                                                                                 >  1:6323277/1‑150 (MQ=255)
                                                                     aCCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCg                                                                                >  2:5127009/1‑150 (MQ=255)
                                                                      ccAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGc                                                                               <  2:2387837/150‑1 (MQ=255)
                                                                           ttGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATACAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATg                                                                          >  2:651489/1‑150 (MQ=255)
                                                                             ggCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCg                                                                        <  1:248800/150‑1 (MQ=255)
                                                                                 aGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAg                                                                    <  2:1692505/150‑1 (MQ=255)
                                                                                  ggTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGa                                                                   <  2:2195127/150‑1 (MQ=255)
                                                                                    tAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAAc                                                                 <  1:1346446/150‑1 (MQ=255)
                                                                                        cTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCa                                                             >  2:2927752/1‑150 (MQ=255)
                                                                                         ttCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCCGGGTGATATCCCGCCGTTTCGCGGCCTTTTTCGCATCTTCGGAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAg                                                            >  1:3818741/1‑150 (MQ=255)
                                                                                         ttCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAg                                                            >  2:1419730/1‑150 (MQ=255)
                                                                                             aGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCGGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTgg                                                        <  2:3843580/150‑1 (MQ=255)
                                                                                                 aaTGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCg                                                    >  1:3921038/1‑150 (MQ=255)
                                                                                                  aTGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCgg                                                   <  1:3760494/150‑1 (MQ=255)
                                                                                                   tGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCggg                                                  <  2:727958/150‑1 (MQ=255)
                                                                                                    ggAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGt                                                 <  1:1099604/150‑1 (MQ=255)
                                                                                                    ggAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGt                                                 <  2:831651/150‑1 (MQ=255)
                                                                                                    ggAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGt                                                 <  1:3510235/150‑1 (MQ=255)
                                                                                                    ggAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGt                                                 <  1:5402921/150‑1 (MQ=255)
                                                                                                       aCGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCt                                              <  2:5846362/150‑1 (MQ=255)
                                                                                                        cGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTg                                             <  1:2893420/150‑1 (MQ=255)
                                                                                                        cGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTg                                             <  1:3299863/150‑1 (MQ=255)
                                                                                                        cGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTg                                             <  1:4488990/150‑1 (MQ=255)
                                                                                                         ggAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTgg                                            >  1:6717103/1‑150 (MQ=255)
                                                                                                         ggAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCCCGGGCGGTCCGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCAAGCGGTAGAACATCAGCGGGTCCGGGTGCTgg                                            >  1:920444/1‑150 (MQ=255)
                                                                                                             tCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGtttt                                        <  1:693887/150‑1 (MQ=255)
                                                                                                               ccGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTc                                      <  1:5616780/150‑1 (MQ=255)
                                                                                                                cGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCa                                     <  2:856426/150‑1 (MQ=255)
                                                                                                                   acaTGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCt                                  >  2:4799346/1‑150 (MQ=255)
                                                                                                                   acaTGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATACAGCAGTTTCGCGGCCATTTTAGCAACTTAGTAGAATATAACGTAGAAGTCGCCCCAGCCGTAGAACAACAGCAGGTCCCGGTGCGGGTTTTTAAGCt                                  >  1:2270896/1‑150 (MQ=255)
                                                                                                                       ggggATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCTGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCTTGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCa                              <  1:1029343/150‑1 (MQ=255)
                                                                                                                       ggggATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCa                              <  1:2446942/150‑1 (MQ=255)
                                                                                                                       ggggATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTAGCATATTCGTAGAAGATATCGTAGAAGTCGCCCATGCGGTAGAACAAAAGCTGGTCCGGGTGCTGGTTTTTAAGCATCCa                              >  1:3126271/1‑150 (MQ=255)
                                                                                                                              gACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTActgc                       <  2:4216098/150‑1 (MQ=255)
                                                                                                                              gACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTActgc                       <  2:2578229/150‑1 (MQ=255)
                                                                                                                              gACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTActgc                       <  2:1304799/150‑1 (MQ=255)
                                                                                                                                      atcacccTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTACTGCTGcatcat               >  2:3732780/7‑150 (MQ=255)
                                                                                                                                      atcacccTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTACTGCTGcatcat               >  2:5913202/7‑150 (MQ=255)
                                                                                                                                        ccgACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTACTGCTGCATCATCg             <  2:5999886/150‑1 (MQ=255)
                                                                                                                                         cgACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGAtctc                                                                                          <  2:5953074/72‑1 (MQ=255)
                                                                                                                                         cgACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGAtctc                                                                                          >  1:3655081/1‑72 (MQ=255)
                                                                                                                                          cccTGACCGCCCGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGt                                                                                                    >  2:5235469/3‑61 (MQ=38)
                                                                                                                                          cccTGACCGCCCGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTACTGCTGCATCATCggg           >  2:4454640/3‑150 (MQ=255)
                                                                                                                                           ccTGACCGCCCGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTACTGCTGCATCATCgggg          >  1:568046/2‑150 (MQ=255)
                                                                                                                                           ccTGACCGCCCGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTACTGCTGCATCATCgggg          >  1:5874801/2‑150 (MQ=255)
                                                                                                                                             tGACCGCGGGCGGTCAGGGTGATATCCAGCCGTTTCGCGGCCTTGTTCGCATCTTCGTAGAAGATCTCGTCGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTACTGCTGCAGCAGCgggggg        >  2:3289824/1‑148 (MQ=255)
                                                                                                                                               aCCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTACTGCTGCATCATCGGGgtgtg      >  1:5537851/1‑150 (MQ=255)
                                                                                                                                                  gcgggcggTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTACTGCTGCATCATCGGGGTGTGTGc   >  2:1584824/1‑150 (MQ=255)
                                                                                                                                                   cgggcggTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTACTGCTGCATCATCGGGGTGTGTGCg  <  2:1153108/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TCCACCGGGCCCTTGCTGGTGGCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTACTGCTGCATCATCGGGGTGTGTGCG  >  NC_002947/1829934‑1830230

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: