Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,273,931 T→A 5.3% I523N (ATC→AAC)  yejA → microcin C ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,273,9310TA 5.3% 242.9 / 3.2 76I523N (ATC→AAC) yejAmicrocin C ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base T (37/35);  new base A (2/2);  total (39/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.10e-01

TCGCGACTATGACATGATGCCGCGCGTATGGCGGGCGATGCCGTGGCCCAGTTCCGATTTACAGATTTCCTGGTCATCGGAATATATCAATTCCACTTATAATGCCCCCGGCGTGCAAAGCCCGGTTATCGACTCGCTGATCAACCAAATTATTGCCGCGCAGGGAAATAAAGAAAAATTACTGCCGTTGGGGCGAGCACTGGATCGCGTATTAACGTGGAATTATTACATGCTGCCAATGTGGTACATGGCGGAAGACCGTCTCGCCTGGTGGGATAAATTCTCC  >  NC_000913/2273791‑2274076
                                                                                                                                            |                                                                                                                                                 
tCGCGACTATGACATGATGCCGCGCGTATGGCGGGCGATGCCGTGGCCCAGTTCCGATTTACAGATTTCCTGGTCATCGGAATATATCAATTCCACTTATAATGCCCCCGGCGTGCAAAGCCCGGTTATCGACTCGCTGATCAACCaaa                                                                                                                                           >  2:28700/1‑149 (MQ=255)
  gcgACTATGACATGATGCCGCGCGTATGGCGGGCGATGCCGTGGCCCAGTTCCGATTTACAGATTTCCTGGTCATCGGAATATATCAATTCCACTTATAATGCCCCCGGCGTGCAAAGCCCGGTTATCGACTCGCTGATCAACCAAat                                                                                                                                          >  7:203620/1‑148 (MQ=255)
  gcgACTATGACATGATGCCGCGCGTATGGCGGGCGATGCCGTGGCCCAGTTCCGATTTACAGATTTCCTGGTCATCGGAATATATCAATTCCACTTATAATGCCCCCGGCGTGCAAAGCCCGGTTATCGACTCGCTGATCAACCAAat                                                                                                                                          <  1:274873/148‑1 (MQ=255)
  gcgACTATGACATGATGCCGCGCGTATGGCGGGCGATGCCGTGGCCCAGTTCCGATTTACAGATTTCCTGGTCATCGGAATATATCAATTCCACTTATAATGCCCCCGGCGTGCAAAGCCCGGTTATCGACTCGCTGATCAACCAAat                                                                                                                                          >  4:317046/1‑148 (MQ=255)
   cgACTATGACATGATGCCGCGCGTATGGCGGGCGATGCCGTGGCCCAGTTCCGATTTACAGATTTCCTGGTCATCGGAATATATCAATTCCACTTATAATGCCCCCGGCGTGCAAAGCCCGGTTATCGACTCGCTGATCAACCAAatta                                                                                                                                        >  4:288495/1‑149 (MQ=255)
    gACTATGACATGATGCCGCGCGTATGGCGGGCGATGCCGTGGCCCAGTTCCGATTTACAGATTTCCTGGTCATCGGAATATATCAATTCCACTTATAATGCCCCCGGCGTGCAAAGCCCGGTTATCGACTCGCTGATCAACCAAattat                                                                                                                                       >  1:321754/1‑149 (MQ=255)
       tATGACATGATGCCGCGCGTATGGCGGGCGATGCCGTGGCCCAGTTCCGATTTACAGATTTCCTGGTCATCGGAATATATCAATTCCACTTATAATGCCCCCGGCGTGCAAAGCCCGGTTATCGACTCGCTGATCAACCAAATTATTGc                                                                                                                                    >  8:331890/1‑149 (MQ=255)
       tATGACATGATGCCGCGCGTATGGCGGGCGATGCCGTGGCCCAGTTCCGATTTACAGATTTCCTGGTCATCGGAATATATCAATTCCACTTATAATGCCCCCGGCGTGCAAAGCCCGGTTATCGACTCGCTGATCAACCAAATTATTGc                                                                                                                                    >  4:121334/1‑149 (MQ=255)
       tATGACATGATGCCGCGCGTATGGCGGGCGATGCCGTGGCCCAGTTCCGATTTACAGATTTCCTGGTCATCGGAATATATCAATTCCACTTATAATGCCCCCGGCGTGCAAAGACCGGTTATCGACTCGCTGATCAACCAAATTATTGc                                                                                                                                    >  5:227329/1‑149 (MQ=255)
         tGACATGATGCCGCGCGTATGGCGGGCGATGCCGTGGCCCAGTTCCGATTTACAGATTTCCTGGTCATCGGAATATATCAATTCCACTTATAATGCCCCCGGCGTGCAAAGCCCGGTTATCGACTCGCTGATCAACCAAATTATTGcc                                                                                                                                   <  1:81073/148‑1 (MQ=255)
         tGACATGATGCCGCGCGTATGGCGGGCGATGCCGTGGCCCAGTTCCGATTTACAGATTTCCTGGTCATCGGAATATATCAATTCCACTTATAATGCCCCCGGCGTGCAAAGCCCGGTTATCGACTCGCTGATCAACCAAATTATTGcc                                                                                                                                   <  6:141402/148‑1 (MQ=255)
            cATGATGCCGCGCGTATGGCGGGCGATGCCGTGGCCCAGTTCCGATTTACAGATTTCCTGGTCATCGGAATATATCAATTCCACTTATAATGCCCCCGGCGTGCAAAGCCCGGTTATCGACTCGCTGATCAACCAAATTATTGCcgcg                                                                                                                                >  3:284756/1‑148 (MQ=255)
                     gcgcgTATGGCGGGCGATGCCGTGGCCCAGTTCCGATTTACAGATTTCCTGGTCATCGGAATATATCAATTCCACTTATAATGCCCCCGGCGTGCAAAGCCCGGTTATCGACTCGCTGATCAACCAAATTATTGCCGCGCAGGGAAAt                                                                                                                       >  8:278408/1‑148 (MQ=255)
                     gcgcgTATGGCGGGCGATGCCGTGGCCCAGTTCCGATTTACAGATTTCCTGGTCATCGGAATATATCAATTCCACTTATAATGCCCCCGGCGTGCAAAGCCCGGTTATCGACTCGCTGATCAACCAAATTATTGCCGCGCAGGGAAATa                                                                                                                      <  1:318460/149‑1 (MQ=255)
                       gcgTATGGCGGGCGATGCCGTGGCCCAGTTCCGATTTACAGATTTCCTGGTCATCGGAATATATCAATTCCACTTATAATGCCCCCGGCGTGCAAAGCCCGGTTATCGACTCGCTGATCAACCAAATTATTGCCGCGCAGGGAAATaaa                                                                                                                    <  3:127598/149‑1 (MQ=255)
                         gTATGGCGGGCGATGCCGTGGCCCAGTTCCGATTTACAGATTTCCTGGTCATCGGAATATATCAATTCCACTTATAATGCCCCCGGCGTGCAAAGCCCGGTTATCGACTCGCTGATCAACCAAATTATTGCCGCGCAGGGAAATAAAGa                                                                                                                  <  4:4968/149‑1 (MQ=255)
                           aTGGCGGGCGATGCCGTGGCCCAGTTCCGATTTACAGATTTCCTGGTCATCGGAATATATCAATTCCACTTATAATGCCCCCGGCGTGCAAAGCCCGGTTATCGACTCGCTGATCAACCAAATTATTGCAGCGCAGGGAAATAAAGaa                                                                                                                 >  6:285909/1‑148 (MQ=255)
                            tGGCGGGCGATGCCGTGGCCCAGTTCCGATTTACAGATTTCCTGGTCATCGGAATATATCAATTCCACTTATAATGCCCCCGGCGTGCAAAGCCCGGTTATCGACTCGCTGATCAAc                                                                                                                                               >  8:272983/1‑117 (MQ=255)
                            tGGCGGGCGATGCCGTGGCCCAGTTCCGATTTACAGATTTCCTGGTCATCGGAATATATCAATTCCACTTATAATGCCCCCGGCGTGCAAAGCCCGGTTATCGACTCGCTGATCAAc                                                                                                                                               <  7:272983/117‑1 (MQ=255)
                            tGGCGGGCGATGCCGTGGCCCAGTTCCGATTTACAGATTTCCTGGGCATCGGAATATATCAATTCCACTTATAATGCCCCCGGCGTGCAAAGCCCGGTTATCGACTCGCTGATCAACCAAATTATTGCCGCGCAGGGAAATAAAGaaaa                                                                                                               <  3:100613/149‑1 (MQ=255)
                              gcgggcgATGCCGTGGCCCAGTTCCGATTTACAGATTTCCTGGTCATCGGAATATATCAATTCCACTTATAATGCCCCCGGCGTGCAAAGCCCGGTTATCGACTCGCTGATCAACCAAATTATTGCCGCGCAGGGAAATAAAGaaaaa                                                                                                              <  8:203620/148‑1 (MQ=255)
                                      tGCCGTGGCCCAGTTCTGATTTACAGATTTCCTGGTCATCGGAATATATCAATTCCACTTATAATGCCCCCGGCGTGCAAAGCCCGGTTATCGACTCGCTGATCAACCAAATTATTGCCGCGCAGGGAAATAAAGAAAAATTACTGCCg                                                                                                     >  4:328002/1‑149 (MQ=255)
                                      tGCCGTGGCCCAGTTCCGATTTACAGATTTCCTGGTCATCGGAATATATCAATTCCACTTATAATGCCCCCGGCGTGCAAAGCCCGGTTATCGACTCGCTGATCAACCAAATTATTGCCGCGCAGGGAAATAAAGAAAAATTACTGCCg                                                                                                     >  1:198704/1‑149 (MQ=255)
                                        ccGTGGCCCAGTTCCGATTTACAGATTTCCTGGTCATCGGAATATATCAATTCCACTTATAATGCCCCCGGCGTGCAAAGCCCGGTTATCGACTCGCTGATCAACCAAATTATTGCCGCGCAGGGAAATAAAGAAAAATTACTGCCGtt                                                                                                   <  8:34097/149‑1 (MQ=255)
                                           tGGCCCAGTTCCGATTTACAGATTTCCTGGTCATCGGAATATATCAATTCCACTTATAATGCCCCCGGCGTGCAAAGCCCGGTTATCGACTCGCTGATCAACCAAATTATTGCCGCGCAGGGAAATAAAGAAAAATTACTGCCGTTggg                                                                                                >  1:272883/1‑149 (MQ=255)
                                             gCCCAGGTCCGATTTACAGAATTCCTGGTCATCGGAATAAATCAATTCCACTTATAATGCCCCCGGCGTGCAAAGCCCGGTTATCGACTCGCTGATCAACCAAATTATTGCCGCGCAGGGAAATAAAGAAAAATTACTGCCGTTGGGGc                                                                                              <  8:342888/149‑1 (MQ=255)
                                               ccAGTTCCGATTTACAGATTTCCTGGTCATCGGAATATATCAATTCCACTTATAATGCCCCCGGCGTGCAAAGCCCGGTTATCGACTCGCTGATCAACCAAATTATTGCCGCGCAGGGAAATAAAGAAAAATTACTGCCGTTGGGGCGa                                                                                            <  2:321754/149‑1 (MQ=255)
                                                 aGTTCCGATTTACAGATTTCCTGGTCATCGGAATATATCAATTCCACTTATAATGCCCCCGGCGTGCAAAGCCCGGTTATCGACTCGCTGATCAACCAAATTATTGCCGCGCAGGGAAATAAAGAAAAATTa                                                                                                           >  5:82773/1‑132 (MQ=255)
                                                 aGTTCCGATTTACAGATTTCCTGGTCATCGGAATATATCAATTCCACTTATAATGCCCCCGGCGTGCAAAGCCCGGTTATCGACTCGCTGATCAACCAAATTATTGCCGCGCAGGGAAATAAAGAAAAATTa                                                                                                           <  6:82773/132‑1 (MQ=255)
                                                   ttCCGATTTACAGATTTCCTGGTCATCGGAATATATCAATTCCACTTATAATGCCCCCGGCGTGCAAAGCCCGGTTATCGACTCGCTGATCAACCAAATTATTGCCGCGCAGGGAAATAAAGAAAAATTACTGCCGTTGGGGCGAGCAc                                                                                        <  3:121334/149‑1 (MQ=255)
                                                    tCCGATTTACAGATTTCCTGGTCATCGGAATATATCAATTCCACTTATAATGCCCCCGGCGTGCAAAGCCCGGTTATCGACTCGCTGATCAACCAAATTATTGCCGCGCAGGGAAATAAAGAAAAATTACTGCCGTTGGGGCGAGCAc                                                                                        <  1:28700/148‑1 (MQ=255)
                                                            aCAGATTTCCTGGTCATCGGAATATATCAATTCCACTTATAATGCCCCCGGCGTGCAAAGCCCGGTTATCGACTCGCTGATCAACCAAATTATTGCCGCGCAGGGAAATAAAGAAAAATTACTGCCGTTGGGGCGAGCACTGGATcgcg                                                                               >  7:319149/1‑149 (MQ=255)
                                                              aGATTTCCTGGTCATCGGAATATATCAATTCCACTTATAATGCCCCCGGCGTGCAAAGCCCGGTTATCGACTCGCTGATCAACCAAATTATTGCCGCGCAGGGAAATAAAGAAAAATTACTGCCGTTGGGGCGAGCACTGGATCGCGTa                                                                             <  2:50836/149‑1 (MQ=255)
                                                               gATTTCCTGGTCATCGGAATATATCAATTCCACTTATAATGCCCCCGGCGTGCAAAGCCCGGTTATCGACTCGCTGATCAACCAAATTATTGCCGCGCAGGGAAATAAAGAAAAATTACTGCCGTTGGGGCGAGCACTGGATCGCGTAt                                                                            >  4:70458/1‑149 (MQ=255)
                                                               gATTTCCTGGTCATCGGAATATATCAATTCCACTTATAATGCCCCCGGCGTGCAAAGCCCGGTTATCGACTCGCTGATCAACCAAATTATTGCCGCGCAGGGAAATAAAGAAAAATTACTGCCGTTGGGGCGAGCACTGGATCGCGTAt                                                                            <  3:70458/149‑1 (MQ=255)
                                                                 tttCCTGGTCATCGGAATATATCAATTCCACTTATAATGCCCCCGGCGTGCAAAGCCCGGTTATCGACTCGCTGATCCCCCAAATTATTGCCGCGCAGGGAAATAAAGAAAAATTACTGCCGTTGGGGCGAGCACGGGATCGCGTAtt                                                                           >  6:103303/1‑148 (MQ=255)
                                                                   tCCTGGTCATCGGAATATATCAATTCCACTTATAATGCCCCCGGCGTGCAAAGCCCGGTTATCGACTCGCTGATCAACCATATTATTGCCGCGCAGGGAAATAAAGAAAAATTACTGCCGTTGGGGCGAGCACTGGATCGCGTATTaa                                                                         <  8:338399/148‑1 (MQ=255)
                                                                      tGGTCATCGGAATATATCAATTCCACTTATAATGCCCCCGGCGTGCAAAGCCCGGTTATCGACTCGCTGATCAACCAAATTATTGCCGCGCAGGGAAATAAAGAAAAATTACTGCCGTTGGGGCGAGCACTGGATCGCGTATTAACGTg                                                                     >  4:174668/1‑149 (MQ=255)
                                                                        gTCATCGGAATATATCAATTCCACTTATAATGCCCCCGGCGTGCAAAGCCCGGTTATCGACTCGCTGATCAACCAAATTATTGCCGCGCAGGGAAATAAAGAAAAATTACTGCCGTTGGGGCGAGCACTGGATCGCGTATTAACGTgg                                                                    <  2:198704/148‑1 (MQ=255)
                                                                           aTCGGAATATATCAATTCCACTTATAATGCCCCCGGCGTGCAAAGCCCGGTTATCGACTCGCTGATCAACCAAATTATTGCCGCGCAGGGAAATAAAGAAAAATTACTGCCGTTGGGGCGAGCACTGGATCGCGTATTAACGTGGAat                                                                 >  7:3998/1‑148 (MQ=255)
                                                                            tCGGAATATATCAATTCCACTTATAATGCCCCCGGCGTGCAAAGCCCGGTTATCGACTCGCTGAACAACCAAATTATTGCCGCGCAGGGAAATAAAGAAAAATTACTGCCGTTGGGGCGAg                                                                                           >  7:26106/1‑121 (MQ=255)
                                                                            tCGGAATATATCAATTCCACTTATAATGCCCCCGGCGTGCAAAGCCCGGTTATCGACTCGCTGAACAACCAAATTATTGCCGCGCAGGGAAATAAAGAAAAATTACTGCCGTTGGGGCGAg                                                                                           <  8:31816/121‑1 (MQ=255)
                                                                            tCGGAATATATCAATTCCACTTATAATGCCCCCGGCGTGCAAAGCCCGGTTATCGACTCGCTGAACAACCAAATTATTGCCGCGCAGGGAAATAAAGAAAAATTACTGCCGTTGGGGCGAg                                                                                           >  7:31816/1‑121 (MQ=255)
                                                                            tCGGAATATATCAATTCCACTTATAATGCCCCCGGCGTGCAAAGCCCGGTTATCGACTCGCTGAACAACCAAATTATTGCCGCGCAGGGAAATAAAGAAAAATTACTGCCGTTGGGGCGAg                                                                                           <  8:26106/121‑1 (MQ=255)
                                                                              ggAATATATCAATTCCACTTATAATGCCCCCGGCGTGCAAAGCCCGGTTATCGACTCGCTGATCAACCAAATTATTGCCGCGCAGGGAAATAAAGAAAAATTACTGCCGTTGGGGCGAGCACTGGATCGCGTATTAACGTGGAattatt                                                             <  3:317046/149‑1 (MQ=255)
                                                                               gAATATATCAATTCCACTTATAATGCCCCCGGCGTGCAAAGCCCGGTTATCGACTCGCTGATCAACCAAATTATTGCCGCGCAGGGAAATAAAGAAAAATTACTGCCGTTGGGGCGAGCACTGGATCGCGTATTAACGTGGAAttatta                                                            <  8:3998/149‑1 (MQ=255)
                                                                                aaTATATCAATTCCACTTATAATGCCCCCGGCGTGCAAAGCCCGGTTATCGACTCGCTGATCAACCAAATTATTGCCGCGCAGGGAAATAAAGAAAAATTACTGCCGTTGGGGCGAGCACTGGATCGCGTATTAACGTGGAAttatta                                                            >  2:265862/1‑148 (MQ=255)
                                                                                      tCAATTCCACTTATAATGCCCCCGGCGTGCAAAGCCCGGTTATCGACTCGCTGATCAACCAAATTATTGCCGCGCAGGGAAATAAAGAAAAATTACTGCCGTTGGGGCGAGCACTGGATCGCGTATTAACGTGGAATTATTACAtgctg                                                     >  1:310928/1‑149 (MQ=255)
                                                                                      tCAATTCCACTTATAATGCCCCCGGCGTGCAAAGCCCGGTTATCGACTCGCTGATCAACCAAATTATTGCCGCGCAGGGAAATAAAGAAAAATTACTGCCGTTGGGGCGAGCACTGGATCGCGTATTAACGTGGAATTATTACAtgctg                                                     >  2:350065/1‑149 (MQ=255)
                                                                                      tCAATTCCACTTATAATGCCCCCGGCGTGCAAAGCCCGGTTATCGACTCGCTGATCAACCAAATTATTGCCGCGCAGGGAAATAAAGAAAAATTACTGCCGTTGGGGCGAGCACTGGATCGCGTATTAACGTGGAATTATTACAtgct                                                      <  7:278408/148‑1 (MQ=255)
                                                                                        aaTTCCACTTATAATGCCCCCGGCGTGCAAAGCCCGGTTATCGACTCGCTGATCAACCAAATTATTGCCGCGCAGGGAAATAAAGAAAAATTACTGCCGTTGGGGCGAGCACTGGATCGCGTATTAACGTGGAATTATTACAtgctgc                                                    >  7:342079/1‑148 (MQ=255)
                                                                                             cACTTATAATGCCCCCGGCGTGCAAAGCCCGGTTATCGACTCGCTGATCAACCAAATTATTGCCGCGCAGGGAAATAAAGAAAAATTACTGCCGTTGGGGCGAGCACTGGATCGCGTATTAACGTGGAATTATTACATGCTGCCAAtgt                                              >  7:90752/1‑149 (MQ=255)
                                                                                              aCTTATAATGCCCCCGGCGTGCAAAGCCCGGTTATCGACTCGCTGATCAACCAAATTATTGCCGCGCAGGGAAATAAAGAAAAATTACTGCCGTTGGGGCGAGCACTGGATCGCGTATTAACGTGGAATTATTACATGCTGCCAAtgtg                                             >  1:17141/1‑149 (MQ=255)
                                                                                              aCTTATAATGCCCCCGGCGTGCAAAGCCCGGTTATCGACTCGCTGATCAACCAAATTATTGCCGCGCAGGGAAATAAAGAAAAATTACTGCCGTTGGGGCGAGCACTGGATCGCGTATTAACGTGGAATTATTACATGCTGCCAAtgtg                                             >  5:184740/1‑149 (MQ=255)
                                                                                              aCTTATAATGCCCCCGGCGTGCAAAGCCCGGTTATCGACTCGCTGATCAACCAAATTATTGCCGCGCAGGGAAATAAAGAAAAATTACTGCCGTTGGGGCGAGCACTGGATCGCGTATTAACGTGGAATTATTACATGCTGCCAAtgtg                                             <  4:284756/149‑1 (MQ=255)
                                                                                              aCTTATAATGCCCCCGGCGTGCAAAGCCCGGTTATCGACTCGCTGATCAACCAAATTATTGCCGCGCAGGGAAATAAAGAAAAATTACTGCCGTTGGGGCGAACACTGGATCGCGTATTAACGTGGAATTATTACATGCTGCCAAtgtg                                             <  2:272883/149‑1 (MQ=255)
                                                                                                  ataATGCCCCCGGCGTGCAAAGCCCTGTTATCGACTCGCTGATCAACCAAATTATTGCCGCGCAGGGAAATAAAGAAAAATTACTGCCGTTGGGGCGAGCACTGGATCGCGTATTAACGTGGAATTATTACATGCTGCCAATGTGGTa                                          <  5:285909/148‑1 (MQ=255)
                                                                                                  ataATGCCCCCGGCGTGCAAAGCCCGGTTATCGACTCGCTGATCAACCAAATTATTGCCGCGCAGGGAAATAAAGAAAAATTACTGCCGTTGGGGCGAGCACTGGATCGCGTATTAAc                                                                        <  1:36704/118‑1 (MQ=255)
                                                                                                  ataATGCCCCCGGCGTGCAAAGCCCGGTTATCGACTCGCTGATCAACCAAATTATTGCCGCGCAGGGAAATAAAGAAAAATTACTGCCGTTGGGGCGAGCACTGGATCGCGTATTAAc                                                                        >  2:36704/1‑118 (MQ=255)
                                                                                                     aTGCCCCCGGCGTGCAAAGCCCGGTTATCGACTCGCTGATCAACCAAATTATTGCCGCGCAGGGAAATAAAGAAAAATTACTGCCGTTGGGGCGAGCACTGGATCGCGTATTAACGTGGAATTATTACATGCTGCCAATGTGGTACATg                                      <  6:184740/149‑1 (MQ=255)
                                                                                                      tGCCCCCGGCGTGCAAAGCCCGGTTATCGACTCGCTGATCAACCAAATTATTGCCGCGCAGGGAAATAAAGAAAAATTACTGCCGTTGGGGCGAGCACTGGATCGCGTATTAACGTGGAATTATTACATGCTGCCAATGTGGTACATgg                                     <  1:99364/149‑1 (MQ=255)
                                                                                                        cccccGGCGTGCAAAGCCCGGTTATCGACTCGCTGATCAACCAAATTATTGCCGCGCAGGGAAATAAAGAAAAATTACTGCCGTTGGGGCGAGCACTGGATCGCGTATTAACGTGGAATTATTACATGCTGCCAATGTGGTACATGGc                                    <  2:17141/148‑1 (MQ=255)
                                                                                                         ccccGGCGTGCAAAGCCCGGTTATCGACTCGCTGATCAACCAAATTATTGCCGCGCAGGGAAATAAAGAAAAATTACTGCCGTTGGGGCGAGCACTGGATCGCGTATTAACGTGGAATTATTACATGCTGCCAATGTGGTACATGGCg                                   >  2:242061/1‑148 (MQ=255)
                                                                                                                 tGCAAAGCCCGGTTATCGACTCGCTGATCAACCAAATTATTGCCGCGCAGGGAAATAAAGAAAAATTACTGCCGTTGGGGCGAGCACTGGATCGCGTATTAACGTGGAATTATTACATGCTGCCGATGTGGTACATGGCGGAAGACCGt                          >  3:149115/1‑149 (MQ=255)
                                                                                                                 tGCAAAGCCCGGTTATCGACTCGCTGATCAACCAAATTATTGCCGCGCAGGGAAATAAAGAAAAATTACTGCCGTTGGGGCGAGCACTGGATCGCGTATTAACGTGGAATTATTACATGCTGCCAATGTGGTACATGGCGGAAGACCGt                          >  2:81533/1‑149 (MQ=255)
                                                                                                                   cAAAGCCCGGTTATCGACTCGCTGATCAACCAAATTATTGCCGCGCAGGGAAATAAAGAAAAATTACTGCCGTTGGGGCGAGCACTGGATCGCGTATTAACGTGGAATTATTACATGCTGCCAATGTGGTACATGGCGGAAGACCGtct                        <  1:350065/149‑1 (MQ=255)
                                                                                                                        cccGGTTATCGACTCGCTGATCAACCAAATTATTGCCGCGCAGGGAAATAAAGAAAAATTACTGCCGTTGGGGCGAGCACTGGATCGCGTATTAACGTGGAATTATTACATGCTGCCAATGTGGTACATGGCGGAAGACCGTCTCGCCt                   >  5:344386/1‑149 (MQ=255)
                                                                                                                           ggTTATCGACTCGCTGATCAACCAAATTATTGCCGCGCAGGGAAATAAAGAAAAATTACTGCCGTTGGGGCGAGCACTGGATCGCGTATTAACGTGGAATTATTACATGCTGCCAATGTGGTACATGGCGGAAGACCGTCTCGCCtggt                >  1:39003/1‑149 (MQ=255)
                                                                                                                            gTTATCGACTCGCTGATCAACCAAATTATTGCCGCGCAGGGAAATAAAGAAAAATTACTGCCGTTGGGGCGAGCACTGGATCGCGTATTAACGTGGAATTATTACATGCTGCCAATGTGGTACATGGCGGAAGACCGTCTCGCCtggtg               <  8:319149/149‑1 (MQ=255)
                                                                                                                             ttATCGACTCGCTGATCAACCAAATTATTGCCGCGCAGGGAAATAAAGAAAAATTACTGCCGTTGGGGCGAGCACTGGATCGCGTATTAACGTGGAATTATTACATGCTGCCAATGTGGTACATGGCGGAAGACCGTCTCGCCtggtgg              >  7:275587/1‑149 (MQ=255)
                                                                                                                                  gACTCGCTGATCAACCAAATTATTGCCGCGCAGGGAAATAAAGAAAAATTACTGCCGTTGGGGCGAGCACTGGATCGCGTATTAACGt                                                                      >  2:28875/1‑88 (MQ=255)
                                                                                                                                  gACTCGCTGATCAACCAAATTATTGCCGCGCAGGGAAATAAAGAAAAATTACTGCCGTTGGGGCGAGCACTGGATCGCGTATTAACGt                                                                      <  1:28875/88‑1 (MQ=255)
                                                                                                                                   aCTCGCTGATCAACCAAATTATTGCCGCGCAGGGAAATAAAGAAAAATTACTGCCGTTGGGGCGAGCACTGGATCGCGTATTAACGTGGAATTATTACATGCTGCCAATGTGGTACATGGCGGAAGACCGTCTCGCCTGGTGGGATaaa        <  1:323538/149‑1 (MQ=255)
                                                                                                                                   aCTCGCTGATCAACCAAATTATTGCCGAGCAGGGAAATAAAGAAAAATTACTGCCGTTGGGGCGAGCACTGGATCGCGTATTAACGTGGAATTATTACATGCTGCCAATGTGGTACATGGCGGAAGACCGTCTCGCCTGGTGGGATaaa        <  6:227329/149‑1 (MQ=255)
                                                                                                                                        cTGATCAACCAAATTATTGCCGCGCAGGGAAATAAAGAAAAATTACTGCCGTTGGGGCGAGCACTGGATCGCGTATTAACGTGGAATTATTACATGCTGCCAATGTGGTACATGGCGGAAGACCGTCTCGCCTGGTGGGATAAATtctc   <  1:81533/149‑1 (MQ=255)
                                                                                                                                          gATCAACCAAATTATTGCCGCGCAGGGAAATAAAGAAAAATTACTGCCGTTGGGGCGAGCACTGGATCGCGTATTAACGTGGAATTATTACATGCTGCCAATGTGGTACATGGCGGAAGACCGTCTCGCCTGGTGGGATAAATTCTcc  >  2:281929/1‑148 (MQ=255)
                                                                                                                                            |                                                                                                                                                 
TCGCGACTATGACATGATGCCGCGCGTATGGCGGGCGATGCCGTGGCCCAGTTCCGATTTACAGATTTCCTGGTCATCGGAATATATCAATTCCACTTATAATGCCCCCGGCGTGCAAAGCCCGGTTATCGACTCGCTGATCAACCAAATTATTGCCGCGCAGGGAAATAAAGAAAAATTACTGCCGTTGGGGCGAGCACTGGATCGCGTATTAACGTGGAATTATTACATGCTGCCAATGTGGTACATGGCGGAAGACCGTCTCGCCTGGTGGGATAAATTCTCC  >  NC_000913/2273791‑2274076

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