Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,629,003 T→C D179G (GAC→GGC)  rbbA ← ribosome‑associated ATPase: ATP‑binding protein/ATP‑binding membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,629,0030TC96.7% 198.3 / NA 60D179G (GAC→GGC) rbbAribosome‑associated ATPase: ATP‑binding protein/ATP‑binding membrane protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (29/29);  minor base G (2/0);  total (31/29)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGCACTTCTCCGGCATTCATCGCTACCAGCCAGTCGAAGCGTTCGGCCTCTTCCATATAGGCGGTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGTCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCTTCATCCCGCCGGAGAGTTTCCCTGCCGGACGATCGCGAAACGGTGCTAACCCGGTGCTGGTCAGCA  >  NC_000913/3628858‑3629147
                                                                                                                                                 |                                                                                                                                                
cAGCACTTCTCCGGCATTCATCGCTACCAGCCAGTCGAAGCGTTCGGCCTCTTCCATATAGGCGGTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAAcc                                                                                                                                             >  2:383890/1‑151 (MQ=255)
          ccGGCATTCATCGCTACCAGCCAGTCGAAGCGTTCGGCCTCTTCCATATAGGCGGTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGttgttg                                                                                                                                   <  1:543080/151‑1 (MQ=255)
            ggCATTCATCGATACCAGCCAGTCGAAGCGTTCGGCCTCTTCCATATAGGCGGTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGc                                                                                                                                 <  1:631018/151‑1 (MQ=255)
             gCATTCATCGCTACCAGCCAGTCGAAGCGTTCGGCCTCTTCCATATAGGCGGTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGc                                                                                                                                 >  1:79293/1‑150 (MQ=255)
                       cTACCAGCCAGTCGAAGCGTTCGGCCTCTTCCATATAGGCGGTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTTGCTCATCAAGGa                                                                                                                       >  1:657657/1‑150 (MQ=255)
                       cTACCAGCCAGTCGAAGCGTTCGGCCTCTTCCATATAGGCGGTGGCGACCAGCACGCTCATATTGCTCTGCCGCCGGCGAATACTGTCGATCAAAACCCAAAAATGGGGGCGGGAGAGCGGGGCAACCCCCGTTGTTGGCTCATCAAGGAt                                                                                                                      >  2:92917/1‑151 (MQ=255)
                          ccagccagTCGAAGCGTTCGGCCTCTTCCATATAGGCGGTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCa                                                                                                                    <  2:572441/150‑1 (MQ=255)
                           cagccagTCGAAGCGTTCGGCCTCTTCCATATAGGCGGTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCaa                                                                                                                   >  2:143612/1‑150 (MQ=255)
                           cagccagTCGAAGCGTTCGGCCTCTTCCATATAGGCGGTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAAc                                                                                                                  <  1:176761/151‑1 (MQ=255)
                            agccagTCGAAGCGTTCGGCCTCTTCCATATAGGCGGTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAAc                                                                                                                  >  2:361223/1‑150 (MQ=255)
                             gccagTCGAAGCGTTCGGCCTCTTCCATATAGGCGGTGGCGGCCCACACGCTCATATTGCTCTGCCGCTTGCGAAAACTGTCGATCAGATCCCCGAACTGGGAGCGGGAGAGCGGGGCAAACCCCGTTGTTGGCTCATCAAGGATCAAc                                                                                                                  >  2:234885/1‑149 (MQ=255)
                               cagTCGAAGCGTTCGGCCTCTTCCATATAGGCGGTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGtt                                                                                                              >  2:124586/1‑151 (MQ=255)
                                agTCGAAGCGTTCGGCCTCTTCCATATAGGCGGTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGtt                                                                                                              >  2:521033/1‑150 (MQ=255)
                                agTCGAAGCGTTCGGCCTCTTCCATATAGGCGGTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTc                                                                                                             <  1:620309/151‑1 (MQ=255)
                                 gTCGAAGCGTTCGGCCTCTTCCATATAGGCGGTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTcc                                                                                                            <  1:42474/151‑1 (MQ=255)
                                  tCGAAGCGTTCGGCCTCTTCCATATAGGCGGTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTcc                                                                                                            >  1:134809/1‑150 (MQ=255)
                                   cGAAGCGTTCGGCCTCTTCCATATAGGCGGTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCg                                                                                                           <  2:32620/150‑1 (MQ=255)
                                   cGAAGCGTTCGGCCTCTTCCATATAGGCGGTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCg                                                                                                           <  2:85338/150‑1 (MQ=255)
                                        cGTTCGGCCTCTTCCATATAGGCGGTGGCGACCAGCACGCCCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGt                                                                                                     >  1:133924/1‑151 (MQ=255)
                                                   ttCCATATAGGCGGTGGCGACCAGCACGCTCCTATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAAcgc                                                                                           >  1:173715/1‑150 (MQ=255)
                                                    tCCATATAGGCGGTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAAcgc                                                                                           <  1:412330/149‑1 (MQ=255)
                                                    tCCATATAGGCGGTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAAcg                                                                                            <  2:145616/148‑1 (MQ=255)
                                                         ataGGCGGTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGcac                                                                                    <  1:442087/151‑1 (MQ=255)
                                                           aGGCGGTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGcaca                                                                                   <  1:251298/150‑1 (MQ=255)
                                                             gCGGTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGcc                                                                                >  2:381158/1‑151 (MQ=255)
                                                               ggTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCa                                                                              <  1:630952/151‑1 (MQ=255)
                                                                   gCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGtt                                                                          <  2:170176/151‑1 (MQ=255)
                                                                    cGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGttt                                                                         <  2:433662/151‑1 (MQ=255)
                                                                        cAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGc                                                                     >  1:538571/1‑151 (MQ=255)
                                                                        cAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGc                                                                     >  1:38581/1‑151 (MQ=255)
                                                                          gCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCtt                                                                   >  1:440018/1‑151 (MQ=255)
                                                                                    atTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCTTCATCccgccg                                                         <  1:167292/151‑1 (MQ=255)
                                                                                     ttGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCTTCATCCCGCCgg                                                        >  2:561587/1‑151 (MQ=255)
                                                                                         tctGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCTTCATCCCGCCGgag                                                      >  2:493606/1‑149 (MQ=255)
                                                                                          ctGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCTTCATCCCGCCGgaga                                                     <  1:124589/149‑1 (MQ=255)
                                                                                            gCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCTTCATCCCGCCGGAGAGt                                                   <  1:646807/149‑1 (MQ=255)
                                                                                                  ggCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCTTCATCCCGCCGGAGAGTTTCCCt                                             <  2:555309/149‑1 (MQ=255)
                                                                                                  ggCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCTTCATCCCGCCGGAGAGTTTCCCt                                             <  1:521416/149‑1 (MQ=255)
                                                                                                     gAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCTTCATCCCGCCGGAGAGTTTCCCTGCCgg                                        <  1:326444/151‑1 (MQ=255)
                                                                                                     gAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCTTCATCCCGCCGGAGAGTTTCCCTGCCgg                                        >  2:60182/1‑151 (MQ=255)
                                                                                                        tACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCTTCATCCCGCCGGAGAGTTTCCCTGCCGGACg                                     >  1:538427/1‑151 (MQ=255)
                                                                                                         aCTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCTTCATCCCGCCGGAGAGTTTCCCTGCCGGACGa                                    <  2:442797/151‑1 (MQ=255)
                                                                                                          cTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCTTCATCCCGCCGGAGAGTTTCCCTGCCGGACGa                                    <  1:97694/150‑1 (MQ=255)
                                                                                                            gTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCTTCATCCCGCCGGAGAGTTTCCCTGCCGGACGATcg                                 >  1:545028/1‑151 (MQ=255)
                                                                                                            gTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCTTCATCCCGCCGGAGAGTTTCCCTGCCGGACGATcg                                 <  1:678823/151‑1 (MQ=255)
                                                                                                             tCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCTTCATCCCGCCGGAGAGTTTCCCTGCCGGACGATcgc                                >  1:506674/1‑151 (MQ=255)
                                                                                                              cGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCTTCATCCCGCCGGAGAGTTTCCCTGCCGGACGATcgcg                               <  2:660006/151‑1 (MQ=255)
                                                                                                                 tCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCTTCATCCCGCCGGAGAGTTTCCCTGCCGGACGATCGCGaa                             >  2:362293/1‑150 (MQ=255)
                                                                                                                    gATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCTTCATCCCGCCGGAGAGTTTCCCTGCCGGACGATCGCGAAACgg                         <  2:356395/151‑1 (MQ=255)
                                                                                                                      tCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCTTCATCCCGCCGGAGAGTTTCCCTGCCGGACGATCGCGAAACGGTg                       >  2:489789/1‑151 (MQ=255)
                                                                                                                        ccAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCTTCATCCCGCCGGAGAGTTTCCCTGCCGGACGATCGCGAAACGGTGCt                     >  1:404232/1‑151 (MQ=255)
                                                                                                                         cAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCTTCATCCCGCCGGAGAGTTTCCCTGCCGGACGATCGCGAAACGGTGCt                     >  1:542858/1‑150 (MQ=255)
                                                                                                                          aGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCTTCATCCCGCCGGAGAGTTTCCCTGCCGGACGATCGCGAAACGGTGCTa                    >  1:103694/1‑150 (MQ=255)
                                                                                                                            aaCTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCTTCATCCCGCCGGAGAGTTTCCCTGCCGGACGATCGCGAAACGGTGCTAAcc                 >  2:193423/1‑151 (MQ=255)
                                                                                                                                gggAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCTTCATCCCGCCGGAGAGTTTCCCTGCCGGACGATCGCGAAACGGTGCTAACCCgg              >  2:453149/1‑150 (MQ=255)
                                                                                                                                  gAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCTTCATCCCGCCGGAGAGTTTCCCTGCCGGACGATCGCGAAACGGTGCTAACCCGGTGc           >  2:335072/1‑151 (MQ=255)
                                                                                                                                   aGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCTTCATCCCGCCGGAGAGTTTCCCTGCCGGACGATCGCGAAACGGTGCTAACCCGGTGCt          >  1:677254/1‑151 (MQ=255)
                                                                                                                                       ggAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCTTCATCCCGCCGGAGAGTTTCCCTGCCGGACGATCGCGAAACGGTGCTAACCCGGTGCTGGTc      <  2:53785/151‑1 (MQ=255)
                                                                                                                                       ggAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGACCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTGTTGCTTCATCCCGCCGGAGAGTTTCCCTGCCGGACGATCGCGAAACGGTGCTAACCCGGTGCTGGTc      <  2:535912/151‑1 (MQ=255)
                                                                                                                                           agCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCTTCATCCCGCCGGAGAGTTTCCCTGCCGGACGATCGCGAAACGGTGCTAACCCGGTGCTGGTcagca  <  2:44876/151‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
CAGCACTTCTCCGGCATTCATCGCTACCAGCCAGTCGAAGCGTTCGGCCTCTTCCATATAGGCGGTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGTCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCTTCATCCCGCCGGAGAGTTTCCCTGCCGGACGATCGCGAAACGGTGCTAACCCGGTGCTGGTCAGCA  >  NC_000913/3628858‑3629147

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑