Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,367,147 T→G K256Q (AAA→CAA)  pspF ← psp operon transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,367,1470TG98.6% 258.4 / NA 69K256Q (AAA→CAA) pspFpsp operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (30/38);  total (31/38)
Fisher's exact test for biased strand distribution p-value = 4.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCTGGTTATATTTCCCTTGTTGCAAACTGAGTTGCAGCAACTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACC  >  NC_000913/1367005‑1367290
                                                                                                                                              |                                                                                                                                               
ttCTGGTTATATTTCCCTTGTTGCAAACTGAGTTGCAGCAACTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTGAAAGGGAt                                                                                                                                         >  1:380713/1‑151 (MQ=255)
    ggTTATATTTCCCTTGTTGCAAACTGAGTTGCAGCAACTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTGAAAGGGATc                                                                                                                                        <  2:165592/148‑1 (MQ=255)
    ggTTATATTTCCCTTGTTGCAAACTGAGTTGCAGCAACTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTGAAAGGGATCaat                                                                                                                                     <  1:745838/151‑1 (MQ=255)
     gTTATATTTCCCTTGTTGCAAACTGAGTTGCAGCAACTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAAAGGCGATAGCGTCTTCAGGCGGACGCGGTTGAAAGGGATCaat                                                                                                                                     >  1:178980/1‑150 (MQ=255)
        atatTTCCCTTGTTGCAAACTGAGTTGCAGCAACTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTGAAAGGGATCAATAatg                                                                                                                                 >  2:392349/1‑151 (MQ=255)
        atatTTCCCTTGTTGCAAACTGAGTTGCAGCAACTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTGAAAGGGATCAATAatg                                                                                                                                 <  2:985705/151‑1 (MQ=255)
         tatTTCCCTTGTTGCAAACTGAGTTGCAGCAACTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTGAAAGGGATCAATAatg                                                                                                                                 >  1:775477/1‑150 (MQ=255)
              cccTTGTTGCAAACTGAGTTGCAGCAACTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTGAAAGGGATCAATAATGATGtc                                                                                                                            <  2:993785/150‑1 (MQ=255)
               ccTTGTTGCAAACTGAGTTGCAGCAACTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTGAAAGGGATCAATAATGATGtcat                                                                                                                          <  1:219242/151‑1 (MQ=255)
               ccTTGTTGCAAACTGAGTTGCAGCAACTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTGAAAGGGATCAATAATGATGtcat                                                                                                                          <  2:1103609/151‑1 (MQ=255)
                       cAAACTGAGTTGCAGCAACTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGCGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTGACAGGGATCAATAATGATGTCATCAAGCgg                                                                                                                   >  2:413302/1‑150 (MQ=255)
                           cTGAGTTGCAGCAACTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTGAAAGGGATCAATAATGATGTCATCAAGCGGATAAt                                                                                                              >  2:116708/1‑151 (MQ=255)
                              aGTTGCAGCAACTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTGAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGc                                                                                                           <  1:93823/151‑1 (MQ=255)
                                   cagcaACTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTGAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTg                                                                                                      <  2:523343/151‑1 (MQ=255)
                                    agcaACTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCATTGAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGc                                                                                                     >  1:24650/1‑151 (MQ=255)
                                      caACTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTGAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCg                                                                                                   <  2:405137/151‑1 (MQ=255)
                                       aaCTCTTTTTCCTGCTGCATCAGAAACTCACGTAAATCCAGCGGCAGTGTTGGAATCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTCAAGGGATCAATAATGAAGTCATCAAGCGGATAATCGCTGGGGCCGt                                                                                                  >  2:938861/1‑151 (MQ=255)
                                         ctctTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTGAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTgg                                                                                                >  1:551712/1‑151 (MQ=255)
                                          tctTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTGAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGc                                                                                               >  2:842087/1‑151 (MQ=255)
                                              tttCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTGAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGAt                                                                                            >  2:315980/1‑150 (MQ=255)
                                                 ccTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTGAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATaca                                                                                         <  1:1107309/150‑1 (MQ=255)
                                                 ccTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTGAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATaca                                                                                         >  1:96132/1‑150 (MQ=255)
                                                    gctgcATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTGAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACt                                                                                       <  2:784985/149‑1 (MQ=255)
                                                    gctgcATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTGAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTg                                                                                      >  2:608857/1‑150 (MQ=255)
                                                      tgcATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTGAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGa                                                                                     >  1:1151942/1‑149 (MQ=255)
                                                      tgcATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTGAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAAc                                                                                   >  2:225949/1‑151 (MQ=255)
                                                      tgcATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTGAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAAc                                                                                   >  1:937258/1‑151 (MQ=255)
                                                       gcATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTGAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGaa                                                                                    <  2:1199323/149‑1 (MQ=255)
                                                         aTCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTGAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGtt                                                                                >  1:1133141/1‑151 (MQ=255)
                                                          tCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTGAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGt                                                                                 <  2:525419/149‑1 (MQ=255)
                                                             gAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTGAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTccac                                                                            >  1:122843/1‑151 (MQ=255)
                                                             gAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTGAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTcca                                                                             <  2:692497/150‑1 (MQ=255)
                                                                   cACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTGAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGtt                                                                      >  2:434789/1‑151 (MQ=255)
                                                                   cACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTGAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGtt                                                                      >  2:819405/1‑151 (MQ=255)
                                                                   cACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTGAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGtt                                                                      <  1:446655/151‑1 (MQ=255)
                                                                     cGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTGAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGGGGCGATACACTGAACGCTCCACCACGtttt                                                                    >  2:811202/1‑151 (MQ=255)
                                                                      gTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTGAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGttttt                                                                   <  1:1204298/151‑1 (MQ=255)
                                                                      gTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTGAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGttttt                                                                   <  2:970520/151‑1 (MQ=255)
                                                                        aaaTCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTGAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCa                                                                 <  2:973622/151‑1 (MQ=255)
                                                                         aaTCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTGAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCaa                                                                <  1:472662/151‑1 (MQ=255)
                                                                           tCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTGAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAAtt                                                              >  1:1196412/1‑151 (MQ=255)
                                                                            ccAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTGAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTc                                                             >  1:865331/1‑151 (MQ=255)
                                                                             cAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGAGGCCGTTGAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGAGGCGATACACGGAACGTTCCACCACGTTTTTCAAtt                                                              >  2:12190/1‑149 (MQ=255)
                                                                                   cAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTGAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAata                                                      <  1:1128333/151‑1 (MQ=255)
                                                                                   cAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTGAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAat                                                       >  2:978159/1‑150 (MQ=255)
                                                                                    aGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTGAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAatat                                                     <  2:750515/151‑1 (MQ=255)
                                                                                       gtTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTGAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATAttt                                                   >  2:624547/1‑150 (MQ=255)
                                                                                         tGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTGAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTcc                                                 <  2:566018/150‑1 (MQ=255)
                                                                                             aGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTGAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGcc                                            >  2:820417/1‑151 (MQ=255)
                                                                                              gCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTGAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGcc                                            <  2:1012434/150‑1 (MQ=255)
                                                                                                            gAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTGAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCa                                <  2:1007056/148‑1 (MQ=255)
                                                                                                            gAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTGAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGc                              <  1:281131/150‑1 (MQ=255)
                                                                                                            gAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTGAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGc                              <  2:165326/150‑1 (MQ=255)
                                                                                                             aaaCGGCGATAGCGTCTTCAGGCGGACGCCGTTGAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCaa                            <  1:308091/151‑1 (MQ=255)
                                                                                                             aaaCGGCGATAGCGTCTTCAGGCGGACGCCGTTGAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCaa                            <  2:191851/151‑1 (MQ=255)
                                                                                                                cGGCGATAGCGTCTTCAGGCGGACGCCGTTGAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGt                         >  1:1192401/1‑151 (MQ=255)
                                                                                                                cGGCGATAGCGTCTTCAGGCGGACGCCGTTGAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGt                         <  1:290388/151‑1 (MQ=255)
                                                                                                                   cGATAGCGTCTTCAGGCGGACGCCGTTGAAAGGGATCAATAATGATGACATCAAGCGGATAATCGTTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTGTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGttt                       <  2:584876/150‑1 (MQ=255)
                                                                                                                     aTAGCGTCTTCAGGCGGACGCCGTTGAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTtctc                    <  1:492838/151‑1 (MQ=255)
                                                                                                                       aGCGTCTTCAGGCGGACGCCGTTGAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTctct                   <  1:661713/150‑1 (MQ=255)
                                                                                                                        gCGTCTTCAGGCGGACGCCGTTGAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTg                  <  2:178972/150‑1 (MQ=255)
                                                                                                                          gTCTTCAGGCGGACGCCGTTGAAATGGATCAATACTGATGTCATCTAGCGGATAATCGCTGGTGCCATGGCGATACGCTGAACATTCCACCACGTTTTTCAATTCCCGAATATTTCCCGTCCAACGATAATTCAGCAATGTTTctct                   <  2:779786/147‑1 (MQ=255)
                                                                                                                          gTCTTCAGGCGGACGCCGTTGAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGgcg               <  2:562564/151‑1 (MQ=255)
                                                                                                                               cAGGCGGACGCCGTTGAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGTGCt           <  2:934767/150‑1 (MQ=255)
                                                                                                                                  gCGGACGCCGTTGAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGt       <  1:791100/151‑1 (MQ=255)
                                                                                                                                    ggACGCCGTTGAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTaa     <  2:1057805/151‑1 (MQ=255)
                                                                                                                                     gACGCCGTTGAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTa      >  1:413021/1‑149 (MQ=255)
                                                                                                                                       cGCCGTTGAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAAcc  >  2:506369/1‑151 (MQ=255)
                                                                                                                                       cGCCGTTGAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGAGCTCCGTAAAcc  >  2:512597/1‑151 (MQ=255)
                                                                                                                                              |                                                                                                                                               
TTCTGGTTATATTTCCCTTGTTGCAAACTGAGTTGCAGCAACTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACC  >  NC_000913/1367005‑1367290

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