Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,273,242 C→T 5.5% Y293Y (TAC→TAT 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,2420CT 5.5% 236.9 / 3.3 72Y293Y (TAC→TATyejAmicrocin C ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base C (31/37);  new base T (2/2);  total (33/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.58e-01

GACACCATTCGCTACGATTATTACCTCGATGATAATGTCGCCTTTGAAGCGTTTAAAGCAGGTGCCTTTGATTTGCGTATGGAAAACGACGCCAAAAACTGGGCCACGCGTTATACCGGTAAAAATTTCGATAAAAAATACATCATCAAAGATGAGCAAAAGAACGAATCAGCCCAGGATACGCGTTGGCTGGCGTTTAATATCCAACGTCCGGTATTCAGCGATCGCCGGGTCCGGGAAGCTATCACTCTCGCCTTTGACTTTGAATGGATGAACAAGGCGTTGTT  >  NC_000913/2273102‑2273388
                                                                                                                                            |                                                                                                                                                  
gACACCATTCGCTACGATTATTACCTCGATGATAATGTCGCCTTTGAAGCGTTTAAAGCAGGTGCCTTTGATTTGCGTATGGAAAACGACGCCAAAAACTGGGCCACGCGTTATACCGGTAAAAATTTCGATAAAAAATACATCATCaa                                                                                                                                            <  1:339243/149‑1 (MQ=255)
gACACCATTCGCTACGATTATTACCTCGATGATAATGTCGCCTTTGAAGCGTTTAAAGCAGGTGCCTTTGATTTGCGTATGGAAAACGACGCCAAAAACTGGGCCACGCGTTATACCGGTAAAAATTTCGATAAAAAATACATCATCaa                                                                                                                                            >  4:501755/1‑149 (MQ=255)
 acacCATTCGCTACGATTATTACCTCGATGATAATGTCGCCTTTGAAGCGTTTAAAGCAGGTGCCTTTGATTTGCGTATGGAAAACGACGCCAAAAACTGGGCCACGCGTTATACCGGTAAAAATTTCGATAAAAAATACATCATCaaa                                                                                                                                           >  3:285125/1‑149 (MQ=255)
  cacCATTCGCTACGATTATTACCTCGATGATAATGTCGCCTTTGAAGCGTTTAAAGCAGGTGCCTTTGATTTGCGTATGGAAAACGACGCCAAAAACTGGGCCACGCGTTATACCGGTAAAAATTTCGATAAAAAATACATCATCAAAg                                                                                                                                          >  3:444700/1‑149 (MQ=255)
          gCTACGATTATTACCTCGATGATAATGTCGCCTTTGAAGCGTTTAAAGCAGGTGCCTTTGATTTGCGTATGGAAAACGACGCCAAAAACTGGGCCACGCGTTATACCGGTAAAAATTTCGATAAAAAATATATCATCAAAGATGAGCaa                                                                                                                                  >  2:676603/1‑149 (MQ=255)
          gCTACGATTATTACCTCGATGATAATGTCGCCTTTGAAGCGTTTAAAGCAGGTGCCTTTGATTTGCGTATGGAAAACGACGCCAAAAACTGGGCCACGCGTTATACCGGTAAAAATTTCGATAAAAAATATATCATCAAAGATGAGCaa                                                                                                                                  >  2:383612/1‑149 (MQ=255)
            tACGATTATTACCTCGATGATAATGTCGCCTTTGAAGCGTTTAAAGCAGGTGCCTTTGATTTGCGTATGGAAAACGACGCCAAAAACTGGGCCACGCGTTATACCGGTAAAAATTTCGATAAAAAATACATCATCAAAGATGAGCaaaa                                                                                                                                <  3:667888/149‑1 (MQ=255)
            tACGATTATTACCTCGATGATAATGTCGCCTTTGAAGCGTTTAAAGCAGGTGCCTTTGATTTGCGTATGGAAAACGACGCCAAAAACTGGGCCACGCGTTATACCGGTAAAAATTTCGATAAAAAATACATCATCAAAGATGAGCaaaa                                                                                                                                <  3:707432/149‑1 (MQ=255)
            tACGATTATTACCTCGATGATAATGTCGCCTTTGAAGCGTTTAAAGCAGGTGCCTTTGATTTGCGTATGGAAAACGACGCCAAAAACTGGGCCACGCGTTATACCGGTAAAAATTTCGATAAAAAATACATCATCAAAGATGAGCaaaa                                                                                                                                >  4:543400/1‑149 (MQ=255)
               gATTATTACCTCGATGATAATGTCGCCTTTGAAGCGTTTAAAGCAGGTGCCTTTGATTTGCGTATGGAAAACGACGCCAAAAACTGGGCCACGCGTTATACCGGTAAAAATTTCGATAAAAAATACATCATCAAAGATGAGCAAAAGaa                                                                                                                             <  2:199019/149‑1 (MQ=255)
               gATTATTACCTCGATGATAATGTCGCCTTTGAAGCGTTTAAAGCAGGTGCCTTTGATTTGCGTATGGAAAACGACGCCAAAAACTGGGCCACGCGTTATACCGGTAAAAATTTCGATAAAAAATACATCATCAAAGACGAGCAAAAGaa                                                                                                                             <  4:84154/149‑1 (MQ=255)
                      aCCTCGATGATAATGTCGCCTTTGAAGCGTTTAAAGCAGGTGCCTTTGATTTGCGTATGGAAAACGACGCCAAAAACTGGGCCACGCGTTATACCGGTAAAAATTTCGATAAAAAATACATCATCAAAGATGAGCAAAAGAACGAATCa                                                                                                                      <  4:444700/149‑1 (MQ=255)
                        cTCGATGATAATGTCGCCTTTGAAGCGTTTAAAGCAGGTGCCTTTGATTTGCGTATGGAAAACGACGCCAAAAACTGGGCCACGCGTTATACCGGTAAAAATTTCGATAAAAAATACATCATCAAAGATGAGCAAAAGAACGAATCAGc                                                                                                                    >  4:684102/1‑149 (MQ=255)
                         tCGATGATAATGTCGCCTTTGAAGCGTTTAAAGCAGGTGCCTTTGATTTGCGTATGGAAAACGACGCCAAAAACTGGGCCACGCGTTATACCGGTAAAAATTTCGATAAAAAATACATCATCAAAGATGAGCAAAGGAACGAATCAGcc                                                                                                                   >  3:35882/1‑149 (MQ=255)
                                 aaTGTCGCCTTTGAAGCGTTTAAAGCAGGTGCCTTTGATTTGCGTATGGAAAACGACGCCAAAAACTGGGCCACGCGTTATACCGGTAAAAATTTCGATAAAAAATACATCATCAAAGATGAGCAAAAGAACGCATCAGCCCAGGATAc                                                                                                           <  3:468720/149‑1 (MQ=255)
                                    gTCGCCTTTGAAGCGTTTAAAGCAGGTGCCTTTGATTTGCGTATGGAAAACGACGCCAAAAACTGGGCCACGCGTTATACCGGTAAAAATTTCGATAAAAAATAca                                                                                                                                                   >  4:538377/1‑106 (MQ=255)
                                    gTCGCCTTTGAAGCGTTTAAAGCAGGTGCCTTTGATTTGCGTATGGAAAACGACGCCAAAAACTGGGCCACGCGTTATACCGGTAAAAATTTCGATAAAAAATAca                                                                                                                                                   <  3:538377/106‑1 (MQ=255)
                                       gCCTTTGAAGCGTTTAAAGCAGGTGCCTTTGATTTGCGTATGGAAAACGACGCCAAAAACTGGGCCACGCGTTATACCGGTAAAAATTTCGATAAAAAATACATCATCAAAGATGAGCAAAAGAACGAATCAGCCCAGGATACGCGTtg                                                                                                     >  1:406649/1‑149 (MQ=255)
                                        ccTTTGAAGCGTTTAAAGCAGGTGCCTTTGATTTGCGTATGGAAAACGACGCCAAAAACTGGGCCACGCGTTATACCGGTAAAAATTTCGATAAAAAATACATCATCAAAGATGAGCAAAAGAACGAATCAGCCCAGGATACGCGTtgg                                                                                                    <  3:298117/149‑1 (MQ=255)
                                          tttGAAGCGTTTAAAGCAGGTGCCTTTGATTTGCGTATGGAAAACGACGCCAAAAACTGGGCCACGCGTTATACCGGTAAAAATTTCGATAAAAAATACATCATCAAAGATGAGCAAAAGAACGAATCAGCCCAGGATACGCGTtggct                                                                                                  <  1:19088/149‑1 (MQ=255)
                                             gAAGCGTTTAAAGCAGGTGCCTTTGATTTGCGTATGGAAAACGACGCCAAAAACTGGGCCACGCGTTATACCGGTAAAAATTTCGATAAAAAATACATCATCAAAGATGAGCAAAAGAACGAATCAGCCCAGGATACGCGTtggctggc                                                                                               >  3:309415/1‑149 (MQ=255)
                                               aGCGTTTAAAGCAGGTGCCTTTGATTTGCGTATGGAAAACGACGCCAAAAACTGGGCCACGCGTTATACCGGTAAAAATTTCGATAAAAAATACATCATCAAAGATGAGCAAAAGAACGAATCAGCCCAGGATACGCGTTGGCTGGCGt                                                                                             >  4:512727/1‑149 (MQ=255)
                                                 cGTTTAAAGCAGGTGCCTTTGATTTGCGTATGGAAAACGACGCCAAAAACTGGGCCACGCGTTATACCGGTAAAAATTTCGATAAAAAATACATCATCAAAGATGAGCAAAAGAACGAATCAGCCCAGGATACGCGTTGGCTGGCGttt                                                                                           <  3:698462/149‑1 (MQ=255)
                                                     tAAAGCAGGTGCCTTTGATTTGCGTATGGAAAACGACGCCAAAAACTGGGCCACGCGTTATACCGGTAAAAATTTCGATAAAAAATACATCATCAAAGATGAGCAAAAGAACGAATCAGCCCAGGATACGCGTTGGCTGGCGTTTAata                                                                                       <  3:114129/149‑1 (MQ=255)
                                                         gCAGGTGCCTTTGATTTGCGTATGGAAAACGACGCCAAAAACTGGGCCACGCGTTATACCGGTAAAAATTTCGATAAAAAATACATCATCAAAGATGAGCAAAAGAACGAATCAGCCCAGGATACGCGTTGGCTGGCGTTTAATATCCa                                                                                   >  1:136210/1‑149 (MQ=255)
                                                          cAGGTGCCTTTGATTTGCGTATGGAAAACGACGCCAAAAACTGGGCCACGCGTTATACCGGTAAAAATTTCGATAAAAAATACATCATCAAAGATGAGCAAAAGAACGAATCAGCCCAGGATACGCGTTGGCTGGCGTTTAATATCCaa                                                                                  <  2:363041/149‑1 (MQ=255)
                                                           aGGTGCCTTTGATTTGCGTATGGAAAACGACGCCAAAAACTGGGCCACGCGTTATACCGGTAAAAATTTCGATAAAAAATACATCATCAAAGATGAGCAAAAGAACGAATCAGCCCAGGATACGCGTtggctggc                                                                                               <  4:164187/135‑1 (MQ=255)
                                                           aGGTGCCTTTGACTGGCGTATGGAAAACGACGCCAAAAACTGGGCCACGCGTTATACCGGTAAAAATTTCGATAAAAAATACATCATCAAAGATGAGCAAAAGAACGAATCAGCCCAGGATACGCGTtggctggc                                                                                               >  3:164187/1‑135 (MQ=255)
                                                             gTGCCTTTGATTTGCGTATGGAAAACGACGCCAAAAACTGGGCCACGCGTTATACCGGTAAAAATTTCGATAAAAAATACATCATCAAAGATGAGCAAAAGAACGAATCAGCCCAGGATACGCGTTGGCTGGCGTTTAATATCCAACGt                                                                               >  1:44015/1‑149 (MQ=255)
                                                                ccTTTGATTTGCGTATGGAAAACGACGCCAAAAACTGGGCCACGCGTTATACCGGTAAAAATTTCGATAAAAAATATATCATCAAAGATGAGCAAAAGAACGAATCAGCCCAGGATACGCGTTGGCTGGCGTTTAATATCCAACGTCCg                                                                            <  1:676603/149‑1 (MQ=255)
                                                                ccTTTGATTTGCGTATGGAAAACGACGCCAAAAACTGGGCCACGCGTTATACCGGTAAAAATTTCGATAAAAAATATATCATCAAAGATGAGCAAAAGAACGAATCAGCCCAGGATACGCGTTGGCTGGCGTTTAATATCCAACGTCCg                                                                            <  1:383612/149‑1 (MQ=255)
                                                                ccTTTGATTTGCGTATGGAAAACGACGCCAAAAACTGGGCCACGCGTTATACCGGTAAAAATTTCGATAAAAAATACATCATCAAAGATGAGCAAAAGAACGAATCAGCCCAGGATACGCGTTGGCTGGCGTTTAATATCCAACGTCCg                                                                            <  2:25601/149‑1 (MQ=255)
                                                                ccTTTGATTTGCGTATGGAAAACGACGCCAAAAACTGGGCCACGCGTTATACCGGTAAAAATTTCGATAAAAAATACATCATCAAAGATGAGCAAAAGAACGAATCAGCCCAGGATACGCGTTGGCTGGCGTTTAATATCCAACGTCCg                                                                            >  4:548617/1‑149 (MQ=255)
                                                                 cTTTGATTTGCGTATGGAAAACGACGCCAAAAACTGGGCCACGCGTTATACCGGTAAAAATTTCGATAAAAAATACATCATCAAAGATGAGCAAAAGAACGAATCAGCCCAGGATACGCGTTGGCTGGCGTTTAATATCCAACGTCCgg                                                                           >  2:113649/1‑149 (MQ=255)
                                                                  tttGATTTGCGTATGGAAAACGACGCCAAAAACTGGGCCACGCGTTATACCGGTAAAAATTTCGATAAAAAATACATCATCAAAGATGAGCAAAAGAACGAATCAGCCCAGGATACGCGTTGGCTGGCGTTTAATATCCAACGTCCGGt                                                                          >  1:478037/1‑149 (MQ=255)
                                                                  tttGATTTGCGTATGGAAAACGACGCCAAAAACTGGGCCACGCGTTATACCGGTAAAAATTTCGATAAAAAATACATCATCAAAGATGAGCAAAAGAACGAATCAGCCCAGGATACGCGTTGGCTGGCGTTTAATATCCAACGTCCGGt                                                                          >  4:209797/1‑149 (MQ=255)
                                                                     gATTTGCGTATGGAAAACGACGCCAAAAACTGGGCCACGCGTTATACCGGTAAAAATTTCGATAAAAAATACATCATCAAAGATGAGCAAAAGAACGAATCAGCCCAGGATACGCGTTGGCTGGCGTTTAATATCCAACGTCCGGTAtt                                                                       >  2:489520/1‑149 (MQ=255)
                                                                       tttGCGTATGGAAAACGACGCCAAAAACTGGGCCACGCGTTATACCGGTAAAAATTTCGATAAAAAATACATCATCAAAGATGAGCAAAAGAACGAATCAGCCCAGGa                                                                                                              <  2:152384/108‑1 (MQ=255)
                                                                       tttGCGTATGGAAAACGACGCCAAAAACTGGGCCACGCGTTATACCGGTAAAAATTTCGATAAAAAATACATCATCAAAGATGAGCAAAAGAACGAATCAGCCCAGGa                                                                                                              >  1:152384/1‑108 (MQ=255)
                                                                       tttGCGTATGGAAAACGACGCCAAAAACTGGGCCACGCGTTATACCGGTAAAAATTTCGATAAAAAATACATCATCAAAGATGAGCAAAAGAACGAATCAGCCCAGGATACGCGTTGGCTGGCGTTTAATATCCAACGTCCGGTATTCa                                                                     <  4:188671/149‑1 (MQ=255)
                                                                           cGTATGGAAAACGACGCCAAAAACTGGGCCACGCGTTATACCGGTAAAAATTTCGATAAAAAATACATCATCAAAGATGAGCAAAAGAACGAATCAGCCCAGGATACGCGTTGGCTGGCGTTTAATATCCAACGTCCGGTATTCAGCGa                                                                 <  4:309415/149‑1 (MQ=255)
                                                                            gTATGGAAAACGACGCCAAAAACTGGGCCACGCGTTATACCGGTAAAAATTTCGATAAAAAATACATCATCAAAGATGAGCAAAAGAACGAATCAGCCCAGGATACGCGTTGGCTGGCGTTTAATATCCAACGTCCGGTATTCAGCGAt                                                                <  4:20824/149‑1 (MQ=255)
                                                                               tGGAAAACGACGCCAAAAACTGGGCCACGCGTTATACCGGTAAAAATTTCGATAAAAAATACATCATCAAAGATGAGCAAAAGAAc                                                                                                                            <  3:695114/86‑1 (MQ=255)
                                                                               tGGAAAACGACGCCAAAAACTGGGCCACGCGTTATACCGGTAAAAATTTCGATAAAAAATACATCATCAAAGATGAGCAAAAGAAc                                                                                                                            >  4:695114/1‑86 (MQ=255)
                                                                                 gAAAACGACGCCAAAAACTGGGCCACGCGTTATACCGGTAAAAATTTCGATAAAAAATACATCATCAAAGATGAGCAAAAGAACGAATCAGCCCAGGATACGCGTTGGCTGGCGTTTAATATCCAACGTCCGGTATTCAGCGATCGCCg                                                           <  4:285125/149‑1 (MQ=255)
                                                                                  aaaaCGACGCCAAAAACTGGGCCACGCGTTATACCGGTAAAAATTTCGATAAAAAATACATCATCAAAGATGAGCAAAAGAACGAATCAGCCCAGGATACGCGTTGGCTGGCGTTTAATATCCAACGTCCGGTATTCAGCGATCGCCgg                                                          >  3:125854/1‑149 (MQ=255)
                                                                                  aaaaCGACGCCAAAAACTGGGCCACGCGTTATACCGGTAAAAATTTCGATAAAAAATACATCATCAAAGATGAGCAAAAGAACGAATCAGCCCAGGATACGCGTTGGCTGGCGTTTAATATCCAACGTCCGGTATTCAGCGATCGCCgg                                                          >  2:709772/1‑149 (MQ=255)
                                                                                    aaCGACGCCAAAAACTGGGCCACGCGTTATACCGGTAAAAATTTCGATAAAAAATACATCATCAAAGATGAGCAAAAGAACGAATCAGCCCAGGATACGCGTTGGCTGGCGTTTAATATCCAACGTCCGGTATTCAGCGATCGCCGGGt                                                        >  3:478869/1‑149 (MQ=255)
                                                                                    aaCGACGCCAAAAACTGGGCCACGCGTTATACCGGTAAAAATTTCGATAAAAAATACATCATCAAAGATGAGCAAAAGAACGAATCAGCCCAGGATACGCGTTGGCTGGCGTTTAATATCCAACGTCCGGTATTCAGCGATCGCCGGGt                                                        <  2:44707/149‑1 (MQ=255)
                                                                                    aaCGACGCCAAAAACTGGGCCACGCGTTATACCGGTAAAAATTTCGATAAAAAATACATCATCAAAGATGAGCAAAAGAACGAATCAGCCCAGGATACGCGTTGGCTGGCGTTTAATATCCAACGTCCGGTATTCAGCGATCGCCGGGt                                                        >  2:359383/1‑149 (MQ=255)
                                                                                            cAAAAACTGGGCCACGCGTTATACCGGTAAAAATTTCGATAAAAAATACATCATCAAAGATGAGCAAAAGAACGAATCAGCCCAGGATACGCGTTGGCTGGCGTTTAATATCCAACGTCCGGTATTCAGCGATCGCCGGGTCCGGGAAg                                                >  3:152332/1‑149 (MQ=255)
                                                                                               aaaCTGGGCCACGCGTTATACCGGTAAAAATTTCGATAAAAAATACATCATCAAAGATGAGCAAAAGAACGAATCAGCCCAGGATACGCGTTGGCTGGCGTTTAATATCCAACGTCCGGTATTCAGCGATCGCCGGGTCCGGGAAGCTa                                             <  4:152332/149‑1 (MQ=255)
                                                                                               aaaCTGGGCCACGCGTTATACCGGTAAAAATTTCGATAAAAAATACATCATCAAAGATGAGCAAAAGAACGAATCAGCCCAGGATACGCGTTGGCTGGCGTTTAATATCCAACGTCCGGTATTCAGCGATCGCCGGGTCCGGGAAGCTa                                             <  1:709772/149‑1 (MQ=255)
                                                                                                            cgTTATACCGGTAAAAATTTCGATAAAAAATACATCATCAAAGATGAGCAAAAGAACGAATCAGCCCAGGATACGCGTTGGCTGGCGTTTAATATCCAACGTCCGGTATTCAGCGATCGCCGGGTCCGGGAAGCTATCACTCTCGCctt                                <  2:74327/149‑1 (MQ=255)
                                                                                                             gTTATACCGGTAAAAATTTCGATAAAAAATACATCATCAAAGATGAGCAAAAGAACGAATCAGCCCAGGATACGCGTTGGCTGGCGTTTAATATCCAACGTCCGGTATTCAGCGATCGCCGGGTCCGGGAAGCTATCACTCTCGCcttt                               <  1:191195/149‑1 (MQ=255)
                                                                                                             gTTATACCGGTAAAAATTTCGATAAAAAATACATCATCAAAGATGAGCAAAAGAACGAATCAGCCCAGGATACGCGTTGGCTGGCGTTTAATATCCAACGTCCGGTATTCAGCGATCGCCAGGTCCGGGAAGCTATCTCTCTCGCcttt                               <  4:275964/149‑1 (MQ=255)
                                                                                                                  aCCGGTAAAAATTTCGATAAAAAATACATCATCAAAGATGAGCAAAAGAACGAATCAGCCCAGGATACGCGTTGGCTGGCGTTTAATATCCAACGTCCGGTATTCAGCGATCGCCGGGTCCGGGAAGCTATCACTCTCGCctttgactt                          <  1:113649/149‑1 (MQ=255)
                                                                                                                  aCCGGTAAAAATTTCGATAAAAAATACATCATCAAAGATGAGCAAAAGAACGAATCAGCCCAGGATACGCGTTGGCTGGCGTTTAATATCCAACGTCCGGTATTCAGCGATCGCCGGGTCCGGGAAGCTATCACTCTCGCctttgactt                          >  1:717837/1‑149 (MQ=255)
                                                                                                                   ccGGTAAAAATTTCGATAAAAAATACATCATCAAAGATGAGCAAAAGAACGAATCAGCCCAGGATACGCGTTGGCTGGCGTTTAATATCCAACGTCCGGTATTCAGCGATCGCCGGGTCCGGGAAGCTATCACTCTCGCctttgacttt                         <  2:136210/149‑1 (MQ=255)
                                                                                                                   ccGGTAAAAATTTCGATAAAAAATACATCATCAAAGATGAGCAAAAGAACGAATCAGCCCAGGATACGCGTTGGCTGGCGTTTAATATCCAACGTCCGGTATTCAGCGATCGCCGGGTCCGGGAAGCTATCACTCTCGCctttgacttt                         <  2:402570/149‑1 (MQ=255)
                                                                                                                       tAAAAATTTCGATAAAAAATACATCATCAAAGATGAGCAAAAGAACGAATCAGCCCAGGATACGCGTTGGCTGGCGTTTAATATCCAACGTCCGGTATTCAGCGATCGCCGGGTCCGGGAAGCTATCACTCTCGCCTTTGACTTTGAAt                     <  2:44015/149‑1 (MQ=255)
                                                                                                                       tAAAAATTTCGATAAAAAATACATCATCAAAGATGAGCAAAAGAACGAATCAGCCCAGGATACGCGTTGGCTGGCGTTTAATATCCAACGTCCGGTATTCAGCGATCGCCGGGTCCGGGAAGCTATCACTCTCGCCTTTGACTTTGAAt                     <  3:548617/149‑1 (MQ=255)
                                                                                                                       tAAAAATTTCGATAAAAAATACATCATCAAAGATGAGCAAAAGAACGAATCAGCCCAGGATACGCGTTGGCTGGCGTTTAATATCCAACGTCCGGTATTCAGCGATCGCCGGGTCCGGGAAGCTATCACTCTCGCCTTTGACTTTGAAt                     >  4:15714/1‑149 (MQ=255)
                                                                                                                       tAAAAATTTCGATAAAAAATACATCATCAAAGATGAGCAAAAGAACGAATCAGCCCAGGATACGCGTTGGCTGGCGTTTAATATCCAACGTCCGGTATTCAGCGATCGCCGGGTCCGGGAAGCTATCACTCTCGCCTTTGACTTTGAAt                     <  4:125854/149‑1 (MQ=255)
                                                                                                                          aaaTTTCGATAAAAAATACATCATCAAAGATGAGCAAAAGAACGAATCAGCCCAGGATACGCGTTGGCTGGCGTTTAATATCCAACGTCCGGTATTCAGCGATCGCCGGGTCCGGGAAGCTATCACTCTCGCCTTTGACTTTGAATGGa                  >  4:274824/1‑149 (MQ=255)
                                                                                                                          aaaTTTCGATAAAAAATACATCATCAAAGATGAGCAAAAGAACGAATCAGCCCAGGATACGCGTTGGCTGGCGTTTAATATCCAACGTCCGGTATTCAGCGATCGCCGGGTCCGGGAAGCTATCACTCTCGCCTTTGACTTTGAATGGa                  <  2:502677/149‑1 (MQ=255)
                                                                                                                          aaaTTTCGATAAAAAATACATCATCAAAGATGAGCAAAAGAACGAATCAGCCCAGGATACGCGTTGGCTGGCGTTTAATATCCAACGTCCGGTATTCAGCGATCGCCGGGTCCGGGAAGCTATCACTCTCGCCTTTGACTTTGAATGGa                  >  1:32839/1‑149 (MQ=255)
                                                                                                                            aTTTCGATAAAAAATACATCATCAAAGATGAGCAAAAGAACGAATCAGCCCAGGATACGCGTTGGCTGGCGTTTAATATCCAACGTCCGGTATTCAGCGATCGCCGGGTCCGGGAAGCTATCACTCTCGCCTTTGACTTTGAATGGATg                >  4:292157/1‑149 (MQ=255)
                                                                                                                             tttCGATAAAAAATACATCATCAAAGATGAGCAAAAGAACGAATCAGCCCAGGATACGCGTTGGCTGGCGTTTAATATCCAACGTCCGGTATTCAGCGATCGCCGGGTCCGGGAAGCTATCACTCTCGCCTTTGACTTTGAATGGATGa               >  1:502709/1‑149 (MQ=255)
                                                                                                                                   tAAAAAATACATCATCAAAGATGAGCAAAAGAACGAATCAGCCCAGGATACGCGTTGGCTGGCGTTTAATATCCAACGTCCGGTATTCAGCGATCGCCGGGTCCGGGAAGCTATCACTCTCGCCTTTGACTTTGAATGGATGAACAAg          >  2:255613/1‑148 (MQ=255)
                                                                                                                                   tAAAAAATACATCATCAAAGATGAGCAAAAGAACGAATCAGCCCAGGATACGCGTTGGCTGGCGTTTAATATCCAACGTCCGGTATTCAGCGATCGCCGGGTCCGGGAAGCTATCACTCTCGCCTTTGACTTTGAATGGATGAACAAg          <  1:255613/148‑1 (MQ=255)
                                                                                                                                    aaaaaaTACATCATCAAAGATGAGCAAAAGAACGAATCAGCCCAGGATACGCGTTGGCTGGCGTTTAATATCCAACGTCCGGTATTCAGCGATCGCCGGGTCCGGGAAGCTATCACTCTCGCCTTTGACTTTGAATGGATGAACAAGGc        <  1:111003/149‑1 (MQ=255)
                                                                                                                                        aaTACATCATCAAAGATGAGCAAAAGAACGAATCAGCCCAGGATACGCGTTGGCTGGCGTTTAATATCCAACGTCCGGTATTCAGCGATCGCCGGGTCCGGGAAGCTATCACTCTCGCCTTTGACTTTGAATGGATGAACAAGGCgttg    <  3:274824/149‑1 (MQ=255)
                                                                                                                                          tACATCATCAAAGATGAGCAAAAGAACGAAACAGCCCAGGATACGCGTTGGCTGGCGTTTAATATCCAACGTCCGGTATTCAGCGATCGCCGGGTCCGGGAAGCTATCACTCTCGCCTTTGACTTTGAATGGATGAACAAGGCgttgtt  <  2:751008/149‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                  
GACACCATTCGCTACGATTATTACCTCGATGATAATGTCGCCTTTGAAGCGTTTAAAGCAGGTGCCTTTGATTTGCGTATGGAAAACGACGCCAAAAACTGGGCCACGCGTTATACCGGTAAAAATTTCGATAAAAAATACATCATCAAAGATGAGCAAAAGAACGAATCAGCCCAGGATACGCGTTGGCTGGCGTTTAATATCCAACGTCCGGTATTCAGCGATCGCCGGGTCCGGGAAGCTATCACTCTCGCCTTTGACTTTGAATGGATGAACAAGGCGTTGTT  >  NC_000913/2273102‑2273388

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