Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,062,371 A→T 6.6% N42I (AAT→ATT)  yihN → MFS transporter family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,062,3710AT 6.6% 422.7 / 14.3 122N42I (AAT→ATT) yihNMFS transporter family protein
Reads supporting (aligned to +/- strand):  ref base A (57/56);  new base T (4/4);  total (62/60)
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.92e-01

Only 100 of 126 total aligned reads displayed.

AAAAATCGCTTAGGGGATGTTATGCTCACGAAAAAGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCGACAGGATTACTGGGTGTTTATTTGAC  >  NC_000913/4062226‑4062518
                                                                                                                                                 |                                                                                                                                                   
aaaaaTCGCTTAGGGGATGTTATGCTCACGAAAAAGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATgg                                                                                                                                                  <  2:132877/149‑1 (MQ=255)
  aaaTCGCTTAGGGGATGTTATGCTCACGAAAAAGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTc                                                                                                                                                <  4:1114283/149‑1 (MQ=255)
     tCGCTTAGGGGATGTTATGCTCACGAAAAAGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCaaa                                                                                                                                             <  2:78752/149‑1 (MQ=255)
      cGCTTAGGGGATGTTATGCTCACGAAAAAGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCATTGGTCAAAt                                                                                                                                            <  4:181926/149‑1 (MQ=255)
        cTTAGGGGATGTTATGCTCACGAAAAAGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTg                                                                                                                                          <  4:22780/149‑1 (MQ=255)
         ttAGGGGATGTTATGCTCACGAAAAAGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCATTGGTCAAAtgga                                                                                                                                         <  1:971048/149‑4 (MQ=255)
          tAGGGGATGTTATGCTCACGAAAAAGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGt                                                                                                                                        <  2:201269/149‑1 (MQ=255)
          tAGGGGATGTTATGCTCACGAAAAAGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGt                                                                                                                                        <  4:249271/149‑1 (MQ=255)
              ggATGTTATGCTCACGAAAAAGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATg                                                                                                                                    <  2:993166/149‑1 (MQ=255)
               gATGTTATGCTCACGAAAAAGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGc                                                                                                                                   >  1:438750/1‑149 (MQ=255)
               gATGTTATGCTCACGAAAAAGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGc                                                                                                                                   <  4:547054/149‑1 (MQ=255)
                  gTTATGCTCACGAAAAAGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTAt                                                                                                                                >  3:542172/1‑149 (MQ=255)
                       gCTCACGAAAAAGAAATGGGCGTTATTTAGTTTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATGTTTCCATTTGAACAATGGTCAAATTGGTAATGCTATGTCgg                                                                                                                           >  3:707293/1‑149 (MQ=255)
                        cTCACGAAAAAGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGt                                                                                                                          <  4:542172/149‑1 (MQ=255)
                         tCACGAAAAAGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTa                                                                                                                         >  3:453421/1‑149 (MQ=255)
                         tCACGAAAAAGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTa                                                                                                                         >  2:856133/1‑149 (MQ=255)
                         tCACGAAAAAGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTa                                                                                                                         >  4:294571/1‑149 (MQ=255)
                         tCACGAAAAAGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTa                                                                                                                         >  1:164776/1‑149 (MQ=255)
                               aaaaGATATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCa                                                                                                                   >  2:816677/1‑149 (MQ=255)
                                aaaGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCAt                                                                                                                  <  4:1109681/149‑1 (MQ=255)
                                aaaGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCAt                                                                                                                  >  1:179974/1‑149 (MQ=255)
                                 aaGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGTCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCAtt                                                                                                                 <  1:772490/149‑1 (MQ=255)
                                 aaGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCAtt                                                                                                                 <  1:795569/149‑1 (MQ=255)
                                 aaGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCAtt                                                                                                                 <  2:361683/149‑1 (MQ=255)
                                  aGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCAttt                                                                                                                <  1:140989/149‑1 (MQ=255)
                                      aTGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTca                                                                                                            <  2:79708/149‑1 (MQ=255)
                                        gggCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGt                                                                                                                                        >  4:397801/1‑119 (MQ=255)
                                        gggCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTcacc                                                                                                          >  3:355709/1‑149 (MQ=255)
                                            gTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAg                                                                                                      <  4:862167/149‑1 (MQ=255)
                                             ttATTTAGACTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGt                                                                                                     <  1:856133/149‑1 (MQ=255)
                                               aTTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTgg                                                                                                   <  2:169568/149‑1 (MQ=255)
                                               aTTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTgg                                                                                                   >  3:33430/1‑149 (MQ=255)
                                                   aGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCtt                                                                                               >  2:545122/1‑149 (MQ=255)
                                                    gTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCttt                                                                                              <  2:160839/149‑1 (MQ=255)
                                                       tATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCtttttt                                                                                           >  2:1120619/1‑149 (MQ=255)
                                                         ttAACACTGTGTGGCGGTGCAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCt                                                                                         <  1:1082392/149‑1 (MQ=255)
                                                         ttAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCt                                                                                         >  2:1023366/1‑149 (MQ=255)
                                                           aaCACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGt                                                                                       <  3:294571/149‑1 (MQ=255)
                                                            acacTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTcaccac                                                                                                        >  1:93744/1‑131 (MQ=255)
                                                              acTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAg                                                                                                      <  2:972622/131‑1 (MQ=255)
                                                              acTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAg                                                                                                      >  1:972622/1‑131 (MQ=255)
                                                              acTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCta                                                                                    <  2:23864/149‑1 (MQ=255)
                                                               cTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTc                                                                                                             <  1:1017875/123‑1 (MQ=255)
                                                               cTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCtat                                                                                   <  2:449300/149‑1 (MQ=255)
                                                                         gTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGcc                                                                         <  1:1023366/149‑1 (MQ=255)
                                                                           aCAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGa                                                                       >  2:18858/1‑149 (MQ=255)
                                                                            cAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGAt                                                                      <  1:131969/149‑1 (MQ=255)
                                                                             aaTTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTca                                                                                                            <  3:77423/110‑1 (MQ=255)
                                                                                   taAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTAc                                                               >  2:631354/1‑149 (MQ=255)
                                                                                   taAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTAc                                                               <  1:385704/149‑1 (MQ=255)
                                                                                       ttACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACcgcg                                                           >  4:929040/1‑149 (MQ=255)
                                                                                        tACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCtat                                                                                   <  3:313469/124‑1 (MQ=255)
                                                                                           cGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTAt                                                                                                                                <  2:372600/76‑1 (MQ=255)
                                                                                           cGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTAt                                                                                                                                >  1:372600/1‑76 (MQ=255)
                                                                                                cTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATacac                                                  >  1:612265/1‑149 (MQ=255)
                                                                                                     aGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGt                                             <  4:453421/149‑1 (MQ=255)
                                                                                                        tGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCAt                                          >  2:681133/1‑149 (MQ=255)
                                                                                                         gCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCAtt                                         <  2:164776/149‑1 (MQ=255)
                                                                                                          cGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATtc                                        <  1:545122/149‑1 (MQ=255)
                                                                                                          cGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATtc                                        <  4:741029/149‑1 (MQ=255)
                                                                                                           gctttATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGt                                                                                                     >  1:715387/3‑87 (MQ=255)
                                                                                                           gTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGt                                                                                                     <  2:715387/87‑1 (MQ=255)
                                                                                                           gTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATtct                                       >  4:369324/1‑149 (MQ=255)
                                                                                                           gTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATtct                                       >  3:753225/1‑149 (MQ=255)
                                                                                                               tataTCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTctcact                                   >  1:1070570/1‑149 (MQ=255)
                                                                                                                  atCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCAt                                >  3:1115506/1‑149 (MQ=255)
                                                                                                                   tCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTctca                                     <  4:753043/143‑1 (MQ=255)
                                                                                                                   tCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTctca                                     >  3:753043/1‑143 (MQ=255)
                                                                                                                     ccGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGc                             >  4:114151/1‑149 (MQ=255)
                                                                                                                      cGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCg                            <  1:375469/149‑1 (MQ=255)
                                                                                                                            aGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCtt                                                                                               >  2:1059284/1‑76 (MQ=255)
                                                                                                                            aGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCtt                                                                                               <  1:1059284/76‑1 (MQ=255)
                                                                                                                            aGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCGACAGGa                      >  4:627311/1‑149 (MQ=255)
                                                                                                                             ggAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGc                                                                          <  2:187844/96‑1 (MQ=255)
                                                                                                                             ggAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGc                                                                          >  1:187844/1‑96 (MQ=255)
                                                                                                                             ggAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCGACAGGAt                     >  1:146041/1‑149 (MQ=255)
                                                                                                                              gAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCa                                                                                                                   <  4:96700/54‑1 (MQ=255)
                                                                                                                              gAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCa                                                                                                                   >  3:96700/1‑54 (MQ=255)
                                                                                                                               aaTATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCGACAGGATTa                   <  4:6578/149‑1 (MQ=255)
                                                                                                                                  atTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCtt                                                                                               >  1:803584/1‑70 (MQ=255)
                                                                                                                                   tttCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTAtttt                                                                            >  3:1060181/1‑88 (MQ=255)
                                                                                                                                   tttCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTAtttt                                                                            <  4:1060181/88‑1 (MQ=255)
                                                                                                                                   tttCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCGACAGGATTACTgg               <  4:753225/149‑1 (MQ=255)
                                                                                                                                   tttCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCGACAGGATTACTgg               >  4:1125855/1‑149 (MQ=255)
                                                                                                                                   ttaccATTTGACCAATGGTCAAATTGGTAATGCTCTGTCGGTAAACTCATTTGTCAGCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGTAGATACACCATGTCATTCTCACTCATTGCGACAGGATTACTgg               <  4:707293/146‑1 (MQ=255)
                                                                                                                                       aaTTTGACCATTGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCAt                                               <  4:1077162/112‑1 (MQ=255)
                                                                                                                                       aaTTTGACCATTGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCGACAGGATTACTGGgtgt           <  3:512696/148‑1 (MQ=255)
                                                                                                                                       aaTTTGACCATTGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCGACAGGATTACTGGgtga           >  3:179385/2‑148 (MQ=255)
                                                                                                                                        aTTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCGACAGGATTACTGGGTGtt          >  3:170405/1‑149 (MQ=255)
                                                                                                                                        aTTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCGACAGGATTACTGGGTGtt          >  2:926465/1‑149 (MQ=255)
                                                                                                                                         tttGACCAGTGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCGACAGGATTACTGGGTGttt         >  2:407272/1‑149 (MQ=255)
                                                                                                                                         tttGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCGACAGGATTACTGGGTGttt         >  2:993070/1‑149 (MQ=255)
                                                                                                                                               cAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCGACa                         <  1:705668/127‑1 (MQ=255)
                                                                                                                                               cAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCGACAGGATTACTGGGTGTTTATTTga   <  2:916031/149‑1 (MQ=255)
                                                                                                                                               cAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCGACAGGATTACTGGGTGTTTATTTga   <  2:612265/149‑1 (MQ=255)
                                                                                                                                               cAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCGACAGGATTACTGGGTGTTTATTTga   >  2:120034/1‑149 (MQ=255)
                                                                                                                                                aaTGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCGACAGGATTACTGGGTGTTTATTTgac  >  3:128539/1‑149 (MQ=255)
                                                                                                                                                aaTGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCGACAGGATTACTGGGTGTTTATTTgac  >  3:589133/1‑149 (MQ=255)
                                                                                                                                                aaTGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCGACAGGATTACTGGGTGTTTATTTgac  <  4:33430/149‑1 (MQ=255)
                                                                                                                                                aaTGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCGACAGGATTACTGGGTGTTTATTTgac  <  1:263504/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
AAAAATCGCTTAGGGGATGTTATGCTCACGAAAAAGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCGACAGGATTACTGGGTGTTTATTTGAC  >  NC_000913/4062226‑4062518

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