Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 309,846 G→A 6.3% P61S (CCG→TCG)  ecpB ← ECP production pilus chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913309,8460GA 6.3% 197.2 / 3.0 63P61S (CCG→TCG) ecpBECP production pilus chaperone
Reads supporting (aligned to +/- strand):  ref base G (30/29);  new base A (2/2);  total (32/31)
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.06e-01

CGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCG  >  NC_000913/309701‑309988
                                                                                                                                                 |                                                                                                                                              
cGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTg                                                                                                                                              <  8:102088/148‑1 (MQ=255)
cGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTg                                                                                                                                              <  7:312329/148‑1 (MQ=255)
cGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGa                                                                                                                                             <  5:329827/149‑1 (MQ=255)
 ggTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTg                                                                                                                                              >  5:24165/1‑147 (MQ=255)
 ggTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGAGCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGaa                                                                                                                                            <  4:166006/149‑1 (MQ=255)
   tagtagCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGTATCACTTCGCTGGCTTGGGCGGGGAGCAGCAAGCTGGCGGGAGTGAGTAGCAACTCATCCGGATTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAaga                                                                                                                                          <  3:201467/149‑1 (MQ=255)
          gcTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGTCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCAACTCCTCCGGCTTGTCCATTGAGCTAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGt                                                                                                                                   >  3:245646/1‑149 (MQ=255)
              tttttCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTc                                                                                                                                 >  6:112403/1‑147 (MQ=255)
                  tCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGAc                                                                                                                            <  5:122755/148‑1 (MQ=255)
                    aTCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGa                                                                                                                         <  8:46202/149‑1 (MQ=255)
                     tCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGa                                                                                                                         >  4:136376/1‑148 (MQ=255)
                       tGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGatat                                                                                                                      >  3:200824/1‑149 (MQ=255)
                         cGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATAtt                                                                                                                     >  1:190588/1‑148 (MQ=255)
                         cGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATAtt                                                                                                                     <  5:112403/148‑1 (MQ=255)
                                 ttATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGCTAAGGc                                                                                                             >  5:103346/1‑148 (MQ=255)
                                      gaagaaGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCa                                                                                                        <  3:19139/148‑1 (MQ=255)
                                      gaagaaGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGATGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCAc                                                                                                       >  2:231206/1‑149 (MQ=255)
                                      gaagaaGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGATGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCAc                                                                                                       >  2:231207/1‑149 (MQ=255)
                                              ggATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCAGTGAGATACCCTGCCCGGCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGCCCACTGGCGGTa                                                                                               >  5:227583/1‑149 (MQ=255)
                                               gATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTAc                                                                                              >  8:105547/1‑149 (MQ=255)
                                               gATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCAGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTAc                                                                                              >  4:117699/1‑149 (MQ=255)
                                               gATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGATGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCATTGATATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTAc                                                                                              >  4:158056/1‑149 (MQ=255)
                                                 tCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTg                                                                                            <  8:219335/149‑1 (MQ=255)
                                                  cACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCAACTCATCCGGCGTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTTATAAGGCGACCACTGTCGGTACTg                                                                                            <  8:304197/148‑1 (MQ=255)
                                                    cTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTGCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGCGCCGTTCGAGACGGATATTGATAAGGCGAGCACTGTCGGTACGGttt                                                                                         >  8:84585/1‑149 (MQ=255)
                                                          tGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATc                                                                                   <  4:200824/149‑1 (MQ=255)
                                                              ttGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGttt                                                                               >  7:216686/1‑149 (MQ=255)
                                                                gggCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGttt                                                                               <  3:136376/147‑1 (MQ=255)
                                                                 ggCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGc                                                                            <  8:57133/149‑1 (MQ=255)
                                                                    gggTAGCAGCAAGCTGGCGGGAGTTTGTAGCAACTCCTCCTGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCa                                                                         <  4:257748/149‑1 (MQ=255)
                                                                    gggTAGCAGCAAGCTGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCt                                                                           >  2:178535/1‑147 (MQ=255)
                                                                      gTAGCAGCAAGCTGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGc                                                                       >  8:42924/1‑149 (MQ=255)
                                                                      gTAGCAGCAAGCTGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGc                                                                       <  6:8348/149‑1 (MQ=255)
                                                                         gcagcaAGCTGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGtg                                                                    <  6:24165/149‑1 (MQ=255)
                                                                                    gCGGGAGTGAGTAGCAACTCATCCGGCTTGGCCATTGCGAGAACCTGCCCGTCGGCCAGCGGTGAAGAGAGCCGTTCGAGACGGATCTTGATAAGGCGACCACTGTCGGTACTGGTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCa                                                           <  6:227583/147‑1 (MQ=255)
                                                                                      gggAGTGAGTAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGt                                                       <  1:178535/149‑1 (MQ=255)
                                                                                              gtAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGttt                                                                                         <  4:301383/107‑1 (MQ=255)
                                                                                              gtAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGttt                                                                                         >  3:301383/1‑107 (MQ=255)
                                                                                              gtAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAAc                                               >  7:34453/1‑149 (MQ=255)
                                                                                                 gCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACga                                             <  3:198184/148‑1 (MQ=255)
                                                                                                 gCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGGGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACga                                             <  4:108021/148‑1 (MQ=255)
                                                                                                 gCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGATGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACga                                             <  1:231207/148‑1 (MQ=255)
                                                                                                 gCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGATGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACga                                             <  1:231206/148‑1 (MQ=255)
                                                                                                  cAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACgatg                                           <  5:262832/149‑1 (MQ=255)
                                                                                                  cAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACgat                                            >  6:262832/1‑148 (MQ=255)
                                                                                                    aCTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGGGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACgctgat                                         >  5:292614/1‑149 (MQ=255)
                                                                                                        aTCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGatatat                                     >  6:188027/1‑149 (MQ=255)
                                                                                                        aTCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGatatat                                     <  6:136165/149‑1 (MQ=255)
                                                                                                        aTCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGatatat                                     >  3:297494/1‑149 (MQ=255)
                                                                                                                 gTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGACGTTCGAGACGGATATTGATAAGGCGTCCACTGTCGGTACTGTGTTTGAGCGTTTTGCTCATCGTGCTGCTGTCACTGTTCATACACGATGATATATCGCCGACAt                            >  7:110782/1‑149 (MQ=255)
                                                                                                                  tCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATc                           <  7:171770/149‑1 (MQ=255)
                                                                                                                    cATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCa                         <  3:158056/149‑1 (MQ=255)
                                                                                                                     aTTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGATCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTTCTCAGCGTGCTGATGTCACTGTTCATAAACGATGATATATCGCCGACATCCAg                        >  4:20817/1‑149 (MQ=255)
                                                                                                                      ttGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAgc                       >  6:313326/1‑149 (MQ=255)
                                                                                                                             aaCCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTggg                >  4:65807/1‑149 (MQ=255)
                                                                                                                                cTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCg             >  1:137310/1‑149 (MQ=255)
                                                                                                                                   cccGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGc               <  4:290904/144‑1 (MQ=255)
                                                                                                                                   cccGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGc               >  3:290904/1‑144 (MQ=255)
                                                                                                                                   cccGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGga           <  4:2141/148‑1 (MQ=255)
                                                                                                                                    ccGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCg             >  3:95212/1‑145 (MQ=255)
                                                                                                                                       tcgtcAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGatat      >  1:125324/1‑149 (MQ=255)
                                                                                                                                          tcAAGCGGTGAAGCGAGCCGTTCGAGACGGATATTGATAAGGCGCCCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTcc   >  6:276366/1‑149 (MQ=255)
                                                                                                                                            aaGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCg  <  7:84585/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
CGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCG  >  NC_000913/309701‑309988

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