Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 546,008 A→T 14.8% intergenic (+48/‑25) rpoC → / → PP_mr08 DNA‑directed RNA polymerase subunit beta'/rpsL_psuedo ncRNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947546,0080AT14.8% 151.2 / 15.0 60intergenic (+48/‑25)rpoC/PP_mr08DNA‑directed RNA polymerase subunit beta'/rpsL_psuedo ncRNA
Reads supporting (aligned to +/- strand):  ref base A (26/25);  new base T (5/4);  total (31/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 = 1.00e+00

GCAAGCGTCGCCGTGATGCCGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTGCAAGACAATAGAGTCGCAAGACAA  >  NC_002947/545866‑546155
                                                                                                                                              |                                                                                                                                                   
gcAAGCGTCGCCGTGATGCCGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGtt                                                                                                                                               <  3:63000/149‑1 (MQ=255)
gcAAGCGTCGCCGTGATGCCGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGtt                                                                                                                                               >  4:63000/1‑149 (MQ=255)
     cGTCGCCGTGATGCCGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCga                                                                                                                                          <  8:40913/149‑1 (MQ=255)
         gCCGTGATGCCGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCgatgat                                                                                                                                      >  8:110174/1‑149 (MQ=255)
          ccGTGATGCCGACAAACCGCTGCGAGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATc                                                                                                                                     <  6:358010/149‑1 (MQ=255)
            gTGATGCCGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCCCTCCCCGCCCGATCATCGGCAACATTTCATGTCGCCGATGATcgc                                                                                                                                   <  2:56231/149‑2 (MQ=25)
                    gACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCAACATTTCATGTCGCCga                                                                                                                                          <  1:337291/134‑1 (MQ=255)
                    gACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCAACATTTCATGTCGCCga                                                                                                                                          >  2:337291/1‑134 (MQ=255)
                     aCAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAgg                                                                                                                          <  8:199751/149‑1 (MQ=255)
                       aaaCCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGc                                                                                                                        >  8:177341/1‑149 (MQ=255)
                            gCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCCCTCCCCGCCCGATCATCGGCAACATTTCATGTCGCCgat                                                                                                                                         <  3:214706/127‑1 (MQ=21)
                            gCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCCCTCCCCGCCCGATCATCGGCAACATTTCATGTCGCCgat                                                                                                                                         >  4:214706/1‑127 (MQ=21)
                                  tgtAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGt                                                                                                             >  4:302427/1‑149 (MQ=255)
                                   gtAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCggg                                                                                                                    <  5:183300/141‑1 (MQ=255)
                                   gtAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCggg                                                                                                                    >  6:183300/1‑141 (MQ=255)
                                      aGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCtt                                                                                                          <  7:110174/148‑1 (MQ=255)
                                           cAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCg                                                                                                                      >  3:3655/1‑131 (MQ=255)
                                                  gTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTc                                                                                                            <  3:226647/134‑1 (MQ=255)
                                                  gTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTc                                                                                                            >  4:226647/1‑134 (MQ=255)
                                                   tGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGc                                                                                                >  4:258052/1‑145 (MQ=255)
                                                   tGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGa                                                                                            <  7:177341/149‑1 (MQ=255)
                                                          tGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCtt                                                                                     >  2:106101/1‑149 (MQ=255)
                                                           gCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCtt                                                                                     >  3:25222/1‑148 (MQ=255)
                                                            cACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTa                                                                                   <  1:106101/149‑1 (MQ=255)
                                                             aCTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCtt                                                                                     >  1:213565/1‑146 (MQ=255)
                                                                         gcTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTAcc                                                                       >  5:297114/1‑148 (MQ=255)
                                                                         gcTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTAcc                                                                       <  6:297114/148‑1 (MQ=255)
                                                                            gAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCttt                                                                                    <  7:10855/132‑1 (MQ=255)
                                                                            gAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCttt                                                                                    >  8:10855/1‑132 (MQ=255)
                                                                                 ccAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAAttt                                                              >  4:67772/1‑149 (MQ=255)
                                                                                  cAGCGGTAATTACGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTg                                                             <  8:302053/149‑1 (MQ=255)
                                                                                  cAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTg                                                             >  1:319407/1‑149 (MQ=255)
                                                                                    gCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGt                                                           <  3:67772/149‑1 (MQ=255)
                                                                                    gCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGt                                                           >  8:202025/1‑149 (MQ=255)
                                                                                       gTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCtt                                                                                     >  7:329949/1‑120 (MQ=255)
                                                                                       gTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCtt                                                                                     <  8:329949/120‑1 (MQ=255)
                                                                                        tAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAgg                                                       <  6:173598/149‑1 (MQ=255)
                                                                                           ttAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATACTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTc                                                    >  4:89992/1‑149 (MQ=255)
                                                                                                tACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGtct                                               >  2:307686/1‑149 (MQ=255)
                                                                                                     gacgaGGCCCCGGCCCTCCCCGCCCGATCATCGGCAACATTTAATGTCGCCGATGATCGCCCGAGGCTTGACGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGcc                                           >  4:322464/5‑148 (MQ=14)
                                                                                                           gCCGCGGCAAGCCTCGGGCTATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAAttttt                                    <  7:202025/149‑1 (MQ=255)
                                                                                                            ccccGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCaa                                         <  2:7566/143‑1 (MQ=255)
                                                                                                            ccccGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCaa                                         >  1:7566/1‑143 (MQ=255)
                                                                                                              ccGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGc                                 >  4:162227/1‑149 (MQ=255)
                                                                                                                   aaGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCt                            >  1:57313/1‑149 (MQ=255)
                                                                                                                    aGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCtt                           <  3:162227/149‑1 (MQ=255)
                                                                                                                       gtCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGTGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGGGAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTGCa                        <  2:213565/148‑1 (MQ=255)
                                                                                                                            gCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTGCAAGa                     <  2:98836/147‑1 (MQ=255)
                                                                                                                            gCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTGCAAGa                     >  1:98836/1‑147 (MQ=255)
                                                                                                                            gCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTGCAAGa                     <  2:102807/147‑1 (MQ=255)
                                                                                                                            gCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTGCAAGa                     >  1:102807/1‑147 (MQ=255)
                                                                                                                            ccGATCATCGGCAACATTTCATGTCGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAAttttt                                    >  6:211163/2‑132 (MQ=38)
                                                                                                                            ccGATCATCGGCAACATTTCATGTCGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAAttttt                                    <  5:211163/131‑1 (MQ=38)
                                                                                                                              gATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTgg                                  <  8:175502/132‑1 (MQ=255)
                                                                                                                               atcatcGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTGCAAGACAATa                <  1:307686/149‑1 (MQ=255)
                                                                                                                                   tcGGCAACATTTCATGTCGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTGCAAGACAATAGAGt            >  6:137106/1‑149 (MQ=255)
                                                                                                                                     ggCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTGCAAGACAATAGAGTc           <  2:319407/148‑1 (MQ=255)
                                                                                                                                          cATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTGCAAGACAATAGAGTCGCAAGa     <  3:258052/149‑1 (MQ=255)
                                                                                                                                            ttaaaTGTTCCCGATGTTGGGGGGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTGCAAGACAATAGAGTCGCAAGAc    <  4:33902/146‑1 (MQ=255)
                                                                                                                                             gAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTGCAAGACAATAGAGTCGCAAGAc    >  3:33902/1‑147 (MQ=255)
                                                                                                                                             gAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTGCAAGACAATAGAGTCGCAAGACaa  >  1:344860/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
GCAAGCGTCGCCGTGATGCCGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTGCAAGACAATAGAGTCGCAAGACAA  >  NC_002947/545866‑546155

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: