Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,045,126 C→A 7.7% intergenic (+25/+32) 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,1260CA 7.7% 342.6 / 15.0 118intergenic (+25/+32)bglA/ygfF6‑phospho‑beta‑glucosidase A/putative NAD(P)‑dependent oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (57/51);  new base A (5/4);  total (63/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.98e-01

Only 100 of 119 total aligned reads displayed.

GCAAACGCTACGGCTTTATCTATGTGAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTCCACCACGCTGCATGGGGA  >  NC_000913/3044977‑3045273
                                                                                                                                                     |                                                                                                                                                   
gCAAACGCTACGGCTTTATCTATGTGAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGTTGATTGCCAGCAACGGCGAGAAGCTTTAAGTAGATGAAGTACCGGATGCAATAc                                                                                                                                                     >  1:2208566/1‑150 (MQ=255)
 cAAACGCTACGGCTTTATCTATGTGAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACt                                                                                                                                                    <  2:693440/150‑1 (MQ=255)
  aaaCGCTACGGCTTTATCTATGTGAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACttg                                                                                                                                                  >  2:2053703/1‑150 (MQ=255)
     cGCTACGGCTTTATCTATGTGAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACttgtt                                                                                                                                                >  1:1029406/1‑150 (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)
              tttATCTATGTGAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAACAAGTATTGCATCCGGt                                                                                                                                       <  2:1145788/150‑1 (MQ=255)
                  tCTATGTGAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGa                                                                                                                                   >  1:150383/1‑150 (MQ=255)
                    tATGTGAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACa                                                                                                                                 >  1:1335516/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)
                                cATGACGACGGTACTGGCGATATGTCGCGTTCACGAAAGAAGAGCTTTAACTGGTACAAAGAGGCGATAGCCAGCAACGGCGAGAAGCTTTAAGTCGTTGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGttt                                                                                                                     <  1:1809200/150‑1 (MQ=255)
                                    acgacgGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTAtt                                                                                                                 <  2:1537369/150‑1 (MQ=255)
                                       acgGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTcc                                                                                                              <  2:1029406/150‑1 (MQ=255)
                                       acgGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTcc                                                                                                              >  2:1232572/1‑150 (MQ=255)
                                        cgGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTccc                                                                                                             >  1:932873/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)
                                               ggCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGATGCGTTAAGTCGATGAAGTACCGGAAGCAATAGTTGTTGCATCCCGTCAGACAACTTCCCTGTTGTATTTCCcgcccgcc                                                                                                      >  2:425443/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)
                                                    tatGTCGCGGGCACGTAAGATGAGCTTTAACTGGTACAAAGAGGTGAGTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATc                                                                                                 <  2:377145/150‑1 (MQ=255)
                                                      tGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGa                                                                                               <  1:2189654/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)
                                                                   tAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGAAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGt                                                                                  <  1:388582/150‑1 (MQ=255)
                                                                       aagaGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACg                                                                              <  1:1982655/150‑1 (MQ=255)
                                                                             tttAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAagag                                                                        >  1:2014494/1‑150 (MQ=255)
                                                                              ttAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAgg                                                                       >  1:1090748/1‑150 (MQ=255)
                                                                               tAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGc                                                                      >  2:166360/1‑150 (MQ=255)
                                                                               tAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGc                                                                      <  1:1866580/150‑1 (MQ=255)
                                                                                 aCTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCCTCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCGCGTGACGGAAGAGGCtt                                                                    >  2:230526/1‑150 (MQ=255)
                                                                                 aCTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCtt                                                                    >  1:531671/1‑150 (MQ=255)
                                                                                 aCTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCAGATGCAATACTTGTTGCATCCGTTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCtt                                                                    <  1:370203/150‑1 (MQ=255)
                                                                                    ggTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGAAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTAt                                                                 <  1:1330717/150‑1 (MQ=255)
                                                                                     gTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATc                                                                <  1:1681729/150‑1 (MQ=255)
                                                                                       aCAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGTTTTATCAc                                                              <  1:1467625/150‑1 (MQ=255)
                                                                                        cAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACt                                                             >  1:601541/1‑150 (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)
                                                                                                   ttGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGa                                                  >  2:1593509/1‑150 (MQ=255)
                                                                                                     gCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACATCTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGAGGTAAGAGGCTTTATCACTTAGTAGCCAGACa                                                >  2:1337789/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)
                                                                                                                                       tACCGGATGCAATACTTTTTTCATACTGTCTTACAAATTCCCTGTTATAATTCAGGCCCGCCAAATCGATCAAACTTCCCTTGACGTCAGAGGCTTTAACACTTAGTAGACAGACAATGGCCTGCGCGACCTCTTCTGCCTGTCcaccac              <  2:1337241/150‑1 (MQ=255)
                                                                                                                                         ccGGATGCAATACTTGTTTCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCGCGCCAAATAGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTGAGTAGCGAGACAATGGCCTTCGCGACCTCTTCTGCCTGTGAAACACGc            >  2:350105/1‑150 (MQ=255)
                                                                                                                                         ccGGATGCAACAAGTATTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTCCACCACGc            >  1:431139/1‑150 (MQ=255)
                                                                                                                                         ccGGATGCAACAAGTATTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGATTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTCCACCACgg            >  2:954573/1‑149 (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)
                                                                                                                                                   tACTTGTTGCATCCGGTCAGACAACTACCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTCCACCACGCTGCATGGGGa  <  2:1853513/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GCAAACGCTACGGCTTTATCTATGTGAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTCCACCACGCTGCATGGGGA  >  NC_000913/3044977‑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: