Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,515,527 C→T 5.8% L256L (CTA→TTA)  ydcV → putative ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,515,5270CT 5.8% 220.8 / 2.6 69L256L (CTA→TTA) ydcVputative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (34/31);  new base T (2/2);  total (36/33)
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.37e-01

TTACGACCTTTACGGCAGGTCATGAACGAACGTTACCGTTGTGGTTGCTCAATCAGCTTGGGCGACCGCGTGATGTACCGGTAACTAACGTGGTGGCACTGCTGGTTATGTTGGTAACAACCTTGCCGATCCTGGGGGCCTGGTGGCTAACCCGCGAAGGCGACAATGGTCAATAACCACTGATACAGGAATATGCTATGCAACATAAGTTACTGATTAACGGAGAACTGGTTAGCGGCGAAGGGGAAAAACAGCCTGTCTATAATCCGGCAACGGGGGACGTTTTACTGGA  >  NC_000913/1515381‑1515672
                                                                                                                                                  |                                                                                                                                                 
ttACGACCTTTACGGCAGGTCATGAACGAACGTTACCGTTGTGGTTGCTCAATCAGCTTGGGCGACCGCGTGATGTCCCGGTAACTAACGTGGTGGCACTGCTGGTTATGTTGGTAACACCCTTGCCGATCCTGGGGGCCTGGTGGCt                                                                                                                                                  >  6:204148/1‑148 (MQ=255)
          tACGGCAGGTCATGAACGAACGTTACCGTTGTGGTTGCTCAATCAGCTTGGGCGACCGCGTGATGTACCGGTAACTAACGTGGTGGCACTGCTGGTTATGTTGGTAACAACCTTGCCGATCCTGGGGGCCTGGTGGCTAACCCGCGaa                                                                                                                                        <  1:14993/148‑1 (MQ=255)
             ggCAGGTCATGAACGAACGTTACCGTTGTGGTTGCTCAATCAGCTTGGGCGACCGCGTGATGTACCGGTAACTAACGTGGTGGCACTGCTGGTTATGTTGGTAACAACCTTGCCGATCCTGGGGGCCTGGTGGCTAACCCGCGAAg                                                                                                                                       <  2:222411/146‑1 (MQ=255)
             ggCAGGTCATGAACGAACGTTACCGTTGTGGTTGCTCAATCACCTTGGGCGACCGCGTGATGTACCGGTAACTAACGTGGTGGCACTGCTGGTTATGTTGGTAACAACCTTGCCGATCCTGGGGGCCTGGTGGCTAACCCGCGAAg                                                                                                                                       >  1:222411/1‑146 (MQ=255)
                   tCATGAACGAACGTTACCGTTGTGGTTGCTCAATCAGCTTGGGCGACCGCGTGATGTACCGGTAACTAACGTGGTGGCACTGCTGGTTATGTTGGTAACAACCTTGCCGATCCTGGGGGCCTGGTGGCTAACCCGCGAAGGCGACAATg                                                                                                                              <  4:290852/149‑1 (MQ=255)
                     aTGAACGAACGTTACCGTTGTGGTTGCTCAATCAGCTTGGGCGACCGCGTGATGTACCGGTAACTAACGTGGGGGCACTGCTGGTTATGTTGGTAACAACCTTGCCGATCCTGGGGGCCTGGTGGCTAACCCGCGAAGGCGACAATGGt                                                                                                                            <  1:286989/149‑1 (MQ=255)
                        aacgaacgTTACCGTTGTGGTTGCTCAATCAGCTTGGGCGACCGCGTGATGTACCGGTAACTAACGTGGTGGCACTGCTGGTTATGTTGGTAACAACCTTGCCGATCCTGGGGGCCTGGTGGCTAACCCGCGAAGGCGACAATGGTCa                                                                                                                          <  8:292728/148‑1 (MQ=255)
                         acgaacgTTACCGTTGTGGTTGCTCAATCAGCTTGGGCGACCGCGTGATGTACCGGTAACTAACGTGGTGGCACTGCTGGTTATGTTGGTAACAACCTTGCCGATCCTGGGGGCCTGGTGGCTAACCCGCGAAGGCGACAATGGTCAAt                                                                                                                        <  8:205116/149‑1 (MQ=255)
                                     gTTGTGGTTGCTCAATCAGCTTGGGCGACCGCGTGAGGTACCGGTAACTAACGTGGTGGCACTGCTGGTGATGTTGGTAACAACCTTGCCGATCCTGGGGGCCTGGTGGCTAACCCGCGAAGGCGCCAATGGTCAATAACCTCTGATa                                                                                                             >  6:205586/1‑148 (MQ=255)
                                      ttGTGGTTGCTCAATCAGCTTGGGCGACCGCGTGATGTACCGGTAACTAACGTGGTGGCACTGCTGGTTATGTTGGTAACAACCTTGCCGATCCTGGGGGCCTGGTGGCTAACCCGCGAAGGCGACAATGGTCAATAACCACTGATa                                                                                                             >  8:75141/1‑147 (MQ=255)
                                            ttGCTCAATCAGCTTGGGCGACCGCGTGATGTACCGGTAACTAACGTGGTGGCACTGCTGGTTATGTTGGTAACAACCTTGCCGATCCTGGGGGCCTGGTGGCTAACCCGCGAAGGCGACAATGGTCAATAACCACTGATACAGGAata                                                                                                     >  6:256895/1‑149 (MQ=255)
                                               cTCAATCAGCTTGGGCGACCGCGTGATGTACCGGTAACTAACGTGGTGGCACTGCTGGTTATGTTGGTAACACCCTTGCCGATCCTGGGGGCCTGGTGGCTAACCCGCGAAGGCGACAATGGTCAATAACCACTGATACAGGAAtatgc                                                                                                  >  2:360209/1‑149 (MQ=255)
                                               cTCAATCAGCTTGGGCGACCGCGTGATGTACCGGTAACTAACGTGGTGGCACTGCTGGTTATGTTGGTAACAACCTTGCCGATCCTGGGGGCCTGGTGGCTAACCCGCGAAGGCGACAATGGTCAATAACCACTGATACAGGAAtatg                                                                                                   <  4:271599/148‑1 (MQ=255)
                                               cTCAATCAGCTTGGGCGACCGCGTGATGTACAGGTAACTAACGTGGTGGCACTGCTGGTTATGTTGGTAACAACCTTGCCGATCCTGGGGGCCTGGTGGCTAACCCGCGAAGGCGACAATGGTCAAGAACCACTGATACAGGAAtatgc                                                                                                  >  7:336007/1‑149 (MQ=255)
                                                tCAATCAGCTTGGGCGACCGCGTGATGTACCGGTAACTAACGTGGTGGCACTGCTGGTTATGTTGGTAACAACCTTGCCGATCCTGGGGGCCTGGTGGCTAACCCGCGAAGGCGACAATGGTCAATAACCACTGATACAGGAAtatgc                                                                                                  >  2:323425/1‑148 (MQ=255)
                                                  aaTCAGCTTGGGCGACCGCGTGATGTACCGGTAACTAACGTGGTGGCACTGCTGGTTATGTTGGTAACAACCTTGCCGATCCTGGGGGCCTGGTGGTTAACCCGCGAAGGCGACAATGGTCAATAACCACTGATACAGGAAtatgctat                                                                                               >  4:296270/1‑149 (MQ=255)
                                                  aaTCAGCTTGGGCGACCGCGTGATGTACCGGTAACTAACGTGGTGGCACTGCTGGTTATGTTGGTAACAACCTTGCCGATCCTGGGGGCCTGGTGGTTAACCCGCGAAGGCGACAATGGTCAATAACCACTGATACAGGAAtatgctat                                                                                               >  4:291344/1‑149 (MQ=255)
                                                  aaTCAGCTTGGGCGACCGCGTGATGTACCGGTAACTAACGTGGTGGCACTGCTGGTTATGGTGGTAACAACCATGCCGATCCTGGGGGCCTGGTGGCTAACCCGCGAAGGCGACAATGGTCAATAACCACTGATACCGCAAtatgctat                                                                                               >  2:254254/1‑149 (MQ=255)
                                                    tCAACTTGGGCGACCGCGTGATGTACCGGTAACTAACGTGGTGGCACTGCTGGTTATGTTGGTAACAACCTTGCCGATCCTGGGGGCCTGGTGGCTAACCCGCGAAGGCGACAATGGTCAATAACCACTGATACAGGAAtatgctatg                                                                                              >  3:192580/1‑148 (MQ=255)
                                                       gCTTGGGCGACCGCGTGATGTACCGGTAACTAACGTGGTGGCACTGCTGGTTATGTTGGTAACAACCTTGCCGATCCTGGGGGCCTGGTGGCTAACCCGCGAAGGCGACAATGGTCAATAACCACTGATACAGGAATATGCTATGCaa                                                                                           >  6:79211/1‑148 (MQ=255)
                                                        cTTGGGCGACCGCGTGATGTACCGGTAACTAACGTGGTGGCACTGCTGGTTATGTTGGTAACAACCTTGCCGATCCTGGGGGCCTGGTGGCTAACCCGCGAAGGCGACAATGGTCAATAACCACTGATACAGGAATATGCTATGCAACa                                                                                         <  5:356233/149‑1 (MQ=255)
                                                               gACCGCGTGATGTACCGGTAACTAACGTGGTGGCGCTGCTGGTTATGTTGGTAACAACCTTGCCGATCCTGGGGGCCTGGTGGCTAACCCGCGAAGGCGACAATGGTCAATAACCACTGATACAGGAATATGCTATGCAACATAAGTTa                                                                                  >  8:211095/1‑149 (MQ=255)
                                                               gACCGCGTGATGTACCGGTAACTAACGTGGTGGCACTGCTGGTTATGTTGGTAACAACCTTGCCGATCCTGGGGGCCTGGTGGCTAACCCGCGAAGGCGACaa                                                                                                                                >  6:218735/1‑103 (MQ=255)
                                                               gACCGCGTGATGTACCGGTAACTAACGTGGTGGCACTGCTGGTTATGTTGGTAACAACCTTGCCGATCCTGGGGGCCTGGTGGCTAACCCGCGAAGGCGACaa                                                                                                                                <  5:218735/103‑1 (MQ=255)
                                                                  cgcgTGATGTACCGGTAACTAACGTGGTGGCACTGCTGGTTATGTTGGTAACAACCTTGCCGATCCTGGGGGCCTGGTGGCTAACCCGCGAAGGCGACAATGGTCAATAACCACTGATACAGGAATATGCTATGCAACATAAGTTACTg                                                                               <  4:316791/149‑1 (MQ=255)
                                                                   gcgTGATGTACCGGTAACTAACGTGGTGGCGCTGCTGGTTATGTTGGTAACAACCTTGCCGATCCTGGGGGCCTGGTGGCTAACCCGCGAAGGCGACAATGGTCAATAACCACTGATACAGGAATATGCTATGCAACATAAGTTACTGa                                                                              <  7:211095/149‑1 (MQ=255)
                                                                   gcgTGATGTACCGGTAACTAACGTGGTGGCACTGCTGGTTATGTTGGTAACAACCTTGCCGATCCTGGGGGCCTGGTGGCTAACCCGCGAAGGCGACAATGGTCAATAACCACTGATACAGGAATATGCTATGCAACATAAGTTACTGa                                                                              <  4:192580/149‑1 (MQ=255)
                                                                   gcgTGATGTACCGGTAACTAACGTGGTGGCACTGCTGGTTATGTTGGTAACAACCTTGCCGATCCTGGGGGCCTGGTGGCTAACCCGCGAAGGCGACAATGGTCAATAACCACTGATACAGGAATATGCTATGCAACATAAGTTACTGa                                                                              <  7:75141/149‑1 (MQ=255)
                                                                     gTGATGTACCGGTAACTAACGTGGTGGCACTGCTGGTTATGTTGGTAACAACCTTGCCGATCCTGGGGGCCTGGTGGCTAACCCGCGAAGGCGACAATGGTCAATAACCACTGATACAGGAATATGCTATGCAACATAAGTTACTGAtg                                                                            <  5:256895/149‑2 (MQ=255)
                                                                        aTGTACCGGTAACTACCGTGGTGGCACTGCTGGTTATGTTGGTAACAACCTTGCCGATCCTGGGGGCCTGGTGGCTAACCCGCGAAGGCGACAATGGGCAATAACCACTGATACAGGAATATGCTATGCAACATAAGTTACTGATTAAc                                                                         <  5:73692/149‑1 (MQ=255)
                                                                        aTGTACCGGTAACTAACGTGGTGGCACTGCTGGTTATGTTGGTAACAACCTTGCCGATCCTGGGGGCCTGGTGGCTAACCCGCGAAGGCGACAATGGTCAATAACCACTGATACAGGAATATGCTATGCAACATAAGTTACTGATTAAc                                                                         <  1:323425/149‑1 (MQ=255)
                                                                                     taacGTGGTGGCACTGCTGGTTATGTTGGTAACAACCTTGCCGATCCTGGGGGCCTGGTGGCTAACCCGCGAAGGCGACAATGGTCAATAACCACTGATACAGGAATATGCTATGCAACATAAGTTACTGATTAACGGAGAACTGGt                                                              >  4:227697/1‑147 (MQ=255)
                                                                                       acGTGGTGGCACTGCTGGTTATGTTGGTAACAACCTTGCCGATCCTGGGGGCCTGGTGGCTAACCCGCGAAGGCGACAATGGTCAATAACCACTGATACAGGAATATGCTATGCAACATAAGt                                                                                    >  4:312849/1‑123 (MQ=255)
                                                                                       acGTGGTGGCACTGCTGGTTATGTTGGTAACAACCTTGCCGATCCTGGGGGCCTGGTGGCTAACCCGCGAAGGCGACAATGGTCAATAACCACTGATACAGGAATATGCTATGCAACATAAGTTACTGATTAACGGAGAACTGGTTAgc                                                          <  6:202588/149‑1 (MQ=255)
                                                                                         gtggtgGCACTGCTGGTTATGTTGGTAACAACCTTGCCGATCCTGGGGGCCTGGTGGCTAACCCGCGAAGGCGACAATGGTCAATAACCACTGATACAGGAATATGCTATGCAACATAAGTTACTGATTAACGGAGAACTGGTTAgcgg                                                        >  5:66519/1‑149 (MQ=255)
                                                                                             tggCACTGCTGGTTATGTTGGTAACAACCTTGCCGATCCTGGGGGCCTGGTGGCTAACCCGCGAAGGCGACAATGGTCAATAACCACTGATACAGGAATATGCTATGCAACATAAGTTACTGATTAACGGAGAACTGGTTAGCGGCGaa                                                    <  3:227697/149‑1 (MQ=255)
                                                                                               gCACTGCTGGTTATGTTGGTAACAACCTTGCCGATCCTGGGGGCCTGGTGGTTAACCCGCGAAGGCGACAATGGTCAATAACCACTGATACAGGAATATGCTATGCAACATAAGTTACTGATTAACGGAGAACTGGTTAGCGGCGAAgg                                                  <  3:291344/149‑1 (MQ=255)
                                                                                               gCACTGCTGGTTATGTTGGTAACAACCTTGCCGATCCTGGGGGCCTGGTGGTTAACCCGCGAAGGCGACAATGGTCAATAACCACTGATACAGGAATATGCTATGCAACATAAGTTACTGATTAACGGAGAACTGGTTAGCGGCGAAgg                                                  <  3:296270/149‑1 (MQ=255)
                                                                                                cACTGCTGGTTATGTTGGTAACAACCTTGCCGATCCTGGGGGCCTGGTGGCTAACCCGCGAAGGCGACAATGGTCAATAACCACTGATACAGGAATATGCTATGCAACATAAGTTACTGATTAACGGAGAACTGGTTAGCGGCGAAggg                                                 <  5:79211/149‑1 (MQ=255)
                                                                                                      tgGTTATGTTGGTAACAACCTTGCCGATCCTGGGGGCCTGGTGGCTAACCCGCGAAGGCGACAATGGTCAATAACCACTGATACAGGAATATGCTATGCAACATAAGTTACTGATTAACGGAGAACTGGTTAGCGGCGAAGGGGaaaaa                                           >  8:115001/1‑149 (MQ=255)
                                                                                                        gTTATGTTGGTAACAACCTTGCCGATCCTGGGGGCCTGGTGGCTAACCCGCGAAGGCGACAATGGTCAATAACCACTGATACAGGAATATGCTATGCAACATAAGTTACTGATTAACCGAGAACTGGTTAGCGGCGAAGGGGaaaaa                                           >  2:346116/1‑147 (MQ=255)
                                                                                                         ttATGTTGGTAACAACCTTGCCGATCCTGGGGGCCTGGTGGCTAACCCGCGAAGGCGACAATGGTCAAAAACCACTGATACAGGAATATGCTATGCAACATAAGTTACTGATTTACGGAGAACTGTTTAGCGGCGAAGGGGAAAAACAg                                        >  6:216148/1‑149 (MQ=255)
                                                                                                          tATGTTGGTAACAACCTTGCCGATCCTGGGGGCCTGGTGGCTAACCCGCGAAGGCGACAATGGTCAATAACCACTGATACAGGAATATGCTATGCAACATAAGTTACTGATTAACGGAGAACTGGTTAGCGGCGAAGGGGAAAAACAGc                                       >  5:323729/1‑149 (MQ=255)
                                                                                                             gTTGGTAACAACCTTGCCGATCCTGGGGGCCTGGTGGCTAACCCGCGAAGGCGACAATGGTCAATAACCACTGATACAGGAATATGCTATGCAACATAAGTTACTGATTAACGGAGAACTGGTTAGCGGCGaa                                                    <  7:255171/133‑1 (MQ=255)
                                                                                                             gTTGGTAACAACCTTGCCGATCCTGGGGGCCTGGTGGCTAACCCGCGAAGGCGACAATGGTCAATAACCACTGATACAGGAATATGCTATGCAACATAAGTTACTGATTAACGGAGAACTGGTTAGCGGCGaa                                                    >  8:255171/1‑133 (MQ=255)
                                                                                                              ttGGTAACAACCTTGCCGATCCTGGGGGCCTGGTGGCTAACCCGCGAAGGCGACAATGGTCAATAACCAc                                                                                                                  >  5:117288/1‑70 (MQ=255)
                                                                                                              ttGGTAACAACCTTGCCGATCCTGGGGGCCTGGTGGCTAACCCGCGAAGGCGACAATGGTCAATAACCAc                                                                                                                  <  6:117288/70‑1 (MQ=255)
                                                                                                              ttGGTAACAACCTTGCCGATCCTGGGGGCCTGGTGGCTAACCCGCGAAGGCGACAATGGTCAATAACCACTGATACAGGAATATGCTATGCAACATAAGTTACt                                                                                <  3:285012/104‑1 (MQ=255)
                                                                                                              ttGGTAACAACCTTGCCGATCCTGGGGGCCTGGTGGCTAACCCGCGAAGGCGACAATGGTCAATAACCACTGATACAGGAATATGCTATGCAACATAAGTTACt                                                                                >  4:285012/1‑104 (MQ=255)
                                                                                                              ttGGTAACAACCTTGCCGATCCTGGGGGCCTGGTGGCTAACCCGCGAAGGCGACAATGGTCAATAACCACTGATACAGGAATATGCTATGCAACATAAGTTACTGATTAACGGAGAACTGGTTAGCGGCGAAGGGGAAAAACAGCCTGt                                   >  1:355391/1‑149 (MQ=255)
                                                                                                                        ccTTGCCGATCCTGGGGGCCTGGTGGCTAACCCGCGAAGGCGACAATGGTCAATAACCACTGATACAGGAATATGCTATGCAACATAAGTTACTGATTAACGGAGAACTGGTTAGCGGCGAAGGGGAAAAACAGCCTGTCTATAATCCg                         >  8:223951/1‑149 (MQ=255)
                                                                                                                            gCCGATCCTGGGGGCCTGGTGGCTAACCCGCGAAGGCGACAATGGTCAATAACCACTGATACAGGAATATGCTATGCAACATAAGTTACTGATTAACGGAGAACTGGTTAGCGGCGAAGGGGAAAAACAGCCTGTCTATAATc                           <  4:188156/143‑1 (MQ=255)
                                                                                                                            gCCGATCCTGGGGGCCTGGTGGCTAACCCGCGAAGGCGACAATGGTCAATAACCACTGATACAGGAATATGCTATGCAACATAAGTTACTGATTAACGGAGAACTGGTTAGCGGCGAAGGGGAAAAACAGCCTGTCTATAATc                           >  3:188156/1‑143 (MQ=255)
                                                                                                                              cGATCCTGGGGGCCTGGTGGCTAACCCGCGAAGGCGACAATGGTCAATAACCACTGATACAGGAATATGCTATGCAACATAAGTTACTGATTAACGGAGAACTGGTTAGCGGCGAAGGGGAAAAACAGCCTGTCTATAATCCgg                        <  2:49290/144‑1 (MQ=255)
                                                                                                                              cGATCCTGGGGGCCTGGTGGCTAACCCGCGAAGGCGACAATGGTCAATAACCACTGATACAGGAATATGCTATGCAACATAAGTTACTGATTAACGGAGAACTGGTTAGCGGCGAAGGGGAAAAACAGCCTGTCTATAATCCgg                        >  1:49290/1‑144 (MQ=255)
                                                                                                                                aTCCTGGGGGCCTGGTGGCTAACCCGCGAAGGCGACAATGGTCAATAACCACTGATACAGGAATATGCTATGCAACATAAGTTACTGATTAACGGAGAACTGGTTAGCGGCGAAGGGGAAAAACAGCCTGTCTATAATCCGGCAACggg                 >  5:91804/1‑149 (MQ=255)
                                                                                                                                   cTGGGGGCCTGGTGGCTAACCCGCGAAGGCGACAATGGTCAATAACCACTGATACAGGAATATGCTATGCAACATAAGTTACTGATTAACGGAGAACTGGTTAGCGGCGAAGGGGAAAAACAGCCTGTCTATAATCCGGCAACGGGGGa              <  6:323729/149‑1 (MQ=255)
                                                                                                                                    tGGGGGCCTGGTGGCTAACCCGCGAAGGCGACAATGGTCAATAACCACTGATACAGGAATATGCTATGCAACATAAGTTACTGATTAACGGAGAACTGGTTAGCGGCGAAGGGGAAAAACAGCCTGTCTATaa                             <  1:347547/133‑1 (MQ=255)
                                                                                                                                    tGGGGGCCTGGTGGCTAACCCGCGAAGGCGACAATGGTCAATAACCACTGATACAGGAATATGCTATGCAACATAAGTTACTGATTAACGGAGAACTGGTTAGCGGCGAAGGGGAAAAACAGCCTGTCTATaa                             >  2:347547/1‑133 (MQ=255)
                                                                                                                                        ggCCTGGTGGCTAACCCGCGAAGGCGACAATGGTCAATAACCACTGATACAGGAATATGCTATGCAACATAAGTTACTGATTAACGGAGAACTGGTTAGCGGCGAAGGGGAAAAACAGCCTGTCTATAATCCGGCAACGGGGGACGttt         <  6:66519/149‑1 (MQ=255)
                                                                                                                                            tggtggCTAACCCGCGAAGGCGACAATGGTCAATAACCACTGATACAGGAATATGCTATGCAACATAAGTTACTGATTAACGga                                                                      >  2:269767/1‑84 (MQ=255)
                                                                                                                                            tggtggCTAACCCGCGAAGGCGACAATGGTCAATAACCACTGATACAGGAATATGCTATGCAACATAAGTTACTGATTAACGga                                                                      <  1:269767/84‑1 (MQ=255)
                                                                                                                                            tggtggCTAACCCGCGAAGGCGACAATGGTCAATAACCACTGATACAGGAATATGCTATGCAACATAAGTTACTGATTAACGGAGAACTGGTTAGCGGCGAAGGGGAAAAACAGCCTGTCTATAATCCGGCAACGGGGGACGTTTTACt     >  1:31551/1‑149 (MQ=255)
                                                                                                                                            tggtggCTAACCCGCGAAGGCGACAATGGTCAATAACCACTGATACAGGAATATGCTATGCAACATAAGTTACTGATTAACGGAGAACTGGTTAGCGGCGAAGGGGAAAAACAGCCTGTCTATAATCCGGCAACGGGGGACGTTTTACt     >  5:229956/1‑149 (MQ=255)
                                                                                                                                             ggtggCTAACCCGCGAAGGCGACAATGGTCAATAACCACTGATACAGGAATATGCTATGCAACATAAGTTACTGATTAACGGAGAACTGGTTAGCGGCGAAGGGGAAAAACAGCCTGTCTATAATCCGGCAACGGGGGACGTTTTACTg    >  1:244792/1‑149 (MQ=255)
                                                                                                                                             ggtggCTAACCCGCGAAGGCGACAATGGTCAATAACCACTGATACAGGAATATGCTATGCAACATAAGTTACTGATTAACGGAGAACTGGTTAGCGGCGAAGGGGAAAAACAGCCTGTCTATAATCCGGCAACGGGGGACGTTTTACTg    >  8:79520/1‑149 (MQ=255)
                                                                                                                                              gtggCTAACCCGCGAAGGCGACAATGGTCAATAACCACTGATACAGGAATATGCTATGCAACATAAGTTACTGATTAACGGAGAACTGGTTAGCGGCGAAGGGGAAAAACAGCCTGTCTATAATCCGGCAACGGGGGACGTTTTACTgg   <  2:31551/149‑1 (MQ=255)
                                                                                                                                              gtggCTAACCCGCGAAGGCGACAATGGTCAATAACCACTGATACAGGAATATGCTATGCAACATAAGTTACTGATTAACGGAGAACTGGTTAGCGGCGAAGGGGAAAAACAGCCTGTCTATAATCCGGCAACGGGGGACGTTTTACTg    <  6:229956/148‑1 (MQ=255)
                                                                                                                                               tggCTAACCCGCGAACGCGACAATGGTCAATAACCACTGATACAGGAATATGCTATGCAACATAAGTTACTGATTAACGGAGAACTGGTTAGCGGCGAAGGGGAAAAACAGCCTGTCTATAATCCGGCAACGGGGGACGTTTTACTGGa  <  5:204148/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TTACGACCTTTACGGCAGGTCATGAACGAACGTTACCGTTGTGGTTGCTCAATCAGCTTGGGCGACCGCGTGATGTACCGGTAACTAACGTGGTGGCACTGCTGGTTATGTTGGTAACAACCTTGCCGATCCTGGGGGCCTGGTGGCTAACCCGCGAAGGCGACAATGGTCAATAACCACTGATACAGGAATATGCTATGCAACATAAGTTACTGATTAACGGAGAACTGGTTAGCGGCGAAGGGGAAAAACAGCCTGTCTATAATCCGGCAACGGGGGACGTTTTACTGGA  >  NC_000913/1515381‑1515672

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: ‑