Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,876,797 G→T 5.4% S451* (TCG→TAG)  ilvA‑II ← threonine deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,876,7970GT 5.4% 593.7 / 13.7 178S451* (TCG→TAG) ilvA‑IIthreonine deaminase
Reads supporting (aligned to +/- strand):  ref base G (84/84);  new base T (8/2);  total (92/86)
Fisher's exact test for biased strand distribution p-value = 1.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 189 total aligned reads displayed.

GTTTGTAGGCCGGGTTTTCGGTCTCGTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGACGAATATGCAGCTTGGCCAATTCGTTG  >  NC_002947/5876648‑5876946
                                                                                                                                                     |                                                                                                                                                     
gTTTGTAGGCCGGGTTTTCGGTCTCGTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCg                                                                                                                                                       >  2:2261062/1‑150 (MQ=255)
gTTTGTAGGCCGGGTTTTCGGTCTCGTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCg                                                                                                                                                       <  2:4138091/150‑1 (MQ=255)
 tttGTAGGCCGGGTTTTCGGTCTCGTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCga                                                                                                                                                      <  2:3958818/150‑1 (MQ=255)
 tttGTAGGCCGGGTTTTCGGTCTCGTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCga                                                                                                                                                      >  2:5044396/1‑150 (MQ=255)
  ttGTAGGCCGGGTTTTCGGTCTCGTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCgag                                                                                                                                                     <  1:1587112/150‑1 (MQ=255)
    gTAGGCCGGGTTTTCGGTCTCGTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGAt                                                                                                                                                   >  1:3374710/1‑150 (MQ=255)
       ggCCGGGTTTTCGGTCTCGTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGtt                                                                                                                                                >  2:3589700/1‑150 (MQ=255)
        gCCGGGTTTTCGGTCTCGTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTc                                                                                                                                               <  1:1088544/150‑1 (MQ=255)
        gCCGGGTTTTCGGTCTCGTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTc                                                                                                                                               <  1:2103836/150‑1 (MQ=255)
          cGGGTTTTCGGTCTCGTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCa                                                                                                                                             >  1:2626884/1‑150 (MQ=255)
          cGGGTTTTCGGTCTCGTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCa                                                                                                                                             >  1:1987719/1‑150 (MQ=255)
                   ggTCTCGTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGcgcccgcc                                                                                                                                    <  1:1126020/150‑1 (MQ=255)
                     tctcGTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCTCCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCa                                                                                                                                  >  1:1551687/1‑150 (MQ=255)
                        cGTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCt                                                                                                                               <  1:3209602/150‑1 (MQ=255)
                         gTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCtt                                                                                                                              >  2:1004658/1‑150 (MQ=255)
                         gTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCtt                                                                                                                              >  2:819633/1‑150 (MQ=255)
                          tCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCttg                                                                                                                             >  2:2033964/1‑150 (MQ=255)
                            ccAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGTTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCttgtt                                                                                                                           <  1:1945565/150‑1 (MQ=255)
                            ccAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCttgtt                                                                                                                           <  2:3257285/150‑1 (MQ=255)
                            ccAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGGTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCttgtt                                                                                                                           <  2:371477/150‑1 (MQ=255)
                              aGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGa                                                                                                                         <  2:3000517/150‑1 (MQ=255)
                              aGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGa                                                                                                                         >  1:567332/1‑150 (MQ=255)
                              aGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGGTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGa                                                                                                                         <  1:1332618/150‑1 (MQ=255)
                              aGTACGGGTAACCGAGCTGGGCCAGCGCTGCCGGCACCAGGTGACGGTCGTCCTCCGGCACCTGCAGCCCCGCCACCCCGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGa                                                                                                                         <  2:2220325/150‑1 (MQ=255)
                               gTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAg                                                                                                                        <  1:1605144/150‑1 (MQ=255)
                                tACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAgg                                                                                                                       <  2:1214775/150‑1 (MQ=255)
                                  cGGGTAACCGATCTTGGCCAGCGCGGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGaa                                                                                                                     <  1:4698290/150‑1 (MQ=255)
                                    ggTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGt                                                                                                                   <  1:2035664/150‑1 (MQ=255)
                                     gTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGACGtt                                                                                                                  >  1:4478335/1‑150 (MQ=255)
                                          cGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCga                                                                                                                                                      >  1:905295/1‑109 (MQ=255)
                                          cGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCga                                                                                                                                                      <  2:905295/109‑1 (MQ=255)
                                            aTCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGt                                                                                                           <  1:4249860/150‑1 (MQ=255)
                                            aTCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGt                                                                                                           <  2:2706136/150‑1 (MQ=255)
                                             tCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTg                                                                                                          >  2:3475266/1‑150 (MQ=255)
                                                  gCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCt                                                                                                     >  1:3217091/1‑150 (MQ=255)
                                                    cAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTgg                                                                                                   <  2:180421/150‑1 (MQ=255)
                                                     aGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGc                                                                                                  <  2:4597566/150‑1 (MQ=255)
                                                          tGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCCGCTTGTTGTGGAAGTTGAACAGTGCCCCTGGCCGCTc                                                                                             >  1:4921032/1‑150 (MQ=255)
                                                          tGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTc                                                                                             <  2:1166856/150‑1 (MQ=255)
                                                             cGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCTGCTCggg                                                                                          <  2:2506035/150‑1 (MQ=255)
                                                               gCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGa                                                                                        >  1:3104876/1‑150 (MQ=255)
                                                                  ccAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAAtt                                                                                     <  2:4168421/150‑1 (MQ=255)
                                                                  ccAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAAtt                                                                                     <  2:4483941/150‑1 (MQ=255)
                                                                      gTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGaa                                                                                 <  1:4558837/150‑1 (MQ=255)
                                                                      gTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGaa                                                                                 <  1:3792014/150‑1 (MQ=255)
                                                                      gTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGaa                                                                                 >  2:4759449/1‑150 (MQ=255)
                                                                       tGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAg                                                                                <  2:1591357/150‑1 (MQ=255)
                                                                           gTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCgca                                                                            >  2:1685646/1‑150 (MQ=255)
                                                                           gTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCgca                                                                            >  2:2348770/1‑150 (MQ=255)
                                                                           gTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCgca                                                                            >  2:944040/1‑150 (MQ=255)
                                                                           gTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCgca                                                                            >  1:4184110/1‑150 (MQ=255)
                                                                           gTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCgca                                                                            <  1:1433213/150‑1 (MQ=255)
                                                                           gTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGGTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCgca                                                                            >  2:1896729/1‑150 (MQ=255)
                                                                            ttCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCgcag                                                                           >  2:4902516/1‑150 (MQ=255)
                                                                             tCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCgcagc                                                                          <  1:2492118/150‑1 (MQ=255)
                                                                             tCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCgcagc                                                                          >  1:4296006/1‑150 (MQ=255)
                                                                                tcctccGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCAcc                                                                       >  2:1517816/1‑150 (MQ=255)
                                                                                  ctccGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTCGTGGACCAGCGAGATGTTCCAGCGCCCGCCCCGCGGGTTGAGGCAGTTGAACAGTGCCCCTGGCCGCTCGGGGACTTCGACGCGCAGCACCAt                                                                     >  2:4916036/1‑150 (MQ=255)
                                                                                  ctccGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCAt                                                                     >  2:3331405/1‑150 (MQ=255)
                                                                                  ctccGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCAt                                                                     <  1:3901075/150‑1 (MQ=255)
                                                                                   tccGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCAtt                                                                    >  2:750705/1‑150 (MQ=255)
                                                                                     cGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTtc                                                                  <  1:2177855/150‑1 (MQ=255)
                                                                                        cACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCAATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTtcgtc                                                               <  2:3256326/150‑1 (MQ=255)
                                                                                         aCCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTtcgtcg                                                              >  1:3710004/1‑150 (MQ=255)
                                                                                            tGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTg                                                           >  1:2177306/1‑150 (MQ=255)
                                                                                             gCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTgg                                                          <  1:4105599/150‑1 (MQ=255)
                                                                                             gCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTgg                                                          >  1:5027917/1‑150 (MQ=255)
                                                                                               aGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCa                                                        >  2:2107029/1‑150 (MQ=255)
                                                                                                  cccGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCg                                                     >  1:2449525/1‑150 (MQ=255)
                                                                                                  cccGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCg                                                     >  1:2449605/1‑150 (MQ=255)
                                                                                                   ccGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCgg                                                    >  2:4800384/1‑150 (MQ=255)
                                                                                                   ccGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCgg                                                    >  1:464861/1‑150 (MQ=255)
                                                                                                   ccGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCgg                                                    >  1:1652563/1‑150 (MQ=255)
                                                                                                    cGCCACCACGCGCCCATCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGc                                                   >  2:1315725/1‑150 (MQ=255)
                                                                                                     gCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGcc                                                  >  1:2475396/1‑150 (MQ=255)
                                                                                                       caccacGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGa                                                <  2:1060761/150‑1 (MQ=255)
                                                                                                        accacGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCCGCCGAg                                               >  1:5007501/1‑150 (MQ=255)
                                                                                                        accacGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTACAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCAATTCGTCGCTGGCACCGGCCGAg                                               >  1:5007527/1‑150 (MQ=255)
                                                                                                          cacGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTg                                             <  2:65004/150‑1 (MQ=255)
                                                                                                               gcCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGgccgc                                        >  2:2621198/1‑150 (MQ=255)
                                                                                                                   gTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTACCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATAGCGTCGCTGGCACCGGCCGAGAGGCCGCCGAc                                    >  1:3781998/1‑150 (MQ=255)
                                                                                                                         cgcgCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATg                              >  1:781100/1‑150 (MQ=255)
                                                                                                                          gcgcCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGa                             <  2:368213/150‑1 (MQ=255)
                                                                                                                           cgcCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGAc                            <  2:2166385/150‑1 (MQ=255)
                                                                                                                              cATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGACGaa                         >  2:462728/1‑150 (MQ=255)
                                                                                                                                    tGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGACGAATATGCa                   <  1:1733323/150‑1 (MQ=255)
                                                                                                                                    tGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGACGAATATGCa                   <  1:2159357/150‑1 (MQ=255)
                                                                                                                                          gctGGAACATCTCGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGACGAATATGCAGCTTgg             >  1:3165607/3‑150 (MQ=255)
                                                                                                                                          gctGGAACATCTCGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGACGAATATGCAGCTTgg             >  2:1643749/3‑150 (MQ=255)
                                                                                                                                           gaggAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGACGAATATGCAGCTTGGc            <  2:1199619/148‑1 (MQ=255)
                                                                                                                                           gTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGACGAATATGCAGCTTGGc            <  1:720942/150‑1 (MQ=255)
                                                                                                                                            tGGAACATCTCGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGa                                                <  1:414952/113‑1 (MQ=255)
                                                                                                                                            tGGAACATCTCGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGACGAATATGCAGCTTGGcc           >  1:2726355/1‑150 (MQ=255)
                                                                                                                                            tGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGACGAATATGCAGCTTGGcc           >  2:195999/1‑150 (MQ=255)
                                                                                                                                             ggAACATCTCGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGACGAATATGCAGCTTGGCCa          >  1:279083/1‑150 (MQ=255)
                                                                                                                                              gAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGACGAATATGCAGCTTGGCCaa         >  1:2418602/1‑150 (MQ=255)
                                                                                                                                              gAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGACGAATATGCAGCTTGGCCaa         >  2:2793707/1‑150 (MQ=255)
                                                                                                                                                aCATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGACGAATATGCAGCTTGGCCAAtt       >  1:3614114/1‑150 (MQ=255)
                                                                                                                                                  aTCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGACGAATATGCAGCTTGGCCAATTCg     >  2:2394654/1‑150 (MQ=255)
                                                                                                                                                    cGAGATGGTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGACGAATATGCAGCTTGGCCAATTCGtt   <  1:1850210/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTTTGTAGGCCGGGTTTTCGGTCTCGTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGACGAATATGCAGCTTGGCCAATTCGTTG  >  NC_002947/5876648‑5876946

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