Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,643,261 G→A 8.3% intergenic (+46/‑77) dbpA → / → yhiN 23S rRNA specific ATP‑dependent RNA helicase/FAD/NAD(P)‑binding domain‑containing oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,643,2610GA 8.3% 171.8 / 4.5 60intergenic (+46/‑77)dbpA/yhiN23S rRNA specific ATP‑dependent RNA helicase/FAD/NAD(P)‑binding domain‑containing oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (27/28);  new base A (2/3);  total (29/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 = 6.93e-01

CAGGCGTTTGTGGCCGTGGAGCGGGCGCTGGCCAAGCAGGCCATGCAGCGGCTGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCT  >  NC_002947/5643114‑5643399
                                                                                                                                                   |                                                                                                                                          
cAGGCGTTTGTGGCCGTGGAGCGGGCGCTGGCCAAGCAGGCCATGCAGCGGCTGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCgg                                                                                                                                            >  2:182875/1‑148 (MQ=255)
 aGGCGTTTGTGGCCGTGGAGCGGGCGCTGGCCAAGCAGGCCATGCAGCGGCTGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCa                                                                                                                                          <  3:169123/149‑1 (MQ=255)
  ggCGTTTGTGGCCGTGGAGCGGGCGCTGGCCAAGCAGGCCATGCAGCGGCTGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAATCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAg                                                                                                                                         >  3:165299/1‑149 (MQ=255)
   gCGTTTGTGTCCGTGGAGCGGGCGCTGGCCAAGCAGGCCATGCAGCGGCTGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTAACCGCGAATGCGGCAGt                                                                                                                                        <  8:223092/149‑1 (MQ=255)
   gCGTTTGTGGCCGTGGAGCGGGCGCTGGCCAAGCAGGCCATGCAGCGGCTGAACAGCGGGAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGt                                                                                                                                        <  1:218721/149‑1 (MQ=255)
      tttGTGGCCGTGGAGCGGGCGCTGGCCAAGCAGGCCATGCAGCGGCTGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGaa                                                                                                                                     >  3:105651/1‑149 (MQ=255)
         gtgGCCGTGGAGCGGGCGCTGGCCAAGCAGGCCATGCAGCGGCTGAACAGCGGCAAGATCAAGGGTCGCTCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATg                                                                                                                                   <  2:239538/148‑1 (MQ=255)
          tgGCCGTGGAGCGGGCGCTGGCCAAGCAGGCCATGCAGCGGCTGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCa                                                                                                                                 >  1:147221/1‑149 (MQ=255)
                  gAGCGGGCGCTGGCCAAGCAGGCCATGCAGCGGCTGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCa                                                                                                                         <  5:6144/149‑1 (MQ=255)
                       ggCGCTGGCCAAGCAGGCCATGCAGCGGCTGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTcgc                                                                                                                    <  2:6945/149‑1 (MQ=255)
                         cgcTGGCCAAGCAGGCCATGCAGCGGCTGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCgg                                                                                                                  <  2:197669/149‑1 (MQ=255)
                         cgcTGGCCAAGCAGGCCATGCAGCGGCTGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCgg                                                                                                                  >  3:99523/1‑149 (MQ=255)
                          gcTGGCCAAGCAGGCCATGCAGCGGCTGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCggg                                                                                                                 >  7:67885/1‑149 (MQ=255)
                          gcTGGCCAAGCAGGCCATGCAGCGGCTGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCggg                                                                                                                 <  7:119273/149‑1 (MQ=255)
                               ccAAGCAGGCCATGCAGCGGCTGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAAc                                                                                                            <  1:182875/149‑1 (MQ=255)
                                   gCAGGCCATGCAGCGGCTGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGc                                                                                                        >  2:72108/1‑149 (MQ=255)
                                             cAGCGGCTGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTGGGAGCGGGTTCACCCGCGAATGCGATGGTGCATTCACTGCCGCATTCGCGGGTAAACCCGCTCCTACAggg                                                                                              <  2:874/149‑1 (MQ=18)
                                             cAGCGGCTGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAATCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAgg                                                                                               <  4:165299/148‑1 (MQ=255)
                                                  gCTGACCAGCGCCAAGATCAAGGGCCGCGCCCTCAAAGTCCGCATTATCTAAACCCAGGCAACCCCCCCTGTAGCAGCGGGTTTACCCGCTAATGCGGCAGTTAATGCACCATCCCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAg                                                                                         <  2:147221/149‑1 (MQ=37)
                                                   cTGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCCCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGCACCCGCTCCCACAGGGTTCAg                                                                                         >  8:271633/1‑148 (MQ=255)
                                                    tGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCa                                                                                       >  3:256324/1‑149 (MQ=255)
                                                        cAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCAttt                                                                                    >  1:163632/1‑148 (MQ=255)
                                                          gCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCaa                                                                                 <  4:105651/149‑1 (MQ=255)
                                                            ggCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCt                                                                                                       >  7:217828/1‑125 (MQ=255)
                                                            ggCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCt                                                                                                       <  8:217828/125‑1 (MQ=255)
                                                            ggCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAAt                                                                               <  1:72108/149‑1 (MQ=255)
                                                              cAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCt                                                                             <  8:67885/149‑1 (MQ=255)
                                                               aaGTTCAAGGGTCGCGCCCTTAATGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTc                                                                                   <  2:101772/142‑1 (MQ=255)
                                                               aaGATCAAGGGTCGCGCCCTGAAAGTCAGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACACGCCAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCAATTc                                                                                   >  1:101772/1‑142 (MQ=255)
                                                                aGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGa                                                                           <  4:256324/149‑1 (MQ=255)
                                                                 gATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCTACAGGGTTCAGCATTTCAAATCTGaa                                                                          <  1:9625/149‑1 (MQ=255)
                                                                   tCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGTAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAg                                                                         <  5:319802/148‑1 (MQ=255)
                                                                   tCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAgg                                                                        >  5:21126/1‑149 (MQ=255)
                                                                      aGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGtt                                                                      >  5:71748/1‑148 (MQ=255)
                                                                         gTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCACGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGAACCATCGCATTAGCGGGTGAACCAGCTCCCCCAGGGTTCAGCATTTCAAATATGAAGGTTCAt                                                                   >  2:277244/1‑148 (MQ=38)
                                                                          tCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCGCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCTTTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATAc                                                                 <  6:21126/149‑1 (MQ=255)
                                                                                 cTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCAc                                                          <  4:66650/149‑1 (MQ=255)
                                                                                    aaaGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGTTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTc                                                        >  5:303773/1‑148 (MQ=255)
                                                                                     aaGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCa                                                      <  8:177432/149‑1 (MQ=255)
                                                                                        tCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCg                                                   >  1:127228/1‑149 (MQ=255)
                                                                                                 tCTAATCCCAGGCAAAGCCCCCTGTCGTAGCGGGTTTACCCGCTAATGCGGCAGTGAATCCACCATCGCATTGGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGatc                                           <  7:271633/148‑1 (MQ=255)
                                                                                                   tAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAg                                                                         >  1:279226/1‑116 (MQ=255)
                                                                                                   tAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAg                                                                         <  2:279226/116‑1 (MQ=255)
                                                                                                      aGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCtt                                     >  7:220249/1‑149 (MQ=255)
                                                                                                        ccAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTgg                                   <  4:99523/149‑1 (MQ=255)
                                                                                                         cAGGCAAAGCCCCGTGTAGGAGTGGTTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTTAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTggc                                  <  2:163632/149‑1 (MQ=255)
                                                                                                            gCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTggcgcc                               >  7:300566/1‑149 (MQ=255)
                                                                                                            gCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTggcgcc                               >  8:156696/1‑149 (MQ=255)
                                                                                                              aaaGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTggcgccgg                             >  4:291316/1‑149 (MQ=255)
                                                                                                                    cccTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGcgccgc                       >  2:171096/1‑149 (MQ=255)
                                                                                                                    cccTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGcgccgc                       >  5:19941/1‑149 (MQ=255)
                                                                                                                              gCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGAtg             >  6:82031/1‑149 (MQ=255)
                                                                                                                              gCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGAtg             >  1:13350/1‑149 (MQ=255)
                                                                                                                               cGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTggcgccggcggcgccggcgccaagg            <  3:291316/149‑8 (MQ=255)
                                                                                                                                   tttACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTgcgc        <  6:19941/149‑1 (MQ=255)
                                                                                                                                    ttACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGcc       >  2:174913/1‑149 (MQ=255)
                                                                                                                                         cgcgAATGCGATGGTGCATTCACTGCCGCATTCGCGGGTAAACCCGCTCCTACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCa     >  3:246504/1‑146 (MQ=25)
                                                                                                                                         cgcgAATGCGATGGTGCATTCACTGCCGCATTCGCGGGTAAACCCGCTCCTACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCa     <  2:276284/146‑1 (MQ=25)
                                                                                                                                         cgcgAATGCGATGGTGCATTCACTGCCGCATTCGCGGGTAAACCCGCTCCTACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCa     <  4:246504/146‑1 (MQ=25)
                                                                                                                                         cgcgAATGCGATGGTGCATTCACTGCCGCATTCGCGGGTAAACCCGCTCCTACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCa     >  1:276284/1‑146 (MQ=25)
                                                                                                                                          gcgAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCt  >  8:112116/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                          
CAGGCGTTTGTGGCCGTGGAGCGGGCGCTGGCCAAGCAGGCCATGCAGCGGCTGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCT  >  NC_002947/5643114‑5643399

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: ‑