Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,156,591 T→G 6.2% intergenic (+73/‑32) yqhD → / → dkgA aldehyde reductase, NADPH‑dependent/2,5‑diketo‑D‑gluconate reductase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,156,5910TG 6.2% 228.6 / 2.8 102intergenic (+73/‑32)yqhD/dkgAaldehyde reductase, NADPH‑dependent/2,5‑diketo‑D‑gluconate reductase A
Reads supporting (aligned to +/- strand):  ref base T (47/47);  new base G (2/6);  total (49/53)
Fisher's exact test for biased strand distribution p-value = 2.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 102 total aligned reads displayed.

GGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCG  >  NC_000913/3156444‑3156733
                                                                                                                                                   |                                                                                                                                              
ggCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCTCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGATAACAACACCTCACCGGAGCCTg                                                                                                                                               >  1:405296/1‑149 (MQ=255)
 gCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGc                                                                                                                                              <  2:20520/149‑1 (MQ=255)
  cATGACCCAACTGGGCGAAAATCAAGACATTACGTTGGATGGCAGCCCCCGTATATACGAAGCCGCCCGATAAGCTTTTTACGCATCAAACTGTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGTAGCCTgat                                                                                                                                             >  1:306058/1‑147 (MQ=255)
      aCCCAACTGGGCGAAAATCCAGACATTACGTGGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATCTCGGCCAGACTTAAGTTCACAACCCCTCACCGGCGCCTGCTCCgg                                                                                                                                         <  1:867128/149‑1 (MQ=255)
      aCCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCAGGCTCCgg                                                                                                                                         <  1:1206087/149‑1 (MQ=255)
       cccAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGt                                                                                                                                        >  2:1131984/1‑149 (MQ=255)
       cccAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGt                                                                                                                                        <  1:574629/149‑1 (MQ=255)
       cccAACTGGGCGAAAATCATGACATTACGTTGGATGGCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACTCCTCACCGGAGCCTGCTCCGGt                                                                                                                                        <  1:1223315/149‑1 (MQ=255)
        ccAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTg                                                                                                                                       >  1:838078/1‑149 (MQ=255)
         cAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGTGCTTTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGa                                                                                                                                      >  2:1394044/1‑149 (MQ=255)
         cAACTGGGCGAAAATCATGACATAACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCAGGCTCCGGTGa                                                                                                                                      <  1:715834/149‑1 (MQ=255)
          aaCTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAg                                                                                                                                     <  1:989183/149‑1 (MQ=255)
           aCTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCTCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGt                                                                                                                                    <  1:451705/149‑1 (MQ=255)
             tGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTc                                                                                                                                  >  1:1198932/1‑149 (MQ=255)
              gggCGAAAATCATGACATTACGTTGGATGGCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCACAAACATTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCa                                                                                                                                 <  2:385637/149‑1 (MQ=255)
               ggCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCat                                                                                                                                <  1:141452/149‑1 (MQ=255)
                gCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCata                                                                                                                               <  1:1204736/149‑1 (MQ=255)
                  gAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCataaa                                                                                                                             >  2:494426/1‑147 (MQ=255)
                  gAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCAtata                                                                                                                             >  1:1306951/1‑149 (MQ=255)
                    aaaTCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCACCGGTGAGTTCATATaaa                                                                                                                           >  2:1007100/1‑149 (MQ=255)
                       tCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGGCGGCTCCGGTGAGTTCATATaattga                                                                                                                        >  2:186394/1‑145 (MQ=255)
                       tCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAagga                                                                                                                        <  1:435355/149‑1 (MQ=255)
                         aTGACATTACGTTGGATGTCAGCCCCAGTATATACGAAGCCGCCCGCTAAGATTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTACGTTCACAACACCTCACCGCAGCCTGATCCGTTGAGATCATATAAaggacg                                                                                                                      >  1:845917/1‑147 (MQ=255)
                          tGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGTGTTCATATAAAggagga                                                                                                                     >  1:1056243/1‑149 (MQ=255)
                            aCATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAAc                                                                                                                   <  1:222875/149‑1 (MQ=255)
                              aTTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGt                                                                                                                 <  2:974610/149‑1 (MQ=255)
                                tACGTTGGCTGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAACCTTTCCTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTGCTTATAAAGGAGGAACGtct                                                                                                               >  2:435597/1‑147 (MQ=255)
                                 aCGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATg                                                                                                              >  1:2373/1‑149 (MQ=255)
                                  cGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATgg                                                                                                             <  1:1152758/149‑1 (MQ=255)
                                   gTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGc                                                                                                            >  1:75249/1‑149 (MQ=255)
                                     tGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTATCGTTTTCGGGCATTTCGGCCAGACTTAAGTTCACAACACCTCACCCGAGCCTGCTCCGGGGAGTTCATATAAAGGAGGAACGTATGGCTa                                                                                                          >  1:404209/1‑149 (MQ=255)
                                      ggATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTaa                                                                                                         <  2:1236004/149‑1 (MQ=255)
                                       gATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCTCCGTAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAAt                                                                                                        <  2:348179/149‑1 (MQ=255)
                                       gATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAAt                                                                                                        >  1:736039/1‑149 (MQ=255)
                                        aTGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCATCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATc                                                                                                       >  2:302112/1‑149 (MQ=255)
                                           tCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGATCACAACACCTCACCGCAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCac                                                                                                    <  1:398577/149‑2 (MQ=255)
                                            cAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAAc                                                                                                   <  2:761860/149‑1 (MQ=255)
                                             aGCCGCCGTATATACGAAGCCGCCCGCGAAGCTTTTTACGCCTCAAACTTTCGGTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTATTCCCAcc                                                                                                  >  2:597508/1‑149 (MQ=255)
                                               ccgccgTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCCACCGt                                                                                                >  2:995018/1‑149 (MQ=255)
                                                      atataCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATTAAGGAGGAACGGATGGCGAATCCTACGGTTATTAAg                                                                                         >  2:1160034/1‑149 (MQ=255)
                                                      atataCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCATCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAg                                                                                         <  1:302112/149‑1 (MQ=255)
                                                        ataCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGGGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCt                                                                                       >  1:691054/1‑149 (MQ=255)
                                                            gAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGGTCAGATTAAGGCGGCCCGTATGGCTAATCCAACCGTTATTAAGCTTCAg                                                                                   >  1:1426304/1‑149 (MQ=255)
                                                               gccgccAGCTAAGCTGTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGAATTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTAGGGCTAATCCAACCGTTATTAAGCTACAGGAt                                                                                <  1:1342822/149‑1 (MQ=255)
                                                                 cgcccgcTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATgg                                                                              <  2:609656/149‑1 (MQ=255)
                                                                  gccctcTAAGGTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGc                                                                             <  1:1503881/149‑1 (MQ=255)
                                                                  gcccgcTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGc                                                                             <  1:1175196/149‑1 (MQ=255)
                                                                     cgcTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAAt                                                                          <  1:681246/149‑1 (MQ=255)
                                                                          aGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCATACTTAAGTTCAAAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTTTGGCTAATACAACCGTTATTAAGCTACAGGTTGTCAATGTCAt                                                                     >  2:1246387/1‑149 (MQ=255)
                                                                            cTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCTGTGAGTTCATATAATGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGGAATGTCATGc                                                                   >  2:1511691/1‑149 (MQ=255)
                                                                              ttttACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGccc                                                                 <  2:742464/149‑1 (MQ=255)
                                                                              ttttACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGccc                                                                 >  1:1521819/1‑149 (MQ=255)
                                                                                    gCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACCCCACCTCACCGGAGCCTGCTCCGGGGAGGTCAGATAAAGGAGGAACGTATGGCTTCTCCACCCGGTATTAAGCTCCAGGATGGCAATGTCATGCCCCAGCTg                                                           >  2:264588/1‑149 (MQ=255)
                                                                                     ccccAAACTTGCGATTTCGGGCATTTCGGCCAGACTTAAGTTCACAACACCTCACCGCAGCCGGCTCCGGTGAGTTCATATAACCGAGGCACGTATGGCTACTCCAACCGTTATTAAGCTACAGCCTGGCAATGGCATGCCCCAGCTgg                                                          <  2:1448039/146‑1 (MQ=255)
                                                                                     ccTCAACCTTTCGTTTTGGGAGATTTCGTCCATACTTAAGTTCACAACACCTCTCCGGAGCCTGCTCAGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGGTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTgg                                                          <  2:805180/149‑1 (MQ=255)
                                                                                     ccTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTgg                                                          <  1:937885/149‑1 (MQ=255)
                                                                                     ccTCAAACTTTCGATTTCGGGCATTTCCTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATTAAGGAGCAACGTATGGCTAATACAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTgg                                                          <  2:230808/149‑1 (MQ=255)
                                                                                     cagcAAACTTTCCTTTTCGTGCATTTCGTCCAGAATTAAGTTCACAACACCTCCCCGTCGCCGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAACATACAGGATTGCAATTTCATGCCCCAGCTgg                                                          <  2:1013319/146‑1 (MQ=255)
                                                                                      cTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTggg                                                         >  1:376290/1‑149 (MQ=255)
                                                                                        cAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCAAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGAc                                                       >  1:1549627/1‑149 (MQ=255)
                                                                                         aaaCTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCGAATCCCCCCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGagt                                                      >  2:1441179/1‑147 (MQ=255)
                                                                                         aaaCCTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACt                                                      <  2:425521/149‑1 (MQ=255)
                                                                                          aaCTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTg                                                     <  2:801351/149‑1 (MQ=255)
                                                                                           aCTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTgg                                                    <  1:580511/149‑1 (MQ=255)
                                                                                             tttCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACATCAACTCACCGGAGCCTGCTCCGGTGAGTTCAAATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGc                                                  <  2:404209/149‑1 (MQ=255)
                                                                                             tttCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCCACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGc                                                  >  2:1228103/1‑149 (MQ=255)
                                                                                             tttCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGc                                                  >  2:863731/1‑149 (MQ=255)
                                                                                             tttCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGc                                                  <  1:443379/149‑1 (MQ=255)
                                                                                              ttCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCg                                                 >  1:576507/1‑149 (MQ=255)
                                                                                                 gTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCt                                              <  2:415089/149‑1 (MQ=255)
                                                                                                 cttttCGGGACTTTCGTCAAGAGATATGTTCCCAACCCCTCACCGGAGCCTGCTCCGGTGAGTTCATAGACAGGCCGACCGTATGGCTCCTCCCACCGTTAATCCGCTACAGGATGGAACTGTCATGCCCCAGCGGGGACTGGGCGTcg                                              <  2:252315/148‑2 (MQ=255)
                                                                                                  ttttCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGATTTCATATAAAGGAGTAACGTATGGCTAATCAAACCGTTATTAAGCTACAGTATGGCAATGTCATCACCCAGCTGGGACTGGCCGACTg                                             >  2:1564375/1‑149 (MQ=255)
                                                                                                  ttttCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTg                                             <  2:2373/149‑1 (MQ=255)
                                                                                                   tttCGGGCATTTCGTCCAGACTAAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTgg                                            <  2:1567835/149‑1 (MQ=255)
                                                                                                    ttCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGgc                                           <  2:1134561/149‑1 (MQ=255)
                                                                                                      cGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGCATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGgcat                                         >  1:622632/1‑148 (MQ=255)
                                                                                                           aTTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCt                                              <  2:298988/139‑1 (MQ=255)
                                                                                                           aTTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCt                                              >  1:298988/1‑139 (MQ=255)
                                                                                                            tttCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGATCGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGt                                   <  1:388861/149‑1 (MQ=255)
                                                                                                               cGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAAt                                <  1:1177143/149‑1 (MQ=255)
                                                                                                                gTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGCGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCCGCTGGGACTGGGCGTCTGGCAAGCAAGTAATg                               >  2:1069061/1‑149 (MQ=255)
                                                                                                                 tCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGa                              <  1:824955/149‑1 (MQ=255)
                                                                                                                    aGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTGATATAAAGGAGGAAGGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGGTGGGACTGGGCGTCTGGCATGCTAGTAATGAggt                           >  2:247523/1‑148 (MQ=255)
                                                                                                                    aGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGa                           >  1:760847/1‑149 (MQ=255)
                                                                                                                     gACTTAATGTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAATCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGGCTGGCAAGCAAGTAATGAGGaa                          <  2:1548198/149‑1 (MQ=255)
                                                                                                                     gACTTAAGTTCACAACACCTCACCGCAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGaa                          <  2:736039/149‑1 (MQ=255)
                                                                                                                        ttAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTa                       <  2:838078/149‑1 (MQ=255)
                                                                                                                        ttAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTa                       >  1:666691/1‑149 (MQ=255)
                                                                                                                        ttAAGTTCACAACACCTCACCGGAGCCGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACAGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTATGGCAAGCAAGTAATGAGGAAGTa                       <  1:1246387/149‑1 (MQ=255)
                                                                                                                            gTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCGAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCa                   >  1:421172/1‑149 (MQ=255)
                                                                                                                             ttCACAACACCTCACCGGAGCCTGCTCCGCTGAGTTCATATAAAGGAGGAACGTAAGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCAGGGACTGGGCCTCTGGCAAGCAAGAAATGAGGAAGTAATCAc                  >  1:24573/1‑149 (MQ=255)
                                                                                                                                acaacaCCTCACCGGAGACTGCTCCGGTGAGTTCATATAAAGTAGGAACGTATGGCTAATCCAACCGTTATTAAGCTCCAGTATGGCAATGGCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGc               <  2:405296/149‑1 (MQ=255)
                                                                                                                                  aacaCCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCACAGCTGGTACTGGGCGTCTGGCATGCAAGTAATGAGGAAGTATTCACCGCCa             >  1:403910/1‑149 (MQ=255)
                                                                                                                                   acacCTCACCGGAGCCTGCTCCGGTGAGTTCATTTAAAGGAGGAACGTATGGCTAATCCAACCGTTATTTAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCAt            >  1:1052768/1‑149 (MQ=255)
                                                                                                                                       cTCACCGGCGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGAAAGTAATGAGGAAGTAATCACCGCCATTCaa        >  2:1156987/1‑149 (MQ=255)
                                                                                                                                       cTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCaa        >  1:1075553/1‑149 (MQ=255)
                                                                                                                                        tCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAAAGTATGGCTAATCCAACCGGTATTAAGATACAGGATGGCAATGTCATGCCCCAGCAGGGACAGGGCGTCTGGCAAGCAAGTAATAAGGAAGTAATCCCCGCCAGTCaaa       >  1:980320/1‑149 (MQ=255)
                                                                                                                                        tCACCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCaaa       >  2:1207235/1‑149 (MQ=255)
                                                                                                                                          aCCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCaaaaa     >  2:644743/1‑149 (MQ=255)
                                                                                                                                             ggAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCg  <  1:1007100/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
GGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCG  >  NC_000913/3156444‑3156733

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: