Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,167,281 T→C 5.4% intergenic (+81/+21) cyaR → / ← ogrK sRNA antisense regulator of ompX mRNA instability, Hfq‑dependent, cAMP‑induced/orphan Ogr protein, positive regulator of P2 growth

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,167,2810TC 5.4% 826.1 / 23.8 243intergenic (+81/+21)cyaR/ogrKsRNA antisense regulator of ompX mRNA instability, Hfq‑dependent, cAMP‑induced/orphan Ogr protein, positive regulator of P2 growth
Reads supporting (aligned to +/- strand):  ref base T (125/104);  new base C (5/8);  total (131/112)
Fisher's exact test for biased strand distribution p-value = 2.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

Only 100 of 246 total aligned reads displayed.

AATGCTAGCTGTACCAGGAACCACCTCCTTAGCCTGTGTAATCTCCCTTACACGGGCTTATTTTTTACGCGTAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTTGCTGCCCTGACGGCAACGGGTGCGGCCTTACGACGTGGACTTCTCCCGGCTTCACAATGTATCGCTGTACCGACTCATAAGTGATGAACGTGGCGCTGCAATTCACGTTCTG  >  NC_000913/2167135‑2167430
                                                                                                                                                  |                                                                                                                                                     
aaTGCTAGCTGTACCAGGAACCACCTCCTTAGCCTGTTTAATCTCCCTTACACGGGCTTATTTTTTACGCTTAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAg                                                                                                                                                    <  1:2159273/150‑1 (MQ=255)
   gCTAGCTGTACCAGGAACCACCTCCTTAGCCTGTGTAATCTCCCTTACACGGGCTTATTTTTTACGCGTAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAgggg                                                                                                                                                 <  1:2325114/150‑1 (MQ=255)
     tAGCTGTACCAGGAACCACCTCCTTAGCCTGTGTAATCTCCCTTACACGGGCTTATTTTTTACGCTTAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCAGTCGAGGGGCt                                                                                                                                               <  1:130424/150‑1 (MQ=255)
     tAGCTGTACCAGGAACCACCTCCTTAGCCTGTGTAATCTCCCTTACACGGGCTTATTTTTTACGCGTAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCt                                                                                                                                               >  2:2077873/1‑150 (MQ=255)
       gCTGTACCCGGAACCACCTCCTTAGCCTGTGTAATCTCCCTTACACGGGATTATTTTTTACGCGTAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCAGTCGAGGGGCttt                                                                                                                                             <  1:2684411/150‑2 (MQ=255)
       gCTGTACCAGGAACCCCCTCCTTAGCCTGTTTAATCTCCCTTACACGGGCTTATTTTTTCCGCGTAATACAATGAAATAAAAGGATTTATTTCTGGACACGCCCACACATTGACCACATCGACAAAAAAGCCCCTCAGTCGAGGGGCttt                                                                                                                                             <  1:1399109/150‑2 (MQ=255)
       gCTGTACCAGGAACCACCTCCTTAGCCTGTGTAATCTCCCTTACACGGGCTTATTTTTTACGCGTAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTc                                                                                                                                             >  2:3968054/1‑150 (MQ=255)
       gCTGTACCAGGAACCACCTCCTTAGCCTGTGTAATCTCCCTTACACGGGCTTATTTTTTACGCGTAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCAGTCGAGGGGCttt                                                                                                                                             <  2:2396425/150‑2 (MQ=255)
         tGTACCAGGAACCACCTCCTTCGCCTGTGTAATCTCCCTTACACGGGCTTATTTTTTACGCGTAATACAATTAAATAAAAGTATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCt                                                                                                                                           <  1:3786675/150‑1 (MQ=255)
         tGTACCAGGAACCACCTCCTTAGCCTGTGTAATCTCCCTTACACGGGCTTATTTTTTACGCGTAATACAATGAAAAAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCt                                                                                                                                           <  1:1218314/150‑1 (MQ=255)
           tACCAGGAACCACCTCCTTAGCCTGTGTAATCTCCCTTACACGGGCTTATTTTTTACGTGTAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGt                                                                                                                                         >  2:3532974/1‑150 (MQ=255)
             ccAGGACCCACCTCCTTAGCCTGTGTAATCTCCCTTACACGGGATTATTTTTTACGCGTAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGttt                                                                                                                                       <  1:1039889/150‑1 (MQ=255)
              cAGGAACCACCTCCTTAGCCTGTGTAATCTCCCTTACACGGGCTTATTTTTTACGCGTAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTc                                                                                                                                      >  1:1910098/1‑150 (MQ=255)
              cAGGAACCACCTCCTTAGCCTGTGTAATCTCCCTTACACGGGCTTATTTTTTACGCGTAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTc                                                                                                                                      >  1:3037418/1‑150 (MQ=255)
                 gAACCACCTCCTTAGCCTGTGTAATCTCCCTTACACGGGCTTATTTTTTACGCGGAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAACAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTaa                                                                                                                                   <  2:3448401/150‑1 (MQ=255)
                   accaccTCCTTAGCCTGTGTAATCTCCCTTACACGGGCTTATTTTTTACGCGTAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAAtt                                                                                                                                 <  1:1138073/150‑1 (MQ=255)
                     caccTCCTTAGCCTGTGTAATCTCCCTTACACGGGCTTATTTTTTACGCGTAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTAc                                                                                                                               <  2:83041/150‑1 (MQ=255)
                        ctcctTAGCCTGTGTAATCTCCCTTACACGGGCTTATTTTTTACGCGTAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATc                                                                                                                            >  1:1767122/1‑150 (MQ=255)
                          cctTAGCCTGTGTAATCTCCCTTACACGGGCTTATTTTTTACGCGTAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCTCATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCca                                                                                                                          <  2:3474961/150‑1 (MQ=255)
                            ttAGCCTGTGTAATCTCCCTTACACGGGCTTATTTTTTACGCGTAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCcaca                                                                                                                        >  1:2073461/1‑150 (MQ=255)
                            ttAGCCTGTGTAATCTCCCTTACACGGGCTTATTTTTTACGCGTAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCcaca                                                                                                                        <  1:2728083/150‑1 (MQ=255)
                                   gtgtAATCTCCCTTACACGGGCTTATTTATTACGCGTAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTtg                                                                                                                 <  1:1430843/150‑1 (MQ=255)
                                   ctgtAATCGCCCTTCCACGGGCTTATTTTTTACGCGTAAGACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACAGTGACCACATCGACAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAAGTACATCCACATAATTtgc                                                                                                                <  2:2382715/149‑1 (MQ=255)
                                    tgtAATCTCCCTTACACGGGCTTATTTTTTACGTGTAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTtgc                                                                                                                <  1:3532974/150‑1 (MQ=255)
                                    tgtAATCTCCCTTACACGGGCTTATTTTTTACGCGTAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAGTGGCTTCCTGTTTCTAATTACATCCACATAATTtgc                                                                                                                >  1:2512841/1‑150 (MQ=255)
                                    tgtAATCTCCCTTACACGGGCTTATTTTTTACGCGTAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTtgc                                                                                                                <  2:770388/150‑1 (MQ=255)
                                        aTCTCCCTTACACGGGCTTATTTTTTACGCGTAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTTGCTGcc                                                                                                            <  2:3952532/150‑1 (MQ=255)
                                           tcCCTTACACGGGCTTATTTTTTACGCGTAATACAATGAAATAAAAGGATTTATTTCTGGTCACGGCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATACTTTGCTGCCCTg                                                                                                         >  1:3883016/1‑150 (MQ=255)
                                              cTTACACGGGCTTATTTTTTACGCGTAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTTGCTGCCCTGACg                                                                                                      >  1:722578/1‑150 (MQ=255)
                                              cTTACACGGGCTTATTTTTTACGCGTAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTTGCTGCCCTGACg                                                                                                      <  2:2685874/150‑1 (MQ=255)
                                               ttACACGGGCTTATTTTTTACGCGTAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTTGCTGCCCTGACgg                                                                                                     >  2:1012331/1‑150 (MQ=255)
                                               ttACACGGGCTTATTTTTTACGCGTAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTTGCTGCCCTGACgg                                                                                                     >  2:3330792/1‑150 (MQ=255)
                                                 acacGGGCTTATTTTTTACGCGTAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTTGCTGCCCTGACGGCa                                                                                                   >  2:315089/1‑150 (MQ=255)
                                                   acGGGCTTATTTTTTACGCGTAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTTGCTGCCCTGACGGCAAc                                                                                                 >  2:810820/1‑150 (MQ=255)
                                                    cGGGCTTATTTTTTACGCGTAATACAATGAAATAACAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTTGCTGCCCTGACGGCAACg                                                                                                <  1:3601329/150‑1 (MQ=255)
                                                     gggCTTATTTTTTACGCGTAATCCAATGAAATAACAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTTGCTGCCCTGACGGCAACgg                                                                                               <  2:1837693/150‑1 (MQ=255)
                                                     gggCTTATTTTTTACGCGTAAGACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATACTTTGCTGCCCTGACGGCAACgg                                                                                               <  1:908279/150‑1 (MQ=255)
                                                       gCTTATTTTTTACGCGTAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTTGCTGCCCTGACGGCAACGGGt                                                                                             <  1:1937321/150‑1 (MQ=255)
                                                         ttATTTTTTACGCGTAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTTGCTGCCCTGACGGCACCGGGTGc                                                                                           >  1:3120411/1‑150 (MQ=255)
                                                         ttATTTTTTACGCGTAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTTGCTGCCCTGACGGCAACGGGTGc                                                                                           >  1:3324072/1‑150 (MQ=255)
                                                         ttATTTTTTACGCGTAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTTGCTGCCCTGACGGCAACGGGTGc                                                                                           >  1:347654/1‑150 (MQ=255)
                                                          tATTTTTTACGCGTAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTTGCTGCCCTGACGGCACCGGGTGCg                                                                                          >  2:4066091/1‑150 (MQ=255)
                                                           aTTTTTTACGCGTAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTTGCTGCCCTGACGGCAACGGGTGCgg                                                                                         >  2:272562/1‑150 (MQ=255)
                                                             tttttACGCGTAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTTGCTGCCCTGACGGCAACGGGTGCGGcc                                                                                       >  2:3993966/1‑150 (MQ=255)
                                                               tttACGCGTAATCCAATGAAATAAAATGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTTGCTGCCCTGACGGCAACGGGTGCGGCCtt                                                                                     <  2:858702/150‑1 (MQ=255)
                                                                ttACGCGTAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACCTTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTTGCTGCCCTGACGGCAACGGGTGCGGCCTTa                                                                                    >  2:3222801/1‑150 (MQ=255)
                                                                   cgcgTAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTTGCTGCCCTGACGGCAACGGGTGCGGCCTTacga                                                                                 <  1:280283/150‑1 (MQ=255)
                                                                      gTAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTTGCTGCCCTGACGGCAACGGGTGCGGCCTTACGACGt                                                                              <  1:1784689/150‑1 (MQ=255)
                                                                       tAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTTGCTGCCCTGACGGCAACGGGTGCGGCCTTACGACGTg                                                                             <  2:1796228/150‑1 (MQ=255)
                                                                       tAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTTGCTGCCCTGACGGCAACGGGTGCGGCCTTACGACGTg                                                                             >  2:1160594/1‑150 (MQ=255)
                                                                       tAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTTGCTGCCCTGACGGCAACGGGTGCGGCCTTACGACGTg                                                                             <  2:277170/150‑1 (MQ=255)
                                                                       tAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTTGCTGCCCTGACGGCAACGGGTGCGGCCTTACGACGTg                                                                             <  1:3022056/150‑1 (MQ=255)
                                                                       tAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTTGCTGCCCTGACGGCAACGGGTGCGGCCTTACGACGTg                                                                             <  1:1970733/150‑1 (MQ=255)
                                                                       tAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTTGCTGCCCTGACGGCAACGGGTGCGGCCTTACGACGTg                                                                             <  1:1769330/150‑1 (MQ=255)
                                                                        aaTACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTTGCTGCCCTGACGGCAACGGGTGCGGCCTTACGACGTgg                                                                            <  1:272552/150‑1 (MQ=255)
                                                                             aaTGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTTGCTGCCCTGACGGCAACGGGTGCGGCCTTACGACGTGGACTtc                                                                       >  2:1499519/1‑150 (MQ=255)
                                                                              aTGAAATAAAATGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTTGCTGCCCTGACGGCAACGGGTGCGGCCTTACGACGTGGACTtct                                                                      <  1:15474/150‑1 (MQ=255)
                                                                               tGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTTGCTGCCCTGACGGCAACGGGTGCGGCCTTACGACGTGGACTtctc                                                                     >  1:417024/1‑150 (MQ=255)
                                                                               tGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTTGCTGCCCTGACGGCAACGGGTGCGGCCTTACGACGTGGACTtctc                                                                     <  1:1584914/150‑1 (MQ=255)
                                                                                gAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTTGCTGCCCTGACGGCAACGGGTGCGGCCTTACGACGTGGACTTCTcc                                                                    >  2:2179878/1‑150 (MQ=255)
                                                                                 aaaTAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTTGCTGCCCTGACGGCAACGGGTGCGGCCTTACGACGTGGACTTCTccc                                                                   <  1:1110354/150‑1 (MQ=255)
                                                                                 aaaTAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTTGCTGCCCTGACGGCAACGGGTGCGGCCTTACGACGTGGACTTCTccc                                                                   >  1:3422250/1‑150 (MQ=255)
                                                                                   aTAAAAGGATTTATGTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTTGCTGCCCTGACGGCAACGGGTGCGGCCTTACGACGTGGACTTCTCCCgg                                                                 <  1:3646403/150‑1 (MQ=255)
                                                                                      aaaGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTTGCTGCCCTGACGGCAACGGGTGCGGCCTTACGACGTGGACTTCTCCCGGCtt                                                              <  2:3037371/150‑1 (MQ=255)
                                                                                         ggATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTTGCTGCCCTGACGGCAACGGGTGCGGCCTTACGACGTGGACTTCTCCCGGCTTcac                                                           >  1:1186873/1‑150 (MQ=255)
                                                                                          gATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTTGCTGCCCTGACGGCAACGGGTGCGGCCTTACGACGTGGACTTCTCCCGGCTTcaca                                                          <  2:469546/150‑1 (MQ=255)
                                                                                          gATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTTGCTGCCCTGACGGCAACGGGTGCGGCCTTACGACGTGGACTTCTCCCGGCTTcaca                                                          >  1:594199/1‑150 (MQ=255)
                                                                                                    tGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACTTCCACATAATTTGCTGCCCTGACGGCAACGGGTGCGGCCTTACGACATGGACTTCTCCCGGCTTCACAATGTATCGCt                                                >  1:3657260/1‑150 (MQ=255)
                                                                                                    tGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTTGCTGCCCTGACGGCAACGGGTGCGGCCTTACGACGTGGACTTCTCCCGGCTTCACAATGTATCGCt                                                <  2:1555043/150‑1 (MQ=255)
                                                                                                              cacacaTTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTTGCTGCCCTGACGGCAACGGGTGCGGCCTTACGACGTGGACTTCTCCCGGCTTCACAATGTATCGCTGTACCGACTc                                      >  2:89431/1‑150 (MQ=255)
                                                                                                              cacacaTTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTTGCTGCCCTGACGGCAACGGGTGCGGCCTTACGACGTGGACTTCTCCCGGCTTCACAATGTATCGCTGTACCGACTc                                      <  2:1767122/150‑1 (MQ=255)
                                                                                                               acacaTTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTTGCTGCCCTGACGGCAACGGGTGCGGCCTTACGACGTGGACTTCTCCCGGCTTCACAATGTATCGCTGTACCGACTCa                                     >  2:2234281/1‑150 (MQ=255)
                                                                                                                   aTTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTTGCTGCCCTGACGGCAACGGGTGCGGCCTTACGACGTGGACTTCTCCCGGCTTCACAATGTATCGCTGTACCGACTCATAAg                                 >  2:1699214/1‑150 (MQ=255)
                                                                                                                      gACCACATCGACAAAATAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTTGCTGCCCTGACGGCAACGGGTGCGGCCTTACGACGTGGACTTCTCCCGGCATCACAATGTATCGCTGTACCGACTCATAAGtga                              >  2:3410825/1‑150 (MQ=255)
                                                                                                                      gACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACCTAATTTGCTGCCCTGACGGCAACGGGTGCGGCCTTACGACGTGGACTTCTCCCGGCTTCACACTGTCTCGCGGTACCGACTCATAAGtga                              >  2:3501769/1‑150 (MQ=255)
                                                                                                                       aCCACATCGACAAAATAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTTGCTGCCCTGACGGCAACGGGTGCGGCCTTACGACGTGGACTTCTCCCGGCTTCACACTGTCTCGCTGTACCGACTCATAAGtgat                             >  2:4030963/1‑150 (MQ=255)
                                                                                                                         cacaTCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTTGCTGCCCTGACGGCAACGGGTGCGGCCTTACGACGTGGACTTCTCCCGGCTTCACAATGTATCGCTGTACCGACTCATAAGtgatga                           <  2:1381028/150‑1 (MQ=255)
                                                                                                                          acaTCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTTGCTGCCCTGACGGCAACGGGTGCGGCCTTACGACGTGGACTTCTCCCGGCTTCACAATGTATCGCTGTACCGACTCATAAGTGATGaa                          >  2:1253144/1‑150 (MQ=255)
                                                                                                                          acaTCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTTGCTGCCCTGACGGCAACGGGTGCGGCCTTACGACGTGGACTTCTCCCGGCTTCACAATGTATCGCTGTACCGACTCATAAGTGATGaa                          <  1:2008282/150‑1 (MQ=255)
                                                                                                                           caTCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTTGCTGCCCTGACGGCAACGGGTGCGGCCTTACGACGTGGACTTCTCCCGGCTTCACAATGTATCGCTGTACCGACTCATAAGTGATGAAc                         >  1:1347557/1‑150 (MQ=255)
                                                                                                                           caTCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTTGCTGCCCTGACGGCAACGGGTGCGGCCTTACGACGTGGACTTCTCCCGGCTTCACAATGTATCGCTGTACCGACTCATAAGTGATGAAc                         >  1:3025627/1‑150 (MQ=255)
                                                                                                                              cGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTTGCTGCCCTGACGGCAACGGGTGCGGCCTTACGACGTGGACTTCTCCCGGCTTCACAATGTATCGCTGTACCGACTCATAAGTGATGAACGTg                      <  1:3968054/150‑1 (MQ=255)
                                                                                                                               gACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTTGCTGCCCTGACGGCAACGGGTGCGGCCTTACGACGTGGACTTCTCCCGGCTTCACAATGTATCGCTGTACCGACTCATAAGTGATGAACGTgg                     >  2:2660611/1‑150 (MQ=255)
                                                                                                                               gACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTTGCTGCCCTGACGGCAACGGGTGCGGCCTTACGACGTGGACTTCTCCCGGCTTCACAATGTATCGCTGTACCGACTCATAAGTGATGAACGTgg                     <  2:2419839/150‑1 (MQ=255)
                                                                                                                                  aaaaaaGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATCATTTGCTGCCCTGACGGCAACGGGTGCGGCCTTACGACGTGGACTTCTCCCGGCTTCACAATGTATCGCTGTACCGACTCATAAGTGATGAACGTGgcgc                  <  1:2151969/150‑1 (MQ=255)
                                                                                                                                  aaaaaaGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTTGCTGCCCTGACGGCAACGGGTGCGGCCTTACGACGTGGACTTCTCCCGGCTTCACAATGTATCGCTGTACCGACTCATAAGTGATGAACGTGgcgc                  >  2:3704312/1‑150 (MQ=255)
                                                                                                                                  aaaaaaGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTTGCTGCCCTGACGGCAACGGGTGCGGCCTTACGACGTGGACTTCTCCCGGCTTCACAATGTATCGCTGTACCGACTCATAAGTGATGAACGTGgcgc                  >  2:1652267/1‑150 (MQ=255)
                                                                                                                                  aaaaaaGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTTGCTGCCCTGACGGCAACGGGTGCGGCCTTACGACGTGGACTTCTCCCGGCTTCACAATGTATCGCTGTACCGACTCATAAGTGATGAACGTGgcgc                  >  2:4007673/1‑150 (MQ=255)
                                                                                                                                    ggaaGCCCCTCAGTCGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTTGCTGCCCTGACGGCAACGGGTGCGGCCTTACGACGTGGACTTCTCCCGGCTTCACAATGTATCGCTGTACCGACTCATAAGTGATGAACGTGGCGCTg                >  2:422178/3‑150 (MQ=255)
                                                                                                                                     aaaGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTTGCTGCCCTGACGGCAACGGGTGCGGCCTTACGACGTGGACTTCTCCCGGCTTCACAATGTATCGCTGTACCGACTCATAAGTGATGAACGTGGCGCTGc               <  2:651262/150‑1 (MQ=255)
                                                                                                                                     aaaGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTTGCTGCCCTGACGGCAACGGGTGCGGCCTTACGACGTGGACTTCTCCCGGCTTCACAATGTATAGCTGTACCGACTCATAAGTGATGAACGTGGCGCTGc               >  1:3947507/1‑150 (MQ=255)
                                                                                                                                          cccTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTTGCTGCCCTGACGGCAACGGGTGCGGCCTTACGACGTGGACTTCTCCCGGCTTCACAATGTATCGCTGTACCGACTCATAAGTGATGAACGTGGCGCTGCAATTc          <  1:690343/150‑1 (MQ=255)
                                                                                                                                           ccTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTTGCTGCCCTGACGGCAACGGGTGCGGCCTTACGACGTGGACTTCTCCCGGCTTCACAATGTATCGCTGTACCGACTCATAAGTGATGAACGTGGCGCTGCAATTCa         >  1:2637844/1‑150 (MQ=255)
                                                                                                                                           ccTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTTGCTGCCCTGACGGCAACGGGTGCGGCCTTACGACGTGGACTTCTCCCGGCTTCACAATGTATCGCTGTACCGACTCATAAGTGATGAACGTGGCGCTGCAATTCa         <  2:3639060/150‑1 (MQ=255)
                                                                                                                                           ccTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTTGCTGCCCTGACGGCAACGGGTGCGGCCTTACGACGTGGACTTCTCCCGGCTTCACAATGTATCGCTGTACCGACTCATAAGTGATGAACGTGGCGCTGCAATTCa         >  1:1679498/1‑150 (MQ=255)
                                                                                                                                            cTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTTGCTGCCCTGACGGCAACGGGTGCGGCCTTACGACGTGGACTTCTCCCGGCTTCACAATGTATCGCTGTACCGACTCATAAGTGATGAACGTGGCGCTGCAATTCAc        >  2:483748/1‑150 (MQ=255)
                                                                                                                                            cTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTTGCTGCCCTGACGGCAACGGGTGCGGCCTTACGACGTGGACTTCTCCCGGCTTCACAATGTATCGCTGTACCGACTCATAAGTGATGAACGTGGCGCTGCAATTCAc        >  1:367329/1‑150 (MQ=255)
                                                                                                                                               gACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTTGCTGCCCTGACGGCAACGGGTGCGGCCTTACGACGTGGACTTCTCCCGGCTTCACAATGTATCGCTGTACCGACTCATAAGTGATGAACGTGGCGCTGCAATTCACGtt     >  2:2213185/1‑150 (MQ=255)
                                                                                                                                                 cTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTTGCTGCCCTGACGGCAACGGGTGCGGCCTTACGACGTGGACTTCTCCCGGCTTCACAATGTATCGCTGTACCGACTCATAAGTGATGAACGTGGCGCTGCAATTCACGTTCt   <  2:1985221/150‑1 (MQ=255)
                                                                                                                                                 cTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTTGCTGCCCTGACGGCAACGGGTGCGGCCTTACGACGTGGACTTCTCCCGGCTTCACAATGTATCGCTGTACCGACTCATAAGTGATGAACGTGGCGCTGCAATTCACGTTCt   <  1:3260626/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
AATGCTAGCTGTACCAGGAACCACCTCCTTAGCCTGTGTAATCTCCCTTACACGGGCTTATTTTTTACGCGTAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTTGCTGCCCTGACGGCAACGGGTGCGGCCTTACGACGTGGACTTCTCCCGGCTTCACAATGTATCGCTGTACCGACTCATAAGTGATGAACGTGGCGCTGCAATTCACGTTCTG  >  NC_000913/2167135‑2167430

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: