Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,959,482 A→C 7.9% 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 7.9% 964.9 / 15.7 348intergenic (+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 (128/179);  new base C (27/8);  total (161/187)
Fisher's exact test for biased strand distribution p-value = 1.07e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 402 total aligned reads displayed.

CGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAAC  >  NC_000913/3959334‑3959631
                                                                                                                                                    |                                                                                                                                                     
   gCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCCCACCCACcgggg                                                                                                                                                   >  2:6337208/1‑149 (MQ=255)
    cAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACagggag                                                                                                                                                  <  1:1363799/150‑1 (MQ=255)
     aaCTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCGCACCCCCACCCAAACCCACacggcga                                                                                                                                                 >  2:5058687/1‑150 (MQ=255)
     aaCTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGAGGCCCTCACCCCGACCCTCTCCCACagggaga                                                                                                                                                 >  1:7001277/1‑150 (MQ=255)
       cTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCTGtgggagagg                                                                                                                                               <  2:2155797/150‑1 (MQ=255)
       cTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACagggagagg                                                                                                                                               <  1:5837585/150‑1 (MQ=255)
       cTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACagggagagg                                                                                                                                               <  1:6029286/150‑1 (MQ=255)
       cTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCCCACCCCGACCCACTCCCACagggcgagg                                                                                                                                               >  1:3464203/1‑150 (MQ=255)
         gCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGAGGCCCTCACCCCGACCCACTCCCACacggcgggggc                                                                                                                                             >  2:5171909/1‑149 (MQ=255)
             gATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGaaa                                                                                                                                         <  2:2833421/150‑1 (MQ=255)
             gATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGaaa                                                                                                                                         <  1:5214038/150‑1 (MQ=255)
              aTGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCACACCCCCACCCACACCCACACGGAGAGGGGGaaaa                                                                                                                                        >  2:2068352/1‑150 (MQ=255)
                gtgAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCCCACCCCGACCCTCTCCCACCGGGAGAGGGAGAAAaca                                                                                                                                      <  1:6128950/150‑1 (MQ=255)
                gtgAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGAGTCCCTCACCCCGACCCACTCCCCCAGGGGGCGGGGGaaaaaa                                                                                                                                      >  1:5619926/1‑148 (MQ=255)
                 tgAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTACGTGCGCGCTGATGCCCTCACCCCGACCCACACCCCCCGGGCGGGGGCGAAAacac                                                                                                                                     >  1:6247509/1‑150 (MQ=255)
                 tgAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGTGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAacac                                                                                                                                     <  2:1854465/150‑1 (MQ=255)
                   aaCTAAGCGATTCGCATCCTTGCGATTGATCAGGTAGGTAAGAAACTGCGCGGATATATGAATGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCTACCCTCTCCCACAGGGAGAGGGAGAAAACACTc                                                                                                                                   <  1:6382373/150‑1 (MQ=37)
                   aaCGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTc                                                                                                                                   <  1:2792473/150‑1 (MQ=255)
                   aaCGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGAGGCCCACACCCCCACCCACGACCCCCGGGCGGGGGCGAAAacccaa                                                                                                                                   >  1:3233593/1‑146 (MQ=255)
                    aCGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCACTCCCACAGGGAGAGGGAGAAAAACCTCa                                                                                                                                  >  2:4099512/1‑150 (MQ=255)
                    aCGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGAGGCCCGCACCCCCCCCCCCTCCCCCCGGGCGCGGGGGACAAAAcaca                                                                                                                                  >  2:698895/1‑150 (MQ=255)
                     cGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCaa                                                                                                                                 <  2:1254923/150‑1 (MQ=255)
                     cGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCaa                                                                                                                                 <  1:23236/150‑1 (MQ=255)
                      gAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTCAGTGCGCGCTGAGGCCCACACCCCCCCCCCCTCCCCCCGGGAGCGGGAGAAAACCCACAAg                                                                                                                                >  1:1179531/1‑150 (MQ=255)
                      gAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAg                                                                                                                                <  2:2345997/150‑1 (MQ=255)
                      gAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCACGCCCACAGGGAGAGGGAGAAAACACTCAAg                                                                                                                                >  1:3667569/1‑150 (MQ=255)
                      gAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCGCACCCCGACCCTCTCCCCCAGGGAGAGGGAGAAAACACGCAAg                                                                                                                                >  2:5428738/1‑150 (MQ=255)
                      gAAGCGATTCGCAGCCATGCGATTGAGCAGGGAGGTCAGAAACTGCGCGGCTATATGACAGCTCTGACCCGGATTGCTGGGGCGGGGGAAGTGCGCGCTGAGGCCCTCCCCCCGCCCCTCGCCCCCCGGGCGCGGGCGACCACACGCAAg                                                                                                                                >  2:5003728/1‑150 (MQ=255)
                        aGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGACACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGc                                                                                                                              <  2:6538781/150‑1 (MQ=255)
                         gCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGAGGCCCTCACCCCCACCCACTCCCCCAGGGCGAGGGGGAAAACCCTCAAGGcc                                                                                                                             >  2:6430805/1‑150 (MQ=255)
                              tCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTtctc                                                                                                                        <  2:298292/150‑1 (MQ=255)
                              tCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGAGGCCCTCACCCCCACCCACAACCACCCGGCGGGGGAGTAAAAACACAAGGCCGtcta                                                                                                                        >  1:6811817/1‑149 (MQ=255)
                                gCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGGCCTCACCCCGACCCTCTACCACATGGAGTGGGAGAAAACACTCAAGGCCttgtgct                                                                                                                      >  1:1047826/1‑145 (MQ=255)
                                  aGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTgg                                                                                                                    <  1:5011328/150‑1 (MQ=255)
                                  aGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCACACCCCCACCCACTACCACAGGGAGAGGGAGAAAACACGCAAGGCCATCTCCTgg                                                                                                                    >  1:4595840/1‑150 (MQ=39)
                                   gCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGga                                                                                                                   <  1:4762751/150‑1 (MQ=255)
                                      aTGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTGAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACCGGGAGAGGGAGAAAACACTCATGGCCTTCTCCTGGAGaa                                                                                                                >  2:3054154/1‑150 (MQ=255)
                                       tGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGAGGCCCTCACCCCGCCCCACACCCCCAGGGAGAGGGGGACAACACACAATGCCATCTACGGgagacg                                                                                                               >  1:3967322/1‑148 (MQ=255)
                                       tGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGAGGCCCGCACCCCCCCCCGCACCCACACCGGGGGGGAGAAAAAACACAAGGCCATCTCCCGGAGCAg                                                                                                               >  2:2689980/1‑150 (MQ=255)
                                        gCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGGCCGCACCCCCACCCACACCCACCCGGCGCGGGGGAAAACACACAAGGCCGTAACCGGGAGAAgg                                                                                                              >  2:360606/1‑150 (MQ=255)
                                        gCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACGCAAGGCCTTCTCCAGGAGAAgg                                                                                                              <  1:1120989/150‑1 (MQ=255)
                                          gATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCAAACCCACCGGGGGAGGGCGAAAAAACACAAAGCCGTCTCCGGGAGAAGGcc                                                                                                            >  1:1453255/1‑150 (MQ=255)
                                            ttGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCtt                                                                                                          <  2:1571253/150‑1 (MQ=255)
                                                 cAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCCCAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTAt                                                                                                     <  2:5503918/150‑1 (MQ=255)
                                                   ggtaggtaAGTATCTGCGCGGCTATATGACAGATATGAATCGTATTGCTGTTGCGTGTTTAGTGCGCGCTGTTGCCCTCACACCGACCCTCTCCCACATGGAGATGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTTCTATTa                                                                                                   <  2:4200399/150‑1 (MQ=255)
                                                    gtaggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCACTCCCACAGGGAGAGGGAGAAAAAACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTat                                                                                                  >  2:5289932/1‑149 (MQ=255)
                                                         taAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCACACCCCCACCCACACCCACAGGGGGTGGGGGGAAAAACAAAAGGCCGTATACAGGGGGAGGCCTTGCTATTATTTTCg                                                                                             >  2:1305935/1‑150 (MQ=255)
                                                            gAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGGTAAGTGCGCGCTGAGGCCCTCCCCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTGCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTa                                                                                          <  1:4373605/150‑1 (MQ=255)
                                                               aCTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAg                                                                                       <  1:2547152/150‑1 (MQ=255)
                                                               aCTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAg                                                                                       <  1:3557777/150‑1 (MQ=255)
                                                                cTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGc                                                                                      <  1:4775991/150‑1 (MQ=255)
                                                                 tGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCGCACCCCCCCCCACACCCCCCCGGAGCGGGAGAAAACACACAAGGCCCTCACCGCGCGAAGGCCGGGCGATGACTTGCgggccaccc                                                                                     >  1:6637700/1‑143 (MQ=255)
                                                                  gcgcGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCGCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCAc                                                                                    <  2:1651628/150‑1 (MQ=255)
                                                                   cgcgGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCAcc                                                                                   <  2:6112/150‑1 (MQ=255)
                                                                    gcgGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTGAGGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCCCCGGGAGAGGGAGAAAACAATCAAGGCCTTCTCCGGGCGAAGGCCTTGCGATTAGTTGCGGTCCAGCACCt                                                                                  >  2:4696369/1‑150 (MQ=37)
                                                                        cTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAat                                                                              <  1:552452/150‑1 (MQ=255)
                                                                             tGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCTCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGTAGGCCTTGCTATTAGTTGCGGTACATCACCTTTAtgttgt                                                                         <  1:1940650/150‑1 (MQ=38)
                                                                             tGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGAtgt                                                                         <  1:3088133/150‑1 (MQ=255)
                                                                                  gATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACCGGGAGAGGGAGAAAACACTCTTGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTCAATGATGTGatat                                                                    <  2:1452787/150‑1 (MQ=255)
                                                                                  gATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTCAATGATGTGatat                                                                    <  1:5356796/150‑1 (MQ=255)
                                                                                        aaaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCCACCCTCTCCCACAGGGGGCGGGGGAAAAAACACAAGGCCTTCTCCTGGAGAAGGCCATGCTATTAGTTGCGGGACAGCACCTTAATGATGTGATATACGAAc                                                              >  2:3318669/1‑150 (MQ=38)
                                                                                                     ttGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACGCAAGGCCTTCTCCGGGAGAAGGCCGTGCTATTAGTTGCGGTACAGCACCTTAATGTTGTGATATCCGCACGGGGTGGGCAGCg                                                 >  2:3045070/1‑150 (MQ=37)
                                                                                                       gCGGGTTAAGTGCGCGCTGATGCCCACACCCCCACCCACACCCACCCGGGGCGGGAGAAAAAACACAAGCCCGTATCCGGGAGAAGGCCGTGCTATTAGTTGGGGGACAGCACCATAATTAAGTGATATACGAACGGGGGGTGCAGCggg                                               >  2:3722938/1‑150 (MQ=255)
                                                                                                         gggTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGcc                                             <  1:335812/150‑1 (MQ=255)
                                                                                                         gggTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGcc                                             <  1:1601553/150‑1 (MQ=255)
                                                                                                         gggTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGcc                                             <  2:2582193/150‑1 (MQ=255)
                                                                                                          ggTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCg                                            <  1:2367313/150‑1 (MQ=255)
                                                                                                          ggTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCg                                            <  2:6142158/150‑1 (MQ=255)
                                                                                                            ttAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGt                                          <  1:2701537/150‑1 (MQ=255)
                                                                                                             tAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACCGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTc                                         <  2:4859202/150‑1 (MQ=255)
                                                                                                             tAAGTGCGCGCTGATGCCCACACCCCCACCCACTCCCACAGGGGGAGGGAGAAAACACACAAGGCCATATACAGGAGAAAGCCATGCTATTAGTTGCCGTACAGCACCATAATTAAGTGATAAACCAACAGGGTGTGCAGCGGGCCCGTc                                         >  2:4459696/1‑150 (MQ=255)
                                                                                                               aGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTTTGCAGCGGGCCGGTCgg                                       <  2:5712789/150‑1 (MQ=255)
                                                                                                                 tGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATGAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCt                                     <  2:2225682/150‑1 (MQ=255)
                                                                                                                  gcgcgcTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCTCCTTAATGATGTGATATCCGACCTGGGTGTGCAGCGGGCCGGTCGGCTc                                    <  2:7216475/150‑1 (MQ=255)
                                                                                                                  gcgcgcTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTc                                    <  1:388203/150‑1 (MQ=255)
                                                                                                                    gcgcTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCa                                  <  1:7078267/150‑1 (MQ=255)
                                                                                                                    gcgcTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCGTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCa                                  <  2:6801474/150‑1 (MQ=255)
                                                                                                                    gcgcTGATGCCCACACCCCCACCCACTCCCCAAGGGAGAGGGAGAAAACACTCAAGGCCATATCCGGGGGAAGGGCGTGCTATTATTTGCGGTACAGCACCGTAATTATGTTATATACCAACAGGGTGTGCAGCGGGCCGGGCGGCTCCa                                  >  2:7100782/1‑150 (MQ=255)
                                                                                                                     cgcTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAg                                 <  2:5973307/150‑1 (MQ=255)
                                                                                                                         gATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTAcc                             <  1:4130850/150‑1 (MQ=255)
                                                                                                                          aTGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCg                            <  1:4257939/150‑1 (MQ=255)
                                                                                                                          aTGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCg                            <  1:6083446/150‑1 (MQ=255)
                                                                                                                            gCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGa                          <  1:1766347/150‑1 (MQ=255)
                                                                                                                            gCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGa                          <  2:1962175/150‑1 (MQ=255)
                                                                                                                               cTCACCCCGACCCTCTCCCCCAGGGAGAGGGAGAAAACACTCAAGGCCTTCGCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACaa                       <  2:3954708/150‑1 (MQ=255)
                                                                                                                               cTCACCCCGACCCTCTCCCACAGGGAGAGGAAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTATTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACaa                       >  1:3036387/1‑150 (MQ=255)
                                                                                                                                tCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCCGTACCGGACAag                      >  2:4349885/1‑150 (MQ=255)
                                                                                                                                   ccccGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAgaga                   <  2:3397503/150‑1 (MQ=255)
                                                                                                                                   ccccGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTACTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAgaga                   <  1:2210189/150‑1 (MQ=255)
                                                                                                                                    cccGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGaa                  >  2:1076948/1‑150 (MQ=255)
                                                                                                                                    cccGACCCTCTCCCACAGGGAGAGGGAGAAAACACACAAGGCCATCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCAACTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGaa                  >  1:4059795/1‑150 (MQ=255)
                                                                                                                                          ccTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCAc            >  1:756780/1‑150 (MQ=255)
                                                                                                                                          ccTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTGGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCAc            >  1:431193/1‑150 (MQ=255)
                                                                                                                                            tctcCCACAGGGAGAGGGAGAATACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACtt          >  2:3209072/1‑150 (MQ=255)
                                                                                                                                            tctcCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCTGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACtt          >  1:1694302/1‑150 (MQ=255)
                                                                                                                                              tcCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGACCTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTa        >  2:1731158/1‑150 (MQ=255)
                                                                                                                                                ccACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCAGACAAGAGAAAACCACTTTATc      >  2:669678/1‑150 (MQ=255)
                                                                                                                                                 cacaGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCg     >  1:2106601/1‑150 (MQ=255)
                                                                                                                                                 cacaGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCg     >  2:3479134/1‑150 (MQ=255)
                                                                                                                                                  acaGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGa    >  1:5749564/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAAC  >  NC_000913/3959334‑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: