Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,693,326 C→G R637R (CGG→CGC bcsA ← cellulose synthase, catalytic subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,693,3260CG100.0% 327.3 / NA 91R637R (CGG→CGCbcsAcellulose synthase, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (51/40);  total (51/40)

CTGCGCCGCCAAGAACAATCAGGTTGTAGAACACCCACACCATACTGACGACCACGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGCCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCGTAGATTTCACTCC  >  NC_000913/3693181‑3693471
                                                                                                                                                 |                                                                                                                                                 
cTGCGCCGCCAAGAACAATCAGGTTGTAGAACACCCACACCATACTGACGACCACGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGgcgc                                                                                                                                                <  1:121492/149‑1 (MQ=255)
     ccgccAAGAACAATCAGGTTGTAGAACACCCACACCATACTGACGACCACGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGAt                                                                                                                                           <  7:297976/149‑1 (MQ=255)
     ccgccAAGAACAATCAGGTTGTAGAACACCCACACCATACTGACGACCACGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGAt                                                                                                                                           <  3:192860/149‑1 (MQ=255)
      cgccAAGAACAATCAGGTTGTAGAACACCCACACCATACTGACGACCACGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATc                                                                                                                                          <  7:347631/149‑1 (MQ=255)
        ccAAGAACAATCAGGTTGTAGAACACCCACACCATACTGACGACCACGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCAc                                                                                                                                        <  4:169525/149‑1 (MQ=255)
        ccAAGAACAATCAGGTTGTAGAACACCCACACCATACTGACGACCACGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCAc                                                                                                                                        >  5:286297/1‑149 (MQ=255)
         cAAGAACAATCAGGTTGTAGAACACCCACACCATACTGACGACCACGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCAcc                                                                                                                                       <  6:91023/149‑1 (MQ=255)
           aGAACAATCAGGTTGTAGAACACCCACACCATACTGACGACCACGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCAcc                                                                                                                                       >  4:245270/1‑147 (MQ=255)
                aaTCAGGTTGTAGAACACCCACACCATACTGACGCCCACGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGa                                                                                                                                >  7:68194/1‑149 (MQ=255)
                 aTCAGGTTGTAGAACACCCACACCATACTGACGACCACGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGAc                                                                                                                               >  8:315276/1‑149 (MQ=255)
                 aTCAGGTTGTAGAACACCCACACCATACTGACGACCACGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGAc                                                                                                                               <  5:234051/149‑1 (MQ=255)
                  tCAGGTTGTAGAACACCCACACCATACTGACGACCACGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACg                                                                                                                              >  3:254277/1‑149 (MQ=255)
                  tCAGGTTGTAGAACACCCACACCATACTGACGACCACGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACg                                                                                                                              >  2:115829/1‑149 (MQ=255)
                    aGGTTGTAGAACACCCACACCATACTGACGACCACGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGTGCGAGATCACCCAGTCGACGTa                                                                                                                            <  8:291168/149‑1 (MQ=255)
                      gTTGTAGAACACCCACACCATACTGACGACCACGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTAct                                                                                                                          >  8:234921/1‑149 (MQ=255)
                       ttGTAGAACACCCACACCATACTGACGACCACGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTActc                                                                                                                         <  3:355779/149‑1 (MQ=255)
                         gTAGAACACCCACACCATACTGCCGACCACGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCCGTCGACGTACtctt                                                                                                                       >  3:106220/1‑149 (MQ=255)
                         gTAGAACACCCACACCATACTGACGACCACGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACtctt                                                                                                                       >  6:158663/1‑149 (MQ=255)
                          tAGAACACCCACACCATACTGACGACCACGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGGTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACtcttc                                                                                                                      <  5:247523/149‑1 (MQ=255)
                            gAACACCGACACCATACTGACGAACACGGTTAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTcttctt                                                                                                                    <  6:286297/149‑1 (MQ=255)
                              acacCCACACCATACTGACGACCACGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCttcttc                                                                                                                   <  4:333515/148‑1 (MQ=255)
                              acacCCACACCATACTGACGACCACGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTcc                                                                                                                  <  4:66976/149‑1 (MQ=255)
                                  ccccacCATCCTGACGACCACGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTAACGCCACGCCCCCCAGGTTGAGCAGCACAAGGAAGATGTAGTGGCGCGAGAGCCCCCAGTCGACGTACGCTTCGTccac                                                                                                                >  6:103098/4‑147 (MQ=255)
                                  ccACACCATACTGACGACCACGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCgaga                                                                                                                                            >  2:277898/1‑119 (MQ=255)
                                  ccACACCATACTGACGACCACGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCgaga                                                                                                                                            <  1:277898/119‑1 (MQ=255)
                                  ccACACCATACTGACGACCACGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTccacca                                                                                                              <  7:315276/149‑1 (MQ=255)
                                  ccACACCATACTGACGACCACGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTcca                                                                                                                 <  7:279798/146‑1 (MQ=255)
                                  ccACACCATACTGACGACCACGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTcca                                                                                                                 >  8:279798/1‑146 (MQ=255)
                                  ccACACCATACTGACGACCACGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTAGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTccac                                                                                                                >  6:63636/1‑147 (MQ=255)
                                   cacaCCATACTGACGACCACGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCACCGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAg                                                                                                             >  2:105210/1‑149 (MQ=255)
                                          tACTGACGACCACGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCt                                                                                                      >  1:69481/1‑149 (MQ=255)
                                             tGACGACCACGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTg                                                                                                   >  7:370651/1‑149 (MQ=255)
                                             tGACGACCACGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTg                                                                                                   >  4:191958/1‑149 (MQ=255)
                                                cgacCACGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCg                                                                                                <  1:115829/149‑1 (MQ=255)
                                                    cACGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGa                                                                                            >  7:89879/1‑149 (MQ=255)
                                                     aCGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGCTGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGAc                                                                                           >  5:183247/1‑149 (MQ=255)
                                                         tGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGATTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGtt                                                                                        >  6:178444/1‑148 (MQ=255)
                                                         tGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCACACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTAGTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGtt                                                                                        >  7:299103/1‑148 (MQ=255)
                                                           aGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGCCGTACTCTTCTTCCCCCAGTCCACCTTTGGCGGTGACGTTaaa                                                                                     >  1:110128/1‑149 (MQ=255)
                                                           aGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTaaa                                                                                     <  3:245270/149‑1 (MQ=255)
                                                              aTCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAAttt                                                                                  <  8:89879/149‑1 (MQ=255)
                                                               tctcGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCATCGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGTGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCCCCTGTCCACCTTTGGCGGTGACGTTAAATTTg                                                                                 >  1:208397/1‑149 (MQ=255)
                                                                   ggTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGc                                                                                <  5:158663/146‑1 (MQ=255)
                                                                    gTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCttt                                                                            <  2:306647/149‑1 (MQ=255)
                                                                        gCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAAttt                                                                                  >  6:83366/1‑139 (MQ=255)
                                                                        gCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAAttt                                                                                  <  5:83366/139‑1 (MQ=255)
                                                                         cGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCg                                                                       <  3:246485/149‑1 (MQ=255)
                                                                         cGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCg                                                                       <  7:201176/149‑1 (MQ=255)
                                                                             ccATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGtt                                                                   <  7:234921/149‑1 (MQ=255)
                                                                              cATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGtt                                                                   >  8:38268/1‑148 (MQ=255)
                                                                               aTAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTcc                                                                                                                  <  6:146532/100‑1 (MQ=255)
                                                                               aTAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTaa                                                                 >  6:227490/1‑149 (MQ=255)
                                                                                 aGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCATTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAACTTTGACGTTGTGCGGGTTAAt                                                                >  6:296163/1‑148 (MQ=255)
                                                                                    aGTAGCGCCAGATGCCTACCGCAACGCCCACCAGATTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAgc                                                            <  5:178444/149‑1 (MQ=255)
                                                                                        gcgcCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGcca                                                        >  7:20779/1‑149 (MQ=255)
                                                                                           ccAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGccac                                                       >  4:112134/1‑147 (MQ=255)
                                                                                             aGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGc                                                   >  1:46053/1‑149 (MQ=255)
                                                                                              gATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGc                                                                                                 >  5:42979/1‑102 (MQ=255)
                                                                                              gATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGc                                                                                                 <  6:42979/102‑1 (MQ=255)
                                                                                               aTGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGt                                                 <  7:307231/149‑1 (MQ=255)
                                                                                                 gCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCACCGTCg                                               >  4:17259/1‑149 (MQ=255)
                                                                                                   cTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCggt                                             >  8:157605/1‑149 (MQ=255)
                                                                                                   cTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCggt                                             >  6:52918/1‑149 (MQ=255)
                                                                                                    tACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACCTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCggtg                                            >  2:5790/1‑149 (MQ=255)
                                                                                                     aCCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCggtgg                                           >  3:269990/1‑149 (MQ=255)
                                                                                                       cGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGgtggtg                                         <  2:69481/149‑1 (MQ=255)
                                                                                                        gCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGc                                        >  1:253517/1‑149 (MQ=255)
                                                                                                        gCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGc                                        >  2:21677/1‑149 (MQ=255)
                                                                                                         cAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGc                                        >  1:6384/1‑148 (MQ=255)
                                                                                                         cAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCg                                       <  2:46053/149‑1 (MQ=255)
                                                                                                          aaCGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGa                                      >  6:71171/1‑149 (MQ=255)
                                                                                                           aCGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCATTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGat                                     <  6:183247/149‑1 (MQ=255)
                                                                                                               ccaccaGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCGCCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATAc                                 >  2:243042/1‑149 (MQ=255)
                                                                                                               ccaccaGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATAc                                 <  4:269990/149‑1 (MQ=255)
                                                                                                               ccaccaGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATAc                                 >  4:205727/1‑149 (MQ=255)
                                                                                                                caccaGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGTCGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATAc                                 >  8:222047/1‑148 (MQ=255)
                                                                                                                 accaGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCa                               >  2:237795/1‑149 (MQ=255)
                                                                                                                    aGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGc                            <  1:34246/149‑1 (MQ=255)
                                                                                                                    aGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGc                            <  1:40036/149‑1 (MQ=255)
                                                                                                                      gTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCa                          >  4:285398/1‑149 (MQ=255)
                                                                                                                          aGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCa                       <  1:21677/148‑1 (MQ=255)
                                                                                                                              ggACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGtt                  <  5:227490/149‑1 (MQ=255)
                                                                                                                                aCAAGGAAGATGTAGGTGCGCGAGATCACCCAGTCGACGGACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCAACGtttt                <  5:296163/149‑2 (MQ=255)
                                                                                                                                   aGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCGTa             <  5:195989/149‑1 (MQ=255)
                                                                                                                                    ggAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGata                                    >  7:28022/1‑125 (MQ=255)
                                                                                                                                    ggAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGata                                    <  8:28022/125‑1 (MQ=255)
                                                                                                                                    ggAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCg               >  2:356964/1‑146 (MQ=255)
                                                                                                                                    ggAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCg               >  3:110634/1‑146 (MQ=255)
                                                                                                                                      aaGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCGTAGAt          <  3:285398/149‑1 (MQ=255)
                                                                                                                                           gTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCGTAGATTTCAc     >  2:8753/1‑149 (MQ=255)
                                                                                                                                           gTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCGTAGATTTCAc     >  2:23852/1‑149 (MQ=255)
                                                                                                                                           gTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGCTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCGTAGATTTCa      >  6:241599/1‑148 (MQ=255)
                                                                                                                                              ggggCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCATGCCAGCACCGTTTCGTAGATTTCACTcc  <  2:174593/149‑1 (MQ=255)
                                                                                                                                              ggggCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCGTAGATTTCACTcc  <  2:253517/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
CTGCGCCGCCAAGAACAATCAGGTTGTAGAACACCCACACCATACTGACGACCACGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGCCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCGTAGATTTCACTCC  >  NC_000913/3693181‑3693471

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: