Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,959,482 A→C 6.3% 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 6.3% 1270.4 / 14.4 459intergenic (+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 (144/272);  new base C (27/12);  total (173/286)
Fisher's exact test for biased strand distribution p-value = 4.04e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 544 total aligned reads displayed.

ACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAAC  >  NC_000913/3959333‑3959631
                                                                                                                                                     |                                                                                                                                                     
aCGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCcaca                                                                                                                                                       <  1:7290938/150‑1 (MQ=255)
 cGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCACACCCCCACCCACACCCACag                                                                                                                                                      >  2:540423/1‑150 (MQ=255)
     cAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACagggag                                                                                                                                                  <  2:4357942/150‑1 (MQ=255)
      aaCTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGTGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACagggaga                                                                                                                                                 <  1:6462743/150‑1 (MQ=255)
      aaCTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACagggaga                                                                                                                                                 <  2:9170323/150‑1 (MQ=255)
      aaCTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACagggaga                                                                                                                                                 <  1:5030236/150‑1 (MQ=255)
        cTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCTGtgggagagg                                                                                                                                               <  1:307237/150‑1 (MQ=255)
        cTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACagggagagg                                                                                                                                               <  1:8611300/150‑1 (MQ=255)
          gCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACagggagaggga                                                                                                                                             <  1:9733217/150‑1 (MQ=255)
           cGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGAGGCCCACACCCCCACCCACACCCCCacggggagggag                                                                                                                                            >  2:2061252/1‑150 (MQ=255)
            gTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCACACCCCCAACCACAACCACAcggcgagggaga                                                                                                                                           >  2:9006908/1‑150 (MQ=255)
              gATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGaaa                                                                                                                                         <  2:5713255/150‑1 (MQ=255)
              attaagAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGaaa                                                                                                                                         >  2:8296879/6‑150 (MQ=255)
                  tgAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAacac                                                                                                                                     <  1:10685100/150‑1 (MQ=255)
                  tgAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCACCCACCCGGAGAGGGAGAAAacac                                                                                                                                     >  2:10099194/1‑150 (MQ=255)
                    aaCGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTc                                                                                                                                   <  2:2909955/150‑1 (MQ=255)
                    aaCGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTc                                                                                                                                   <  1:10422231/150‑1 (MQ=255)
                      cGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCaa                                                                                                                                 <  2:7573155/150‑1 (MQ=255)
                        aaGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAgg                                                                                                                               <  2:1830679/150‑1 (MQ=255)
                        aaGCGAGTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAgg                                                                                                                               <  2:2755927/150‑1 (MQ=255)
                         aGCGATTCGCAGCCAGGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTTATGCCCTCACCCCGACCCTCTCCCACATGGAGAGGGAGAAAACACTCAAGGc                                                                                                                              <  1:5112990/150‑1 (MQ=255)
                          gCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGcc                                                                                                                             <  1:843234/150‑1 (MQ=255)
                            gATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGGTAAGTGCGCGCTGAGGCCCGCACCCCCCCCCACGCCCCCCGGGCGCGGGCGCACACACACCAGGccag                                                                                                                           >  2:4879225/1‑148 (MQ=255)
                              ttCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCGCACCCCCACCCTCAACCACAGGGGGAGGGCGCAAACACACAAGGCCAtcg                                                                                                                         >  2:5258504/1‑149 (MQ=255)
                                 gCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCt                                                                                                                      <  1:9896130/150‑1 (MQ=255)
                                 gCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCGCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTccg                                                                                                                      >  1:1445094/1‑149 (MQ=255)
                                   aGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCACTCCCACAGGGGGAGGGAGAAAACACTCAAGGCCTTCTACTgg                                                                                                                    >  1:2804129/1‑150 (MQ=255)
                                   aGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCACACCCACAGGGAGAGGGAGAAAACACACAAGGCCTTCTCCggg                                                                                                                    >  2:4612811/1‑150 (MQ=255)
                                     ccATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACACGGAGAGGGCGAAAAAAATCAAGGCCATAACCTGgag                                                                                                                  >  1:6662773/1‑150 (MQ=38)
                                       aTGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTTATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTACTGGAGaa                                                                                                                <  2:6216516/150‑1 (MQ=255)
                                       aTGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGaa                                                                                                                <  1:3938081/150‑1 (MQ=255)
                                         gCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGACACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAgg                                                                                                              <  2:9817878/150‑1 (MQ=255)
                                          cGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCTCAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGc                                                                                                             <  2:9325630/150‑1 (MQ=255)
                                           gATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGcc                                                                                                            <  1:3556283/150‑1 (MQ=255)
                                            aTTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCAGGAGAAGGCCt                                                                                                           <  2:11014437/150‑1 (MQ=255)
                                            aTTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCAGGAGAAGGCCt                                                                                                           <  1:2035226/150‑1 (MQ=255)
                                             ttGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCtt                                                                                                          <  1:10960104/150‑1 (MQ=255)
                                                   aggtaggtAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCAGGAGAAGGCCTTgagtgt                                                                                                    <  2:10276933/150‑6 (MQ=255)
                                                   aggtaggtAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACGCAAGGCCTTCTCCAGGAGAAGGCCTTgagtgt                                                                                                    <  2:2404470/150‑6 (MQ=255)
                                                    ggtaggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTa                                                                                                   <  1:5391621/150‑1 (MQ=255)
                                                    ggtaggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTa                                                                                                   <  2:6984827/150‑1 (MQ=255)
                                                     gtaggtaAGAAACTGTGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGAGGCCCGCACCCCCCCCCACACCCACCGGGCGAGGGAGACAACACCCAAGGCCATCACCGGGGGAAGGCCGTGCGAGTAg                                                                                                  >  1:2432076/1‑150 (MQ=255)
                                                     gtaggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAg                                                                                                  <  1:6008964/150‑1 (MQ=255)
                                                     gtaggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCCACCCTCACCCACAGGGAGAGGGAGAAAAAACTCAAGGCCGTCTCCTGGAGAAGGCCATGCTATTAg                                                                                                  >  2:6790440/1‑150 (MQ=39)
                                                     gtaggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCGCACCCCCACCCACTCCCCCAGGGAGGGGGAGAAAACACACAAGGCCATCTCCGGGAGAAGGCCGTGCGATTAg                                                                                                  >  1:7408577/1‑150 (MQ=255)
                                                     gtaggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCGCACCCCCACCCAAGCCCAAACGGGGTGGGAGAAAACACACAATGCCTTCTCCTGGAGAAGGGCGTGCTAtttt                                                                                                  >  1:1133242/1‑148 (MQ=255)
                                                       aggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGtt                                                                                                <  1:2262244/150‑1 (MQ=255)
                                                       aggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGtt                                                                                                <  1:1699773/150‑1 (MQ=255)
                                                         gtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCCACCCTCACCCACAGGGAGAGGGAGAAAACACTCAAGGCCATCTCCTGGAGAAGGCCATGCTATTAGTTGc                                                                                              >  2:7112957/1‑150 (MQ=255)
                                                             gAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCACACCCCCACCCACACCCCCCCGGAGAGGGGGGAAACACACAAGGCCATAGACAGGGGAAGGCCCGGCTAATAGTTggggga                                                                                          >  1:6243445/1‑148 (MQ=255)
                                                                   gcgcGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCGCAACCCCACCCTCTCCCACAGGGAGAGGGCGAAAAAACTCAAAGCCTTATACTGGAGGAGGCCATGCTATTAGTTGCGGttaatcac                                                                                    >  1:7729102/1‑143 (MQ=37)
                                                                     gcgGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCCACCCACACCCACACGGAGAGGGAGAAAACACACAATGCCTTCTCCAGGCGAAGGCCGTGCTATTAGTTGCGGTACAGCACCt                                                                                  >  1:10002923/1‑150 (MQ=38)
                                                                        gCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGTGAGGGAGAAAACACTCAAGGCCTTCTCCTTGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTaa                                                                               <  2:7537013/150‑1 (MQ=255)
                                                                          tataTGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAatg                                                                             <  2:9850094/150‑1 (MQ=255)
                                                                          tataTGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCACACCCCCACCCTCTCCCACACGGAGAGCGAGAAAACACTCAAAGCCAAATACAGGAGAAGGCCGAGCTATAACTAGCCCTACAGCACATTaaag                                                                             >  2:8162924/1‑148 (MQ=255)
                                                                            tatGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAatgat                                                                           <  1:2678595/150‑1 (MQ=255)
                                                                            tatGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAatgat                                                                           <  2:3488843/150‑1 (MQ=255)
                                                                             atGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAatgatg                                                                          <  2:4073730/150‑1 (MQ=255)
                                                                               gACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGAtgtg                                                                        <  1:8252989/150‑1 (MQ=255)
                                                                               gACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCGCACCCCCCGGGAGAGGGAGAAAACACTAAAGGCCTTCTCCAGGGGAAGGCCATGCAAGTAGTTGCGGTACAGCACCATAAAGAtttg                                                                        >  2:7812604/1‑150 (MQ=37)
                                                                                aCAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCGAACCCCGACCCACACCCCCAGGGAGAGGGAGAAAAAAAAAAAGGCCATCTCCTGGAGAAGGCCGTGCTATTAATTGCGGTACAGCACCATAATTATGTGa                                                                       >  2:4472346/1‑150 (MQ=255)
                                                                                  aGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGAGGCCCGCACCCCGACCCAAACCCCCACGGAGAGGGTGAAAAAACACAAGGCCATCTCCAGGTGCAGGCCGGGCTATTAGTTTCCGTTCCGCACCATAATTAAGTGata                                                                     >  1:6791209/1‑150 (MQ=255)
                                                                                    atatGAATCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATc                                                                   <  2:5943586/150‑1 (MQ=255)
                                                                                    atatGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCACACCCCCACCCAAAACCACAGGGGGATGGAGAAAAAAAAAAAGGCCATATACGGGAGGATGCCGTGGGATTATTTTCGGTTCAGCACCATAAATATGTtataaa                                                                   >  1:2414190/1‑148 (MQ=255)
                                                                                         aaaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAAc                                                              <  2:7431512/150‑1 (MQ=255)
                                                                                            cGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTgg                                                           <  1:1780497/150‑1 (MQ=255)
                                                                                            cGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTgg                                                           <  1:4742024/150‑1 (MQ=255)
                                                                                            cGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTgg                                                           <  2:8468205/150‑1 (MQ=255)
                                                                                            cGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTGCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTgg                                                           <  1:1320132/150‑1 (MQ=255)
                                                                                             gTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCACACCCCCACCCACAACCCCAGGGCGAGGGAGAAAAAACACAAGGCCATCTACAGGAGAAGGCCTTGCTATTAGTTTCGGTACACCACCATAATGATTTTATATACGaaaaggg                                                          >  1:7970939/1‑145 (MQ=255)
                                                                                               aTTGCTGTTGCGGGTTAAGTGCGCGCTGAGGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTCGTTGCGGGACAGCACCTTAATGATGTGATATCCGAACTGGgtg                                                        <  2:8668056/150‑1 (MQ=255)
                                                                                                ttGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGgtgt                                                       <  1:5775998/150‑1 (MQ=255)
                                                                                                ttGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCGCACCCCGACCCACTACCACACGGCGAGGGCGCAAACACACAAGGCCGTCTCCAGGCGAACGCCCTGCAATTAGTTGCGGGACCGCACCTTAAAGTTGTGATATACGAAATGGgtgg                                                       >  2:2530484/1‑149 (MQ=255)
                                                                                                 tGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGACAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGtgtg                                                      <  1:4513845/150‑1 (MQ=255)
                                                                                                      ttGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCg                                                 <  2:10527428/150‑1 (MQ=255)
                                                                                                           ggTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCg                                            <  1:4306146/150‑1 (MQ=255)
                                                                                                             ttAAGTGCGCGCTGATGCCCTCACCCCGACCCACAACAACAGGGGGGGGGTGAAAACACTCAAGGGCCTCTGCTTGAGGGTGCCTTGCGGTTTGGTGGGGTCCCGCACCATAAGTGTGTGGTATGCGTACTGGGGGGGCAGGGGGcccgg                                          >  2:390841/1‑147 (MQ=255)
                                                                                                                aGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGATAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCgg                                       <  2:159911/150‑1 (MQ=255)
                                                                                                                   gcgcgcTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTc                                    <  1:556764/150‑1 (MQ=255)
                                                                                                                     gcgcTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCa                                  <  1:8720602/150‑1 (MQ=255)
                                                                                                                        cTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTa                               <  1:1018681/150‑1 (MQ=255)
                                                                                                                           aTGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCg                            >  2:3340072/1‑150 (MQ=255)
                                                                                                                                cTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACaa                       <  2:6092238/150‑1 (MQ=255)
                                                                                                                                cTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACaa                       <  1:9061525/150‑1 (MQ=255)
                                                                                                                                 tCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAag                      <  1:3645535/150‑1 (MQ=255)
                                                                                                                                  cACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAaga                     <  2:1074994/150‑1 (MQ=255)
                                                                                                                                  cACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAaga                     <  1:10070681/150‑1 (MQ=255)
                                                                                                                                   aCCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAagag                    >  2:9452899/1‑150 (MQ=255)
                                                                                                                                     cccGACCCTCTCCCACAGTGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGaa                  <  2:5042352/150‑1 (MQ=255)
                                                                                                                                     cccGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGAGGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGaa                  >  2:1346977/1‑150 (MQ=255)
                                                                                                                                        gACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAAc               <  2:5781408/150‑1 (MQ=255)
                                                                                                                                        gACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAAc               >  2:5516384/1‑150 (MQ=255)
                                                                                                                                          cccTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCa             <  1:1567627/150‑1 (MQ=255)
                                                                                                                                           ccTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCAc            >  2:1066615/1‑150 (MQ=255)
                                                                                                                                                cccACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTAt       <  2:7182928/150‑1 (MQ=255)
                                                                                                                                                 ccACAGGGAGAGGGATAAAACTCTCATGGCCTTCTCCTGGAGAAGGCCTTGATATTAGTTGCGGTACAGCACGTTTATGATGTGATATCCGTACTGGGTGTGCAGCGGGCCTGTCGTCTCCAGTACCGGAAAATAGATAACAACTTTATc      <  2:6583936/150‑1 (MQ=255)
                                                                                                                                                 ccACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCCGTACCGGACAAGAGAAAACCACTGTATc      >  2:5070568/1‑150 (MQ=255)
                                                                                                                                                 ccACAGAGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATc      >  2:7363676/1‑150 (MQ=255)
                                                                                                                                                    caGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGaa   <  2:856929/150‑1 (MQ=255)
                                                                                                                                                    caGGGAGAGGGAGAAAACACTCAAGGCCTGCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGTACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGaa   <  2:8004324/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAAC  >  NC_000913/3959333‑3959631

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: