Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,062,371 A→T 5.9% 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 5.9% 346.6 / 10.7 103N42I (AAT→ATT) yihNMFS transporter family protein
Reads supporting (aligned to +/- strand):  ref base A (50/46);  new base T (2/4);  total (52/51)
Fisher's exact test for biased strand distribution p-value = 4.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.32e-01

Only 100 of 105 total aligned reads displayed.

AAAAATCGCTTAGGGGATGTTATGCTCACGAAAAAGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCGACAGGATTACTGGGTGTTTATTT  >  NC_000913/4062226‑4062515
                                                                                                                                                 |                                                                                                                                                
aaaaaTCGCTTAGGGGATGTTATGCTCACGAAAAAGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCATTgg                                                                                                                                               >  1:757422/1‑149 (MQ=255)
 aaaaTCGCTTAGGGGATGTTATGCTCACGAAAAAGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGt                                                                                                                                              <  3:858550/149‑1 (MQ=255)
 aaaaTCGCTTAGGGGATGTTATGCTCACGAAAAAGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGt                                                                                                                                              <  1:818213/149‑1 (MQ=255)
  aaaTCGCTTAGGGGATGTTATGCTCACGAAAAAGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTc                                                                                                                                             <  4:645756/149‑1 (MQ=255)
      cGCTTAGGGGATGTTATGCTCACGAAAAAGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCATTGGTCAAAt                                                                                                                                         <  1:128069/149‑1 (MQ=255)
      cGCTTAGGGGATGTTATGCTCACGAAAAAGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCATTGGTCAAAt                                                                                                                                         <  3:682193/149‑1 (MQ=255)
      cGCTTAGGGGATGTTATGCTCACGAAAAAGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCATTGGTCAAAt                                                                                                                                         <  2:638876/149‑1 (MQ=255)
       gCTTAGGGGATGTTATGCTCACGAAAAAGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAAtt                                                                                                                                        >  3:637980/1‑149 (MQ=255)
       gCTTAGGGGATGTTATGCTCACGAAAAAGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGAACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAAtt                                                                                                                                        <  2:673933/149‑1 (MQ=255)
        cTTAGGGGATGTTATGCTCACGAAAAAGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCATTGGTCAAAtgg                                                                                                                                       <  2:757422/149‑3 (MQ=255)
        cTTAGGGGATGTTATGCTCACGAAAAAGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTg                                                                                                                                       <  3:64031/149‑1 (MQ=255)
         ttAGGGGATGTTATGCTCACGAAAAAGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTgg                                                                                                                                      <  2:544924/149‑1 (MQ=255)
           aGGGGATGTTATGCTCACGAAAAAGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTa                                                                                                                                    >  4:227914/1‑149 (MQ=255)
           aGGGGATGTTATGCTCACGAAAAAGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTa                                                                                                                                    >  2:587093/1‑149 (MQ=255)
            ggggATGTTATGCTCACGAAAAAGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTaa                                                                                                                                   >  2:122360/1‑149 (MQ=255)
             gggATGTTATGCTCACGAAAAAGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAAt                                                                                                                                  >  1:504083/1‑149 (MQ=255)
             gggATGTTATGCTCACGAAAAAGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGATACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAAt                                                                                                                                  >  3:500771/1‑149 (MQ=255)
               gATGTTATGCTCACGAAAAAGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGc                                                                                                                                >  1:609150/1‑149 (MQ=255)
                  gTTATGCTCACGAAAAAGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTAt                                                                                                                             <  4:715630/149‑1 (MQ=255)
                   ttATGCTCACGAAAAAGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATg                                                                                                                            >  4:193615/1‑149 (MQ=255)
                   ttATGCTCACGAAAAAGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATg                                                                                                                            >  2:597100/1‑149 (MQ=255)
                         tCACGAAAAAGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTa                                                                                                                      >  1:713485/1‑149 (MQ=255)
                           aCGAAAAAGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTata                                                                                                                    <  2:698175/149‑3 (MQ=255)
                              aaaaaGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTc                                                                                                                 >  3:36380/1‑149 (MQ=255)
                              aaaaaGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTc                                                                                                                 >  3:237013/1‑149 (MQ=255)
                                aaaGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCAt                                                                                                               <  1:597100/149‑1 (MQ=255)
                                aaaGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCAt                                                                                                               >  3:44570/1‑149 (MQ=255)
                                 aaGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCAtt                                                                                                              <  2:172038/149‑1 (MQ=255)
                                      aTGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGc                                                                                                                                <  3:279027/126‑1 (MQ=255)
                                      aTGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGc                                                                                                                                >  4:279027/1‑126 (MQ=255)
                                      aTGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTca                                                                                                         >  4:767509/1‑149 (MQ=255)
                                      aTGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTca                                                                                                         >  4:31530/1‑149 (MQ=255)
                                       tGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTcac                                                                                                        <  4:117385/149‑1 (MQ=255)
                                         ggCGTTATTTAGCCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTcacca                                                                                                      >  2:685114/1‑149 (MQ=255)
                                          gagTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGTTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTcaccac                                                                                                     >  4:725442/3‑149 (MQ=255)
                                               aTTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTgg                                                                                                <  3:235540/149‑1 (MQ=255)
                                                tttAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTggg                                                                                               >  3:300737/1‑149 (MQ=255)
                                                 ttAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGc                                                                                              <  2:792824/149‑1 (MQ=255)
                                                     tCTATTAACACTGTGTGGCGGTACAATTTATAACTTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCACTGCTCAAATTGCTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCtttt                                                                                          <  2:61029/149‑1 (MQ=255)
                                                       tATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCtttttt                                                                                        >  4:733436/1‑149 (MQ=255)
                                                       tATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCtttttt                                                                                        >  2:418644/1‑149 (MQ=255)
                                                         ttAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCt                                                                                      >  4:110168/1‑149 (MQ=255)
                                                         ttAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCt                                                                                      <  1:644308/149‑1 (MQ=255)
                                                          tAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTg                                                                                     >  1:163555/1‑149 (MQ=255)
                                                           aaCACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGt                                                                                    >  2:855495/1‑149 (MQ=255)
                                                           aaCACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGt                                                                                    >  4:727365/1‑149 (MQ=255)
                                                            acacTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTc                                                                                   <  1:536516/149‑1 (MQ=255)
                                                            acacTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTc                                                                                   <  1:587093/149‑1 (MQ=255)
                                                             cacTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCt                                                                                  <  2:713485/149‑1 (MQ=255)
                                                                tgtgtgGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCtatt                                                                               <  3:178926/149‑1 (MQ=255)
                                                                     ggCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTattttttt                                                                          <  1:756385/149‑1 (MQ=255)
                                                                        ggTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGc                                                                       <  3:749664/149‑1 (MQ=255)
                                                                         gTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGcc                                                                      >  2:318082/1‑149 (MQ=255)
                                                                           aCAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAATTTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGa                                                                    >  2:761406/1‑149 (MQ=255)
                                                                             aaTTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCaa                                                                                                                                                  >  4:49844/1‑69 (MQ=255)
                                                                             aaTTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCaa                                                                                                                                                  <  3:49844/69‑1 (MQ=255)
                                                                             aaTTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATa                                                                  >  2:59377/1‑149 (MQ=255)
                                                                             aaTTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATa                                                                  <  1:79136/149‑1 (MQ=255)
                                                                             aaTTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATa                                                                  >  4:556002/1‑149 (MQ=255)
                                                                               tttATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATaaa                                                                <  1:122360/149‑1 (MQ=255)
                                                                                    aaaTTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTAcc                                                           <  2:163555/149‑1 (MQ=255)
                                                                                      aTTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACcgc                                                         >  2:251800/1‑149 (MQ=255)
                                                                                       ttACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACcgcg                                                        >  2:728467/1‑149 (MQ=255)
                                                                                                 tGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACAcc                                              >  3:846049/1‑149 (MQ=255)
                                                                                                   aaaGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCAt                                            <  3:13205/149‑1 (MQ=255)
                                                                                                   aaaGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGCCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCAt                                            >  4:521676/1‑149 (MQ=255)
                                                                                                      gATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTc                                         <  1:251800/149‑1 (MQ=255)
                                                                                                        tGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCAt                                       >  4:802960/1‑149 (MQ=255)
                                                                                                          cGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTCTTTTGCCGATAAACTACCGCGCAGATACACCATGTCATtc                                     <  4:500771/149‑1 (MQ=255)
                                                                                                          cGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATtc                                     <  3:733436/149‑1 (MQ=255)
                                                                                                          cGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATtc                                     <  1:626624/149‑1 (MQ=255)
                                                                                                          cGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATtc                                     >  4:349028/1‑149 (MQ=255)
                                                                                                           gTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATtct                                    <  3:286063/149‑1 (MQ=255)
                                                                                                            ttttATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAGTGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATGCACCATGTCATtctc                                   >  2:257224/1‑149 (MQ=255)
                                                                                                             tttATATCCCGATTCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTctca                                  <  1:630394/149‑1 (MQ=255)
                                                                                                             tttATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTctca                                  <  3:53595/149‑1 (MQ=255)
                                                                                                                 tatCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTctcactca                              <  4:44570/149‑1 (MQ=255)
                                                                                                                    cccGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTg                           >  4:132581/1‑149 (MQ=255)
                                                                                                                      cGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCg                         >  2:669827/1‑149 (MQ=255)
                                                                                                                       gATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCGa                        <  4:36380/149‑1 (MQ=255)
                                                                                                                          gCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCGACAg                     >  1:265693/1‑149 (MQ=255)
                                                                                                                          gCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCGACAg                     <  2:775294/149‑1 (MQ=255)
                                                                                                                            aGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGAt                                                                   <  1:163/101‑1 (MQ=255)
                                                                                                                            aGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGAt                                                                   >  2:163/1‑101 (MQ=255)
                                                                                                                            aGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCGACAGGa                   >  2:108383/1‑149 (MQ=255)
                                                                                                                             ggAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCGACAGGAt                  <  3:31530/149‑1 (MQ=255)
                                                                                                                              gAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCGACAGGAtt                 >  2:658544/1‑149 (MQ=255)
                                                                                                                               aaTATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTCCTGTCTATTGATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCGACAGGATTa                >  3:534082/1‑149 (MQ=255)
                                                                                                                                 tatTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCGACAGGATTACt              <  4:390358/149‑1 (MQ=255)
                                                                                                                                  atTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCGACAGGATTACTg             <  1:761406/149‑1 (MQ=255)
                                                                                                                                   tttCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGcc                                                                      <  4:611772/91‑1 (MQ=255)
                                                                                                                                   tttCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGcc                                                                      >  3:611772/1‑91 (MQ=255)
                                                                                                                                   tttCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCGACAGGATTACTgg            >  4:517072/1‑149 (MQ=255)
                                                                                                                                   tttCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCCTTTTTCTGTCTATTTATTTTGCCGATAAACTAc                                                            <  3:170637/101‑1 (MQ=255)
                                                                                                                                   tttCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCCTTTTTCTGTCTATTTATTTTGCCGATAAACTAc                                                            >  4:170637/1‑101 (MQ=255)
                                                                                                                                    ttCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCGACAGGATTACTggg           <  4:237013/149‑1 (MQ=255)
                                                                                                                                       aaTTTGACCATTGGCCAAATTGGTAATGCTATGTCGCTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCGACAGGATTACTGGgtgt        >  4:409925/2‑149 (MQ=255)
                                                                                                                                           tGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCGACAGGATTACTGGGTGTTTAt    <  2:454052/149‑1 (MQ=255)
                                                                                                                                            gACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCGACAGGATTACTGGGTGTTCAtt   >  1:575803/1‑149 (MQ=255)
                                                                                                                                             aCCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCGACAGGATTACTGGGTGTTTAttt  <  3:646686/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
AAAAATCGCTTAGGGGATGTTATGCTCACGAAAAAGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCGACAGGATTACTGGGTGTTTATTT  >  NC_000913/4062226‑4062515

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: