Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,643,301 C→T 7.9% intergenic (+86/‑37) 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,3010CT 7.9% 185.9 / 4.6 63intergenic (+86/‑37)dbpA/yhiN23S rRNA specific ATP‑dependent RNA helicase/FAD/NAD(P)‑binding domain‑containing oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (31/27);  new base T (3/2);  total (34/29)
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.75e-01

TGCAGCGGCTGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCTCAGTGCCCGTAGCGGCCGCCGGGTGCTGCTGCTCGACCACGCCAACAAGC  >  NC_002947/5643157‑5643449
                                                                                                                                                |                                                                                                                                                     
tGCAGCGGCTGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAg                                                                                                                                                   >  6:319430/1‑149 (MQ=255)
   aGCGGCTGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGt                                                                                                                                                >  7:259958/1‑149 (MQ=255)
         tGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACACGCGAATGCGGCAGTGAATTCACCATCGCATTCGCGGGTGAACACTCTCCCACAGGGTTctcca                                                                                                                                          >  8:297665/1‑145 (MQ=38)
          gAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCAt                                                                                                                                         >  4:306615/1‑149 (MQ=255)
            aCAGCGGCAAGATCATGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCAtt                                                                                                                                        >  2:115513/1‑148 (MQ=255)
            aCAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCa                                                                                                                                             <  8:105732/143‑1 (MQ=255)
            aCAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCa                                                                                                                                             >  7:105732/1‑143 (MQ=255)
            aCAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCAttt                                                                                                                                       <  5:78107/149‑1 (MQ=255)
             cAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTc                                                                                                                                      <  8:259958/149‑1 (MQ=255)
                 ggCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAAt                                                                                                                                 <  2:69210/149‑1 (MQ=255)
                  gCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATc                                                                                                                                >  5:364270/1‑149 (MQ=255)
                      gATCAAGGGTCGCGCCCTGAAAGTACTCATTATCTATAGCCAGGCAAAGCCCCCTGTATTAGCGGGTTTACCCGCTAATGCTTCAGTGAATTCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGaa                                                                                                                            <  7:297665/149‑1 (MQ=37)
                       aTCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAggg                                                                                                                                                 >  1:63724/1‑128 (MQ=255)
                       aTCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAggg                                                                                                                                                 <  2:63724/128‑1 (MQ=255)
                        tCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCaaa                                                                                                                                   >  7:100238/1‑141 (MQ=255)
                        tCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCaaa                                                                                                                                   <  8:100238/141‑1 (MQ=255)
                        tCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAgg                                                                                                                          <  5:319430/149‑1 (MQ=255)
                          aaGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGtt                                                                                                                        <  4:246180/149‑1 (MQ=255)
                                 gcgcCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCaaa                                                                                                                                   >  6:280698/1‑132 (MQ=255)
                                 gcgcCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCaaa                                                                                                                                   <  5:280698/132‑1 (MQ=255)
                                     ccTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCa                                                                                                             >  4:234452/1‑149 (MQ=255)
                                       tGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACt                                                                                                           >  7:268689/1‑149 (MQ=255)
                                       tGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACt                                                                                                           >  1:314406/1‑149 (MQ=255)
                                          aaGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCAc                                                                                                            >  8:214613/1‑145 (MQ=255)
                                          aaGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCAc                                                                                                            <  7:214613/145‑1 (MQ=255)
                                               cGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGAc                                                                                                   <  3:306615/149‑1 (MQ=255)
                                                 cATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGt                                                                                                 >  3:189802/1‑149 (MQ=255)
                                                        tAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGatcatc                                                                                          >  1:264107/1‑149 (MQ=255)
                                                            gCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTg                                                                                      <  8:268689/149‑1 (MQ=255)
                                                            gCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTg                                                                                      <  2:264107/149‑1 (MQ=255)
                                                             ccAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGaa                                                                                                                            >  4:117546/1‑110 (MQ=255)
                                                             ccAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGaa                                                                                                                            <  3:117546/110‑1 (MQ=255)
                                                               aGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCaaa                                                                                                                                   >  4:325747/1‑102 (MQ=255)
                                                               aGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCaaa                                                                                                                                   <  3:325747/102‑1 (MQ=255)
                                                                 gCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTggcgcc                                                                                 <  1:115513/149‑1 (MQ=255)
                                                                       cccccTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTggcgccggcgcc                                                                           >  3:124207/1‑149 (MQ=255)
                                                                       cccccTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTggcgccggcgcc                                                                           <  6:364270/149‑1 (MQ=255)
                                                                         cccTGTGGGAGCGGGTTCACCCGCGAATGCGATGGTGCATTCACTGCCGCATTCGCGGGTAAACCCGCTCCTACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGcgccgc                                                                         >  1:9380/1‑149 (MQ=18)
                                                                          ccTGTGGGAGCGGGTTCACCCGCGAATGCGATGGTGCATTCACTGCCGCATTCGCGGGTAAACCCGCTCCTACAGGGTTCAGCATTTCAAACTCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTg                                                                                      <  1:72824/136‑1 (MQ=14)
                                                                          ccTGTGGGAGCGGGTTCACCCGCGAATGCGATGGTGCATTCACTGCCGCATTCGCGGGTAAACCCGCTCCTACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTg                                                                                      >  2:72824/1‑135 (MQ=17)
                                                                          ccTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCgccgcc                                                                        >  2:296790/1‑149 (MQ=255)
                                                                            tGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACtt                                                                                                          >  2:91719/1‑112 (MQ=255)
                                                                                ggcgCGGGTTAACCCGCTAATGCGGAAGTGAATGCACCATCGCTTTCGCGGGTTAACCGGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCt                                                                   <  2:314406/145‑1 (MQ=255)
                                                                                 gAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGa                                                                 >  7:30824/1‑149 (MQ=255)
                                                                                  aGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGAt                                                                <  3:234452/149‑1 (MQ=255)
                                                                                    cGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGAtgt                                                              <  4:189802/149‑1 (MQ=255)
                                                                                    cGGGTTCACCCGCGAATGCGATGGTGCATTCACTGCCGCATTCGCGGGTAAACCCGCTCCTACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGAtgt                                                              <  2:9380/149‑1 (MQ=21)
                                                                                      ggTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTgc                                                            <  1:296790/149‑1 (MQ=255)
                                                                                             ccGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTc                                                                                                                                      >  6:271300/1‑69 (MQ=255)
                                                                                             ccGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTc                                                                                                                                      <  5:271300/69‑1 (MQ=255)
                                                                                              cgcgAATGCGATGGTGCATTCACTGCCGCATTCGCGGGTAAACCCGCTCCTACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCt                                                    >  5:2797/1‑149 (MQ=25)
                                                                                               gcgAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCt                                                    >  5:130527/1‑148 (MQ=255)
                                                                                                             gAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTggc                                                                                    <  6:8305/102‑1 (MQ=255)
                                                                                                             gAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTggc                                                                                    >  5:8305/1‑102 (MQ=255)
                                                                                                                     cATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCTCAGTGCCCGTAGCGGCCGCCggg                             >  7:224460/1‑149 (MQ=255)
                                                                                                                     cATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCTCAGTGCCCGTAGCGGCCGCCgg                              <  8:224460/148‑1 (MQ=255)
                                                                                                                     cATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCACCATATCAAATCTGAAGCTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCTCAGTGCCCCTAGCGGCCGCCggg                             >  7:163755/1‑149 (MQ=255)
                                                                                                                       ccGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTg                                                                                      <  8:34582/89‑1 (MQ=255)
                                                                                                                       ccGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTg                                                                                      >  7:34582/2‑90 (MQ=255)
                                                                                                                          cATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAg                                                      >  5:308152/1‑119 (MQ=255)
                                                                                                                            ttCGCGGGTGAACCCGCTACCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCTCAGTGCCCGTAGCGGCCGCCGGGtgctgct                      <  8:163755/149‑1 (MQ=255)
                                                                                                                                      aaCCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCTCAGTGCCCGTAGCGGCCGCCGGGTGCTGCTGCTCGACCAc            >  8:159271/1‑149 (MQ=255)
                                                                                                                                         ccGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCTCAGTGCCCGTAGCGGCCGCCGGGTGCTGCTGCTCGACCACGc          >  6:52935/1‑148 (MQ=255)
                                                                                                                                                cacaGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCTCAGTGCCCGTAGCGGCCGCCGGGTGCTGCTGCTCGACCAc            <  6:46714/139‑1 (MQ=255)
                                                                                                                                                cacaGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCTCAGTGCCCGTAGCGGCCGCCGGGTGCTGCTGCTCGACCAc            >  5:46714/1‑139 (MQ=255)
                                                                                                                                                cacaGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCTCAGTGCCCGTAGCGGCCGCCGGGTGCTGCTGCTCGACCACGCCAACAAGc  >  6:306877/1‑149 (MQ=255)
                                                                                                                                                cacaGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCTCAGTGCCCGTAGCGGCCGCCGGGTGCTGCTGCTCGACCACGCCAACAAGc  <  4:124207/149‑1 (MQ=255)
                                                                                                                                                cacaGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCTCAGTGCCCGTAGCGGCCGCCGGGTGCTGCTGCTCGACCACGCCAACAAGc  >  2:153291/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
TGCAGCGGCTGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCTCAGTGCCCGTAGCGGCCGCCGGGTGCTGCTGCTCGACCACGCCAACAAGC  >  NC_002947/5643157‑5643449

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: