Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,959,457 T→G 6.5% intergenic (+12/+75) 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,4570TG 6.5% 1474.5 / 25.1 515intergenic (+12/+75)ilvC/ppiCketol‑acid reductoisomerase, NAD(P)‑binding/peptidyl‑prolyl cis‑trans isomerase C (rotamase C)
Reads supporting (aligned to +/- strand):  ref base T (165/304);  new base G (39/7);  total (204/311)
Fisher's exact test for biased strand distribution p-value = 7.44e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 529 total aligned reads displayed.

GCTATTCCGGAAGGCGCGGTAGATAACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGG  >  NC_000913/3959308‑3959606
                                                                                                                                                     |                                                                                                                                                     
gCTATTCCGGAAGGCGCGGTAGATAACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGAt                                                                                                                                                       >  1:1308899/1‑150 (MQ=255)
  tATTCCGGAAGGCGCGGTAGATAACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGc                                                                                                                                                     <  2:6508795/150‑1 (MQ=255)
    ttCCGGAAGGCGCGGTAGATAACGGGCAACTGCGTGATGTGAACGAAGCTATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGccc                                                                                                                                                   >  1:3465757/1‑150 (MQ=255)
       cGGAAGGCGCGGTAGATAACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCa                                                                                                                                                >  1:1637566/1‑150 (MQ=255)
               gcgGTAGATAACGGGCAACTGCGTGTTGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGAcc                                                                                                                                        <  1:2611414/150‑1 (MQ=255)
               gcgGTAGATAACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGAcc                                                                                                                                        <  1:7670068/150‑1 (MQ=255)
               gcgGTAGATAACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGAcc                                                                                                                                        <  1:6891096/150‑1 (MQ=255)
                cgGTAGATAACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTCAGTGCGCGCTGATGCCCACACCCCGAccc                                                                                                                                       >  1:2540784/1‑150 (MQ=255)
                cgGTAGATAACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGAccc                                                                                                                                       >  1:4975379/1‑150 (MQ=255)
                  gTAGATAACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCACACCCCGAcccac                                                                                                                                     >  2:2333528/1‑148 (MQ=255)
                    aGATAACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCtctc                                                                                                                                   <  2:8992903/150‑1 (MQ=255)
                    aGATAACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCACACCCCCAcccgcac                                                                                                                                   >  2:8732226/1‑146 (MQ=255)
                     gATAACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTcc                                                                                                                                  <  2:7772098/150‑1 (MQ=255)
                     gATAACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTcc                                                                                                                                  <  1:5962094/150‑1 (MQ=255)
                     gATAACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGAGGCCCTCACCCCGccccacgac                                                                                                                                  >  2:2819010/1‑145 (MQ=255)
                         aCGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCcaca                                                                                                                              <  1:8958608/150‑1 (MQ=255)
                          cGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTCAGTGAGCGCTGATGCCCTCACCCCGACCCTCTCCCACag                                                                                                                             >  1:6605742/1‑150 (MQ=255)
                          cGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACag                                                                                                                             <  2:500902/150‑1 (MQ=255)
                          cGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCGCCCCCCCACCCACAccccccc                                                                                                                             >  2:6970534/1‑146 (MQ=255)
                          cGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGAGGCCCTCACCCCCACCCACTCCCACag                                                                                                                             >  1:1763349/1‑150 (MQ=255)
                           gggCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTccctgtgg                                                                                                                            <  2:4671824/150‑6 (MQ=255)
                               aaCTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACGGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACagggaga                                                                                                                        <  1:3880323/150‑1 (MQ=255)
                               aaCTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACagggaga                                                                                                                        <  1:9957356/150‑1 (MQ=255)
                               aaCTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCACTCCCACagggaga                                                                                                                        >  1:5844835/1‑150 (MQ=255)
                               aaCTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCACACCCCGACCCACACCCACagggcgg                                                                                                                        >  1:7777567/1‑147 (MQ=255)
                               aaCTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGAGATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCACTCCCACagggaga                                                                                                                        <  2:8492039/150‑1 (MQ=255)
                                  tGCGTGATGTGAACGAAGCGATTCGCTGCCATGCGATTGAGCAGGTAGTTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACagggagaggg                                                                                                                     <  1:7903337/150‑1 (MQ=255)
                                  tGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCTGtgggagaggg                                                                                                                     <  2:3078065/150‑1 (MQ=255)
                                  tGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTAGATGACAGAGATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCGCggggagaggg                                                                                                                     >  1:6096846/1‑150 (MQ=255)
                                      tGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGaa                                                                                                                 <  1:154171/150‑1 (MQ=255)
                                        aTGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGaaaa                                                                                                               <  1:4286985/150‑1 (MQ=255)
                                        aTGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGaaaa                                                                                                               <  2:9451598/150‑1 (MQ=255)
                                             aaCGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCCCCGGGGGAGGGGGAACACACTc                                                                                                          >  1:9057534/1‑150 (MQ=255)
                                             aaCGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCCCAGGGCGAGGGAGTAAacacac                                                                                                          >  2:275308/1‑148 (MQ=255)
                                             aaCGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTc                                                                                                          >  2:5154924/1‑150 (MQ=255)
                                             aaCGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGATGCGGGTTAATTTAGCGCTGATTCCAACAACACCACACACTCCCACAGGCAGCGGGAGAAAACACta                                                                                                          >  2:3235457/1‑149 (MQ=255)
                                               cGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCaa                                                                                                        <  1:8175309/150‑1 (MQ=255)
                                                gAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGGCCTCACCCCGACCCACACCCACAGGGAGAGGGAGAAAACACACAAg                                                                                                       >  1:5516105/1‑150 (MQ=255)
                                                gAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACACAAg                                                                                                       >  2:7453393/1‑150 (MQ=255)
                                                gAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGAGGCCCCCACCCCCACACGCACCCAAAGGGAGGGGGAGAAAacacaaaaa                                                                                                       >  2:1741678/1‑145 (MQ=255)
                                                gAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGGTAAGTGCGCGCTGATTCCCACACCCCCACCCTCAACCACCCGGAGAGGGAGAAAACACACAAg                                                                                                       >  2:9860178/1‑150 (MQ=255)
                                                  aGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGc                                                                                                     <  1:8554020/150‑1 (MQ=255)
                                                  aGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGc                                                                                                     >  2:6345912/1‑150 (MQ=255)
                                                        tCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTtctc                                                                                               <  2:2952901/150‑1 (MQ=255)
                                                        tCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCACAACCCCACCCAATCCCCCACGGGGGGGGTGAAAAAACACAAGGccgtatc                                                                                               >  1:3658297/1‑145 (MQ=255)
                                                         cGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGACACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTcc                                                                                              <  2:2127659/150‑1 (MQ=255)
                                                          gCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCACACCCCGACCCACTACCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTccg                                                                                             >  2:9352107/1‑149 (MQ=255)
                                                           cAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACGCAAGGCCTTCTccag                                                                                            <  2:6272867/150‑3 (MQ=255)
                                                             gCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGga                                                                                          <  1:1844312/150‑1 (MQ=255)
                                                             gCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGga                                                                                          <  2:341914/150‑1 (MQ=255)
                                                              ccATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGGCCGCACCCCGACCCACTCCCACCGGGAGAGGGAGAAAAAAAAAAAGGCCTTCTGCTgggg                                                                                         >  1:7126895/1‑148 (MQ=255)
                                                                 tGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGa                                                                                                                                          <  2:10112965/98‑1 (MQ=255)
                                                                 tGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGAGGCCCACACCCCGACCCACACCCACAGGGAGAGGGCGAAAACACTCAAGGCCTTCTCCGGGAGAAg                                                                                      >  1:7454953/1‑150 (MQ=255)
                                                                  gCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAgg                                                                                     <  1:2935798/150‑1 (MQ=255)
                                                                  gCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCACACCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAgg                                                                                     >  2:10286808/1‑150 (MQ=255)
                                                                     aTTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCt                                                                                  <  2:828606/150‑1 (MQ=255)
                                                                     aTTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCACAACCCGACCCACACCCACCGGGCGAGGGGGAAAACACACAACGCCGTCTACAGGAGAAGGccg                                                                                  >  1:440155/1‑149 (MQ=255)
                                                                       tGAGCAGGTAGGTAAGACACTGCGCGGCTATATGACCGAAATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTg                                                                                <  1:5568823/150‑1 (MQ=255)
                                                                       tGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTAGTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCCAGGCCTTCTCCTGGAGAAGGCCTTg                                                                                <  1:4543250/150‑1 (MQ=255)
                                                                        gAGCAGGTAGGTAAGTAACTGCGCGGCTATATGACAGATATGAAAGGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTGCCACAGGGAGAGGGAGAAAACACTCACGGCGTTCTCCTGGAGAAGGCCTTGc                                                                               <  2:5872835/150‑1 (MQ=38)
                                                                        gAGCAGGTAGGTAAGACACTGCGCGGCTATAGGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCCCAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGc                                                                               <  1:3122183/150‑1 (MQ=255)
                                                                            aggtaggtAAGAAACTGCGCGGCTATATGACAGATAGGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACCGGGAGAGGGAGAAAACACTCAAGGCCTTCGCCAGGAGAAGGCCTTgagtgt                                                                           <  2:990670/150‑6 (MQ=38)
                                                                                 ggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTg                                                                      <  1:324094/150‑1 (MQ=255)
                                                                                   taAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCg                                                                    <  1:6497968/150‑1 (MQ=255)
                                                                                   taAGAAACTGCGCGGCTATATGACAGATATGAAACGTAGTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCg                                                                    <  1:3993775/150‑1 (MQ=255)
                                                                                    aaGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCGCACCCCCACCCACGCCCCCAGGGGGGGGGGGAAAAAACACAAGGCCGTAGCCCGGGGGAGGCCGTGCGAATAGTTTCgg                                                                   >  2:2045894/1‑150 (MQ=255)
                                                                                      gAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTa                                                                 <  2:1624803/150‑1 (MQ=255)
                                                                                      gAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTa                                                                 <  1:7232549/150‑1 (MQ=255)
                                                                                       aaaCTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGAGGCCCTCACCCCGACCCTCTACCACAGGGCGAGGGAGAAAACACTCAAGGCCTTCTCCGGGAGAAGGCCGTGCTGTTATTTgacggac                                                                >  2:5601842/1‑144 (MQ=255)
                                                                                        aaCTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACa                                                               <  2:6676414/150‑1 (MQ=255)
                                                                                           tGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCAAACCCCCCCCCACACCCCCCCGGCGGGGGGGGAAACACGCAAAGCCCTATACGGGGGGAGGGCGTGGGATTATTTGCGGGACAGcc                                                            >  2:3855747/1‑149 (MQ=255)
                                                                                            gcgcGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGTGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCAc                                                           <  2:3753052/150‑1 (MQ=255)
                                                                                              gctgCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGAGGCCCTCACCCCCACCCAATACCACAGGGGGAGGGGGTAAAAACACAAGGCCATAACCTGGTGAATGCCATGCTATTAGTTTGGGTACAGCACCt                                                         >  1:6752842/4‑150 (MQ=255)
                                                                                                ggCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGGCCACACCCCCACCCACACCCACAGGGAGAGGGGGAAAAAACACAAAGCCAAATACGGGAGAAGGCCGTGCTATTTGTTGGGGTACAGCAccatt                                                       >  1:6971890/1‑147 (MQ=255)
                                                                                                 gCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGCGAGGGGGAAAACACACAAGGCCGTCTCCGGGAGAAGGCCGTGCTAATAGGTGCGGGACAGCACCTTaa                                                      >  2:1373986/1‑150 (MQ=37)
                                                                                                  cTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAat                                                     <  1:209618/150‑1 (MQ=255)
                                                                                                       tGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGAtgt                                                <  2:7889426/150‑1 (MQ=255)
                                                                                                         aCAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGGCCACACCCCCACCCACAACCACACGGAGAGGGAGAAAAAACAAAAGGCCATCACCAGGAGAAGGCCGTTCTATTAGTTGCCGTGCAGCACCATAAATATGTGa                                              >  2:915640/1‑150 (MQ=255)
                                                                                                         aCAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGAGGCCCACACCCCCACCCACTCCCCCAGCGGGAGGGAGTAAAAACAAAACCCCCTAGACAGGAGAAGGCCGTGCTATTATTTGGGGGTCAGCACCAAAATGAtgtgt                                              >  1:4740195/1‑149 (MQ=255)
                                                                                                          cAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGat                                             <  1:957030/150‑1 (MQ=255)
                                                                                                             atatGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATc                                          <  1:1211675/150‑1 (MQ=255)
                                                                                                                  aaaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCACACCCCGACCCACACCCACAGGGTGAGGGAGAAAAAACACAATGCCGTCTCCTGGGGCAGGCCGGGCGATTAGTTTCGGTGCAGCACCCGAATGATGTGGTATAAGAAc                                     >  1:2371458/1‑150 (MQ=255)
                                                                                                                     cGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTgg                                  <  1:10006957/150‑1 (MQ=255)
                                                                                                                     cGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTgg                                  <  2:8581593/150‑1 (MQ=255)
                                                                                                                        aTTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGgtg                               <  1:8226660/150‑1 (MQ=255)
                                                                                                                            cTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACCGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCa                           <  1:10215174/150‑1 (MQ=255)
                                                                                                                            cTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCACGGCCGTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCa                           <  2:3992409/150‑1 (MQ=255)
                                                                                                                            cTGTTGCGGGTTAAGTGCGCGCTGATGCCCACACCCCGACCCACTCCCACAGGGCGGGGGAGAAAAAACACAATGCCGAATCCTGGAGAAGGCCTTGCGATTAGTTGCGGTACAGCCCCGTAAAGATGTGATATCCGAACTGGGTGTGCa                           >  2:9760864/1‑150 (MQ=255)
                                                                                                                               ttGCGTGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAATGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCg                        <  1:5387128/150‑1 (MQ=255)
                                                                                                                                tGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCgg                       <  2:7134539/150‑1 (MQ=255)
                                                                                                                                  cGGGTTAAGTGCGCGCTGATGCCCTCACCCCCACCCACTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAATGCCATGCTATTAGTTGCGGTACAGCACCTTAATGAAGTGATATCCGAACAGGGTGTGCAGCGGGc                     >  2:1497118/1‑150 (MQ=38)
                                                                                                                                         aGTGCGCGCTGAGGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCCAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCgg              <  2:7392718/150‑1 (MQ=255)
                                                                                                                                          gTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGc             <  1:2967865/150‑1 (MQ=255)
                                                                                                                                           tGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCt            <  2:8098149/150‑1 (MQ=255)
                                                                                                                                                 cTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTa      <  1:8781499/150‑1 (MQ=255)
                                                                                                                                                  tGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTAc     <  2:683277/150‑1 (MQ=255)
                                                                                                                                                   gATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTAcc    <  2:4435824/150‑1 (MQ=255)
                                                                                                                                                   gATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACGGGGTGTGCAGCGGGCCGGTCGGCTCCAGTAcc    >  2:5416026/1‑150 (MQ=255)
                                                                                                                                                     tGCCCTCACCCCGACCCTCTCCCACCGGGAGAGGGAGAAAACACTCGAGGCCTTCTCCTGGAGAAGGCCGTGCTATTAGGTGCGGTCCAGCACCTTAATGATGTGATAACCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCgg  >  2:762783/1‑150 (MQ=255)
                                                                                                                                                     tGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCgg  >  2:1238299/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCTATTCCGGAAGGCGCGGTAGATAACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGG  >  NC_000913/3959308‑3959606

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: