Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,418,629 T→C 8.6% intergenic (+18/‑413) yhdV → / → yhdW putative outer membrane protein/pseudogene, ABC transporter periplasmic binding protein family;putative transport; Not classified; putative periplasmic binding transport protein; periplasmic‑binding component of ABC superfamily

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,418,6290TC 8.6% 337.2 / 18.2 116intergenic (+18/‑413)yhdV/yhdWputative outer membrane protein/pseudogene, ABC transporter periplasmic binding protein family;putative transport; Not classified; putative periplasmic binding transport protein; periplasmic‑binding component of ABC superfamily
Reads supporting (aligned to +/- strand):  ref base T (62/43);  new base C (3/7);  total (66/50)
Fisher's exact test for biased strand distribution p-value = 9.97e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.83e-01

Only 100 of 119 total aligned reads displayed.

TCGTCTGGTTCACACCAATACCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCAATTTAATTTTCACACTT  >  NC_000913/3418482‑3418773
                                                                                                                                                   |                                                                                                                                                
tcgtcTGGTTCACACCAATACCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGtt                                                                                                                                                 >  2:703621/1‑149 (MQ=255)
  gtcTGGTTCACACCAATACCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTcc                                                                                                                                               <  4:505517/149‑1 (MQ=255)
    cTGGTTCACACCAATACCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCcgc                                                                                                                                             <  1:1009364/149‑1 (MQ=255)
     tGGTTCACACCAATACCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCcgcg                                                                                                                                            >  4:566767/1‑149 (MQ=255)
     tGGTTCACACCAATACCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGGTGcgcg                                                                                                                                            >  3:817025/1‑149 (MQ=255)
       gTTCACACCAATACCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAACTCCgcgcc                                                                                                                                          <  2:990293/149‑2 (MQ=255)
        ttCACACCAATACCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGctc                                                                                                                                         >  3:851436/1‑149 (MQ=255)
        ttCACACCAATACCAGCATGCAAGGGACGACACAGGATAACAGGGAAACGCCTTGCGCTATTGCCGGCGGGGCAGCAACTGTTGCTGGCCTGTCGATGGGGATGATCGCGCTACCGAAATAACTGAACATTGCCCGGGGTTCCGCGctc                                                                                                                                         >  4:191239/1‑149 (MQ=255)
         tCACACCAATACCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAACTCCGCGccct                                                                                                                                        <  4:657115/149‑1 (MQ=255)
         tCACACCAATACCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGAGCGCGGAACTCCGCGccct                                                                                                                                        <  4:102843/149‑1 (MQ=255)
         tCACACCAATACCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGAGCGCGGAACTCCGCGccct                                                                                                                                        <  3:657510/149‑1 (MQ=255)
          cacaCCAATACCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCtt                                                                                                                                       <  4:771763/149‑1 (MQ=255)
           acacCAATACCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCttt                                                                                                                                      >  1:1063739/1‑149 (MQ=255)
             acCAATACCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGTGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCttttt                                                                                                                                    >  1:760989/1‑149 (MQ=255)
             acCAATACCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCttttt                                                                                                                                    <  1:53941/149‑1 (MQ=255)
               cAATACCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATCGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCCCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTaa                                                                                                                                  <  4:791793/149‑1 (MQ=255)
                  tACCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGa                                                                                                                               <  1:556563/149‑1 (MQ=255)
                   aCCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGAt                                                                                                                              <  2:1087417/149‑1 (MQ=255)
                    ccTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATa                                                                                                                             <  2:1063739/149‑1 (MQ=255)
                      tGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACa                                                                                                                           >  3:565117/1‑149 (MQ=255)
                       gTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAg                                                                                                                          >  4:445437/1‑149 (MQ=255)
                       gTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGAAGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAg                                                                                                                          <  4:473175/149‑1 (MQ=255)
                         aTGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTa                                                                                                                        <  3:19034/149‑1 (MQ=255)
                           gAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTAtc                                                                                                                      >  3:918807/1‑149 (MQ=255)
                           gAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTAtc                                                                                                                      >  3:764017/1‑149 (MQ=255)
                           gAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTAtc                                                                                                                      >  2:809723/1‑149 (MQ=255)
                              aGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTAtctct                                                                                                                   <  3:239229/149‑1 (MQ=255)
                              aGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTAtctct                                                                                                                   >  1:1087574/1‑149 (MQ=255)
                                      cTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGa                                                                                                           <  2:997470/149‑1 (MQ=255)
                                       tCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGAt                                                                                                          >  2:545760/1‑149 (MQ=255)
                                         aGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATaa                                                                                                        >  2:681137/1‑149 (MQ=255)
                                             tAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTa                                                                                                                        >  2:1000482/1‑129 (MQ=255)
                                             tAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTa                                                                                                                        <  1:1000482/129‑1 (MQ=255)
                                             tAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAata                                                                                                    >  4:706373/1‑149 (MQ=255)
                                             tAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAACTCCGCGccct                                                                                                                                        >  1:678940/1‑113 (MQ=255)
                                             tAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAACTCCGCGccct                                                                                                                                        <  2:678940/113‑1 (MQ=255)
                                              aaCTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAatat                                                                                                   >  3:648247/1‑149 (MQ=255)
                                              aaCTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAatat                                                                                                   >  3:931025/1‑149 (MQ=255)
                                                cTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTg                                                                                                 <  4:52729/149‑1 (MQ=255)
                                                   ggAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATaaaa                                                                                                      >  3:778407/1‑141 (MQ=255)
                                                   ggAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATaaaa                                                                                                      <  4:778407/141‑1 (MQ=255)
                                                   ggAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCta                                                                                              >  3:432149/1‑149 (MQ=255)
                                                   ggAAACGGCTGAGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCta                                                                                              >  3:614190/1‑149 (MQ=255)
                                                    gAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATaaaa                                                                                                      >  4:302580/1‑140 (MQ=255)
                                                    gAAACGGATGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATaaaa                                                                                                      <  3:302580/140‑1 (MQ=255)
                                                       aCGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATAtt                                                                                          >  3:57911/1‑149 (MQ=255)
                                                       aCGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATAtt                                                                                          >  4:992185/1‑149 (MQ=255)
                                                            tGGGGCTATTGCCGGCGGGGCAGCTACTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGCCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTctgct                                                                                     <  3:191239/149‑1 (MQ=255)
                                                              gggCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTgg                                                                                   >  2:193310/1‑149 (MQ=255)
                                                                  tATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTaa                                                                               >  3:1023793/1‑149 (MQ=255)
                                                                        cggcggGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCct                                                                         <  1:1080959/149‑1 (MQ=255)
                                                                          gcggGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCctgc                                                                       <  1:746672/149‑1 (MQ=255)
                                                                            ggggCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAGCTGAACAGGGCTCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCctgctg                                                                     <  4:1034117/149‑1 (MQ=255)
                                                                            ggggCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCctgctg                                                                     >  1:648634/1‑149 (MQ=255)
                                                                               gCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGAtt                                                                  >  3:1030135/1‑149 (MQ=255)
                                                                                 aGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTa                                                                <  4:1032830/149‑1 (MQ=255)
                                                                                       tgTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGt                                                          >  3:71864/1‑149 (MQ=255)
                                                                                             tGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCAt                                                            >  1:196200/1‑141 (MQ=255)
                                                                                             tGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCAt                                                            <  2:196200/141‑1 (MQ=255)
                                                                                             tGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATc                                                    >  1:294249/1‑149 (MQ=255)
                                                                                                ccTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATcgcg                                                 <  2:530753/149‑1 (MQ=255)
                                                                                                  tGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGtt                                               <  1:681137/149‑1 (MQ=255)
                                                                                                   gACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTa                                              <  1:25900/149‑1 (MQ=255)
                                                                                                   gACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTa                                              >  1:947793/1‑149 (MQ=255)
                                                                                                    aCGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTAt                                             >  1:883163/1‑149 (MQ=255)
                                                                                                      gATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTAttt                                           >  4:285763/1‑149 (MQ=255)
                                                                                                          gggATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCgct                                       >  2:696647/1‑149 (MQ=255)
                                                                                                          gggATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCgct                                       <  4:57911/149‑1 (MQ=255)
                                                                                                            gATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCgctgc                                     <  1:1051918/149‑1 (MQ=255)
                                                                                                              tcatcGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCtt                                   >  2:48168/1‑149 (MQ=255)
                                                                                                                atcGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATaa                                                                                                             <  2:315088/73‑1 (MQ=255)
                                                                                                                atcGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTcc                                                                           <  1:690160/107‑1 (MQ=255)
                                                                                                                atcGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTc                                 >  4:840391/1‑149 (MQ=255)
                                                                                                                atcGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTc                                 >  1:736812/1‑149 (MQ=255)
                                                                                                                atcGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTc                                 <  4:12444/149‑1 (MQ=255)
                                                                                                                atcGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTc                                 >  3:12444/1‑149 (MQ=255)
                                                                                                                    ctctATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATGGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGaaa                             <  2:760989/149‑1 (MQ=255)
                                                                                                                     tctATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAAt                            <  3:202276/149‑1 (MQ=255)
                                                                                                                       tATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATaa                          >  1:353246/1‑149 (MQ=255)
                                                                                                                       tATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATaa                          >  3:458926/1‑149 (MQ=255)
                                                                                                                        aTCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAAc                         >  3:790730/1‑149 (MQ=255)
                                                                                                                         tCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACg                        >  3:726837/1‑149 (MQ=255)
                                                                                                                         tCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACg                        >  3:642877/1‑149 (MQ=255)
                                                                                                                         tCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACg                        >  4:362121/1‑149 (MQ=255)
                                                                                                                         tCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACg                        >  3:122637/1‑149 (MQ=255)
                                                                                                                           gAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCa                      <  4:307128/149‑1 (MQ=255)
                                                                                                                              aTAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCa                   >  1:228918/1‑149 (MQ=255)
                                                                                                                               tAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCaa                  <  1:193310/149‑1 (MQ=255)
                                                                                                                                aaCTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCaat                 >  2:897438/1‑149 (MQ=255)
                                                                                                                                 aCTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCaatt                <  4:842111/149‑1 (MQ=255)
                                                                                                                                    gAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCaatttaa             >  3:89640/1‑149 (MQ=255)
                                                                                                                                    gAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCaatttaa             <  2:769819/149‑1 (MQ=255)
                                                                                                                                      aaagagcgcGGAACTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCaatttaatt           >  1:518772/6‑149 (MQ=255)
                                                                                                                                       cAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCaatttaattt          >  4:952161/1‑149 (MQ=255)
                                                                                                                                        aGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCAATTTAAtttt         >  4:814033/1‑149 (MQ=255)
                                                                                                                                          ggCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGACTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCAATTTAATTTTca       <  3:623542/149‑1 (MQ=255)
                                                                                                                                             gcgGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCAATTTAATTTTCacac    >  1:934578/1‑149 (MQ=255)
                                                                                                                                               ggAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCAATTTAATTTTCACACtt  >  4:556867/1‑149 (MQ=255)
                                                                                                                                                   ttCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGt                                                >  2:971023/1‑99 (MQ=255)
                                                                                                                                                   ttCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGt                                                <  1:971023/99‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
TCGTCTGGTTCACACCAATACCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCAATTTAATTTTCACACTT  >  NC_000913/3418482‑3418773

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