Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,567,728 A→T 8.8% A577A (GCT→GCA mutL ← DNA mismatch repair protein MutL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,567,7280AT 8.8% 449.9 / 26.4 148A577A (GCT→GCAmutLDNA mismatch repair protein MutL
Reads supporting (aligned to +/- strand):  ref base A (64/70);  new base T (9/4);  total (73/75)
Fisher's exact test for biased strand distribution p-value = 1.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.52e-01

Only 100 of 151 total aligned reads displayed.

GTTTGTCCAGATCGTCCAGCCCCATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGGG  >  NC_002947/5567580‑5567877
                                                                                                                                                    |                                                                                                                                                     
 tttGTCCAGATCGTCCAGCCCCATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGc                                                                                                                                                     <  1:1008844/150‑1 (MQ=255)
    gTCCAGATCGTCCAGCCCCATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCAt                                                                                                                                                  >  2:1258770/1‑150 (MQ=255)
     tCCAGATCGTCCAGCCCCATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATg                                                                                                                                                 <  2:4939337/150‑1 (MQ=255)
       cAGATCGTCCAGCCCCATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGt                                                                                                                                               >  1:2763058/1‑150 (MQ=255)
       cAGATCGTCCAGCCCCATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGt                                                                                                                                               >  2:3172760/1‑150 (MQ=255)
        aGATCGTCCAGCCCCATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTg                                                                                                                                              <  2:4553965/150‑1 (MQ=255)
        aGATCGTCCAGCCCCATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTg                                                                                                                                              >  1:1553839/1‑150 (MQ=255)
        aGATCGTCCAGCCCCATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCATGGCTTGCCACGGTg                                                                                                                                              <  1:4214017/150‑1 (MQ=255)
         gATCGTCCAGCCCCATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGc                                                                                                                                             >  2:774104/1‑150 (MQ=255)
         gATCGTCCAGCCCCATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGc                                                                                                                                             >  1:3397226/1‑150 (MQ=255)
          aTCGTCCAGCCCCATCTGGGTCCAGGTAGGACGACGATGGTTGCTCTGGCCGCTGCGCTCGGTGTTGTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGcc                                                                                                                                            <  1:3680391/150‑1 (MQ=255)
          aTCGTCCAGCCCCATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCACTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCATGGCTTGCCACGGTGcc                                                                                                                                            <  1:1927007/150‑1 (MQ=255)
            cGTCCAGCCCCATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGa                                                                                                                                          <  1:4039802/150‑1 (MQ=255)
            cGTCCAGCCCCATCTGGGTCCAGGGAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCCTGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGa                                                                                                                                          <  2:4315784/150‑1 (MQ=255)
             gTCCAGCCCCATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGag                                                                                                                                         >  2:2599265/1‑150 (MQ=255)
             gTCCAGCCCCATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGag                                                                                                                                         >  1:2137794/1‑150 (MQ=255)
                 aGCCCCATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTTGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGagcagc                                                                                                                                     >  2:41188/1‑150 (MQ=255)
                 aGCCCCATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCTCGGCACCATGGCTTGCCACGGTGCCGtgcgcg                                                                                                                                     <  2:2224116/150‑4 (MQ=255)
                 aGCCCCATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGagcagc                                                                                                                                     >  2:1663854/1‑150 (MQ=255)
                    cccATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCa                                                                                                                                  <  1:66786/150‑1 (MQ=255)
                     ccATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCAt                                                                                                                                 >  2:3423848/1‑150 (MQ=255)
                      cATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCAtt                                                                                                                                >  1:2139162/1‑150 (MQ=255)
                      cATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCAtt                                                                                                                                <  2:4320741/150‑1 (MQ=255)
                       aTCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTg                                                                                                                               <  1:3059332/150‑1 (MQ=255)
                         cTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAg                                                                                                                             <  1:4579281/150‑1 (MQ=255)
                         cTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAg                                                                                                                             <  1:1374743/150‑1 (MQ=255)
                          tGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAgg                                                                                                                            >  2:213636/1‑150 (MQ=255)
                          tGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCACCTCAATGAgc                                                                                                                            <  2:1057813/150‑2 (MQ=255)
                           gggTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCACGCCATGGTGCCGAGCAGCTCATTGCGGt                                                                                                                           <  2:3723027/150‑1 (MQ=255)
                                   gTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGa                                                                                                                   >  1:1585871/1‑150 (MQ=255)
                                      ggACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATg                                                                                                                <  1:3940576/150‑1 (MQ=255)
                                      ggACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATg                                                                                                                >  1:1498408/1‑150 (MQ=255)
                                      ggACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCAGGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATg                                                                                                                >  2:3497133/1‑150 (MQ=255)
                                       gacgacCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGc                                                                                                               <  1:2213500/150‑1 (MQ=255)
                                       gacgacCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGc                                                                                                               <  2:1056038/150‑1 (MQ=255)
                                             cctGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCAGGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCAGGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCAc                                                                                                         >  2:1435751/3‑150 (MQ=255)
                                             cATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCAc                                                                                                         >  2:2849777/1‑150 (MQ=255)
                                                  ttGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTg                                                                                                    >  2:2275408/1‑150 (MQ=255)
                                                      aCTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGt                                                                                                >  1:2201463/1‑150 (MQ=255)
                                                      aCTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGt                                                                                                >  2:3172789/1‑150 (MQ=255)
                                                       cTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTa                                                                                               >  2:1616519/1‑150 (MQ=255)
                                                       cTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTa                                                                                               <  2:2655053/150‑1 (MQ=255)
                                                        tGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGTAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTAc                                                                                              <  1:5115522/150‑1 (MQ=255)
                                                        tGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTg                                                                                                       >  2:1407611/1‑141 (MQ=255)
                                                        tGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTg                                                                                                       <  1:1407611/141‑1 (MQ=255)
                                                         ggCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACt                                                                                             <  1:3336493/150‑1 (MQ=255)
                                                           ccGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTcc                                                                                           <  2:495475/150‑1 (MQ=255)
                                                               tGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATaa                                                                                       <  1:3908066/150‑1 (MQ=255)
                                                                gcgcTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAg                                                                                      <  1:4026095/150‑1 (MQ=255)
                                                                    tCGGTGTTTTCCATTTCGAGAAGAAGGGCGTTCTTCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGTATGCGATCACTGGTGCCGTACTCCATAAGGTCg                                                                                  >  2:2541904/1‑150 (MQ=255)
                                                                          ttttCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGc                                                                            <  2:2493772/150‑1 (MQ=255)
                                                                            ttCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCAc                                                                          >  1:534298/1‑150 (MQ=255)
                                                                                   tCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGc                                                                   <  2:3901211/150‑1 (MQ=255)
                                                                                       gcagcaGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCa                                                               <  2:2991744/150‑1 (MQ=255)
                                                                                         agcagGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGTCAAGCACGTCCTGCACCAgg                                                             >  2:267865/1‑150 (MQ=255)
                                                                                           cagGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCg                                                           <  2:2705336/150‑1 (MQ=255)
                                                                                                cGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTgg                                                      <  1:35275/150‑1 (MQ=255)
                                                                                                 gTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCTCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGc                                                     <  1:3084043/150‑1 (MQ=255)
                                                                                                  ttCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCt                                                    >  2:4135545/1‑150 (MQ=255)
                                                                                                  ttCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCt                                                    >  2:1554859/1‑150 (MQ=255)
                                                                                                   tCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCtt                                                   >  2:479108/1‑150 (MQ=255)
                                                                                                    cATCTCGGGAATTGCCAGACGCCGGTTGGGGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTc                                                  <  1:2794557/150‑1 (MQ=255)
                                                                                                     aTCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCa                                                 >  1:1108646/1‑150 (MQ=255)
                                                                                                     aTCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCa                                                 <  1:1746486/150‑1 (MQ=255)
                                                                                                     aTCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCa                                                 <  1:1993420/150‑1 (MQ=255)
                                                                                                       ctcGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGc                                               <  2:1566847/150‑1 (MQ=255)
                                                                                                        tcGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGcc                                              >  1:4318907/1‑150 (MQ=255)
                                                                                                        tcGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGcc                                              >  2:3353861/1‑150 (MQ=255)
                                                                                                        tcGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGcc                                              >  1:896917/1‑150 (MQ=255)
                                                                                                            gAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCt                                          <  2:4529413/150‑1 (MQ=255)
                                                                                                             aaTTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCtt                                         <  1:743613/150‑1 (MQ=255)
                                                                                                               ttGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGa                                       <  2:2471786/150‑1 (MQ=255)
                                                                                                                  ccAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCa                                    >  2:4060774/1‑150 (MQ=255)
                                                                                                                   cAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCaa                                   >  2:3667191/1‑150 (MQ=255)
                                                                                                                         ccGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCgg                             <  1:874470/150‑1 (MQ=255)
                                                                                                                              tGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCt                        <  2:281966/150‑1 (MQ=255)
                                                                                                                               ggCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCTCGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTg                       <  2:3408007/150‑1 (MQ=255)
                                                                                                                                gcgcgcACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTgg                      <  2:2137794/150‑1 (MQ=255)
                                                                                                                                gcgcgcACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCAGTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCAGCTTGAGCAAGGCCGGGATCTgg                      >  2:3050736/1‑150 (MQ=255)
                                                                                                                                 cgcgcACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGc                     <  1:3832050/150‑1 (MQ=255)
                                                                                                                                   tgctcGGCACCATGGCTTGCCACGGTGCCGTGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGa                   >  2:3307539/5‑150 (MQ=255)
                                                                                                                                     ctcGGCACCATGGCTTGCCACGGTGCCGTGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAAt                 >  1:4202812/3‑150 (MQ=255)
                                                                                                                                     ctcGGCACCATGGCTTGCCACGGTGCCGTGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAAt                 >  1:1451377/3‑150 (MQ=255)
                                                                                                                                     ctcGGCACCATGGCTTGCCACGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAAt                 >  2:653502/3‑150 (MQ=255)
                                                                                                                                      tcGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATc                >  1:3180825/2‑150 (MQ=255)
                                                                                                                                      aCGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATc                >  2:2007877/1‑150 (MQ=255)
                                                                                                                                          cACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTTCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCa            <  1:4692913/150‑1 (MQ=255)
                                                                                                                                          cACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCa            >  2:212228/1‑150 (MQ=255)
                                                                                                                                          cACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCa            <  1:1546057/150‑1 (MQ=255)
                                                                                                                                           aCCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGTTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGTGATCTGGCGAATCGCCAg           <  2:4768833/150‑1 (MQ=255)
                                                                                                                                           aCCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAg           <  1:3180863/150‑1 (MQ=255)
                                                                                                                                           aCCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAg           >  1:2541580/1‑150 (MQ=255)
                                                                                                                                           aCCATGGCTTGCCACGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCCAGGCCGGGATCTGGCGAATCGCCAg           >  2:560934/1‑150 (MQ=255)
                                                                                                                                            ccGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAgg          <  1:4394424/150‑1 (MQ=255)
                                                                                                                                            ccGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAgg          <  1:1685389/150‑1 (MQ=255)
                                                                                                                                                ggCAAGCCACGGTGCCGTGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGtct      >  2:252777/1‑150 (MQ=255)
                                                                                                                                                 gCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGtctc     >  2:3760870/1‑150 (MQ=255)
                                                                                                                                                   aaGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCgg   <  1:2599265/150‑1 (MQ=255)
                                                                                                                                                   aaGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCgg   <  1:1003204/150‑1 (MQ=255)
                                                                                                                                                    aGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCggg  <  1:3536076/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GTTTGTCCAGATCGTCCAGCCCCATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGGG  >  NC_002947/5567580‑5567877

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: