Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,933,398 A→T 6.6% H97Q (CAT→CAA 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,3980AT 6.6% 298.5 / 8.1 89H97Q (CAT→CAAfucAL‑fuculose‑1‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base A (34/48);  new base T (2/4);  total (37/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATTTTTGAGAGCCAGCGCAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACAATATGCGACTCCGTCAGTTTT  >  NC_000913/2933250‑2933545
                                                                                                                                                    |                                                                                                                                                   
gATTTTTGAGAGCCAGCGCAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCaa                                                                                                                                                     <  1:742106/149‑1 (MQ=255)
gATTTTTGAGAGCCAGCGCAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCaa                                                                                                                                                     >  1:84779/1‑149 (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:626765/1‑149 (MQ=255)
             cAGCGCAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGAt                                                                                                                                        >  3:699706/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                                                                                                                >  1:592857/1‑149 (MQ=255)
                                     gcgTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCg                                                                                                                <  4:691665/149‑1 (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:558423/149‑1 (MQ=255)
                                             cAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCtt                                                                                                        >  1:649941/1‑149 (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                                                                                                     <  1:689695/128‑1 (MQ=255)
                                                                     gAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGa                                                                                                     >  2:689695/1‑128 (MQ=255)
                                                                        tAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGcc                                                                                            <  1:165487/134‑1 (MQ=255)
                                                                        tAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGcc                                                                                            >  2:165487/1‑134 (MQ=255)
                                                                             ttACCGCCAGCCGCCGCAATCATGTAGTGACTAGCGGGGATCGATCGGTTAGGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCt                                                                        >  2:868905/1‑149 (MQ=255)
                                                                             ttACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCt                                                                        >  2:698113/1‑149 (MQ=255)
                                                                              tACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCtt                                                                       <  1:194537/149‑1 (MQ=255)
                                                                                 cGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAg                                                                    <  4:796021/149‑1 (MQ=255)
                                                                                 cGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAg                                                                    <  1:635409/149‑1 (MQ=255)
                                                                                   ccAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAggg                                                                  <  1:676148/149‑1 (MQ=255)
                                                                                     aGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGa                                                                <  2:554204/149‑1 (MQ=255)
                                                                                     aGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGa                                                                <  3:812674/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                                                 <  3:581023/142‑1 (MQ=255)
                                                                                                           aTAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCAt                                                 >  4:581023/1‑142 (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                                                     <  3:17768/126‑1 (MQ=255)
                                                                                                                       gatcgGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTttccttcc                                                     >  4:17768/1‑126 (MQ=255)
                                                                                                                        atcgGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATa                             >  1:824645/1‑149 (MQ=255)
                                                                                                                        atcgGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATa                             >  1:119176/1‑149 (MQ=255)
                                                                                                                          cgGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATaaa                           <  4:463743/149‑1 (MQ=255)
                                                                                                                          cgGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATaaa                           <  4:875707/149‑1 (MQ=255)
                                                                                                                          cgGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATaaa                           <  2:592857/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                 <  3:356141/149‑1 (MQ=255)
                                                                                                                                    tGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACAATATGc                 <  4:400520/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    <  1:868905/149‑1 (MQ=255)
                                                                                                                                                 gCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACAATATGCGACTCCGTCAGtt    <  4:580292/149‑1 (MQ=255)
                                                                                                                                                   aaTGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACAATATGCGACTCCGTCAGtttt  >  2:12368/1‑149 (MQ=255)
                                                                                                                                                   aaTGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACAATATGCGACTCCGTCAGtttt  >  1:866234/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GATTTTTGAGAGCCAGCGCAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACAATATGCGACTCCGTCAGTTTT  >  NC_000913/2933250‑2933545

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: