Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 546,014 T→C 14.5% intergenic (+54/‑19) 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,0140TC14.5% 143.7 / 15.5 62intergenic (+54/‑19)rpoC/PP_mr08DNA‑directed RNA polymerase subunit beta'/rpsL_psuedo ncRNA
Reads supporting (aligned to +/- strand):  ref base T (26/27);  new base C (5/4);  total (31/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 = 9.48e-01

GCAAGCGTCGCCGTGATGCCGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTGCAAGACAATAGAGTCGCAAGACAATCA  >  NC_002947/545866‑546158
                                                                                                                                                    |                                                                                                                                                
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                                                              >  8:202025/1‑149 (MQ=255)
                                                                                    gCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGt                                                              <  3:67772/149‑1 (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:102807/1‑147 (MQ=255)
                                                                                                                            gCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTGCAAGa                        <  2:102807/147‑1 (MQ=255)
                                                                                                                            gCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTGCAAGa                        >  1:98836/1‑147 (MQ=255)
                                                                                                                            ccGATCATCGGCAACATTTCATGTCGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAAttttt                                       <  5:211163/131‑1 (MQ=38)
                                                                                                                            ccGATCATCGGCAACATTTCATGTCGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAAttttt                                       >  6:211163/2‑132 (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)
                                                                                                                                               aaTGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGAGTGGGTGCAATATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTGCAAGACAATAGAGTCGCAAGACAAt    <  3:302427/148‑1 (MQ=255)
                                                                                                                                               aaTGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTGCAAGACAATAGAGTCGCAAGACAAt    >  2:324989/1‑148 (MQ=255)
                                                                                                                                                aTGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTGCAAGACAATAGAGTCGCAAGACAATCa  <  4:57414/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
GCAAGCGTCGCCGTGATGCCGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTGCAAGACAATAGAGTCGCAAGACAATCA  >  NC_002947/545866‑546158

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: