Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,240,047 A→G R3R (CGT→CGC yjbT ← PF17089 family protein YjbT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,240,0470AG100.0% 1350.0 / NA 308R3R (CGT→CGCyjbTPF17089 family protein YjbT
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (168/140);  total (168/140)

Only 100 of 308 total aligned reads displayed.

TCAGCAGTACCGAGCGGACCAGCATAAGCAGGAAGAACAGAGACTGATAACATCAAAGCAGCAAAATAAGGCTTCATTTTTACCACCTTTATCAGGTTACGTTTCATTTGTTCCAGAGGAACATTGTCGATTTTTCGCGCATTGCTGGTGGCTGGGAATCACCTGAATGGGTGATTTTTGAATTACCGGCTTTGGT  >  NC_000913/4239949‑4240144
                                                                                                  |                                                                                                 
tCAGCAGTACCGAGCGGACCAGCATAAGCAGGAAGAACAGAGACTGATAACATCAAAGCAGCAAAATAAGGCTTCATTTTTACCACCTTTATCAGGTTGCg                                                                                                 <  1:6623006/101‑1 (MQ=255)
  agcagTACCGAGCGGACCAGCATAAGCAGGAAGAACAGAGACTGATAACATCAAAGCAGCAAAATAAGGCTTCATTTTTACCACCTTTATCAGGTTGCGtt                                                                                               <  2:6623399/101‑1 (MQ=255)
  agcagTACCGAGCGGACCAGCATAAGCAGGAAGAACAGAGACTGATAACATCAAAGCAGCAAAATAAGGCTTCATTTTTACCACCTTTATCAGGTTGCGtt                                                                                               <  1:6623244/101‑1 (MQ=255)
   gcagTACCGAGCGGACCAGCATAAGCAGGAAGAACAGAGACTGATAACATCAAAGCAGCAAAATAAGGCTTCATTTTTACCACCTTTATCAGGTTGCGttt                                                                                              <  2:6622977/101‑1 (MQ=255)
   gcagTACCGAGCGGACCAGCATAAGCAGGAAGAACAGAGACTGATAACATCAAAGCAGCAAAATAAGGCTTCATTTTTACCACCTTTATCAGGTTGCGttt                                                                                              >  1:6623532/1‑101 (MQ=255)
     agTACCGAGCGGACCAGCATAAGCAGGAAGAACAGAGACTGATAACATCAAAGCAGCAAAATAAGGCTTCATTTTTACCACCTTTATCAGGTTGCGTTTCa                                                                                            <  1:6622840/101‑1 (MQ=255)
       tACCGAGCGGACCAGCATAAGCAGGAAGAACAGAGACTGATAACATCAAAGCAGCAAAATAAGGCTTCATTTTTACCACCTTTATCAGGTTGCGTTTCAtt                                                                                          <  2:6623148/101‑1 (MQ=255)
           gAGCGGACCAGCATAAGCAGGAAGAACAGAGACTGATAACATCAAAGCAGCAAAATAAGGCTTCATTTTTACCACCTTTATCAGGTTGCGTTTCATTTGtt                                                                                      <  1:6623218/101‑1 (MQ=255)
            aGCGGACCAGCATAAGCAGGAAGAACAGAGACTGATAACATCAAAGCAGCAAAATAAGGCTTCATTTTTACCACCTTTATCAGGTTGCGTTTCATTTGTTc                                                                                     <  2:6623128/101‑1 (MQ=255)
             gCGGACCAGCATAAGCAGGAAGAACAGAGACTGATAACATCAAAGCAGCAAAATAAGGCTTCATTTTTACCACCTTTATCAGGTTGCGTTTCATTTGTTcc                                                                                    <  2:6623175/101‑1 (MQ=255)
             gCGGACCAGCATAAGCAGGAAGAACAGAGACTGATAACATCAAAGCAGCAAAATAAGGCTTCATTTTTACCACCTTTATCAGGTTGCGTTTCATTTGTTcc                                                                                    >  2:6623540/1‑101 (MQ=255)
             gCGGACCAGCATAAGCAGGAAGAACAGAGACTGATAACATCAAAGCAGCAAAATAAGGCTTCATTTTTACCACCTTTATCAGGTTGCGTTTCATTTGTTcc                                                                                    <  1:6623293/101‑1 (MQ=255)
               ggACCAGCATAAGCAGGAAGAACAGAGACTGATAACATCAAAGCAGCAAAATAAGGCTTCATTTTTACCACCTTTATCAGGTTGCGTTTCATTTGTTCCag                                                                                  >  1:6623544/1‑101 (MQ=255)
               ggACCAGCATAAGCAGGAAGAACAGAGACTGATAACATCAAAGCAGCAAAATAAGGCTTCATTTTTACCACCTTTATCAGGTTGCGTTTCATTTGTTCCag                                                                                  >  2:6623543/1‑101 (MQ=255)
                 aCCAGCATAAGCAGGAAGAACAGAGACTGATAACATCAAAGCAGCAAAATAAGGCTTCATTTTTACCACCTTTATCAGGTTGCGTTTCATTTGTTCCagag                                                                                <  2:6623160/101‑1 (MQ=255)
                  ccAGCATAAGCAGGAAGAACAGAGACTGATAACATCAAAGCAGCAAAATAAGGCTTCATTTTTACCACCTTTATCAGGTTGCGgtt                                                                                              >  2:6623547/1‑86 (MQ=255)
                   cAGCATAAGCAGGAAGAACAGAGACTGATAACATCAAAGCAGCAAAATAAGGCTTCATTTTTACCCCCTTTATCAGGTTGCGTTTCATTTGTTCCAGAGGa                                                                              <  1:6623127/101‑1 (MQ=255)
                    aGCATAAGCAGGAAGAACAGAGACTGATAACATCAAAGCAGCAAAATAAGGCTTCATTTTTACCACCTTTATCAGGTTGCGTTTCATTTGTTCCAGAGGaa                                                                             <  2:6623106/101‑1 (MQ=255)
                       aTAAGCAGGAAGAACAGAGACTGATAACATCAAAGCAGCAAAATAAGGCTTCATTTTTACCACCTTTATCAGGTTGCGTTTCATTTGTTCCAGAGGAACAt                                                                          >  1:6623552/1‑101 (MQ=255)
                         aaGCAGGAAGAACAGAGACTGATAACATCAAAGCAGCAAAATAAGGCTTCATTTTTACCACCTTTATCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTg                                                                        >  2:6623553/1‑101 (MQ=255)
                           gCAGGAAGAACAGAGACTGATAACATCAAAGCAGCAAAATAAGGCTTCATTTTTACCACCTTTATCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTc                                                                      <  1:6623137/101‑1 (MQ=255)
                            cAGGAAGAACAGAGACTGATAACATCAAAGCAGCAAAATAAGGCTTCATTTTTACCACCTTTATCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTCg                                                                     >  2:6623557/1‑101 (MQ=255)
                             aGGAAGAACAGAGACTGATAACATCAAAGCAGCAAAATAAGGCTTCATTTTTACCACCTTTATCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTCGa                                                                    >  2:6623561/1‑101 (MQ=255)
                               gaagaaCAGAGACTGATAACATCAAAGCAGCAAAATAAGGCTTCATTTTTACCACCTTTATCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTCGAtt                                                                  <  2:6623079/101‑1 (MQ=255)
                               gaagaaCAGAGACTGATAACATCAAAGCAGCAAAATAAGGCTTCATTTTTACCACCTTTATCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTCGAtt                                                                  <  1:6623193/101‑1 (MQ=255)
                                aagaaCAGAGACTGATAACATCAAAGCAGCAAAATAAGGCTTCATTTTTACCACCTTTATCAGGTTGCGTTTCATTTGTTCCAGAGGa                                                                              <  1:6623158/88‑1 (MQ=255)
                                aagaaCAGAGACTGATAACATCAAAGCAGCAAAATAAGGCTTCATTTTTACCACCTTTATCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTCGAttt                                                                 >  1:6623563/1‑101 (MQ=255)
                                 agaaCAGAGACTGATAACATCAAAGCAGCAAAATAAGGCTTCATTTTTACCACCTTTATCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTCGAtttt                                                                <  1:6623084/101‑1 (MQ=255)
                                 agaaCAGAGACTGATAACATCAAAGCAGCAAAATAAGGCTTCATTTTTACCACCTTTATCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTCGAtttt                                                                >  1:6623568/1‑101 (MQ=255)
                                 agaaCAGAGACTGATAACATCAAAGCAGCAAAATAAGGCTTCATTTTTACCACCTTTATCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTCGAtttt                                                                >  2:6623567/1‑101 (MQ=255)
                                 agaaCAGAGACTGATAACATCAAAGCAGCAAAATAAGGCTTCATTTTTACCACCTTTATCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTCGAtttt                                                                <  2:6623228/101‑1 (MQ=255)
                                  gaaCAGAGACTGATAACATCAAAGCAGCAAAATAAGGCTTCATTTTTACCACCTTTATCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTCGAttttt                                                               <  2:6623236/101‑1 (MQ=255)
                                   aaCAGAGACTGATAACATCAAAGCAGCAAAATAAGGCTTCATTTTTACCACCTTTATCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTCGATTTTTc                                                              <  2:6623028/101‑1 (MQ=255)
                                   aaCAGAGACTGATAACATCAAAGCAGCAAAATAAGGCTTCATTTTTACCACCTTTATCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTCGATTTTTc                                                              >  1:6623573/1‑101 (MQ=255)
                                    aCAGAGACTGATAACATCAAAGCAGCAAAATAAGGCTTCATTTTTACCACCTTTATCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTCGATTTTTcg                                                             >  2:6623576/1‑101 (MQ=255)
                                    aCAGAGACTGATAACATCAAAGCAGCAAAATAAGGCTTCATTTTTACCACCTTTATCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTCGATTTTTcg                                                             <  2:6623348/101‑1 (MQ=255)
                                     cAGAGACTGATAACATCAAAGCAGCAAAATAAGGCTTCATTTTTACCACCTTTATCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTCGATTTTTcgc                                                            >  1:6623577/1‑101 (MQ=255)
                                         gaCTGATAACATCAAAGCAGCAAAATAAGGCTTCATTTTTACCACCTTTATCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTCGATTTTTCGCGCAt                                                        >  2:6623582/1‑101 (MQ=255)
                                         gaCTGATAACATCAAAGCAGCAAAATAAGGCTTCATTTTTACCACCTTTATCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTCGATTTTTCGCGCAt                                                        <  1:6623360/101‑1 (MQ=255)
                                          aCTGATAACATCAAAGCAGCAAAATAAGGCTTCATTTTTACCACCTTTATCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTCGATTTTTCGCGCAtt                                                       <  1:6623449/101‑1 (MQ=255)
                                           cTGATAACATCAAAGCAGCAAAATAAGGCTTCATTTTTACCACCTTTATCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTCGATTTTTCGCGCATTg                                                      <  1:6623311/101‑1 (MQ=255)
                                            tGATAACATCAAAGCAGCAAAATAAGGCTTCATTTTTACCACCTTTATCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTCGATTTTTCGCGCATTGc                                                     >  1:6623585/1‑101 (MQ=255)
                                            tGATAACATCAAAGCAGCAAAATAAGGCTTCATTTTTACCACCTTTATCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTCGATTTTTCGCGCATTGc                                                     >  2:6623586/1‑101 (MQ=255)
                                             gATAACATCAAAGCAGCAAAATAAGGCTTCATTTTTACCACCTTTATCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTCGATTTTTCGCGCATTGCt                                                    >  1:6623590/1‑101 (MQ=255)
                                             gATAACATCAAAGCAGCAAAATAAGGCTTCATTTTTACCACCTTTATCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTCGATTTTTCGCGCATTGCt                                                    <  2:6623098/101‑1 (MQ=255)
                                                aaCATCAAAGCAGCAAAATAAGGCTTCATTTTTACCACCTTTATCAGGTTGCGTTTCGTTTGTTCCAGAGGAACATTGTCGATTTTTCGCGCATTGCtggt                                                 >  1:6623598/1‑101 (MQ=255)
                                                aaCATCAAAGCAGCAAAATAAGGCTTCATTTTTACCACCTTTATCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTCGATTTTTCGCGCATTGCtggt                                                 >  2:6623596/1‑101 (MQ=255)
                                                 aCATCAAAGCAGCAAAATAAGGCTTCATTTTTACCACCTTTATCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTCGATTTTTCGCGCATTGCtggtg                                                >  2:6623601/1‑101 (MQ=255)
                                                 aCATCAAAGCAGCAAAATAAGGCTTCATTTTTACCACCTTTATCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTCGATTTTTCGCGCATTGCtggtg                                                >  1:6623602/1‑101 (MQ=255)
                                                  cATCAAAGCAGCAAAATAAGGCTTCATTTTTACCACCTTTATCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTCGATTTTTCGCGCATTGCtggtgg                                               <  1:6623144/101‑1 (MQ=255)
                                                    tCAAAGCAGCAAAATAAGGCTTCATTTTTACCACCTTTATCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTCg                                                                     >  2:6623606/1‑77 (MQ=255)
                                                     cAAAGCAGCAAAATAAGGCTTCATTTTTACCACCTTTATCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTCGATTTTTCGCGCATTGCTGGTGGCTg                                            <  1:6623271/101‑1 (MQ=255)
                                                     cAAAGCAGCAAAATAAGGCTTCATTTTTACCACCTTTATCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTCGATTTTTCGCGCATTGCTGGTGGCTg                                            >  1:6623608/1‑101 (MQ=255)
                                                      aaaGCAGCAAAATAAGGCTTCATTTTTACCACCTTTATCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTCGATTTTTCGCGCATTGCTGGTGGCTgg                                           <  1:6623535/101‑1 (MQ=255)
                                                      aaaGCAGCAAAATAAGGCTTCATTTTTACCACCTTTATCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTCGATTTTTCGCGCATTGCTGGTGGCTgg                                           >  1:6623614/1‑101 (MQ=255)
                                                       aaGCAGCAAAATAAGGCTTCATTTTTACCACCTTTATCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTCGATTTTTCGCGCATTGCTGGTGGCTggg                                          >  1:6623619/1‑101 (MQ=255)
                                                        agcagcAAAATAAGGCTTCATTTTTACCACCTTTATCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTCGATTTTTCGCGCATTGCTGGTGGCTGGGa                                         >  2:6623620/1‑101 (MQ=255)
                                                         gcagcaAAATAAGGCTTCATTTTTACCACCTTTATCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTCGATTTTTCGCGCATTGCTGGTGGCTGGGaa                                        >  1:6623624/1‑101 (MQ=255)
                                                          cagcaAAATAAGGCTTCATTTTTACCACCTTTATCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTCGATTTTTCGCGCATTGCTGGTGGCTGGGAAt                                       >  1:6623628/1‑101 (MQ=255)
                                                           agcaAAATAAGGCTTCATTTTTACCACCTTTATCAGGTTGCGTTTCAt                                                                                           <  2:6622915/48‑1 (MQ=255)
                                                           agcaAAATAAGGCTTCATTTTTACCACCTTTATCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTCGATTTTTCGCGCATTGCTGGTGGCTGGGAATc                                      <  1:6622747/101‑1 (MQ=255)
                                                           agcaAAATAAGGCTTCATTTTTACCACCTTTATCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTCGATTTTTCGCGCATTGCTGGTGGCTGGGAATc                                      <  1:6623277/101‑1 (MQ=255)
                                                            gcaAAATAAGGCTTCATTTTTACCACCTTTATCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTCGATTTTTCGCGCATTGCTGGTGGCTGGGAATCa                                     >  2:6623632/1‑101 (MQ=255)
                                                             caAAATAAGGCTTCATTTTTACCACCTTTATCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTCGATTTTTCGCGCATTGCTGGTGGCTGGGAATCAc                                    >  1:6623634/1‑101 (MQ=255)
                                                             caAAATAAGGCTTCATTTTTACCACCTTTATCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTCGATTTTTCGCGCATTGCTGGTGGCTGGGAATCAc                                    <  1:6623094/101‑1 (MQ=255)
                                                              aaaaTAAGGCTTCATTTTTACCACCTTTATCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTCGATTTTTCGCGCATTGCTGGTGGCTGGGAATCAcc                                   >  1:6623639/1‑101 (MQ=255)
                                                              aaaaTAAGGCTTCATTTTTACCACCTTTATCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTCGATTTTTCGCGCATTGCTGGTGGCTGGGAATCAcc                                   >  2:6623638/1‑101 (MQ=255)
                                                               aaaTAAGGCTTCATTTTTACCACCTTTATCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTCGATTTTTCGCGCATTGCTGGTGGCTGGGAATCACCt                                  >  1:6623643/1‑101 (MQ=255)
                                                                aaTAAGGCTTCATTTTTACCACCTTTATCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTCGATTTTTCGCGCAt                                                        >  1:6623651/1‑78 (MQ=255)
                                                                aaTAAGGCTTCATTTTTACCACCTTTATCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTCGATTTTTCGCGCATTGCTGGTGGCTGGGAATCACCTg                                 >  1:6623647/1‑101 (MQ=255)
                                                                aaTAAGGCTTCATTTTTACCACCTTTATCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTCGATTTTTCGCGCATTGCTGGTGGCTGGGAATCACCTg                                 >  2:6623650/1‑101 (MQ=255)
                                                                      gCTTCATTTTTACCACCTTTATCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTCGATTTTTCGCGCATTGCTGGTGGCTGGGAATCACCTg                                 >  2:6623655/1‑95 (MQ=255)
                                                                         tCATTTTTACCACCTTTATCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTCGATTTTTCGCGCATTGCTGGTGGCTGGGAATCACCTGAATGGGTGa                        >  1:6623658/1‑101 (MQ=255)
                                                                         tCATTTTTACCACCTTTATCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTCGATTTTTCGCGCATTGCTGGTGGCTGGGAATCACCTGAATGGGTGa                        <  1:6623234/101‑1 (MQ=255)
                                                                          cATTTTTACCACCTTTATCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTCGATTTTTCGCGCATTGCTGGTGGCTGGGAATCACCTGAATGGGTGAt                       >  2:6623661/1‑101 (MQ=255)
                                                                          cATTTTTACCACCTTTATCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTCGATTTTTCGCGCATTGCTGGTGGCTGGGAATCACCTGAATGGGTGAt                       >  1:6623663/1‑101 (MQ=255)
                                                                           aTTTTTACCACCTTTATCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTCGATTTTTCGCGCATTGCTGGTGGCTGGGAATCACCTGAATGGGTGAtt                      <  2:6623017/101‑1 (MQ=255)
                                                                            tttttACCACCTTTATCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTCGATTTTTCGCGCATTGCTGGTGGCTGGGAATCACCTGAATGGGTGAttt                     >  1:6623667/1‑101 (MQ=255)
                                                                             ttttACCCCCTTTATCAGGTTGCGTTTCATTTGTTCCAGAGGAACAt                                                                          <  2:6623047/47‑1 (MQ=255)
                                                                             ttttACCACCTTTATCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTCGATTTTTCGCGCATTGCTGGTGGCTGGGAATCACCTGAATGGGTGAtttt                    >  2:6623669/1‑101 (MQ=255)
                                                                              tttACCACCTTTATCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTCGATTTTTCGCGCATTGCTGGTGGCTGGGAATCACCTGAATGGGTGAttttt                   <  1:6623258/101‑1 (MQ=255)
                                                                              tttACCACCTTTATCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTCGATTTTTCGCGCATTGCTGGTGGCTGGGAATCACCTGAATGGGTGAttttt                   <  1:6623208/101‑1 (MQ=255)
                                                                                tACCACCTTTATCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTCGATTTTTCGCGCATTGCTGGTGGCTGGGAATCACCTGAATGGGTGATTTTTGa                 <  1:6623344/101‑1 (MQ=255)
                                                                                tACCACCTTTATCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTCGATTTTTCGCGCATTGCTGGTGGCTGGGAATCACCTGAATGGGTGATTTTTGa                 >  2:6623673/1‑101 (MQ=255)
                                                                                tACCACCTTTATCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTCGATTTTTCGCGCATTGCTGGTGGCTGGGAATCACCTGAATGGGTGATTTTTGa                 >  2:6623676/1‑101 (MQ=255)
                                                                                 accaccTTTATCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTCGATTTTTCGCGCATTGCTGGTGGCTGGGAATCACCTGAATGGGTGATTTTTGaa                >  2:6623680/1‑101 (MQ=255)
                                                                                  ccaccTTTATCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTCGATTTTTcgc                                                            >  1:6623682/1‑56 (MQ=255)
                                                                                     ccTTTATCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTCGATTTTTCGCGCATTGCTGGTGGCTGGGAATCACCTGAATGGGTGATTTTTGAATTAc            <  2:6623456/101‑1 (MQ=255)
                                                                                     ccTTTATCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTCGATTTTTCGCGCATTGCTGGTGGCTGGGAATCACCTGAATGGGTGATTTTTGAATTAc            >  1:6623688/1‑101 (MQ=255)
                                                                                     ccTTTATCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTCGATTTTTCGCGCATTGCTGGTGGCTGGGAATCACCTGAATGGGTGATTTTTGAATTAc            >  2:6623687/1‑101 (MQ=255)
                                                                                      cTTTATCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTCGATTTTTCGCGCATTGCTGGTGGCTGGGAATCACCTGAATGGGTGATTTTTGAATTAcc           <  1:6623182/101‑1 (MQ=255)
                                                                                       tttATCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTCGATTTTTCGCGCATTGCTGGTGGCTGGGAATCACCTGAATGGGTGATTTTTGAATTACCg          >  2:6623691/1‑101 (MQ=255)
                                                                                        ttATCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTCGATTTTTCGCGCATTGCTGGTGGCTGGGAATCACCTGAATGGGTGATTTTTGAATTACCgg         >  1:6623696/1‑101 (MQ=255)
                                                                                         tATCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTCGATTTTTCGCGCATTGCTGGTGGCTGGGAATCACCTGAATGGGTGATTTTTGAATTACCGGc        <  2:6623256/101‑1 (MQ=255)
                                                                                         tATCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTCGATTTTTCGCGCATTGCTGGTGGCTGGGAATCACCTGAATGGGTGATTTTTGAATTACCGGc        >  1:6623699/1‑101 (MQ=255)
                                                                                           tCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTCGATTTTTCGCGCATTGCTGGTGGCTGGGAATCACCTGAATGGGTGATTTTTGAATTACCGGCtt      <  2:6623289/101‑1 (MQ=255)
                                                                                           tCAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTCGATTTTTCGCGCATTGCTGGTGGCTGGGAATCACCTGAATGGGTGATTTTTGAATTACCGGCtt      <  2:6623214/101‑1 (MQ=255)
                                                                                            cAGGTTGCGTTTCATTTGTTCCAGAGGAACATTGTCGATTTTTCGCGCATTGCTGGTGGCTGGGAATCACCTGAATGGGTGATTTTTGAATTACCGGCttt     <  1:6623197/101‑1 (MQ=255)
                                                                                              ggTTGCGTTTCATTTGTTCCAGAGGAACATTGTCGATTTTTCGCGCATTGCTGTTGGCTGGGAATCACCTGAATGGGTGATTTTTGAATTACCGGCTTTgg   >  1:6623703/1‑101 (MQ=255)
                                                                                               gTTGCGTTTCATTTGTTCCAGAGGAACATTGTCGATTTTTCGCGCATTGCTGGTGGCTGGGAATCACCTGAATGGGTGATTTTTGAATTACCGGCTTTGGt  <  1:6623059/101‑1 (MQ=255)
                                                                                                  |                                                                                                 
TCAGCAGTACCGAGCGGACCAGCATAAGCAGGAAGAACAGAGACTGATAACATCAAAGCAGCAAAATAAGGCTTCATTTTTACCACCTTTATCAGGTTACGTTTCATTTGTTCCAGAGGAACATTGTCGATTTTTCGCGCATTGCTGGTGGCTGGGAATCACCTGAATGGGTGATTTTTGAATTACCGGCTTTGGT  >  NC_000913/4239949‑4240144

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