Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,933,395 G→T 5.8% C98* (TGC→TGA fucA ← L‑fuculose‑1‑phosphate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,933,3950GT 5.8% 293.2 / 6.4 87C98* (TGC→TGAfucAL‑fuculose‑1‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base G (35/47);  new base T (2/3);  total (37/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.44e-01

GATTTTTGAGAGCCAGCGCAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACAATATGCGACTCCGTCAGTT  >  NC_000913/2933250‑2933543
                                                                                                                                                 |                                                                                                                                                    
gATTTTTGAGAGCCAGCGCAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCaa                                                                                                                                                   >  1:84779/1‑149 (MQ=255)
gATTTTTGAGAGCCAGCGCAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCaa                                                                                                                                                   <  1:742106/149‑1 (MQ=255)
 aTTTTTGAGAGCCAGCGCAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAAt                                                                                                                                                  <  2:854132/149‑1 (MQ=255)
    tttGAGAGCCAGCGCAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGaa                                                                                                                                               <  1:600887/149‑1 (MQ=255)
      tGAGAGCCAGCGCAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACg                                                                                                                                             <  3:611411/149‑1 (MQ=255)
        agagCCAGCGCAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGc                                                                                                                                           >  3:755577/1‑149 (MQ=255)
         gagCCAGCGCAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTTCATTGCACGGCa                                                                                                                                          <  1:12605/149‑1 (MQ=255)
         gagCCAGCGCAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTTCATTGCACGGCa                                                                                                                                          <  4:174282/149‑1 (MQ=255)
          agCCAGCGCAAGATGTTCAGAAAGTTCGCGTGTTCCACAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTTCATTGCACGGCag                                                                                                                                         <  3:574132/149‑2 (MQ=255)
             cAGCGCAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGAt                                                                                                                                      >  3:699706/1‑149 (MQ=255)
             cAGCGCAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGAt                                                                                                                                      >  3:626765/1‑149 (MQ=255)
                cgcAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATtgt                                                                                                                                   <  3:718969/149‑1 (MQ=255)
                 gcAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATtgtg                                                                                                                                  <  2:255955/149‑1 (MQ=255)
                   aaCATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGaa                                                                                                                                <  3:172660/149‑1 (MQ=255)
                      aTGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGaacaa                                                                                                                             <  1:350201/149‑1 (MQ=255)
                       tGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGaacaac                                                                                                                            >  2:28185/1‑149 (MQ=255)
                         ttCAGAAAGTTCCCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCAGGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACcg                                                                                                                          >  4:802233/1‑149 (MQ=255)
                           cAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGTATGGAAACTGCGGTGCAATGAACGGCATGATTGTGAACAACcgcg                                                                                                                        >  1:693198/1‑149 (MQ=255)
                           cAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACcgcg                                                                                                                        >  3:368380/1‑149 (MQ=255)
                            aGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGt                                                                                                                       >  3:796021/1‑149 (MQ=255)
                             gAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGtt                                                                                                                      <  1:640644/149‑1 (MQ=255)
                             gAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGtt                                                                                                                      >  3:691665/1‑149 (MQ=255)
                             gAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGtt                                                                                                                      <  1:332204/149‑1 (MQ=255)
                             gAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCAGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGtt                                                                                                                      <  2:638114/149‑1 (MQ=255)
                              aaaGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTg                                                                                                                     <  4:542710/149‑1 (MQ=255)
                                aGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGc                                                                                                                   <  2:651183/149‑1 (MQ=255)
                                  ttCGCGTGTTCCAACGGTCGCATCAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATAGGTTAAGAATGGAAACTGCCGTGCCCTGAACGGCATGATCGTGCACAACCGCGTTGGCAt                                                                                                                 >  2:891143/1‑149 (MQ=255)
                                    cgcgTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATcc                                                                                                               <  4:868230/149‑1 (MQ=255)
                                     gcgTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCg                                                                                                              <  4:691665/149‑1 (MQ=255)
                                     gcgTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCg                                                                                                              >  1:592857/1‑149 (MQ=255)
                                          ttCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTg                                                                                                         <  4:626765/149‑1 (MQ=255)
                                           tCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGc                                                                                                        <  3:39142/149‑1 (MQ=255)
                                           tCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGc                                                                                                        <  3:338435/149‑1 (MQ=255)
                                            ccATAGGTTGCATAAGGCGCGCAAGGGATAGAATTACCACCAGCCGCCGCAAGCATGTAGTGAACAGCTGGGATCGATCGGTTAAGAATGGAAACTGCGGTGCAATGAACGGCATGGCTGTTAACAACCGCGTTGGCATCCGGgctgct                                                                                                       <  2:237619/149‑1 (MQ=255)
                                            ccAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCt                                                                                                       >  4:486641/1‑149 (MQ=255)
                                             cAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCtt                                                                                                      >  1:649941/1‑149 (MQ=255)
                                             cAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCtt                                                                                                      <  1:558423/149‑1 (MQ=255)
                                              aaaGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCttt                                                                                                     <  1:628689/149‑1 (MQ=255)
                                                 ggTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGAt                                                                                                  <  2:883499/149‑1 (MQ=255)
                                                      cATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCt                                                                                             <  2:737297/149‑1 (MQ=255)
                                                       aTAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTg                                                                                            >  3:762858/1‑149 (MQ=255)
                                                           ggCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAGGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCat                                                                                        >  1:639582/1‑149 (MQ=255)
                                                                 cAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGaa                                                                                  >  2:475233/1‑149 (MQ=255)
                                                                   aGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAAc                                                                                >  3:716261/1‑149 (MQ=255)
                                                                     gAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGa                                                                                                   >  2:689695/1‑128 (MQ=255)
                                                                     gAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGa                                                                                                   <  1:689695/128‑1 (MQ=255)
                                                                        tAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGcc                                                                                          >  2:165487/1‑134 (MQ=255)
                                                                        tAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGcc                                                                                          <  1:165487/134‑1 (MQ=255)
                                                                             ttACCGCCAGCCGCCGCAATCATGTAGTGACTAGCGGGGATCGATCGGTTAGGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCt                                                                      >  2:868905/1‑149 (MQ=255)
                                                                             ttACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCt                                                                      >  2:698113/1‑149 (MQ=255)
                                                                              tACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCtt                                                                     <  1:194537/149‑1 (MQ=255)
                                                                                 cGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAg                                                                  <  1:635409/149‑1 (MQ=255)
                                                                                 cGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAg                                                                  <  4:796021/149‑1 (MQ=255)
                                                                                   ccAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAggg                                                                <  1:676148/149‑1 (MQ=255)
                                                                                     aGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGa                                                              <  3:812674/149‑1 (MQ=255)
                                                                                     aGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGa                                                              <  2:554204/149‑1 (MQ=255)
                                                                                         gccgcAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCtt                                                          >  4:356141/1‑149 (MQ=255)
                                                                                              aaTCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTttcctt                                                     <  3:479000/149‑1 (MQ=255)
                                                                                                tCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTttccttcc                                                   <  4:601230/149‑1 (MQ=255)
                                                                                                      aGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGt                                             <  4:638372/149‑1 (MQ=255)
                                                                                                       gTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGtt                                            >  2:725219/1‑149 (MQ=255)
                                                                                                       gTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGtt                                            >  3:552079/1‑149 (MQ=255)
                                                                                                          aaTAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTAc                                         >  3:485683/1‑149 (MQ=255)
                                                                                                           aTAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCAt                                               >  4:581023/1‑142 (MQ=255)
                                                                                                           aTAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCAt                                               <  3:581023/142‑1 (MQ=255)
                                                                                                            tAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCAt                                                                           <  4:677606/113‑1 (MQ=255)
                                                                                                            tAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCAt                                                                           >  3:677606/1‑113 (MQ=255)
                                                                                                              gCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGtt                                     <  3:248677/149‑1 (MQ=255)
                                                                                                                ggggATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGc                                   <  4:552079/149‑1 (MQ=255)
                                                                                                                    atcgatcgGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATc                               >  2:469120/1‑149 (MQ=255)
                                                                                                                       gatcgGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTttccttcc                                                   >  4:17768/1‑126 (MQ=255)
                                                                                                                       gatcgGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTttccttcc                                                   <  3:17768/126‑1 (MQ=255)
                                                                                                                        atcgGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATa                           >  1:824645/1‑149 (MQ=255)
                                                                                                                        atcgGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATa                           >  1:119176/1‑149 (MQ=255)
                                                                                                                          cgGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATaaa                         <  4:875707/149‑1 (MQ=255)
                                                                                                                          cgGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATaaa                         <  2:592857/149‑1 (MQ=255)
                                                                                                                          cgGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATaaa                         <  4:463743/149‑1 (MQ=255)
                                                                                                                            gTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGa                       >  2:202942/1‑149 (MQ=255)
                                                                                                                            gTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGa                       >  3:685183/1‑149 (MQ=255)
                                                                                                                             ttAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGAc                      >  4:418583/1‑149 (MQ=255)
                                                                                                                                 gAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACAata                  <  2:243262/149‑1 (MQ=255)
                                                                                                                                    tGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACAATATGc               <  4:400520/149‑1 (MQ=255)
                                                                                                                                    tGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACAATATGc               <  3:356141/149‑1 (MQ=255)
                                                                                                                                           tGCCGTTCATTGCACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATGGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACAATATGCGACTCCg        >  1:757387/1‑149 (MQ=255)
                                                                                                                                           tGCCGTTCATTGCACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACAATATGCGACTCCg        >  3:610777/1‑149 (MQ=255)
                                                                                                                                            gCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACAATATGCGACTCCGt       <  2:824645/149‑1 (MQ=255)
                                                                                                                                             gcttGCATTGAACGGTTTGAGTGAGAAGAACCGGGTTGGCGTCCGGTCTGGTTTGATAGGCTGCCAGATGCAAACGCCATTCGCTTGAGGGGGGCTTACCTTCCGTATGTTTACCGTTGCCATCAATGAAGACAATATGCGCCTCCGTc      <  2:511764/146‑1 (MQ=255)
                                                                                                                                                 gCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACaa                    <  2:352732/131‑1 (MQ=255)
                                                                                                                                                 gCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACaa                    >  1:352732/1‑131 (MQ=255)
                                                                                                                                                 gCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACAATATGCGACTCCGTCAGtt  <  4:580292/149‑1 (MQ=255)
                                                                                                                                                 gCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACAATATGCGACTCCGTCAGtt  <  1:868905/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GATTTTTGAGAGCCAGCGCAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACAATATGCGACTCCGTCAGTT  >  NC_000913/2933250‑2933543

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