Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,606,480 G→C 5.1% M127I (ATG→ATC)  lsrF → 3‑hydroxy‑2,4‑pentadione 5‑phosphate thiolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,606,4800GC 5.1% 441.0 / 9.5 129M127I (ATG→ATC) lsrF3‑hydroxy‑2,4‑pentadione 5‑phosphate thiolase
Reads supporting (aligned to +/- strand):  ref base G (68/54);  new base C (2/5);  total (70/59)
Fisher's exact test for biased strand distribution p-value = 2.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 138 total aligned reads displayed.

TGTTTGAACATGCCGATGTATTAATGTGTACGCGCGGCATTTTGCGCAGCGTAGTTCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGTTGATGCCGGAATGAAAGTGGGAATGCCGACCATGGCCGTGACTGGCGTG  >  NC_000913/1606332‑1606627
                                                                                                                                                    |                                                                                                                                                   
tGTTTGAACATGCCGATGTATTAATGTGTACGCGCGGCATTTTGCGCAGGGTAGTTCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCgatgg                                                                                                                                                    <  2:1500828/150‑1 (MQ=255)
 gTTTGAACATGCCGATGTATTAATGTGTACGCGCGGCATTTTGCGCAGCGTAGTTCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCgatgga                                                                                                                                                   >  1:18339/1‑150 (MQ=255)
 gTTTGAACATGCCGATGTATTAATGTGTACGCGCGGCATTTTGCGCAGCGTAGTTCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCgatgga                                                                                                                                                   <  1:2293960/150‑1 (MQ=255)
    tGAACATGCCGATGTATTAATGTGTACGCGCGGCATTTTGCGCAGCGTAGTTCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGa                                                                                                                                                >  1:2277934/1‑150 (MQ=255)
    tGAACATGCCGATGTATTAATGTGTACGCGCGGCATTTTGCGCAGCGTAGTTCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGa                                                                                                                                                >  2:751528/1‑150 (MQ=255)
       aCATGCCGATGTATTAATGTGTACGCGCGGCATTTTGCGCAGCGTAGTTCCCCCTGCGACCTATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTCATCCATCGATAACGc                                                                                                                                             <  1:382445/150‑1 (MQ=255)
       aCATGCCGATGTATTAATGTGTACGCGCGGCATTTTGCGCAGCGTAGTTCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGc                                                                                                                                             >  1:1014114/1‑150 (MQ=255)
          tGCCGATGTATTAATGTGTACGCGCGGCATTTTGCGCAGCGTAGTTCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGt                                                                                                                                          >  2:1340939/1‑150 (MQ=255)
           gCCGATGTATTAATGTGTACGCGCGGCATTTTGCGCAGCGTAGTTCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTCATCCATCGATAACGccacg                                                                                                                                         <  1:1795905/150‑4 (MQ=255)
            ccGATGTATTAATGTGTACGCGCGGCATTTTGCGCAGCGTAGTTCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTCATCCATCGATAACGccacgg                                                                                                                                        <  1:1439163/150‑5 (MQ=255)
            ccGATGTATTAATGTGTACGCGCGGCATTTTGCGCAGCGTAGTTCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTCATCCATCGATAACGccacgg                                                                                                                                        <  1:2426407/150‑5 (MQ=255)
             cGATGTATTAATGTGTACGCGCGGCATTTTGCGCAGCGTAGTTCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTgcg                                                                                                                                       >  2:1178993/1‑150 (MQ=255)
             cGATGTATTAATGTGTACGCGCGGCATTTTGCGCAGCGTAGTTCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTgcg                                                                                                                                       <  1:457929/150‑1 (MQ=255)
             cGATGTATTAATGTGTACGCGCGGCATTTTGCGCAGCGTAGTTCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTCATCCATCGATAACGccacggc                                                                                                                                       <  1:166574/150‑6 (MQ=255)
               aTGTATTAATGTGTACGCGCGGCATTTTGCGCAGCGTAGTTCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGcc                                                                                                                                     <  2:1580508/150‑1 (MQ=255)
                tgcaTTAATGTGTACGCGCGGCATTTTGCGCAGCGTAGTTCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCt                                                                                                                                    <  1:2188287/147‑1 (MQ=255)
                 gTATTAATGTGTACGCGCGGCATTTTGCGCAGCGTAGTTCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTg                                                                                                                                   <  1:379540/150‑1 (MQ=255)
                  tATTAATGTGTACGCGCGGCATTTTGCGCAGCGTAGTTCCCCCTGCGACCAATAGGGCGGTGGTACTGCGGGCGTCAGGTGCGAACTATATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGa                                                                                                                                  >  1:1797628/1‑150 (MQ=255)
                  tATTAATGTGTACGCGCGGCATTTTGCGCAGCGTAGTTCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGa                                                                                                                                  >  2:29434/1‑150 (MQ=255)
                  tATTAATGTGTACGCGCGGCATTTTGCGCAGCGTAGTTCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGa                                                                                                                                  >  1:1369245/1‑150 (MQ=255)
                   aTTAATGTGTACGCGCGGCATTTTGCGCAGCGTAGTTCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGaa                                                                                                                                 >  1:2406369/1‑150 (MQ=255)
                    ttAATGTGTACGCGCGGCATTTTGCGCAGCGTAGTTCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAAc                                                                                                                                <  2:1769297/150‑1 (MQ=255)
                      aaTGTGTACGCGCGGCATTTTGCGCAGCGTAGTTCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAg                                                                                                                              >  2:408860/1‑150 (MQ=255)
                        tgtgTACGCGCGGCATTTTGCGCAGCGTAGTTCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGtt                                                                                                                            >  2:2127908/1‑150 (MQ=255)
                         gtgtACGCGCGGCATTTTGCGCAGCGTAGTTCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTg                                                                                                                           >  1:1464897/1‑150 (MQ=255)
                         gtgtACGCGCGGCATTTTGCGCAGCGTAGTTCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTg                                                                                                                           >  2:568637/1‑150 (MQ=255)
                          tgtACGCGCGGCATTTTGCGCAGCGTAGTTCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTgc                                                                                                                          >  2:1197239/1‑150 (MQ=255)
                           gtACGCGCGGCATTTTGCGCAGCGTAGTTCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTgcg                                                                                                                         >  2:1379039/1‑150 (MQ=255)
                            tACGCGCGGCATTTTGCGCAGCGTAGTTCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTTTCGTAATTATGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCTCCTTAACAGTTtcgc                                                                                                                        <  2:1150971/150‑1 (MQ=255)
                             aCGCGCGGCATTTTGCGCAGCGTAGTTCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGcgcg                                                                                                                       >  2:2073939/1‑150 (MQ=255)
                              cgcgcgGCATTTTGCGCAGCGTAGTTCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCgg                                                                                                                      >  1:1502004/1‑150 (MQ=255)
                                cgcgACATTTTGCGCAGCGTAGTTCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTg                                                                                                                    <  1:891999/150‑1 (MQ=255)
                                  cgGCATTTTGCGCAGCGTAGTTCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTggc                                                                                                                  >  2:22498/1‑150 (MQ=255)
                                      aTTTTGCGCAGCGTAGTTCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGgcggcg                                                                                                              >  1:565775/1‑150 (MQ=255)
                                        tttGCGCAGCGTAGTTCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCa                                                                                                            <  2:18339/150‑1 (MQ=255)
                                            cgcAGCGTAGTTCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGtt                                                                                                        <  2:808973/150‑1 (MQ=255)
                                              cAGCGTAGTTCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTta                                                                                                      <  2:171987/150‑1 (MQ=255)
                                              cAGCGTAGTTCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTta                                                                                                      <  2:1327046/150‑1 (MQ=255)
                                              cAGCGTAGTTCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTta                                                                                                      <  2:1951650/150‑1 (MQ=255)
                                              cAGCGTAGTTCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTta                                                                                                      >  1:1824655/1‑150 (MQ=255)
                                              cAGCGTAGTTCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTta                                                                                                      >  2:2478536/1‑150 (MQ=255)
                                                gCGTAGTTCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTtata                                                                                                    >  2:1886822/1‑150 (MQ=255)
                                                gCGTAGTTCCCCCTGCGACCAATAGGCCGGGGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTtata                                                                                                    <  2:860078/150‑1 (MQ=255)
                                                  gTAGTTCCCCCTGCGACCACTAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACGCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATc                                                                                                  <  1:1405522/150‑1 (MQ=255)
                                                   tAGTTCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCg                                                                                                 >  1:2147216/1‑150 (MQ=255)
                                                     gTTCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGc                                                                                               <  2:1222811/150‑1 (MQ=255)
                                                      ttCCCCCTGCGACCAATAGGCCGGTGTTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCa                                                                                              <  2:2205809/150‑1 (MQ=255)
                                                      ttCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATTACGCCGTGCGCCTTAACAGTTGCTCGGTGGCGGCGCAGTTTTATATCGGCa                                                                                              <  2:2166772/150‑1 (MQ=255)
                                                       tCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAg                                                                                             >  1:329000/1‑150 (MQ=255)
                                                         ccccTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCg                                                                                           <  2:2382543/150‑1 (MQ=255)
                                                         ccccTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCg                                                                                           >  1:1215985/1‑150 (MQ=255)
                                                            cTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAat                                                                                        >  2:1927801/1‑150 (MQ=255)
                                                              gCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAatat                                                                                      <  1:2478536/150‑1 (MQ=255)
                                                                    aaTAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGTATGACGCCGTGCGACTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACAt                                                                                >  1:1828323/1‑150 (MQ=255)
                                                                      tAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCa                                                                              <  1:721036/150‑1 (MQ=255)
                                                                         gCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTc                                                                           <  2:1989005/150‑1 (MQ=255)
                                                                          ccGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCg                                                                          >  1:1029479/1‑150 (MQ=255)
                                                                               ggtACTGCGGGCGTCAGGTGCGAACTCTATTCTTGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCaa                                                                     <  2:1194988/150‑1 (MQ=255)
                                                                                  aCTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCaaaaa                                                                  <  1:2467330/150‑1 (MQ=255)
                                                                                  aCTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCaaaaa                                                                  >  1:461807/1‑150 (MQ=255)
                                                                                      cGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAAtatt                                                              <  1:194946/150‑1 (MQ=255)
                                                                                        ggCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAAtattat                                                            >  1:829545/1‑150 (MQ=255)
                                                                                          cGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTc                                                          <  2:1377468/150‑1 (MQ=255)
                                                                                           gTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATAAGAACATCAGTCGATCAAAAATATTATTCa                                                         <  1:29434/150‑1 (MQ=255)
                                                                                            tCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAg                                                        >  2:980723/1‑150 (MQ=255)
                                                                                            tCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAg                                                        >  2:2332611/1‑150 (MQ=255)
                                                                                            tCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAg                                                        <  2:1064613/150‑1 (MQ=255)
                                                                                               ggTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATAGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTg                                                     >  2:308171/1‑150 (MQ=255)
                                                                                                 tGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGt                                                   <  2:1637882/150‑1 (MQ=255)
                                                                                                  gCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGtt                                                  <  2:2277934/150‑1 (MQ=255)
                                                                                                     aaCTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGTTGAt                                               >  2:1776620/1‑150 (MQ=255)
                                                                                                      aCTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGTTGATg                                              >  1:761525/1‑150 (MQ=255)
                                                                                                        tctATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGTTGATGcc                                            >  1:915830/1‑150 (MQ=255)
                                                                                                         ctATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGTTGATGCCg                                           >  2:976353/1‑150 (MQ=255)
                                                                                                         ctATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGTTGATGCCg                                           >  1:737098/1‑150 (MQ=255)
                                                                                                          tATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGTTGATGCCgg                                          <  1:1197293/150‑1 (MQ=255)
                                                                                                            ttCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACCTCAGTCGATCAAAAATATTATTCAGCTGGTTGATGCCGGaa                                        <  2:2147216/150‑1 (MQ=255)
                                                                                                            ttCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGTTGATGCCGGaa                                        <  2:1812430/150‑1 (MQ=255)
                                                                                                                ggCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGTTGATGCCGGAATGaa                                    >  1:1024428/1‑150 (MQ=255)
                                                                                                                 gCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGTTGATGCCGGAATGaaa                                   <  1:1055271/150‑1 (MQ=255)
                                                                                                                  cGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGTTGATGCCGGAATGAAAg                                  <  2:426718/150‑1 (MQ=255)
                                                                                                                   ggAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGTTGATGCCGGAATGAAAGt                                 >  2:1636641/1‑150 (MQ=255)
                                                                                                                    gAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAAAAGTTGCGCGGTGGCGGCGCAGCTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGTTGATGCCGGAATGAAAGTg                                >  1:561925/1‑150 (MQ=255)
                                                                                                                      aTTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGTTGATGCCGGAATGAAAGTggg                              >  2:1589159/1‑150 (MQ=255)
                                                                                                                         aaGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGTTGATGCCGGAATGAAAGTGGGAAt                           >  1:491549/1‑150 (MQ=255)
                                                                                                                         aaGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGTTGATGCCGGAATGAAAGTGGGAAt                           >  1:805532/1‑150 (MQ=255)
                                                                                                                           gTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGTTGATGCCGGAATGAAAGTGGGAATGc                         <  1:16658/150‑1 (MQ=255)
                                                                                                                           gTAATGAAGCCGGGGCGGTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGTCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGTTGATGCCGGAATGAAAGTGGGAATGc                         <  2:1560766/150‑1 (MQ=255)
                                                                                                                                   cgcacggCGTCATCCATCGATAACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGTTGATGCCGGAATGAAAGTGGGAATGCCGACCATg                 >  1:1949730/6‑150 (MQ=255)
                                                                                                                                     cGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGTTGATGCCGGAATGAAAGTGGGAATGCCGACCATGGc               <  2:1464897/150‑1 (MQ=255)
                                                                                                                                      gTGGCGTTATCGATGGATGACGCCCTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGATCATCAGTCGATCAAAAATATTATTAAGCTGGTTGATGCCGGAATGAAAGTGGGAATGCCGACCATGGcc              >  2:83376/1‑150 (MQ=255)
                                                                                                                                           gTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGTTGATGCCGGAATGAAAGTGGGAATGCCGACCATGGCCGTGAc         >  2:636361/1‑150 (MQ=255)
                                                                                                                                            ttATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGTTGATGCCGGAATGAAAGTGGGAATGCCGACCATGGCCGTGACt        >  1:739753/1‑150 (MQ=255)
                                                                                                                                             tATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGTTGATGCCGGAATGAAAGTGGGAATGCCGACCATGGCCGTGACTg       >  2:297828/1‑150 (MQ=255)
                                                                                                                                                 gatggatgACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGTTGATGCCGGAATGAAAGTGGGAATGCCGACCATGGCCGTGACTGGCGt   >  1:1906193/1‑150 (MQ=255)
                                                                                                                                                 gatggatgACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGTTGATGCCGGAATGAAAGTGGGAATGCCGACCATGGCCGTGACTGGCGt   >  2:464512/1‑150 (MQ=255)
                                                                                                                                                 gatggatgACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGTTGATGCCGGAATGAAAGTGGGAATGCCGACCATGGCCGTGACTGGCGt   >  2:473314/1‑150 (MQ=255)
                                                                                                                                                 gatggatgACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGTTGATGCCGGAATGAAAGTGGGAATGCCGACCATGGCCGTGACTGGCGt   >  1:1189037/1‑150 (MQ=255)
                                                                                                                                                  atggatgACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGTTGATGCCGGAATGAAAGTGGGAATGCCGACCATGGCCGTGACTGGCGTg  <  2:1502004/150‑1 (MQ=255)
                                                                                                                                                  atggatgACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGTTGATGCCGGAATGAAAGTGGGAATGCCGACCATGGCCGTGACTGGCGTg  <  1:2332611/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TGTTTGAACATGCCGATGTATTAATGTGTACGCGCGGCATTTTGCGCAGCGTAGTTCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGTTGATGCCGGAATGAAAGTGGGAATGCCGACCATGGCCGTGACTGGCGTG  >  NC_000913/1606332‑1606627

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