Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 518,860 A→G 5.0% F97S (TTC→TCC)  ybbO ← putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913518,8600AG 5.0% 303.2 / 4.8 99F97S (TTC→TCC) ybbOputative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (35/59);  new base G (3/2);  total (38/61)
Fisher's exact test for biased strand distribution p-value = 3.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.94e-01

GATGATGTCATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAACATCATCCGGTTTC  >  NC_000913/518713‑519009
                                                                                                                                                   |                                                                                                                                                     
gatgatGTCATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGAGGAAAGGGGGCCATACATGCCGAAt                                                                                                                                                     >  1:1266079/1‑150 (MQ=255)
    atgTCATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCgg                                                                                                                                                 <  1:431824/150‑1 (MQ=255)
     ggTCATCACAATACGCCCTTCACCGTGCGGTAACAGCGCGGGTAACAGCCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGGGGAAAGGGGGCCATACATGCCGAATCCGGc                                                                                                                                                <  2:91716/149‑1 (MQ=255)
       tcatcaCAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCAt                                                                                                                                              <  2:1550994/150‑1 (MQ=255)
       tcatcaCAATACGCACTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTCTGCGCCGAAAAAGTTGGCGGAAAACAGCTGTTCCATCTGCGCACGCCTGATGGTGGAAAGCGGGCCATACATGCCGAATCCGGAAt                                                                                                                                              >  2:48645/1‑150 (MQ=255)
        catcaCAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGGATTCGGCAtg                                                                                                                                             <  1:151342/150‑2 (MQ=255)
         accacaATCCGCCCTTCACCGAGCGGAAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTTCTCTTCCAGCTGCGCACGTCTGATGGAGGAAAGGGGGCCATACATGCCGAATCCGGCATTg                                                                                                                                            <  1:263278/148‑1 (MQ=255)
         aacacaATACGCCCTTCACCGTGCGGTAACAGCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTg                                                                                                                                            <  2:1589446/148‑1 (MQ=255)
          tcaCAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGt                                                                                                                                           >  2:928890/1‑150 (MQ=255)
           cacaATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGGATTCGGCAtgtat                                                                                                                                          <  2:1522216/150‑5 (MQ=255)
           cacaATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTAGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGtt                                                                                                                                          >  1:642854/1‑150 (MQ=255)
             caATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTaa                                                                                                                                        <  2:1207807/150‑1 (MQ=255)
             caATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGATGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTaa                                                                                                                                        >  1:142275/1‑150 (MQ=255)
               aTACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAg                                                                                                                                      <  2:733054/150‑1 (MQ=255)
                 aCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGAt                                                                                                                                    <  1:1528838/150‑1 (MQ=255)
                  cGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATc                                                                                                                                   >  2:1481932/1‑150 (MQ=255)
                    cccTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCAGGCATTGTTAAAGATccc                                                                                                                                 >  2:71412/1‑150 (MQ=255)
                     ccTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCGGTTCCATCTGCGCACGGCGGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCa                                                                                                                                <  1:101288/150‑1 (MQ=255)
                      cTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATcccta                                                                                                                               >  2:278188/1‑148 (MQ=255)
                      cTTCACCGTGAGGTAACATCGCGGGTAACAGGCTCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCAt                                                                                                                               <  1:843481/150‑1 (MQ=255)
                         cACCGTGCGGAAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACa                                                                                                                            <  1:30904/150‑1 (MQ=255)
                            cGTGCGGTAACATCGCTGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGAc                                                                                                                         <  2:716835/150‑1 (MQ=255)
                            cGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGAc                                                                                                                         <  2:656/150‑1 (MQ=255)
                            cGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGAGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGATCCATCTGCTCACGGCTGATTGTTGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGAc                                                                                                                         <  1:278188/150‑1 (MQ=255)
                             gTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACa                                                                                                                        >  2:446097/1‑150 (MQ=255)
                              ggCGGTAACATCGCGGGTAACAGGCGCATGGTGTGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACaa                                                                                                                       <  2:1778290/149‑1 (MQ=255)
                                  gTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTAt                                                                                                                   >  1:971332/1‑150 (MQ=255)
                                       aTCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAACACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTc                                                                                                              <  2:1099293/150‑1 (MQ=255)
                                       aTCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAATCTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTc                                                                                                              >  1:225020/1‑150 (MQ=255)
                                            gggTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGc                                                                                                         <  1:944394/150‑1 (MQ=255)
                                             ggTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGATTCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCg                                                                                                        >  2:1524244/1‑150 (MQ=255)
                                               tAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGAt                                                                                                      <  2:1742/150‑1 (MQ=255)
                                                aaCAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATc                                                                                                     <  2:1546717/150‑1 (MQ=255)
                                                aaCAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATc                                                                                                     <  1:370754/150‑1 (MQ=255)
                                                  cAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCAc                                                                                                   >  2:1443459/1‑150 (MQ=255)
                                                   agtcgcATGGTGAGCTGGTGTGCTCCTAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCAcc                                                                                                  <  2:1266079/147‑1 (MQ=255)
                                                   aGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCAcc                                                                                                  <  2:997350/150‑1 (MQ=255)
                                                    ggCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCt                                                                                                 >  1:849013/1‑150 (MQ=255)
                                                      cgcATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGTCCATACATTCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCATGGCGATCACCtcg                                                                                               >  2:1716344/1‑150 (MQ=255)
                                                       gcATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATTCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCtcgt                                                                                              >  2:1585529/1‑150 (MQ=255)
                                                         ttGGTGAGCTGGTGTGCGCCGAAATAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACTGCTGATGGTGGAAAGTGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCtcgtcg                                                                                            <  1:1141332/149‑1 (MQ=255)
                                                           ggTGAGCGGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGc                                                                                          <  1:1662487/150‑1 (MQ=255)
                                                            gTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCt                                                                                         >  1:1351705/1‑150 (MQ=255)
                                                             tGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTg                                                                                        <  2:347810/150‑1 (MQ=255)
                                                              gAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTgc                                                                                       >  2:217549/1‑150 (MQ=255)
                                                                gCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGTGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTgcgc                                                                                     <  1:533473/150‑1 (MQ=255)
                                                                 cTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGcgcg                                                                                    <  2:541079/150‑1 (MQ=255)
                                                                  tGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGa                                                                                   <  1:1274651/150‑1 (MQ=255)
                                                                   ggTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGAt                                                                                  <  1:1786986/150‑1 (MQ=255)
                                                                      gtgCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGTGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCaa                                                                               <  1:1047338/150‑1 (MQ=255)
                                                                       tgCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAac                                                                              >  2:1737909/1‑150 (MQ=255)
                                                                           ccGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGAGGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGCTCACCTCGTCGGCTGCGCGATCAACACtt                                                                          >  2:295603/1‑150 (MQ=255)
                                                                            cGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCTGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGTTCAGGGCGATAACCTAGTCGGCTGCGAGATCAACACttt                                                                         >  1:357057/1‑150 (MQ=255)
                                                                            cGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGTTCAGGTCGATCACCTCGTCGGCTGCGCGATCAACACttt                                                                         >  2:1390181/1‑150 (MQ=255)
                                                                                aaaGGTGGCCGAAAACTGCTGATCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGTTCATGGCGATCACCTGGTCGGCTGCGCGATCAACACTTTCTgg                                                                     <  2:135865/150‑1 (MQ=255)
                                                                                 aaGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGt                                                                    <  1:1315403/150‑1 (MQ=255)
                                                                                    ttGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGaa                                                                 >  2:1431786/1‑150 (MQ=255)
                                                                                       gCGGAATACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGATTCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATcc                                                              <  1:393148/150‑1 (MQ=255)
                                                                                        cGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGTGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCa                                                             <  1:216635/150‑1 (MQ=255)
                                                                                        cGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCa                                                             >  2:704095/1‑150 (MQ=255)
                                                                                          gAAAACTGCTGATCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAga                                                           <  1:1590399/150‑1 (MQ=255)
                                                                                            aaaCTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAgatc                                                         <  1:928890/150‑1 (MQ=255)
                                                                                            aaaCTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAgatc                                                         <  1:832326/150‑1 (MQ=255)
                                                                                            aaaCTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAgatc                                                         <  1:1271817/150‑1 (MQ=255)
                                                                                             aaCTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAgatcg                                                        >  1:1500725/1‑150 (MQ=255)
                                                                                             aaCTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAgatcg                                                        >  1:1507758/1‑150 (MQ=255)
                                                                                              aCTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGTTTATTTTTATGATATGGTGTGCTGCGAGATCTACACTTTCTGGTGAATCCAgatcga                                                       <  2:1046236/150‑1 (MQ=255)
                                                                                               ctgctgTT‑CATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGGGAATCCAgatcgatc                                                     <  1:857085/150‑1 (MQ=255)
                                                                                                tgctgTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAgatcgatc                                                     <  1:826478/150‑1 (MQ=255)
                                                                                                tgctgTTCCAACTGCGCACGGCTGATGGTTGAAATGTGGCCATACATGCCGAATCCGGCATTGTTATAGATCCCATACTGTCAATTTTTTTTTTTTTTTATGATATGGTGTGCTGCGAGATCTACACTTTCTGGTGAATCCAgatcgatc                                                     <  2:830805/150‑1 (MQ=255)
                                                                                                tgctgTTACATCTGATCACGGATTATGTTGGAAAGTGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTCTCTGTCAGGGCGAGCACCTCGGCGGCTGCGCGATCAACACTTTCTGGTGAATCCAgatcgatc                                                     <  2:1190197/150‑1 (MQ=255)
                                                                                                   tgTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAac                                                  <  1:1481932/150‑1 (MQ=255)
                                                                                                      tCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACg                                               >  1:1825946/1‑150 (MQ=255)
                                                                                                       ccATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGc                                              >  1:585910/1‑150 (MQ=255)
                                                                                                        cATCTGCGCACGGCGGATGTGGGAAAGGGGGCCCTACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCTGTTTTTTTTATGAGATGGTCTGCTGCGAGATCTACACTTTCTGGTGAATCCAGATCGATCAACACGcc                                             <  2:1333546/150‑1 (MQ=255)
                                                                                                         aTCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCg                                            >  1:1252036/1‑150 (MQ=255)
                                                                                                            tGCCCACGGCGTATGGTGGAACGGGGGCCAAACAGGCCGAATCCGGCATTGTTAAAGAACCCATACAGACAATTTTCGGTATGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTCAATCCAGATCGATCAACACGCCGGTa                                         <  1:1645962/150‑1 (MQ=255)
                                                                                                                 aCGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATct                                    <  1:1551273/150‑2 (MQ=255)
                                                                                                                    gCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGTCATTGTTAAAGATCCCATACAGACAATTATCGGTCATTGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCAt                                 <  1:71412/150‑1 (MQ=255)
                                                                                                                       gATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACCCGCCGGTAAATCCCATGCt                              >  2:339960/1‑150 (MQ=255)
                                                                                                                        aTGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTg                             <  1:451654/150‑1 (MQ=255)
                                                                                                                          ggtggAAAGGGGGCCATACATGCCTATTCCGGCATTGTTAAAGATCCCATAAAGAAAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCTATCAACACTTTCTGGTGATTCCAGATCGATCAACACGCCGGTAAATCCCATGCTGtt                           <  1:1703630/150‑1 (MQ=255)
                                                                                                                            tggAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGATGCGCGATAAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCa                         >  1:1163082/1‑150 (MQ=255)
                                                                                                                              gAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCTGTTATATCACATCGTGTGCTGCGAGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATg                       <  2:225020/150‑1 (MQ=255)
                                                                                                                               aaaGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATgc                      <  1:534712/150‑1 (MQ=255)
                                                                                                                               aaaGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATgc                      >  1:997218/1‑150 (MQ=255)
                                                                                                                                 aGGGCGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCTCGATCAACACTTTCTGGTGAATCCAGTTCGATCAACACGCCGGTAAATCCCATGCTGTTCATgcgc                    <  2:339796/150‑1 (MQ=255)
                                                                                                                                  gggggCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCt                   >  2:1083330/1‑150 (MQ=255)
                                                                                                                                    gggCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCa                 <  1:217549/150‑1 (MQ=255)
                                                                                                                                     ggCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCaa                >  2:864248/1‑150 (MQ=255)
                                                                                                                                     ggCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCaa                <  1:1737909/150‑1 (MQ=255)
                                                                                                                                     ggCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCaa                <  1:534477/150‑1 (MQ=255)
                                                                                                                                      tttaacaaTGCCGGATTCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGAGATCAACACTTTCTTGTGAATCCAGATCGATCACCCCGCCGTTAAATCCCATGCTGTTACTTCGCTCCAc               >  2:482609/8‑150 (MQ=255)
                                                                                                                                        cATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAAcat             <  1:446097/150‑1 (MQ=255)
                                                                                                                                         aacaaTGCCGGATTCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTTCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAAcatc            >  2:1226899/5‑150 (MQ=255)
                                                                                                                                         aacaaTGCCGGATTCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACAATTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAAcatc            >  1:323836/5‑150 (MQ=255)
                                                                                                                                          tACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAAcatca           <  1:733881/150‑1 (MQ=255)
                                                                                                                                               gCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATTCTGTTCATGCGCTCAACATCATCCgg      >  2:467812/1‑150 (MQ=255)
                                                                                                                                                  gAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCATGATCAACACTTTCTGGTGAATCCAGATCGAGCAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAACATCATCCGGttt   <  2:1268866/150‑1 (MQ=255)
                                                                                                                                                   aaTCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAACATCATCCGGTTTc  >  1:964640/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GATGATGTCATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAACATCATCCGGTTTC  >  NC_000913/518713‑519009

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: ‑