Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,959,482 A→C 5.1% intergenic (+37/+50) ilvC → / ← ppiC ketol‑acid reductoisomerase, NAD(P)‑binding/peptidyl‑prolyl cis‑trans isomerase C (rotamase C)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,959,4820AC 5.1% 558.0 / 2.9 202intergenic (+37/+50)ilvC/ppiCketol‑acid reductoisomerase, NAD(P)‑binding/peptidyl‑prolyl cis‑trans isomerase C (rotamase C)
Reads supporting (aligned to +/- strand):  ref base A (73/111);  new base C (11/3);  total (88/114)
Fisher's exact test for biased strand distribution p-value = 9.04e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 229 total aligned reads displayed.

ACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGA  >  NC_000913/3959333‑3959629
                                                                                                                                                      |                                                                                                                                                   
aCGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCcaca                                                                                                                                                     <  2:3268133/150‑1 (MQ=255)
 cGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACag                                                                                                                                                    <  1:4147480/150‑1 (MQ=255)
 cGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACag                                                                                                                                                    <  1:977366/150‑1 (MQ=255)
 cGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACag                                                                                                                                                    <  1:1522556/150‑1 (MQ=255)
  gggCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGTGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGTTATGAAAGGTATTGCTGTTGCGGGTTAAGTGGGCGCTGATGGCCTCACCCCGTGCCTCTACCACagg                                                                                                                                                   >  2:3330811/1‑150 (MQ=255)
  gggCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGAGGGCCGCACCCCCACCCACACCCACagg                                                                                                                                                   >  2:206829/1‑150 (MQ=255)
      aaCTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACagggaga                                                                                                                                               <  1:106506/150‑1 (MQ=255)
      aaCTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACagggaga                                                                                                                                               <  1:3785815/150‑1 (MQ=255)
        cTGCGAGATGTGAACGAAGCGATACGCAGCCAAGCGATTGAGCAGGAAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAAAGTATTGCTGTTCCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACagggagagg                                                                                                                                             <  2:4215492/150‑1 (MQ=38)
        attcGTGATGTGAACTAAGCGATTCGCAGCAATGCGATTGAGCAGGTATGTAAGAAACTTCGCTGATATATGTCAGATATGAAACGTATTTCTGTTGCGGGTTAATTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACagggagagg                                                                                                                                             <  2:2868608/147‑1 (MQ=37)
          gCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACagggagaggga                                                                                                                                           <  1:351252/150‑1 (MQ=255)
               ttgtgAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACTGATTTGAAACGTATTGCTGTTGCGGTTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACCGGGAGAGGGAGaaaa                                                                                                                                      <  1:1319319/149‑1 (MQ=255)
               aTGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGaaaa                                                                                                                                      <  1:3445946/150‑1 (MQ=255)
               aTGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGaaaa                                                                                                                                      <  1:2824636/150‑1 (MQ=255)
                 gtgAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGATAGGGAGAAAaca                                                                                                                                    <  1:777406/150‑1 (MQ=255)
                 gtgAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGAGGCCCTCACCCCCACCCACACCCACAGGGGGTGGGAGaaaaaa                                                                                                                                    >  1:1456617/1‑148 (MQ=255)
                    aaCGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTc                                                                                                                                 <  2:3868407/150‑1 (MQ=255)
                    aaCGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTc                                                                                                                                 <  2:1755397/150‑1 (MQ=255)
                    aaCGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTc                                                                                                                                 <  2:3107776/150‑1 (MQ=255)
                    aaCGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGAGGCCCGCACCCCCCCCCGCGCCCCCCCGGCGCGGGCGAAAacacgc                                                                                                                                 >  2:3035947/1‑148 (MQ=255)
                     aCGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGAGGCCCTCACCCCGACCCACACCCCCAGGGAGGGGGAGAAAAcacaca                                                                                                                                >  1:2323017/1‑150 (MQ=255)
                      cGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCaa                                                                                                                               <  1:186058/150‑1 (MQ=255)
                       gAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAg                                                                                                                              <  1:3131332/150‑1 (MQ=255)
                       gAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCCACCCACACCCACAGGGAGAGGGAGTAAAAACACAAg                                                                                                                              >  1:2989321/1‑150 (MQ=39)
                       gAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGAGGGCCACACCCCCACCCACACCCACACGGAGAGGGAGGAAAAACTaaat                                                                                                                              >  1:3254865/1‑149 (MQ=255)
                       gAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGAGGCCCACACCCCCACCCACTCCCACACGGCGAGGGGGAAAACACACaac                                                                                                                              >  1:44433/1‑149 (MQ=255)
                         agagATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGc                                                                                                                            <  2:3737658/147‑1 (MQ=255)
                         aGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGc                                                                                                                            <  2:597473/150‑1 (MQ=255)
                          gCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCACACCCCCACCCTAACCCCCACCGAGATGGCGAAAAAACACAAGGcc                                                                                                                           >  2:449925/1‑150 (MQ=255)
                          gCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGAGGGCCGCACCCCCCCCCAAAACCACCGGGGGAGGGAGAAAACACTAAAGGcc                                                                                                                           >  2:449850/1‑150 (MQ=255)
                             aTTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTtc                                                                                                                        <  1:1596020/150‑1 (MQ=255)
                                cGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTcc                                                                                                                     <  1:4118941/150‑1 (MQ=255)
                                cGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTcc                                                                                                                     <  1:4034832/150‑1 (MQ=255)
                                 gCATCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGATAAAACACTCAAGGCCTTCTCCt                                                                                                                    <  2:2541144/150‑1 (MQ=255)
                                 gCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCCACCCGCACCCACACGGAGAGGGAGAAAACACACAATGCCGTCTccg                                                                                                                    >  1:3942068/1‑149 (MQ=38)
                                  cAGCCATGCGATTGAGTAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCACACCCCGACCCCCACCCCCACGGCCCGGGAGAAAACACTCAAGGCCTACTccgg                                                                                                                   >  2:3327091/1‑148 (MQ=255)
                                  cAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTg                                                                                                                   <  2:2425018/150‑1 (MQ=255)
                                    gCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGga                                                                                                                 <  1:2014949/150‑1 (MQ=255)
                                       aTGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGaa                                                                                                              <  1:1136879/150‑1 (MQ=255)
                                       aTGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGaa                                                                                                              <  1:2690561/150‑1 (MQ=255)
                                        tGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCGACACCCCCCCCACACCCACACCGCGAGGGCGAAAAAACACAAAGCCTTAACCAGGAGAAg                                                                                                             >  2:2933358/1‑150 (MQ=255)
                                         gCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAgg                                                                                                            <  1:2636121/150‑1 (MQ=255)
                                          cGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGc                                                                                                           <  2:2141224/150‑1 (MQ=255)
                                           gATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCAGGAGAAGGcc                                                                                                          <  1:1697065/150‑1 (MQ=255)
                                           gATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGAGGCCCGCACCCCCCCCCACACCCACCCGGAGAGGGCGGAAAAAAACAAGGCCCTCGACGGGAGAAgggc                                                                                                          >  1:40265/1‑148 (MQ=255)
                                            aTTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCGCACCCACAGGGAGAGGGAGAAAACACACAAGGCCATCTCCGGGAGCAGGccg                                                                                                         >  1:2767794/1‑149 (MQ=38)
                                            aTTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCCACCCTCAACCACAGGGAGAGGGAGAAAACACACAAGGCCATCTCCAGGAGAAGGCCt                                                                                                         >  2:2243330/1‑150 (MQ=38)
                                             ttGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCAGGAGAAGGCCtt                                                                                                        <  2:490661/150‑1 (MQ=255)
                                              tGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTg                                                                                                       <  2:1306830/150‑1 (MQ=255)
                                                aGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGGTAAGTGCGCGCTGAGGCCCTCACCCCGACCCACTCCCACAGGGTGGGGGAGAAAAAACGCAATGCCTTATCCGGGCGACGGCCGTGcg                                                                                                     >  1:2838079/1‑149 (MQ=255)
                                                 gCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTa                                                                                                    <  1:2066539/150‑1 (MQ=255)
                                                   aggtaggtAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCGCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTAtt                                                                                                  <  1:3224725/150‑1 (MQ=255)
                                                     gtaggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACCGGGAGAGGGCGAAAACACTCAAGGCCTTCTCCGGGAGTAGGCCTTTCTATTAg                                                                                                >  2:1362499/1‑150 (MQ=38)
                                                          taAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTAATAGTAGCg                                                                                           >  2:2750980/1‑150 (MQ=255)
                                                          taAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGAGGCCCGCCCCCCCCCCCGCACCCCCACGGGGGGGGGGAAAACACACAAGGCCGTATCCGGGGGCAGGCCGGGCGATTAGGTGCg                                                                                           >  1:1873978/1‑150 (MQ=255)
                                                           aaGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCgg                                                                                          <  2:3594418/150‑1 (MQ=255)
                                                           aaGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCgg                                                                                          <  1:3564995/150‑1 (MQ=255)
                                                            aGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCGCACCCCCACCCACACCCACAGGGAGTGGGCGAAAAAACACAAGGCCAACTCCAGGAGAAGGCCCTGCTAATAGTTGCGGt                                                                                         >  2:3361235/1‑150 (MQ=255)
                                                              aaaCTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCACACCCCCACCCGCTCCCACACGGGGAGGGAGAAAAAACACAAAGCCATCACCAGGAGAAAGCCGTGCGATTAGTTTCGGttc                                                                                       >  2:67035/1‑148 (MQ=255)
                                                              aaaCTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGAGGCCCTCACCCCGACCCGCACCCCCAGGGAGGGGGAGGAAAAAATAAAGGCCTTCAACGGGTGAATGCCGTGGGATTATTTGCGGTAc                                                                                       >  1:491141/1‑150 (MQ=255)
                                                                 cTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGc                                                                                    <  1:245778/150‑1 (MQ=255)
                                                                  tGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCa                                                                                   <  2:4121575/150‑1 (MQ=255)
                                                                     gcgGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGGCCGCACCCCCACCCACGACCACAGGGAGGGGGGGAAAAAACAAAAAGCCGTCTACAGGAGAAGGGCGGTCGATTAGTTGCGGGTCAGCAcca                                                                                >  1:3335039/1‑149 (MQ=255)
                                                                       ggCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCATATCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCttt                                                                              >  1:3065491/1‑149 (MQ=255)
                                                                       ggCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGAGGCCCTCACCCCGACCCACAACCACACGGAGAGGGAGAAAACACACAAGGCCGTCACCTGGAGGAGGCCGTGGGCGTAGTTGGGGGACAGCACCATa                                                                              >  1:2602866/1‑150 (MQ=255)
                                                                        gCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTaa                                                                             <  1:635361/150‑1 (MQ=255)
                                                                         cTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACCCTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAat                                                                            <  2:294250/150‑1 (MQ=255)
                                                                         cTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAATGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAat                                                                            <  2:2040744/150‑1 (MQ=255)
                                                                         cTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAat                                                                            <  1:3595724/150‑1 (MQ=255)
                                                                           atatGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGGCCGCACCCCCACCCACACCCACACCGAGAGGGTGAAAAAAAACAATTCCGAATACCGGGGAATGCCATGCTATTAGTTGGGGTTCAGCAccgacatgt                                                                          >  1:1788606/1‑143 (MQ=255)
                                                                                aCAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGa                                                                     <  2:121241/150‑1 (MQ=255)
                                                                                  aGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTAACCCCGACCCTCTCCCACAGGGCGAGGGAGAAAACACTCAAGGCCATCACCAGGAGAAGGCCGGTCTAATAGTTGACGTACAGCACCATAAAGATGTGaaa                                                                   >  2:4015301/1‑148 (MQ=255)
                                                                                     tatGAAACGTATTGCTGTTGCGGGTTTAGTGCGCGCTGATGCCCTCACCCCGACACTCTCCCACAGGGAGAGGGTGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATcc                                                                <  2:2888076/150‑1 (MQ=255)
                                                                                         aaaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCACACCCCCACCCACAACCACAGGGAGAGGGAGAAAAAACACAAGGCCATCTCCAGGAGAAGGCCGTGCTAATATTTTCGGTACAGCACCATAAAGATTTGATATACCAAc                                                            >  2:179454/1‑150 (MQ=255)
                                                                                         aaaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGAGGCCCTCACCCCCACCCACTCCCACACGGAGTGGGAGAAAACACTCAAGGCCTTCTCCGGGCGAAGGCCGTGCGATTAGTTGCCGTGCAGCACCATAATGAGGTGATAACCgccc                                                            >  2:2351662/1‑147 (MQ=255)
                                                                                            cGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTgg                                                         <  1:880512/150‑1 (MQ=255)
                                                                                              tATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGgt                                                       <  1:2884173/150‑1 (MQ=255)
                                                                                                ttGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCGCACCCCCACCCACACCCACAGGGAGAGGGAGAAAAAACACAAGGCCTTAACCGGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATACGAACTGGgtgt                                                     >  2:2716392/1‑150 (MQ=38)
                                                                                                       tGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCgg                                              <  2:1060450/150‑1 (MQ=255)
                                                                                                             ttAAGTGCGCGCTGATGCCCACACCCCCACCCAAACCCACAGGGAGATGGAGAAAAAAAAAAAGGCCGTAAACATGAGAAGGCCTTTCTATTAGTTTCGGTACAGCACCTTAAATAATTGATATACGAACTGGGTGTGCAGCGGGCCGGt                                        >  2:2976304/1‑150 (MQ=255)
                                                                                                                   gcgcgcTGAGGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTc                                  <  1:2933358/150‑1 (MQ=255)
                                                                                                                         tGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTAc                            <  2:1677645/150‑1 (MQ=255)
                                                                                                                          gATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTAcc                           <  1:2075058/150‑1 (MQ=255)
                                                                                                                             gCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGACAACACTCCAGGCCTTCTCCTGGAGAAGGCCTAGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGa                        <  2:1073041/150‑1 (MQ=255)
                                                                                                                               ccTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACGCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTCCAGCACCTTAAGGATGGGATATCCGACCTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACa                      >  2:2755239/1‑150 (MQ=39)
                                                                                                                                cTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACaa                     >  2:3627221/1‑150 (MQ=255)
                                                                                                                                    aCCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCTGGCCGGTCGGCTCCAGTACCGGACAagag                  >  2:397772/1‑150 (MQ=255)
                                                                                                                                      cccGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGaa                >  2:2674791/1‑150 (MQ=255)
                                                                                                                                         ctcccTCTCCCTGTGGGAGAGGGGGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAAc             >  1:1323904/3‑150 (MQ=255)
                                                                                                                                         ctcccTCTCCCTGTGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAAc             >  1:2400776/3‑150 (MQ=255)
                                                                                                                                           cccTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCa           <  2:1468139/150‑1 (MQ=255)
                                                                                                                                            ccTCTCCCTGTGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCAc          >  1:4182373/1‑150 (MQ=255)
                                                                                                                                            ccTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCAc          >  1:2676911/1‑150 (MQ=255)
                                                                                                                                                tcCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGATGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTa      >  2:1541178/1‑150 (MQ=255)
                                                                                                                                                 cccACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTAt     <  2:1078870/150‑1 (MQ=255)
                                                                                                                                                  ccACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATc    >  2:1095254/1‑150 (MQ=255)
                                                                                                                                                  ccACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATc    >  2:1914651/1‑150 (MQ=255)
                                                                                                                                                  ccACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATc    >  2:1834004/1‑150 (MQ=255)
                                                                                                                                                   cacaGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCg   >  1:3683976/1‑150 (MQ=255)
                                                                                                                                                    acaGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGa  >  1:3025786/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
ACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGA  >  NC_000913/3959333‑3959629

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: