Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,959,501 T→A 14.2% intergenic (+56/+31) 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,5010TA14.2% 563.0 / 16.9 209intergenic (+56/+31)ilvC/ppiCketol‑acid reductoisomerase, NAD(P)‑binding/peptidyl‑prolyl cis‑trans isomerase C (rotamase C)
Reads supporting (aligned to +/- strand):  ref base T (60/110);  new base A (32/2);  total (96/113)
Fisher's exact test for biased strand distribution p-value = 5.90e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 216 total aligned reads displayed.

GAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCT  >  NC_000913/3959352‑3959650
                                                                                                                                                     |                                                                                                                                                     
    gAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAg                                                                                                                                                   <  1:1404652/150‑1 (MQ=255)
    gAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAg                                                                                                                                                   <  1:478020/150‑1 (MQ=255)
    gAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCACCCACAGGGGGAGGGAGAAAAAACACAAg                                                                                                                                                   >  2:613765/1‑150 (MQ=255)
    gAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCACACCCCCACCCGCTCCCACACGGAGAGGGGGAAAACACTCAAg                                                                                                                                                   >  1:1298311/1‑150 (MQ=255)
    gAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGAGGCCCACACCCCGACCCACTCCCCCAGGGAGCGGGAGAAAACACGCAAg                                                                                                                                                   >  1:897187/1‑150 (MQ=255)
     aaGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACACAAgg                                                                                                                                                  >  2:2047089/1‑150 (MQ=255)
      aGCGATTCGCAGCCATGTGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGc                                                                                                                                                 <  1:2070216/150‑1 (MQ=255)
      aGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCACACAGGGAGAGGGAGAAAACACTCAAGGc                                                                                                                                                 <  2:545501/150‑1 (MQ=255)
        cGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCCCAGGGAGAGGGAGAAAACACTCAAGGCCt                                                                                                                                               >  2:2412290/1‑150 (MQ=255)
           ttCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGAGGCCCGCACCCCCCCCCACACCCCCCCGGAGGGGGGGAAAAAACACAAAGCCGtct                                                                                                                                            >  2:3465499/1‑150 (MQ=255)
            tCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTtctc                                                                                                                                           <  2:3992358/150‑1 (MQ=255)
             cGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCACGCCCCCAGGGCGCGGGAGCAAACACACAAGGCCTTCTcc                                                                                                                                          >  2:906396/1‑150 (MQ=255)
               cAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTg                                                                                                                                        <  1:331642/150‑1 (MQ=255)
                aGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTgg                                                                                                                                       <  1:3286601/150‑1 (MQ=255)
                aGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTgg                                                                                                                                       <  1:759689/150‑1 (MQ=255)
                   cATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGGGCTCACCCCGAACCACTACCAAAGGGTGAGGGTGTAAAAACTCAATGGCTTGTACTGgaga                                                                                                                                    >  1:2224932/1‑150 (MQ=255)
                   cATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGgaga                                                                                                                                    <  1:956915/150‑1 (MQ=255)
                    aTGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGaa                                                                                                                                   <  1:2669798/150‑1 (MQ=255)
                      gCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAgg                                                                                                                                 <  2:1567388/150‑1 (MQ=255)
                       cGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCAGGAGAAGGc                                                                                                                                <  2:3428236/150‑1 (MQ=255)
                        gATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGcc                                                                                                                               <  2:2769139/150‑1 (MQ=255)
                        gATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGcc                                                                                                                               <  2:1459522/150‑1 (MQ=255)
                        gATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCACACCCCGACCCACTCCCACAGGGGGAGGGGGAAAACACTCAAGGCCGTCTCCTGGAGAAGGcc                                                                                                                               >  1:227033/1‑150 (MQ=255)
                         aTTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCCACCCAATACCACACGGAGAGGGAGAAAAAAAAAAAAGCCAACACCTGGAGAAGGCCt                                                                                                                              >  1:3216900/1‑150 (MQ=255)
                            gAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGc                                                                                                                           <  1:1410956/150‑1 (MQ=255)
                             aGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGAGATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACCGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCt                                                                                                                          <  1:2548566/150‑1 (MQ=255)
                                 ggtaggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTa                                                                                                                      <  1:1940242/150‑1 (MQ=255)
                                  gtaggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAg                                                                                                                     <  1:3153305/150‑1 (MQ=255)
                                  gtaggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAg                                                                                                                     <  2:3674783/150‑1 (MQ=255)
                                   taggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGt                                                                                                                    <  2:100396/150‑1 (MQ=255)
                                       taAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATTCCCTCACCCCGACCCTCTCCCACAGGGGGGGGGGGAAAAAACACAAGGCCGTTTCCTGGGGAAGGCCATGCTAGTATTTTCg                                                                                                                >  2:883813/1‑150 (MQ=255)
                                       taAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCg                                                                                                                <  1:233938/150‑1 (MQ=255)
                                        aaGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCgg                                                                                                               <  2:2936702/150‑1 (MQ=255)
                                        aaGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACACAAGGCCATCTCCTGGAGAAGGCCTTGCTATTAGTTGCgg                                                                                                               >  2:3022215/1‑150 (MQ=255)
                                        aaGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCAACCCCCCACCCCCACCCACCGGGAGGGGGTGAAAACAAACAAGGCCTTAACCGGGCGTCTGCCGTGCGATTAGTTGCgg                                                                                                               >  2:1203346/1‑150 (MQ=255)
                                         aGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGCCCCTCAACCACAGGGCGGGGGTGTAAACACACAAGGCCGTCTCCAGGGGAAGGGCTTGGGATTAGTTggggg                                                                                                              >  2:1926773/1‑149 (MQ=255)
                                          taaaCTGCGCGGCTATATGACAGATATGAAACTTATTGCTGTTGCGGGTTAAGTGCGAGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGATGGAGAAAACACTCAAGGCCTTCTCCTTTATAAGGCCTTGCTATTAGTTGCGGTa                                                                                                             <  1:2647893/149‑1 (MQ=255)
                                                 cgcgGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCAcc                                                                                                      <  1:4108434/150‑1 (MQ=255)
                                                  gcgGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCACACCCCCACCCACTCCCCCACGGAGAGGGCGAAAACACAAAATGCCTACTACAGGATAAGGCCGTTCTATTAGTTTCGGTACAGCAcca                                                                                                     >  1:1679325/1‑149 (MQ=255)
                                                    ggCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCCCCGGGAGAGGGCGAAAACACTCAAGGCCTTCTCCAGGGGGCGGCCGTGCGATTAGGTGCGGGACAGCAccgtt                                                                                                   >  2:4030243/1‑147 (MQ=37)
                                                    ggCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGTCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTa                                                                                                   >  2:2728908/1‑150 (MQ=255)
                                                     gCTATATGACAGATATGACACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTaa                                                                                                  <  1:2734770/150‑1 (MQ=255)
                                                     gCTAGATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTaa                                                                                                  <  2:3622768/150‑1 (MQ=255)
                                                          atGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAatgatg                                                                                             <  1:3827203/150‑1 (MQ=255)
                                                          atGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAatgatg                                                                                             <  1:18420/150‑1 (MQ=255)
                                                          atGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAatgatg                                                                                             <  1:1058251/150‑1 (MQ=255)
                                                            gccAGCTATGACCCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAACCACTCACGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGAtgtg                                                                                           <  1:2124658/148‑1 (MQ=38)
                                                            gACAGATAAGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACACAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGAtgtg                                                                                           <  1:3942305/150‑1 (MQ=255)
                                                             aCAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCCACCCACGCCCACCCGGAGAGGGAGAAAAAACACAAGGCCTTCAACAGGAGAAGGCCATGCTATTATTTTCGGGAAAGCACCATAATTATGTGa                                                                                          >  1:1648256/1‑150 (MQ=255)
                                                                       aaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACt                                                                                <  2:678118/150‑1 (MQ=255)
                                                                          gTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCACAACCCCACCCACACCCACCGGGAGAGGGAGAAAAAACACAATGCCGTCTCCGGGAGAAGGCCGTGCGATTAGTTGCGGTACAGCACCGTAAATATGTGATATACGAACAggg                                                                             >  2:3293341/1‑150 (MQ=255)
                                                                           tATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGTGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGgt                                                                            <  1:3268922/150‑1 (MQ=255)
                                                                            aTTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGgtg                                                                           <  1:176992/150‑1 (MQ=255)
                                                                                cTGTTGCGGGTTAAGTGCGCGCTGATGCCCGCACCCCCACCCACTCCCACAGGGCGGGGGAGAAAACACACAAGGCCATCTCCGGGCGAAGGCCGTGGTATTAGTTTCGGGACAGCACCTTAAATATGTGTTATACGCACTTGGTGTGCa                                                                       >  2:2867617/1‑150 (MQ=255)
                                                                                cTGTTGCGGGTTAAGTGCGCGCTGAGGCCCACACCCCCCCCCCCACCCCCACGGCCAGGGAGAAAACACACAAGGCCAACGCCACGCGAAGGCCCGGCTAGTAGGTGCGGGACCCCACCCTAAAGAGGGGTAATCCGCACGGGGGGTGCa                                                                       >  2:3500854/1‑150 (MQ=255)
                                                                                 tGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAg                                                                      <  1:1237757/150‑1 (MQ=255)
                                                                                   ttGCGGGTTAAGTGCGCGCTGATGCCCGAACCCCGACCCAAACCCACAGGGCGAGGGAGAAAAAACACAAGGCCATCTACTGGAGAAGGCCTTGCTATTAGTTGCGGTAAAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCg                                                                    >  2:440023/1‑150 (MQ=255)
                                                                                    tGCGGGTTACGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCgg                                                                   <  2:3664564/150‑1 (MQ=255)
                                                                                        ggTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCg                                                               <  1:1847156/150‑1 (MQ=255)
                                                                                          ttAAGTGCGCGCTGATGCCCACACCCCCACCCTCACACACCCGGAGCTGGAGAAAACACACAAGGCCATATCCCGGGGGAGGCCGTGCGATTACTTTGGGGACCGCCCCATAATGATGGGATATACCAAAGGGGGGTTCAGCGGGCCggg                                                             >  2:869358/1‑149 (MQ=255)
                                                                                              gTGCGCGCTGAGGCCCTCACCCCGACCCTCGCCCACCGGGAGAGGGAGACAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGc                                                         <  1:3167654/150‑1 (MQ=255)
                                                                                               tGCGCGCTGATGCCCACACCCCGACCCTCACCCACAGGGAGAGGGAGAAAAAACACAAGGCCTTCTCCGGGAGAAGGCCGTGCTATTAGTTGCGGTACAGCACCATAATGATGTGATATCCGAACAGGGTGTGCAGCGGGCCGGTCGGca                                                        >  2:2513509/1‑149 (MQ=38)
                                                                                                gcgcgcTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTc                                                       <  2:136937/150‑1 (MQ=255)
                                                                                                gcgcgcTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTc                                                       <  2:3307081/150‑1 (MQ=255)
                                                                                                gcgcgcGGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAACCACTCCCGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTc                                                       <  2:1786902/150‑1 (MQ=255)
                                                                                                         tGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTGGCTATAAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCgg                                              >  2:2085225/1‑150 (MQ=255)
                                                                                                          gCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGa                                             <  2:1834669/150‑1 (MQ=255)
                                                                                                           cccTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTCCCGGAc                                            >  2:4056713/1‑150 (MQ=255)
                                                                                                           cccTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGAc                                            >  2:3806450/1‑150 (MQ=255)
                                                                                                            ccTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACa                                           <  2:2472755/150‑1 (MQ=255)
                                                                                                              tCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAag                                         >  2:2377400/1‑150 (MQ=255)
                                                                                                                aCCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAagag                                       >  2:3133009/1‑150 (MQ=255)
                                                                                                                    cGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCACGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGaaaa                                   <  2:1841237/150‑1 (MQ=255)
                                                                                                                       cccTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCa                                >  1:1906684/1‑150 (MQ=255)
                                                                                                                       cccTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCa                                >  1:3065275/1‑150 (MQ=255)
                                                                                                                        ccTCTCCCTGTGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCAc                               >  2:706544/1‑150 (MQ=255)
                                                                                                                          tctcCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACtt                             >  2:2509085/1‑150 (MQ=255)
                                                                                                                          tctcCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGGACCGGACAAGAGAAAACCACtt                             >  1:3993455/1‑150 (MQ=255)
                                                                                                                             cccACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGa                                                                                          <  1:2512117/86‑1 (MQ=255)
                                                                                                                             cccACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGATATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTAt                          <  2:1144800/150‑1 (MQ=255)
                                                                                                                              ccACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATc                         >  2:3444711/1‑150 (MQ=255)
                                                                                                                                    ggagtgggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTATTGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGAAAAGAGAAAACCACTTTATCGAACGc                   <  1:3878457/150‑1 (MQ=255)
                                                                                                                                      tgagggataAAACACTCCCGGCCTTCGCCTGGAGATGGCGTTGCTTTTCGTGGCGGTACTGCACCTTTTTGATGTGTTATCCGACCGGGGTGTGCCGCGGGCGGGTGGGCTCCTGTACCGGACAAGAGAACACCACTTTATCGAACGCCg                 <  1:3118784/149‑1 (MQ=255)
                                                                                                                                      agagggagaAAACACTCAAGGCCTTCTCCTGGAGAAGTCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCg                 >  2:2562675/1‑150 (MQ=255)
                                                                                                                                         gggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGTTGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGaa              <  2:1281403/150‑1 (MQ=255)
                                                                                                                                         gggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGTTGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGaa              <  1:3652847/150‑1 (MQ=255)
                                                                                                                                         gggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGaa              <  2:4022240/150‑1 (MQ=255)
                                                                                                                                          ggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAAc             >  1:624452/1‑150 (MQ=255)
                                                                                                                                          ggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAAc             >  2:322952/1‑150 (MQ=255)
                                                                                                                                           gagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCCGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAccc            >  1:1465717/1‑150 (MQ=255)
                                                                                                                                              aaaaCACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCCGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATc         >  2:3345743/1‑150 (MQ=255)
                                                                                                                                              aaaaCACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATc         >  2:235032/1‑150 (MQ=255)
                                                                                                                                              aaaaCACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATc         >  2:135090/1‑150 (MQ=255)
                                                                                                                                               aaaCACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTGTATCGAACGCCGGAACCATCt        >  1:1740475/1‑150 (MQ=255)
                                                                                                                                                aaCACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTg       >  1:3986358/1‑150 (MQ=255)
                                                                                                                                                   acTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGAcc    >  1:425638/1‑150 (MQ=255)
                                                                                                                                                     tCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCt  >  1:1360087/1‑150 (MQ=255)
                                                                                                                                                     tCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCt  >  1:2837950/1‑150 (MQ=255)
                                                                                                                                                     tCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCt  >  1:3539149/1‑150 (MQ=255)
                                                                                                                                                     tCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCt  >  1:965014/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCT  >  NC_000913/3959352‑3959650

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: