Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,045,127 T→G 8.6% intergenic (+26/+31) bglA → / ← ygfF 6‑phospho‑beta‑glucosidase A/putative NAD(P)‑dependent oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,045,1270TG 8.6% 313.4 / 17.6 113intergenic (+26/+31)bglA/ygfF6‑phospho‑beta‑glucosidase A/putative NAD(P)‑dependent oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (55/48);  new base G (6/4);  total (61/52)
Fisher's exact test for biased strand distribution p-value = 7.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.84e-01

Only 100 of 118 total aligned reads displayed.

CAAACGCTACGGCTTTATCTATGTGAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTCCACCACGCTGCATGGGGA  >  NC_000913/3044978‑3045273
                                                                                                                                                     |                                                                                                                                                  
cAAACGCTACGGCTTTATCTATGTGAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACt                                                                                                                                                    <  2:693440/150‑1 (MQ=255)
    cGCTACGGCTTTATCTATGTGAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACttgtt                                                                                                                                                >  1:1029406/1‑150 (MQ=255)
     gCTACGGCTTTATCTATGTGAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACttgttg                                                                                                                                               <  1:1614613/150‑1 (MQ=255)
     gCTACGGCTTTATCTATGTGAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACttgttg                                                                                                                                               <  1:1443790/150‑1 (MQ=255)
        aCGGCTTTATCTATGTGAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCAt                                                                                                                                            >  1:64182/1‑150 (MQ=255)
         cGGCTTTATCTATGTGAATAAACATGACGACCGTACTGTCGATATGTCGCGTTAACGTAAGAATAGCTTTATCTGGTACAAATAGGTGATTGACAGCAACGGCGAGAAGATTTAAGTCGATGAAATACCGGATTCAATACTAGTTGCAta                                                                                                                                           >  1:1075977/1‑149 (MQ=255)
          ggCTTTATCTATGTTAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGGGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAACAAGTATTGCATcc                                                                                                                                          <  1:1127550/150‑1 (MQ=255)
          ggCTTTATCTATGTGAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAACAAGTATTGCATcc                                                                                                                                          <  2:1771308/150‑1 (MQ=255)
           gCTTTATCTATGTGAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCGTTAACTGGTTAAAAGAGGTGATTGCCAGCAACGGCGAGCAGCTGTAAGTCGATGAAGTACCGGATGCAATACTAGTTGCATCCg                                                                                                                                         >  2:381450/1‑150 (MQ=255)
            cTTTATCTATGTGAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAACAAGTATTGCATCCgg                                                                                                                                        <  1:27084/150‑1 (MQ=255)
               tATCTATGTGAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCa                                                                                                                                     <  2:4515/150‑1 (MQ=255)
               tATCTATGTGAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCa                                                                                                                                     <  1:2127200/150‑1 (MQ=255)
                 tCTATGTGAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGa                                                                                                                                   >  1:150383/1‑150 (MQ=255)
                   tATGTGAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACa                                                                                                                                 >  1:1335516/1‑150 (MQ=255)
                   tATGTGAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTGTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACa                                                                                                                                 >  2:2071161/1‑150 (MQ=255)
                     tgtgAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCTGTCAGACaaa                                                                                                                               >  2:846564/1‑149 (MQ=255)
                     tgtgAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAAc                                                                                                                               <  1:566311/150‑1 (MQ=255)
                         aaTAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTcc                                                                                                                           <  2:2213361/150‑1 (MQ=255)
                          aTAAACATCACGACGTTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTccc                                                                                                                          <  1:230142/150‑1 (MQ=255)
                             aaCATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGt                                                                                                                       >  2:1021529/1‑150 (MQ=255)
                                   acgacgGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTAtt                                                                                                                 <  2:1537369/150‑1 (MQ=255)
                                      acgGTACTGGCGATATGTCGCGTTCAGGTAAGAATAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGTAGCTTTAAGTCGATGAAGTACCGGCTGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTcc                                                                                                              <  1:1462438/150‑1 (MQ=255)
                                      acgGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTcc                                                                                                              <  2:1881771/150‑1 (MQ=255)
                                      acgGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTcc                                                                                                              <  2:1029406/150‑1 (MQ=255)
                                      acgGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTcc                                                                                                              >  2:1232572/1‑150 (MQ=255)
                                            cTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTcccgcccg                                                                                                        <  2:1848982/150‑1 (MQ=255)
                                            cTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTcccgcccg                                                                                                        <  2:1273022/150‑1 (MQ=255)
                                             tGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATACGGTCAGACAACTTACCTGTTTTATTTCcagccagc                                                                                                       >  1:181469/1‑150 (MQ=255)
                                              ggCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCcgcccgcc                                                                                                      >  1:1919095/1‑150 (MQ=255)
                                              ggCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCcgcccgcc                                                                                                      >  1:2102100/1‑150 (MQ=255)
                                              ggCGATATGTCGCGTGCACGTAAGAAGCGCGGTAACTGGTACAAAGAGGTGATTGCCAGCACCGGCGAGAAGCTTTATGTCGATGAAGTACCGGATGCAATACGGGTTGCATCCGGTCAGACAACTTCGCTGTGTTATTTCCcgcccgcc                                                                                                      >  2:767076/1‑150 (MQ=255)
                                                   tatGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATc                                                                                                 >  2:2164632/1‑150 (MQ=255)
                                                   tatGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATc                                                                                                 <  1:554111/150‑1 (MQ=255)
                                                   tatGTCGCGGGCACGTAAGATGAGCTTTAACTGGTACAAAGAGGTGAGTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATc                                                                                                 <  2:377145/150‑1 (MQ=255)
                                                      gTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGAt                                                                                              >  1:1414779/1‑150 (MQ=255)
                                                      gTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGAt                                                                                              >  2:1918488/1‑150 (MQ=255)
                                                      gTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGAt                                                                                              >  1:636047/1‑150 (MQ=255)
                                                       tCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATa                                                                                             <  2:64182/150‑1 (MQ=255)
                                                        cgcgTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATaa                                                                                            >  2:625525/1‑150 (MQ=255)
                                                           gTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAAc                                                                                         <  1:1611718/150‑1 (MQ=255)
                                                              cACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTc                                                                                      >  1:1276646/1‑150 (MQ=255)
                                                                 gTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCg                                                                                   >  2:2142199/1‑150 (MQ=255)
                                                                            tttAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAagag                                                                        >  1:2014494/1‑150 (MQ=255)
                                                                             ttAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAgg                                                                       >  1:1090748/1‑150 (MQ=255)
                                                                              tAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGc                                                                      <  2:135925/150‑1 (MQ=255)
                                                                              tAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGc                                                                      <  1:1866580/150‑1 (MQ=255)
                                                                              tAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGc                                                                      >  2:166360/1‑150 (MQ=255)
                                                                                aCTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCCTCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCGCGTGACGGAAGAGGCtt                                                                    >  2:230526/1‑150 (MQ=255)
                                                                                aCTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCtt                                                                    >  1:531671/1‑150 (MQ=255)
                                                                                aCTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCAGATGCAATACTTGTTGCATCCGTTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCtt                                                                    <  1:370203/150‑1 (MQ=255)
                                                                                 cTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCttt                                                                   >  1:818413/1‑150 (MQ=255)
                                                                                  tGGTACAATGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATTGAATACTTGTTGCATGCGGTCAGACATCTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTGCCGTGACGTAAGAGGCTTTa                                                                  <  2:123012/150‑1 (MQ=255)
                                                                                   ggTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGAAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTAt                                                                 <  1:1330717/150‑1 (MQ=255)
                                                                                       cAAAGACGTGATTGCCAGCAACGGCGAGAAGCTATAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCATGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACt                                                             >  1:1514479/1‑150 (MQ=255)
                                                                                        aaaGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACtt                                                            <  2:907902/150‑1 (MQ=255)
                                                                                        aaaGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACtt                                                            >  2:1978105/1‑150 (MQ=255)
                                                                                         aaGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTta                                                           <  1:16837/150‑1 (MQ=255)
                                                                                           gagGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACACCTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTtagt                                                         >  1:177057/1‑150 (MQ=255)
                                                                                           gagGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGAAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCACGTGACGTAAGAGGCTTTATCACTtagt                                                         >  1:1832598/1‑150 (MQ=255)
                                                                                              gTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGc                                                      >  1:979413/1‑150 (MQ=255)
                                                                                               tGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGcc                                                     <  1:1016697/150‑1 (MQ=255)
                                                                                               ggATTGCCAGCAACGGCGAGAAGCTTTAAGGCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGcc                                                     <  1:1773686/149‑1 (MQ=255)
                                                                                                 aTTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAg                                                   <  2:2208566/150‑1 (MQ=255)
                                                                                                    gCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACATCTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGAGGTAAGAGGCTTTATCACTTAGTAGCCAGACa                                                >  2:1337789/1‑150 (MQ=255)
                                                                                                     ccAGCAACGGCGAGAAGCTTTAAGTCTATGAAGTACCGCATGCAATACTTGATGCATCCGGTCAGACAACGACCCTGTTTTATTTCCCGCCCGCCAACTCGATAAAAATTCCAGTCACGTAAGAGGCTTAATCACTTAGTAGCCAGACaa                                               >  1:345273/1‑150 (MQ=255)
                                                                                                        gCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCTGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAACCTTCCCGTGACGTATGAGGCTTTATCACTTAGTAGCCAGACAATgg                                            <  2:94959/150‑1 (MQ=255)
                                                                                                           aCGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCt                                         <  1:850269/150‑1 (MQ=255)
                                                                                                                gagaAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGcgcg                                    <  2:932873/150‑1 (MQ=255)
                                                                                                                gagaAGCTTTAAGTCGATGAAGTAACGGATGCAATACTTGTTGCATCCGGGCAGACAACTTCCATGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGcgcg                                    <  2:193865/150‑1 (MQ=255)
                                                                                                                 agaAGCTTTAAGTCGATGAAGTACCGGATGCAATACTGGTTGCCTCCGGTCAGACAACTTCCATGTTTTATTTCCCGCCCGCCAAATCGATTAACCTTCCCGGGACGTAAGAGGCTTTATCACTTAGTAGCCCGACAATGTCATGCGCGa                                   >  2:417917/1‑150 (MQ=255)
                                                                                                                   aaGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAACTCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGAcc                                 <  2:1459615/150‑1 (MQ=255)
                                                                                                                   aaGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCCGACAATGGCCTGCGCGAcc                                 <  2:1823425/150‑1 (MQ=255)
                                                                                                                   aaGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGGCAATGGCCTGCGCGAcc                                 <  1:53579/150‑1 (MQ=255)
                                                                                                                   aaGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGAcc                                 <  1:1014704/150‑1 (MQ=255)
                                                                                                                    aGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGAcat                                >  2:408556/1‑148 (MQ=255)
                                                                                                                     gCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGCCAATGGCCTGCGCGACctc                               >  1:2161787/1‑150 (MQ=255)
                                                                                                                         tAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCtcttct                           >  1:1721041/1‑150 (MQ=255)
                                                                                                                          aaGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTg                          >  2:1888953/1‑150 (MQ=255)
                                                                                                                          aaGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTg                          <  2:723597/150‑1 (MQ=255)
                                                                                                                           aGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGc                         >  1:2126128/1‑150 (MQ=255)
                                                                                                                           aGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGc                         >  1:1694896/1‑150 (MQ=255)
                                                                                                                            gTCGATGTAGTACCGGATGCAATTCTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGATTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTgtc                        <  2:1151723/150‑3 (MQ=255)
                                                                                                                             tCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGccg                       >  1:880580/1‑149 (MQ=255)
                                                                                                                             tCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCt                       >  1:948736/1‑150 (MQ=255)
                                                                                                                              cGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTg                      >  1:1568299/1‑150 (MQ=255)
                                                                                                                              cGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTGACCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGGGACGTAAGAGGCTTTTTCGCTGAGTAGCCAGACAATGGCCTGCGAGACCTCTTCTGCCtt                      >  2:1442976/1‑149 (MQ=255)
                                                                                                                                 tGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAACTCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTcc                   <  2:1075157/150‑1 (MQ=255)
                                                                                                                                  gAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAAATTCGCTGATTTATTACCCGCACGCCAAATCGATAAAACTTCCGCTGACGTAAGAGGCTTTATAACTTAGTAGCCAGACAATGGCCTGCGCGACCTCGTCTGCCTTTcca                  >  2:1293663/1‑150 (MQ=255)
                                                                                                                                     tgaCCGGATGCAACAAGTATTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTccaccc               >  1:1881017/3‑149 (MQ=255)
                                                                                                                                     tgaCCGGATGCAACAAGTATTGCATCAGGTCAGACAACTTACCTTTTTTATTTCCCGCCCGCCAAATCGATAAAAATTCACGTGAAGTAAGAGGCTTGATCACTTAGTAGACAGACAATGGCCGGCGCGACCTAGTCTGCCTGTccacca               >  2:1721321/3‑150 (MQ=255)
                                                                                                                                        ccGGATGCAATACTTGTTGCATCCGGTCAGACAACATCCCTGTTTTATTTCCCGCCCGACAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTCCACCACGc            >  2:350170/1‑150 (MQ=255)
                                                                                                                                        ccGGATGCAACAAGTATTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTCCACCACGc            >  1:431139/1‑150 (MQ=255)
                                                                                                                                        ccGGATGCAACAAGTATTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTCCACCACGc            >  1:1789415/1‑150 (MQ=255)
                                                                                                                                        ccGGATGCAACAAGTATTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGATTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTCCACCACgg            >  2:954573/1‑149 (MQ=255)
                                                                                                                                          ggATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTCCACCACGCTg          >  2:1699325/1‑150 (MQ=255)
                                                                                                                                            aTGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGACTGCGCGACCTCTTCTGCCTGTCCACCACGCTGCa        >  1:2096138/1‑150 (MQ=255)
                                                                                                                                              gCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTCCACCACGCTGCATg      >  1:591714/1‑150 (MQ=255)
                                                                                                                                                attACTTGTTGCATCCGTTCAGACTACTTCCCTGTTTTATTTCCCGCCAGCCCCATCGCTCACAGTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGCCAATGGCCTGCGCGACCTCTTCTGCCTGTCCACCACGCTGCATggg    <  1:381450/148‑1 (MQ=255)
                                                                                                                                                aaTACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTCCACCACGCTGCATggg    <  2:1038476/150‑1 (MQ=255)
                                                                                                                                                 aTACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTCCACCACGCTGCATgggg   <  2:1335516/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CAAACGCTACGGCTTTATCTATGTGAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTCCACCACGCTGCATGGGGA  >  NC_000913/3044978‑3045273

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: