Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,728,109 A→C 5.8% D421E (GAT→GAG xylB ← xylulokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,728,1090AC 5.8% 1040.7 / 3.2 375D421E (GAT→GAGxylBxylulokinase
Reads supporting (aligned to +/- strand):  ref base A (126/212);  new base C (27/3);  total (160/215)
Fisher's exact test for biased strand distribution p-value = 1.28e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 381 total aligned reads displayed.

TTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGC  >  NC_000913/3727961‑3728258
                                                                                                                                                    |                                                                                                                                                     
tttCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAt                                                                                                                                                      >  1:3408758/1‑150 (MQ=255)
 ttCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATc                                                                                                                                                     <  1:820809/150‑1 (MQ=255)
 ttCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATc                                                                                                                                                     <  1:3504523/150‑1 (MQ=255)
    tcgtcgTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCTccgc                                                                                                                                                  >  1:3887801/1‑148 (MQ=255)
    tcgtcgTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATcccc                                                                                                                                                  >  2:806249/1‑150 (MQ=255)
    tcgtcgTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATcccc                                                                                                                                                  >  2:2247541/1‑150 (MQ=255)
    tcgtcgTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATcccc                                                                                                                                                  >  1:2817146/1‑150 (MQ=255)
        cgtGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATcccccccc                                                                                                                                              <  1:491000/150‑1 (MQ=255)
        cgtGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATcccccccc                                                                                                                                              <  1:3989400/150‑1 (MQ=255)
         gtGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCg                                                                                                                                             <  2:1540722/150‑1 (MQ=255)
           ggCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTa                                                                                                                                           <  1:358027/150‑1 (MQ=255)
              tGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACGTCGCCCCCCgggcgg                                                                                                                                        >  2:635258/1‑145 (MQ=255)
                  aaGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACGTCCCCCCCCGGACGGTacg                                                                                                                                    >  2:532737/1‑148 (MQ=255)
                   aGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTGCgggactc                                                                                                                                   >  2:2311139/1‑145 (MQ=255)
                    gCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCGCCCCCCCCGTCCGGGACTCg                                                                                                                                  >  2:2384001/1‑150 (MQ=255)
                     cGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGa                                                                                                                                 <  1:762452/150‑1 (MQ=255)
                      ggCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAg                                                                                                                                <  2:162089/150‑1 (MQ=255)
                        cATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAgct                                                                                                                              <  1:1898663/150‑1 (MQ=255)
                          tAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAgctgc                                                                                                                            <  1:2063669/150‑1 (MQ=255)
                            aCGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGctgctg                                                                                                                          <  2:1974408/150‑1 (MQ=255)
                                   gcATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTg                                                                                                                   <  1:1411526/150‑1 (MQ=255)
                                    cATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGa                                                                                                                  <  2:3535753/150‑1 (MQ=255)
                                     aTCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGat                                                                                                                 <  1:1279969/150‑1 (MQ=255)
                                       cTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCTCCCCCGCCGGGCGGGGGTCGGGCTGCTGGCCGCTggggg                                                                                                               >  1:1124926/1‑146 (MQ=255)
                                        tGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGGACGGGCATCGAGCTGCTGGCCGcggggggc                                                                                                              >  2:3653968/1‑143 (MQ=255)
                                           tAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGc                                                                                                           <  1:1355436/150‑1 (MQ=255)
                                                  gACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCTCGCCCCCCGTGCGGGGAGGGGGCTGCGGGCCGCTGAGGTCCGcgcgcctc                                                                                                    >  2:2535863/1‑143 (MQ=255)
                                                   aCTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCt                                                                                                   <  1:3064187/150‑1 (MQ=255)
                                                        tctcACGGTAGGGGCGGCAACAAGTCAATGCGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCGGACCGCTGATATCCGCCAGCATCTGACGc                                                                                              <  2:2720533/150‑1 (MQ=255)
                                                          tAACGGTAGTTGCGGCAACAATTCATTGAGCGTTTTCTCTGTATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCa                                                                                            <  2:485774/150‑1 (MQ=255)
                                                          tAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCa                                                                                            <  1:233257/150‑1 (MQ=255)
                                                          tAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCa                                                                                            <  1:3643182/150‑1 (MQ=255)
                                                           aaCGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCTCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAg                                                                                           <  2:3408758/150‑1 (MQ=255)
                                                              ggTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTCCGGGAGGCGAGCTGCGGGCCGCGGATGTCCGCCCGCCGCTGACGCCCGtcc                                                                                        >  2:3063761/1‑148 (MQ=255)
                                                              ggTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGGGATCGTGCTGCTGGCCGCTGGTGTCCGCCGGCGTCTGGCGCCGGTAc                                                                                        >  2:1237742/1‑150 (MQ=255)
                                                               gTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGt                                                                                          <  2:1952664/147‑1 (MQ=255)
                                                                tAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTc                                                                                      <  1:1790882/150‑1 (MQ=255)
                                                                 aGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCa                                                                                     <  2:3477917/150‑1 (MQ=255)
                                                                   ttGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGGTGTCCGCCCGCATCTGCCGCCCGTACGCACt                                                                                   >  2:3118111/1‑150 (MQ=255)
                                                                    tGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTa                                                                                  <  2:2678353/150‑1 (MQ=255)
                                                                      cGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGGACGGTCATCGAGCTGCTGACCGCTGGGGTCCGCCAGCATCTGGCGCCAGGACTCCCTGcg                                                                                >  2:1650456/1‑150 (MQ=255)
                                                                          aacaaTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTGCGGGCATCGAGCTGCTGACCGCTGATAGCCGCCCGCATCTGACGCCAGTGCTCACTACGcgcc                                                                            >  2:423156/1‑150 (MQ=255)
                                                                          aacaaTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAGCCCCCCCCGTGCGGTGCTCGGGCTGCTGGGCGCGGATCTCCGCCCGCCGCTGGCGCCGGTACTCGCTACGcgcc                                                                            >  1:871006/1‑150 (MQ=255)
                                                                             aaTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGGGCTGCTGACCGCTGATATCCGCCCGCGGCTGACGCCAGGACTCACTGCGcgccccg                                                                         >  2:1801860/1‑150 (MQ=255)
                                                                               ttCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCCGTGCTGGCCCCACCGGCCGCGCCGTACGGTGGTCGTGCTGCGGGCGGGGGGGGGCGGGCAGCGTCGGGCGCCGGGACTCACGACGcgccccgcc                                                                       >  1:533816/1‑150 (MQ=255)
                                                                               ttCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATCTCCGCCTGCATCTGGCGCCCGTCCTCACTACGcgccccgcc                                                                       >  2:3701637/1‑150 (MQ=255)
                                                                               ttCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAGCCCCCCCCGGAGGGGAGTGGGGCTGGTGGCCGCGGATGTGCGGCAGCCGCTGGCGCCCGGGCGGGCGACGcgccccgcc                                                                       >  1:3768913/1‑150 (MQ=255)
                                                                                tCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCGCCCCCCCCGTGGGGGAGGCGGGCGGCGGGGCGCGGGTGTCCGCCCGCGGCGGACGCGCGGCCGCACTGCGcgccccgccc                                                                      >  1:2178685/1‑150 (MQ=255)
                                                                                tCAATGAGCGATTTCTCTGGATGCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGcgccccgccc                                                                      <  1:1960600/150‑1 (MQ=255)
                                                                                    tGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGGACGCGCGCTGCGGACCGCGGAGCTCCGCCCGCAGCGGTCGCCAGGACTCACGGCGCGCCCCGccccccg                                                                  >  2:760517/1‑148 (MQ=255)
                                                                                      aGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTCGGGTAATCGGGCTGCTGGCCGCTGTTATCCGCCAGCGGCTGGCGCCAGTACTCACTGCGCgccccgccccgcctc                                                                >  2:3465300/1‑145 (MQ=255)
                                                                                       gCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTGCGGGACTCGAGCGGCTGGCCGCGGGTCTCCGCCCGCAGCTGGCGCCAGTTCGCTCTACGCGCCCCGCCCCCCATCa                                                               >  2:2017541/1‑150 (MQ=255)
                                                                                       gCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCa                                                               <  2:1303944/150‑1 (MQ=255)
                                                                                         gATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTGCGGTAATCGAGCTGCGGACCGCTGATATCCGCCAGCCTCTGGCGCCAGTACTCACTACGCGCCCCGCCCCCACTCCAc                                                             >  2:1096940/1‑150 (MQ=255)
                                                                                         gATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAAc                                                             <  1:985384/150‑1 (MQ=255)
                                                                                         gATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAAc                                                             <  1:294951/150‑1 (MQ=255)
                                                                                             tctcTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCATCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTaa                                                         <  2:3214376/150‑1 (MQ=255)
                                                                                              ctctGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAAGCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAac                                                        <  1:3312539/150‑1 (MQ=255)
                                                                                                    aTTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCGGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTg                                                  <  2:184715/150‑1 (MQ=255)
                                                                                                      tCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGTGCTGCTGACCGCTGGTGTCCGCCAGCATCGGGCGCCCGTTCTCACTACGCGCCCCGCCCCCGCTCGACGGAACACGTTGCg                                                >  2:1395446/1‑150 (MQ=255)
                                                                                                      tCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCg                                                <  2:1026169/150‑1 (MQ=255)
                                                                                                      tCGCCGCGATCTGCGCCAGCCTTGATGCGCCCAGTGCTGGCCCCACATGCCCCCCCGGGCGGTAATCGGGCTGCTGGCCGCGGATTTCCGCCGGCATCTGACGCGGGTACGGACTTCGCGCCCCGCCCCCAGGCGGCGTAACACGTGGCg                                                >  2:682014/1‑150 (MQ=255)
                                                                                                       cgccgcGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCgg                                               <  1:2721534/150‑1 (MQ=255)
                                                                                                         ccgcGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTCGGGTAATCGGGCTGCTGACCGCGGATATCCGCCAGCGTCTGGCGCCAGTACTCTCTACGCGCCCCGCCCCCCCGCGCCGGAACACTTTGCgggg                                             >  2:3772236/1‑148 (MQ=255)
                                                                                                         ccgcGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGtt                                             <  1:607812/150‑1 (MQ=255)
                                                                                                         ccgcGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGtt                                             <  2:259723/150‑1 (MQ=255)
                                                                                                          cgcgATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTTCGGTAATCGAGCTGCTGGCCGCTGGTGTCCGCCAGCATCTGGCGCCCGTTCTCACTACGCGCCCCGCCCCCAGGCCACGTAACACGGTGCGGttg                                            >  2:278328/1‑149 (MQ=255)
                                                                                                          cgcgATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGttt                                            <  1:384377/150‑1 (MQ=255)
                                                                                                             gATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCCTCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTCACACTGTGCGGTTTAAt                                         >  2:2468500/1‑150 (MQ=255)
                                                                                                              aTCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATa                                        <  2:1568245/150‑1 (MQ=255)
                                                                                                                 tgtgcCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCTCCCCCCCCGGTCGGGAATCGAGCTGCTGACCGCTGATATCCGCCGGCATCTGACGCCCGTACTCACTGCGCGCCCCGCCCCCAATCAACGTAGCACGTTGCGGGTTAGTGCCg                                     >  2:40534/4‑150 (MQ=255)
                                                                                                                   cgcCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCTCCCCCCCCGTACGGTGATCGAGCTGCTGTCCGCTGAGATCCGCCCGCAGCTGACGCCCGTACTCGCTACGCGCCCCGCCCCCAATCGACGTAACACGTTGCGGTTTAATGCCGcg                                   >  2:3780486/1‑149 (MQ=255)
                                                                                                                     ccAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTATTAAACCGCaaa                                 <  1:1091568/150‑3 (MQ=255)
                                                                                                                      cAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGGAATCGAGCTGCTGACCGCTGATATCCGCCAGCCTCTGACGCCAGTACTCACTACGCGCCCCGCCCCCCATCAACGGCACACGTTGCGGGGTAGTACCGCCGGc                                >  2:3183930/1‑150 (MQ=255)
                                                                                                                          cTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGc                            <  1:1657173/150‑1 (MQ=255)
                                                                                                                            tgctgcGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCAc                          <  1:2432245/150‑1 (MQ=255)
                                                                                                                            tgctgcGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCAc                          <  2:1187857/150‑1 (MQ=255)
                                                                                                                              ctgcGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCGATCAACGTAACACTTTGCGGTTTAATAGCGCAGGCATGCAGGa                        >  2:2858661/1‑150 (MQ=255)
                                                                                                                               tgcGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGAc                       <  2:713688/150‑1 (MQ=255)
                                                                                                                                 cgcCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACAt                     >  2:1757711/1‑150 (MQ=255)
                                                                                                                                   cccAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGCGCTGCTGACCGCTGATGTCCGCCAGCGTCTGACGCCAGTACTCACTACGCGCCCCGCCCCCCATCAACGGAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATcc                   >  1:3456016/1‑150 (MQ=255)
                                                                                                                                     cAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCAt                 >  2:2054433/1‑150 (MQ=255)
                                                                                                                                      aGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATg                <  2:3981738/150‑1 (MQ=255)
                                                                                                                                      aGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATg                <  2:3861135/150‑1 (MQ=255)
                                                                                                                                       gTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGc               <  1:1031522/150‑1 (MQ=255)
                                                                                                                                        tGCTGGCCCCACATCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCAGGCACGACATCCATGCCa             >  2:2107095/1‑150 (MQ=255)
                                                                                                                                         gCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCa             <  1:3226371/150‑1 (MQ=255)
                                                                                                                                         gCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCa             <  1:3029598/150‑1 (MQ=255)
                                                                                                                                          cTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACGACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCAt            >  2:324768/1‑150 (MQ=255)
                                                                                                                                           tGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATc           <  1:1698720/150‑1 (MQ=255)
                                                                                                                                             gCCCCACTTACACCCCCCTACGGTAATCGAGCTGCTGACCCCAGATATCCCCCAGAAAATCACCCCAGTACTAACTAAGCGCCCCCCACCCACACAACCTAACACTATGCGGATTAATACAGAAGGCATGCCCGACATCCAAGCCATCTg         >  2:2885769/1‑150 (MQ=255)
                                                                                                                                             gCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGCTATCCGCCAGCATCTGACGCCAGGACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTGTGCGGGTTAATACCGCAGGCATGCACGCTATCCATGCCAGCTg         >  2:2943291/1‑150 (MQ=255)
                                                                                                                                             gCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACCTCCATGCCATCTg         >  1:1514718/1‑150 (MQ=255)
                                                                                                                                              ccccACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGc        <  2:1871401/150‑1 (MQ=255)
                                                                                                                                              ccccACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGc        <  1:2508298/150‑1 (MQ=255)
                                                                                                                                                ccACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCa      <  1:4083718/150‑1 (MQ=255)
                                                                                                                                                ccACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCa      <  2:1485286/150‑1 (MQ=255)
                                                                                                                                                ccACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCa      <  1:3114202/150‑1 (MQ=255)
                                                                                                                                                ccACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCa      <  1:2580146/150‑1 (MQ=255)
                                                                                                                                                ccACATCCCCCCCCGTACGGTAATCGAGCAGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTATGCGGTTTAATACCGCAGGCATGCACGACATCCAAGCCAACAGCCa      >  2:1721006/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGC  >  NC_000913/3727961‑3728258

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: