Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 884,088 Δ1 bp 5.7% coding (416/1233 nt) mdfA → multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913884,0880G. 5.7% 231.1 / 7.8 70coding (416/1233 nt)mdfAmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base G (32/34);  new base . (2/2);  total (34/36)
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.41e-01

TGTGGTTTATCGTCACCTGTCTGGCAATATTGCTGGCGCAAAACATTGAACAATTCACCCTGTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTG  >  NC_000913/883941‑884236
                                                                                                                                                   |                                                                                                                                                    
tgtggtTTATCGTCACCTGTCTGGCAATATTGCTGGCGCAAAACATTGAACAATTCACCCTGTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTg                                                                                                                                                      >  1:245778/1‑148 (MQ=255)
  tggtTTATCGTCACCTGTCTGGCAATATTGCTGGCGCAAAACATTGAACAATTCACCCTGTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTAt                                                                                                                                                   <  1:27436/149‑1 (MQ=255)
         tCGTCACCTGTCTGGCAATATTGCTGGCGCAAAACATTGAACAATTCACCCTGTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGAt                                                                                                                                             >  2:60927/1‑148 (MQ=255)
            tCACCTGTCTGGCAATATTGCTGGCGCAAAACATTGAACAATTCACCCTGTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCAATGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCAcc                                                                                                                                         <  5:181831/149‑1 (MQ=255)
                cTGTCTGGCAATATTGTTGGCGCAAAACATTGAACAATTCACCCTGTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCgcgc                                                                                                                                     >  4:347010/1‑149 (MQ=255)
                   tCTGGCAATATTGCTGGCGCAAAACATTGAACAATTCACCCTGTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGa                                                                                                                                  <  3:86994/149‑1 (MQ=255)
                         aaTATTGCTGGCGCAAAACATTGAACAATTCACCCTGTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCg                                                                                                                             >  2:284330/1‑148 (MQ=255)
                          atatTGCTGGCGCAAAACATTGAACAATTCACCCTGTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGaa                                                                                                                           <  6:234992/149‑1 (MQ=255)
                               gCTGGCGCAAAACATTGAACAATTCACCCTGTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTg                                                                                                                       <  7:326655/148‑1 (MQ=255)
                                 tGGCGCAAAACATTGAACAATTCACCCTGTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGgcg                                                                                                                    >  2:57225/1‑149 (MQ=255)
                                 tGGCGAAAAACAATGAACAATACACCCTGTTGAGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTTGGATACGCCGCAAATCAGGAATACATCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGTCGAACGTGgcg                                                                                                                    >  6:286151/1‑149 (MQ=255)
                                     gcAAAACATTGAACAATTCACCCTGTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGa                                                                                                                                  >  6:177867/1‑131 (MQ=255)
                                            aTTGAACTATTCACCCTGTTGCGCTTCTTGCTGGGCATAAGCCTCTGTTGCATTGGCGCTGTGCGATACGCCGCAATTCAGTAATCCTTCGAAGAGGCGGTTTGTATCAAGCTCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTcc                                                                                                         <  6:254739/149‑1 (MQ=255)
                                            aTTGAACAATTCACCCTGTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTgg                                                                                                                      >  1:348240/1‑136 (MQ=255)
                                            aTTGAACAATTCACCCTGTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTgg                                                                                                                      <  2:348240/136‑1 (MQ=255)
                                              tGAACAATTCACCCTGTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTTTCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGc                                                                                                       >  5:32399/1‑149 (MQ=255)
                                               gAACAATTCACCCTGTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCGGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGc                                                                                                       <  8:119021/148‑1 (MQ=255)
                                                   aaTTCACCCTGTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGCACGTGGCGCGGATTGCTCCGCTACtt                                                                                                  >  2:199300/1‑149 (MQ=255)
                                                   aaTTCACCCTGTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACtt                                                                                                  >  1:75947/1‑149 (MQ=255)
                                                            tGTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCtg                                                                                         >  2:217734/1‑149 (MQ=255)
                                                              ttGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCtggt                                                                                       >  4:123954/1‑149 (MQ=255)
                                                              ttGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCtg                                                                                         <  8:152073/147‑1 (MQ=255)
                                                              ttGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGCGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGGCCGCtggt                                                                                       <  6:354394/149‑1 (MQ=255)
                                                               tGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCtggtg                                                                                      >  6:309241/1‑149 (MQ=255)
                                                                gcgcTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCtggtg                                                                                      <  1:57225/148‑1 (MQ=255)
                                                                   cttcttGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGgcg                                                                                  >  4:252946/1‑149 (MQ=255)
                                                                       ttGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCgcggc                                                                              >  6:239992/1‑149 (MQ=255)
                                                                          cAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGt                                                                            <  1:60927/148‑1 (MQ=255)
                                                                            gggCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGa                                                                         >  2:92422/1‑149 (MQ=255)
                                                                                 tAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCTAt                                                                    >  6:275956/1‑149 (MQ=255)
                                                                                 tAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCAt                                                                    >  5:219282/1‑149 (MQ=255)
                                                                                 tAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCAt                                                                    <  1:109640/149‑1 (MQ=255)
                                                                                   aGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGATGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCAtg                                                                   <  7:140408/148‑1 (MQ=255)
                                                                                   aGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCAtgt                                                                  <  6:219282/149‑1 (MQ=255)
                                                                                    gCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCAtgtg                                                                 <  2:245778/149‑1 (MQ=255)
                                                                                    gCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTT‑TATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGtgc                                                                >  5:174327/1‑149 (MQ=255)
                                                                                    gCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTT‑TATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGtgc                                                                >  5:169079/1‑149 (MQ=255)
                                                                                      ctctGGTTCCTTGGCGCTGGGGGATACGCCGCAATTCAGGACGCCTTCGACGAGGCGGTTTGTAGCAAGCTCACCGCGCGGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGGCCGCTGGTGGGCGCGGCGTGGATCCATGtgct                                                               <  1:199300/149‑1 (MQ=255)
                                                                                       tctGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGtgctg                                                              >  6:357743/1‑149 (MQ=255)
                                                                                         tGTTGCCTTGGCGCTGTGGGATACGCCGCAAGTCAGGAATCCTTGGTCGAGGCGGTTTGTATCAAGATCACCGCGCAGATGGCGAACGTGGCGCTGAATTCTCAGCTACTTCGTCCGCTGGTGGGCGCGGCGGGGATCCATGtgctg                                                              <  2:360718/147‑1 (MQ=255)
                                                                                           tttCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGcc                                                            <  5:239992/147‑1 (MQ=255)
                                                                                           tttCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCt                                                          >  4:349781/1‑149 (MQ=255)
                                                                                               aTTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTgg                                                        >  4:347240/1‑147 (MQ=255)
                                                                                                  ggCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTT‑TATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAggg                                                   <  6:169079/148‑1 (MQ=255)
                                                                                                  ggCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTT‑TATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAggg                                                   <  6:174327/148‑1 (MQ=255)
                                                                                                  ggCGCTGTGGGATACGCCGCAATTCAGGAAGCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAgg                                                    <  1:54648/148‑1 (MQ=255)
                                                                                                       tgtgGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGAtgt                                              >  7:193041/1‑149 (MQ=255)
                                                                                                         tgGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTTTCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGAtgttt                                            <  6:32399/149‑1 (MQ=255)
                                                                                                         tgGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTTATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGAtgttt                                            <  2:25396/149‑1 (MQ=255)
                                                                                                           ggATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTgg                                                        <  3:36688/135‑1 (MQ=255)
                                                                                                           ggATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTgg                                                        >  4:36688/1‑135 (MQ=255)
                                                                                                           ggATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATCTCGAACGTGGCTATGATTGCTCCGCTACTAGGTCCGCTCGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGAtgtttgt                                          >  1:149176/1‑149 (MQ=255)
                                                                                                              tACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGtttgtttt                                       >  1:105575/1‑149 (MQ=255)
                                                                                                                 gccgcAATTCAGTTCTCATTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGAGGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGGGGATCCATGTGCTGCCCTGGGAGGGGAggtttgtt                                         <  5:247679/144‑1 (MQ=255)
                                                                                                                 gccgcAATACAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACTTGGCGCTGATTGCTCCGCTACTAGGTCCGCTGGTGGGCGCGGCGTGGCTCCATGTGCTGCCCTCGGAGCGGAtgtttgtt                                         >  6:247679/1‑144 (MQ=255)
                                                                                                                  ccgcAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAg                                                     <  7:96625/131‑1 (MQ=255)
                                                                                                                  ccgcAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAg                                                     >  8:96625/1‑131 (MQ=255)
                                                                                                                    gcAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGc                                 <  5:357743/149‑1 (MQ=255)
                                                                                                                     cAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGcc                                <  5:309241/149‑1 (MQ=255)
                                                                                                                        ttCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCa                             >  7:56313/1‑149 (MQ=255)
                                                                                                                         tCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGcc                                <  7:204295/145‑1 (MQ=255)
                                                                                                                         tCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGcc                                >  8:204295/1‑145 (MQ=255)
                                                                                                                          cAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCAt                            <  1:92422/148‑1 (MQ=255)
                                                                                                                           aGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTTGTCCGCTGGTGGGCGAGGCGTGGATACATGTGCTGCCCTGGTAGGGGATGTGTGTTGTGTTTGCCGCATTg                          >  6:257288/1‑149 (MQ=255)
                                                                                                                              aaTCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCa                       <  1:348148/149‑1 (MQ=255)
                                                                                                                                    tCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGAtc                 >  4:331979/1‑149 (MQ=255)
                                                                                                                                       aaGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTcc              <  1:284330/149‑1 (MQ=255)
                                                                                                                                           ggCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGA‑GGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCttttt         <  3:349781/149‑1 (MQ=255)
                                                                                                                                           ggCGGTTTGTAGCAAGATCACCGCGCTGATGGCGAACGTGGAGCTGATTGCTCCGGTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCtttt          <  5:286151/149‑1 (MQ=255)
                                                                                                                                               gTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTtgttttgttt                                   >  6:309877/1‑120 (MQ=255)
                                                                                                                                               gTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTtgttttgttt                                   <  5:309877/120‑1 (MQ=255)
                                                                                                                                                   gTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGCTCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTg  >  3:195539/1‑149 (MQ=255)
                                                                                                                                                   gTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTg  >  8:82532/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TGTGGTTTATCGTCACCTGTCTGGCAATATTGCTGGCGCAAAACATTGAACAATTCACCCTGTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTG  >  NC_000913/883941‑884236

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