Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,062,371 A→T 7.3% 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 7.3% 272.0 / 11.4 82N42I (AAT→ATT) yihNMFS transporter family protein
Reads supporting (aligned to +/- strand):  ref base A (40/36);  new base T (2/4);  total (42/40)
Fisher's exact test for biased strand distribution p-value = 4.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.84e-01

TTAAAAATCGCTTAGGGGATGTTATGCTCACGAAAAAGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCGACAGGATTACTGGGTGTTTATTTG  >  NC_000913/4062224‑4062516
                                                                                                                                                   |                                                                                                                                                 
ttAAAAATCGCTTAGGGGATGTTATGCTCACGAAAAAGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAAt                                                                                                                                                  <  4:229144/149‑1 (MQ=255)
 tAAAAATCGCTTAGGGGATGTTATGCTCACGAAAAAGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATg                                                                                                                                                 <  2:496299/149‑1 (MQ=255)
    aaaTCGCTTAGGGGATGTTATGCTCACGAAAAAGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTc                                                                                                                                              >  3:329456/1‑149 (MQ=255)
     aaTCGCTTAGGGGATGTTATGCTCACGAAAAAGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCa                                                                                                                                             >  4:296867/1‑149 (MQ=255)
        cGCTTAGGGGATGTTATGCTCACGAAAAAGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCATTGGTCAAAt                                                                                                                                          <  2:443373/149‑1 (MQ=255)
        cGCTTAGGGGATGTTATGCTCACGAAAAAGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCATTGGTCAAAt                                                                                                                                          <  2:505586/149‑1 (MQ=255)
          cTTAGGGGATGTTATGCTCACGAAAAAGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCATTGGTCAAAtgg                                                                                                                                        <  3:594460/149‑3 (MQ=255)
          cTTAGGGGATGTTATGCTCACGAAAAAGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTg                                                                                                                                        <  4:124226/149‑1 (MQ=255)
           ttAGGGGATGTTATGCTCACGAAAAAGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTgg                                                                                                                                       >  4:37025/1‑149 (MQ=255)
                         gCTCACGAAAAAGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCgg                                                                                                                         <  3:684069/149‑1 (MQ=255)
                           tCACGAAAAAGAAATGGGCGTTATTTAGTCTGTTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTa                                                                                                                       >  1:697030/1‑149 (MQ=255)
                           tCACGAAAAAGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTa                                                                                                                       >  3:78148/1‑149 (MQ=255)
                           tCACGAAAAAGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTa                                                                                                                       >  3:87264/1‑149 (MQ=255)
                           tCACGAAAAAGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTa                                                                                                                       >  1:369975/1‑149 (MQ=255)
                           tCACGAAAAAGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTa                                                                                                                       <  3:323179/149‑1 (MQ=255)
                             aCGAAAAAGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTaaa                                                                                                                     >  2:250434/1‑149 (MQ=255)
                               gAAAAAGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACt                                                                                                                   >  4:96919/1‑149 (MQ=255)
                                 aaaaGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGACAGATGGGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCa                                                                                                                 >  4:506957/1‑149 (MQ=255)
                                 aaaaGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCa                                                                                                                 >  1:223766/1‑149 (MQ=255)
                                  aaaGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCAt                                                                                                                >  2:767424/1‑149 (MQ=255)
                                  aaaGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCAt                                                                                                                >  4:720289/1‑149 (MQ=255)
                                      aaaTGGTCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGt                                                                                                            <  2:581/149‑1 (MQ=255)
                                        aTGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTca                                                                                                          <  1:235423/149‑1 (MQ=255)
                                                  tttAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTggg                                                                                                >  4:189040/1‑149 (MQ=255)
                                                     aGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCtt                                                                                             <  3:478174/149‑1 (MQ=255)
                                                         tATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATgg                                                                                                                                                <  1:313600/94‑1 (MQ=255)
                                                         tATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATgg                                                                                                                                                >  2:313600/1‑94 (MQ=255)
                                                         tATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCtttttt                                                                                         <  3:543700/149‑1 (MQ=255)
                                                          aTTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTc                                                                                                           >  3:591564/1‑130 (MQ=255)
                                                          aTTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTc                                                                                                           <  4:591564/130‑1 (MQ=255)
                                                          aTTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTc                                                                                        >  4:40994/1‑149 (MQ=255)
                                                           ttAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCATTGGTCAAAt                                                                                                                                          >  1:141769/1‑98 (MQ=255)
                                                           ttAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCATTGGTCAAAt                                                                                                                                          <  2:141769/98‑1 (MQ=255)
                                                             aaCACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGt                                                                                     <  4:626921/149‑1 (MQ=255)
                                                              acacTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTc                                                                                    <  2:217625/149‑1 (MQ=255)
                                                              acacTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTc                                                                                    >  1:713581/1‑149 (MQ=255)
                                                                acTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCta                                                                                  <  2:496284/149‑1 (MQ=255)
                                                                 cTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCtat                                                                                 <  3:744391/149‑1 (MQ=255)
                                                                  tgtgtgGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCtatt                                                                                >  1:713586/1‑149 (MQ=255)
                                                                      tgGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAAc                                                                                                                    <  2:361575/109‑1 (MQ=255)
                                                                      tgGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAAc                                                                                                                    >  1:361575/1‑109 (MQ=255)
                                                                          ggTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGc                                                                        <  1:414630/149‑1 (MQ=255)
                                                                          ggTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGc                                                                        <  2:258110/149‑1 (MQ=255)
                                                                                  ttATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAAc                                                                >  2:741530/1‑149 (MQ=255)
                                                                                       aaTTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACcg                                                           >  2:116290/1‑149 (MQ=255)
                                                                                             cGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTAtttt                                                                          <  4:28161/128‑1 (MQ=255)
                                                                                             cGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTAtttt                                                                          >  3:28161/1‑128 (MQ=255)
                                                                                                  cTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATacac                                                <  3:455345/149‑1 (MQ=255)
                                                                                                   tGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACAcc                                               >  4:564186/1‑149 (MQ=255)
                                                                                                         aTGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATATACTACCGCGCAGATACACCATGTCa                                         >  3:286343/1‑149 (MQ=255)
                                                                                                         aTGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCa                                         >  2:641718/1‑149 (MQ=255)
                                                                                                          tGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCAt                                        >  2:146522/1‑149 (MQ=255)
                                                                                                           gCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCAtt                                       >  1:247559/1‑149 (MQ=255)
                                                                                                            cGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATtc                                      >  1:124806/1‑149 (MQ=255)
                                                                                                              ttttATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATtctc                                    >  4:168465/1‑149 (MQ=255)
                                                                                                                  atatCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTctcactc                                <  2:651273/149‑1 (MQ=255)
                                                                                                                     tCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCt                                                                                       <  1:677773/91‑1 (MQ=255)
                                                                                                                     tCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCt                                                                                       >  2:677773/1‑91 (MQ=255)
                                                                                                                     tCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCAtt                             <  1:767424/149‑1 (MQ=255)
                                                                                                                      cccGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTg                            >  1:491474/1‑149 (MQ=255)
                                                                                                                       ccGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGc                           <  1:146522/149‑1 (MQ=255)
                                                                                                                        cGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCg                          >  1:674935/1‑149 (MQ=255)
                                                                                                                           tGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCGACa                       <  3:168465/149‑1 (MQ=255)
                                                                                                                              aGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCGACAGGa                    <  2:246275/149‑1 (MQ=255)
                                                                                                                               ggAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCGACAGGAt                   >  1:50897/1‑149 (MQ=255)
                                                                                                                                gAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCGACAGGAtt                  >  1:231288/1‑149 (MQ=255)
                                                                                                                                     tttCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCGACAGGATTACTgg             <  3:665440/149‑1 (MQ=255)
                                                                                                                                     tttCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCGACAGGATTACTgg             <  1:354866/149‑1 (MQ=255)
                                                                                                                                     tttCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCGACAGGATTACTgg             >  4:109653/1‑149 (MQ=255)
                                                                                                                                      ttCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTctca                                   <  1:174055/126‑1 (MQ=255)
                                                                                                                                      ttCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTctca                                   >  2:174055/1‑126 (MQ=255)
                                                                                                                                       tCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCGACAGGATTACTGGgt           <  1:250434/149‑1 (MQ=255)
                                                                                                                                       tCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCGACAGGATTACTGGgt           <  2:223766/149‑1 (MQ=255)
                                                                                                                                        ccATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCGACAGGATTACTGGgtg          <  4:87264/149‑1 (MQ=255)
                                                                                                                                        caaTTTGACCATTGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCGACAGGATTACTGGgtg          >  3:554237/3‑149 (MQ=255)
                                                                                                                                         cATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCGACAGGATTACTGGgtgt         <  2:247559/149‑1 (MQ=255)
                                                                                                                                          aTTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCGACAGGATTACTGGGTGtt        >  2:677681/1‑149 (MQ=255)
                                                                                                                                          aTTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCGACAGGATTACTGGGTGtt        >  4:77783/1‑149 (MQ=255)
                                                                                                                                          aTTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCGACAGGATTACTGGGTGtt        >  3:228846/1‑149 (MQ=255)
                                                                                                                                               aCCAATGGTCAAATTGGTTATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATATACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCGACAGGATGACTGGGTGTTTAttt   <  4:78148/149‑1 (MQ=255)
                                                                                                                                               aCCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCGACAGGATTACTGGGTGTTTAttt   <  2:713586/149‑1 (MQ=255)
                                                                                                                                                ccAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCGACAGGATTACTGGGTGTTTATTTg  <  4:86098/149‑1 (MQ=255)
                                                                                                                                                  aaTGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATacac                                                >  4:694708/1‑101 (MQ=255)
                                                                                                                                                  aaTGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATacac                                                <  3:694708/101‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TTAAAAATCGCTTAGGGGATGTTATGCTCACGAAAAAGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCGACAGGATTACTGGGTGTTTATTTG  >  NC_000913/4062224‑4062516

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: