Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,900,327 A→C 5.1% intergenic (‑54/‑265) ygcW ← / → yqcE putative SDR family oxidoreductase/putative MFS transporter, inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,900,3270AC 5.1% 277.4 / 3.9 80intergenic (‑54/‑265)ygcW/yqcEputative SDR family oxidoreductase/putative MFS transporter, inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (34/41);  new base C (2/3);  total (36/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.55e-01

CCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTACTCTCACACAAC  >  NC_000913/2900182‑2900474
                                                                                                                                                 |                                                                                                                                                   
ccTAACCCGCTATTCCCACCGGTAAAAATTGCGGTTTTACCTTACAGCGAGAAAAACTCCATTGAAAACGCTTTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAAGTAACTTCCCCCGCGagag                                                                                                                                                  <  4:874501/149‑1 (MQ=255)
ccTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGagag                                                                                                                                                  <  4:567234/149‑1 (MQ=255)
ccTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGagag                                                                                                                                                  <  3:382194/149‑1 (MQ=255)
 cTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAgg                                                                                                                                                 <  3:343397/149‑1 (MQ=255)
 cTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAgg                                                                                                                                                 >  3:198943/1‑149 (MQ=255)
 cTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAgg                                                                                                                                                 >  2:244512/1‑149 (MQ=255)
  tAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAggg                                                                                                                                                >  3:425023/1‑149 (MQ=255)
    aaCCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGaa                                                                                                                                              <  4:67873/149‑1 (MQ=255)
    aaCCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGaa                                                                                                                                              <  3:486689/149‑1 (MQ=255)
    aaCCGCTATTCCCACCGGAAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGaa                                                                                                                                              <  2:544323/149‑1 (MQ=255)
     aCCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCCTTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCACTTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAg                                                                                                                                             <  4:79823/149‑1 (MQ=255)
      ccGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGt                                                                                                                                            <  3:134816/149‑1 (MQ=255)
      ccGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGt                                                                                                                                            <  1:802207/149‑1 (MQ=255)
       cGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCTCTCGCGGGGGAAGtt                                                                                                                                           <  4:459207/149‑1 (MQ=255)
        gCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTa                                                                                                                                          >  4:78841/1‑149 (MQ=255)
          tATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCTCTCGCGGGGGAAGTTAcc                                                                                                                                        <  3:285656/149‑1 (MQ=255)
               ccACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTtct                                                                                                                                   <  1:105186/149‑1 (MQ=255)
               ccACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTtct                                                                                                                                   >  2:837167/1‑149 (MQ=255)
               ccACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTtct                                                                                                                                   >  4:63663/1‑149 (MQ=255)
               ccACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTtct                                                                                                                                   <  4:536579/149‑1 (MQ=255)
                cACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAAGTTCCCCCGCGAGAGGGAAGTTACCTTtctc                                                                                                                                  >  3:718044/1‑149 (MQ=255)
                cACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTtctc                                                                                                                                  >  1:22577/1‑149 (MQ=255)
                  ccGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCAt                                                                                                                                <  2:22577/149‑1 (MQ=255)
                     gTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTAt                                                                                                                             >  3:242492/1‑149 (MQ=255)
                        aCAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTAtctc                                                                                                                          <  4:242492/149‑1 (MQ=255)
                          aaTTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTg                                                                                                                        <  4:418542/149‑1 (MQ=255)
                           aTTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTgg                                                                                                                       <  1:217908/149‑1 (MQ=255)
                                        cTTTAAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGGTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATaa                                                                                                          <  3:842323/149‑1 (MQ=255)
                                         tttCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGCAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCCAGTTACCAGAGATAAt                                                                                                         >  2:434886/1‑149 (MQ=255)
                                           tCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGCAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCCAGTTACCAGAGATAAtga                                                                                                       <  1:434886/149‑3 (MQ=255)
                                            cAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAAtttt                                                                                                      <  3:841581/149‑1 (MQ=255)
                                            cAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAAtttt                                                                                                      <  3:781520/149‑1 (MQ=255)
                                              gggAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTa                                                                                                    >  4:121470/1‑149 (MQ=255)
                                              gggAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTa                                                                                                    <  4:425023/149‑1 (MQ=255)
                                              gggAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTa                                                                                                    <  2:154105/149‑1 (MQ=255)
                                               ggAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAg                                                                                                   >  1:415713/1‑149 (MQ=255)
                                                gagaAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGc                                                                                                  >  2:816874/1‑149 (MQ=255)
                                                   aaaaaaTCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATa                                                                                               >  1:384635/1‑149 (MQ=255)
                                                     aaaaTCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTGACTTACCCCGCGAGAGGGAAGTTACCTTTCTCAtt                                                                                                                               >  3:834265/1‑115 (MQ=255)
                                                     aaaaTCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCAtt                                                                                                                               <  4:834265/115‑1 (MQ=255)
                                                      aaaTCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAgg                                                                                            <  2:415713/149‑1 (MQ=255)
                                                        aTCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCGGTAAGGCAAAAAAAAAGAAAGGAAACTTCCCCCGGGAGAGGGAAGTTCCCTTTCTCATTATCTCTGGAAACTGGAGATAATTTTTAGCATAAGGGc                                                                                          >  1:176263/1‑149 (MQ=255)
                                                        aTCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGc                                                                                          >  2:857592/1‑149 (MQ=255)
                                                             ttGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCTTAAGGGCATAAc                                                                                     <  1:59275/149‑1 (MQ=255)
                                                               gagaACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTg                                                                                   >  3:755854/1‑149 (MQ=255)
                                                                agaACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGc                                                                                  <  1:584581/149‑1 (MQ=255)
                                                                agaACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGc                                                                                  <  1:336501/149‑1 (MQ=255)
                                                                   aCGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTaa                                                                               >  1:370998/1‑149 (MQ=255)
                                                                   aCGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTaa                                                                               <  2:248679/149‑1 (MQ=255)
                                                                    cGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAAt                                                                              <  1:244512/149‑1 (MQ=255)
                                                                      cATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAAtgt                                                                            <  2:525559/149‑1 (MQ=255)
                                                                        ttGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGa                                                                          <  4:856166/149‑1 (MQ=255)
                                                                          gagagaTTTGATTGACCTAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATc                                                                        <  2:176263/149‑1 (MQ=255)
                                                                          gagagaTTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTCTCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATc                                                                        <  1:577297/149‑1 (MQ=255)
                                                                              gaTTCGATTGACATAGTTAGCTCCAATTTCGGTGAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTATCTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTaa                                                                    <  1:579495/149‑1 (MQ=255)
                                                                              gaTTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTaa                                                                    >  1:4339/1‑149 (MQ=255)
                                                                               aTTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAAt                                                                   <  4:596019/149‑1 (MQ=255)
                                                                               aTTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAAt                                                                   >  2:427887/1‑149 (MQ=255)
                                                                                ttCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATc                                                                  >  4:139164/1‑149 (MQ=255)
                                                                                   gATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATcaca                                                               <  4:755854/149‑1 (MQ=255)
                                                                                              ttAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTc                                                    >  3:832554/1‑149 (MQ=255)
                                                                                              ttAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTc                                                    <  1:223787/149‑1 (MQ=255)
                                                                                                   tCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTAtt                                               >  4:225769/1‑149 (MQ=255)
                                                                                                    ccAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCCTAACTGCTAATGTGATCGTAATCACAg                                                              <  2:390624/133‑1 (MQ=255)
                                                                                                    ccAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAg                                                              >  1:390624/1‑133 (MQ=255)
                                                                                                    ccAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTg                                              >  2:59569/1‑149 (MQ=255)
                                                                                                     cAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGt                                             <  4:109049/149‑1 (MQ=255)
                                                                                                                  gCAAAAAAAAAGAAAGGTAACTTACCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTaa                                <  2:177523/149‑1 (MQ=255)
                                                                                                                       aaaaaaGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGc                           <  2:384635/149‑1 (MQ=255)
                                                                                                                        aaaaaGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCa                          <  3:331553/149‑1 (MQ=255)
                                                                                                                                  gTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGAt                >  2:642424/1‑149 (MQ=255)
                                                                                                                                   tAACTTCCCTCTCGCGGGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGAtt               >  4:873265/1‑149 (MQ=255)
                                                                                                                                   tAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGAtt               >  3:420162/1‑149 (MQ=255)
                                                                                                                                   tAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGAtt               >  3:104912/1‑149 (MQ=255)
                                                                                                                                    aaCTTCCCTCTCGCGGGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAg                            <  4:279095/135‑1 (MQ=255)
                                                                                                                                    aaCTTCCCTCTCGCGGGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAg                            >  3:279095/1‑135 (MQ=255)
                                                                                                                                      cTTCCCCCGCGAGAGGGAAGTTACCTTTCTCAGTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAATCGTGATTAct            >  4:564141/1‑149 (MQ=255)
                                                                                                                                       ttCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTActc           >  1:184257/1‑149 (MQ=255)
                                                                                                                                       ttCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGGATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTActc           >  4:865385/1‑149 (MQ=255)
                                                                                                                                        tCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTActct          >  2:499412/1‑149 (MQ=255)
                                                                                                                                         cccccGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGATGCGTGATTACtctc         <  3:63663/149‑1 (MQ=255)
                                                                                                                                            ccGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTACTCTcaca      >  1:96641/1‑149 (MQ=255)
                                                                                                                                            ccGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTACTCTcaca      <  2:4339/149‑1 (MQ=255)
                                                                                                                                                gagaGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTCTCAAATAATTTTAAAAAGCAGAAGCGTGATTACTCTCACacaac  >  1:426000/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTACTCTCACACAAC  >  NC_000913/2900182‑2900474

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: