Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 543,727 A→G 6.7% N656S (AAC→AGC)  rpoC → DNA‑directed RNA polymerase subunit beta'

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947543,7270AG 6.7% 212.4 / 4.8 74N656S (AAC→AGC) rpoCDNA‑directed RNA polymerase subunit beta'
Reads supporting (aligned to +/- strand):  ref base A (32/37);  new base G (3/2);  total (35/39)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.54e-01

CAGTGCTACCGCGTGGTTGGTCTGAAAGAGACCGTTATCTTCGCCGACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGAAGTGTCCAAGGCGATGATGGCCAACCTCTCGAAAGAGAA  >  NC_002947/543579‑543874
                                                                                                                                                    |                                                                                                                                                    
cAGTGCTACCGCGTGGTTGGTCTGAAAGAGACCGTTATCTTCGCCGACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTaa                                                                                                                                                      >  8:276676/1‑149 (MQ=255)
     cTACCGCGTGGTTGGTCTGAAAGAGACCGTTATCTTCGCCGACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCt                                                                                                                                                 >  7:59563/1‑148 (MQ=255)
     cTACCGCGTGGTTGGTCTGAAAGAGACCGTTATCTTCGCCGACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCt                                                                                                                                                 <  8:59563/148‑1 (MQ=255)
      ctgggtgTGGTTGGTCTGAAAGAGACCGTTATCTTAGCCGACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAACGCCCGCATCATCGGTAACGATAc                                                                                                                                               >  2:205680/7‑149 (MQ=255)
        ccGCGTGGTTGGTCTGAAAGAGACCGTTATCTTCGCCGACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCg                                                                                                                                             <  2:321905/149‑1 (MQ=255)
         cgcgTGGTTGGTCTGAAAGAGACCGTTATCTTCGCCGACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGa                                                                                                                                            >  7:285332/1‑149 (MQ=255)
          gcgTGGTTGGTCTGAAAGAGACCGTTATCTTCGCCTATCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCTGCGTTTCCATCGGTGTTAACGAGTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGTTAACGCTACCGAt                                                                                                                                           <  3:95604/149‑1 (MQ=255)
          gcgTGGTTGGTCTGAAAGAGACCGTTATCTTCGCCGACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGAt                                                                                                                                           >  8:190284/1‑149 (MQ=255)
           cgTGGTTGGTCTGAAAGAGACCGTTATCTTCGCCGACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGAtg                                                                                                                                          <  5:3334/149‑1 (MQ=255)
             tggttggtCTGAAAGAGACCGTTATCTTCGCCGACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAGCGTTACCGAtgat                                                                                                                                        <  3:123534/149‑2 (MQ=255)
               gttggtCTGAAAGAGACCGTTATCTTCGCCGACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGAtgaagt                                                                                                                                      <  7:276676/149‑1 (MQ=255)
                ttggtCTGAAAGAGACCGTTATCTTCGCCGACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAGCGTTACCGAtgatgcg                                                                                                                                     <  5:66045/149‑5 (MQ=255)
                   gtCTGAAAGAGACCGTTATCTTCGCCGACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGAtgaagtgaa                                                                                                                                   <  8:104155/148‑1 (MQ=255)
                     cTGAAAGAGACCGTTATCTTCGCCGACCAACTGATCTACACCGCTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCACCATCATCCGTAACGCTACCGATGAAGAGAAgg                                                                                                                                 >  6:90016/1‑148 (MQ=255)
                       gAAAGAGACCGTTATCTTCGCCGACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGaa                                                                                                                               >  8:201239/1‑148 (MQ=255)
                        aaaGAGACCGTTATCTTCGCCGACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGAtgaa                                                                                                                                        <  8:257560/138‑1 (MQ=255)
                        aaaGAGACCGTTATCTTCGCCGACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGAtgaa                                                                                                                                        >  7:257560/1‑138 (MQ=255)
                        aaaGAGACCGTTATCTTCGCCGACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAAt                                                                                                                             >  2:17132/1‑149 (MQ=255)
                        aaaGAGACCGTTATCTTCGCCGACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAAt                                                                                                                             >  1:103259/1‑149 (MQ=255)
                          agagACCGTTATCTTCGCCGACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATc                                                                                                                            >  1:130003/1‑148 (MQ=255)
                             gaCCGTTATCTTCGCCGACCAACTGATGTACACCGTGGTCGCATAATGGACCATTTCCGGCGGTTCCCTCGGTGTTACCGACTTCGTTATCCCGGACGAGAAATCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCgaga                                                                                                                        <  4:27003/149‑1 (MQ=255)
                              aCCGTTATCTTCGCCGACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCCTACCGATGAAGTGAAGGAAATCgaga                                                                                                                        <  1:14306/149‑1 (MQ=255)
                               ccGTTATCTTCGCCGACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGc                                                                                                                      <  8:92740/149‑1 (MQ=255)
                                      cTTCGCCGACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACg                                                                                                               >  6:217280/1‑149 (MQ=255)
                                        tCGCCGACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCgaga                                                                                                                        <  7:293546/138‑1 (MQ=255)
                                        tCGCCGACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCgaga                                                                                                                        >  8:293546/1‑138 (MQ=255)
                                         cGCCGACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGc                                                                                                              <  2:103259/147‑1 (MQ=255)
                                             gACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCgg                                                                                                        >  3:241491/1‑149 (MQ=255)
                                             gACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCgg                                                                                                        >  6:68834/1‑149 (MQ=255)
                                                cAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCga                                                                                                                          >  3:60245/1‑128 (MQ=255)
                                                cAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCga                                                                                                                          <  4:60245/128‑1 (MQ=255)
                                                 aaCTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTg                                                                                                    <  8:285332/149‑1 (MQ=255)
                                                 aaCTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCCGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTg                                                                                                    <  3:306878/149‑1 (MQ=255)
                                                      aTGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAAc                                                                                               >  8:74109/1‑149 (MQ=255)
                                                        gTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGGGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTaccg                                                                                              >  4:14611/1‑145 (MQ=255)
                                                         tACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAAccc                                                                                             <  5:90016/148‑1 (MQ=255)
                                                          acacCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAg                                                                                           >  7:197548/1‑149 (MQ=255)
                                                               ggTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCga                                                                                      <  4:241491/149‑1 (MQ=255)
                                                                gTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCgag                                                                                     <  5:144812/149‑1 (MQ=255)
                                                                   tCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGaa                                                                                   <  3:14611/148‑1 (MQ=255)
                                                                    cGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAg                                                                                  <  7:190284/148‑1 (MQ=255)
                                                                     gCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCa                                                                                                                    >  4:81501/1‑113 (MQ=255)
                                                                     gCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGt                                                                                 <  3:241309/148‑1 (MQ=255)
                                                                         aCTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTacaa                                                                             >  2:8445/1‑148 (MQ=255)
                                                                             gACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAgg                                                                        >  7:274633/1‑149 (MQ=255)
                                                                             gACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAg                                                                         <  1:36244/148‑1 (MQ=255)
                                                                                cATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATAGAGAGCCAGTACGCCTCCGGCCTGGTAACACAGGGCGAGAAGTACAACAAGGTCa                                                                     >  7:285542/1‑149 (MQ=255)
                                                                                  tttCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATc                                                                   <  4:215277/149‑1 (MQ=255)
                                                                                       ggCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACtt                                                              <  5:217280/149‑1 (MQ=255)
                                                                                         cGTTTCCATCGGTGTTACCGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTtgt                                                            <  8:134781/149‑1 (MQ=255)
                                                                                               cATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTtg                                                             <  7:306527/142‑1 (MQ=255)
                                                                                               cATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTtg                                                             >  8:306527/1‑142 (MQ=255)
                                                                                                 tCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAg                                                    >  1:294109/1‑149 (MQ=255)
                                                                                                 tCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAg                                                    <  7:123182/149‑1 (MQ=255)
                                                                                                   ggTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCAGCGACTTGTGGTCGAAGGc                                                  >  6:288591/1‑149 (MQ=255)
                                                                                                       ttAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAac                                              <  7:201239/149‑1 (MQ=255)
                                                                                                         aaCGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAacga                                            >  3:157167/1‑149 (MQ=255)
                                                                                                             aCTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGaa                                        <  8:197548/149‑1 (MQ=255)
                                                                                                             aCTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGaa                                        <  8:274633/149‑1 (MQ=255)
                                                                                                               ttCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATAGACTTGTGGTCGAAGGCTAACGACGAAgt                                      >  6:240115/1‑149 (MQ=255)
                                                                                                                  gTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGAAGTGTc                                   <  5:68834/149‑1 (MQ=255)
                                                                                                                  gTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTCGTGGTCGAAGGCGAACGACGAAgtgt                                    <  1:8445/148‑1 (MQ=255)
                                                                                                                        ccGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTtgtg                                                           >  7:90300/1‑119 (MQ=255)
                                                                                                                        ccGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTtgtg                                                           <  8:90300/119‑1 (MQ=255)
                                                                                                                        ccGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGAAGTGTCCAAGGc                             >  7:50918/1‑149 (MQ=255)
                                                                                                                          ggACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGAAGTGTCCAAGGCga                           <  1:17132/149‑1 (MQ=255)
                                                                                                                               agaAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGAAGTGTCCAAGGCGatgatg                      >  5:232095/1‑149 (MQ=255)
                                                                                                                                     cacGCATCATCCGTAACGCTAACGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGAAGTGTCCAAGGCGATGATGGCCAAc                <  8:285542/147‑1 (MQ=255)
                                                                                                                                      ccGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGAAGTGTCCAAGGCGATGATGGCCAAcc               <  7:74109/149‑1 (MQ=255)
                                                                                                                                       cacttcatcGGTAGCGTTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGAAGTGTCCAAGGCGATGATGGCCAACct              >  8:87389/5‑149 (MQ=255)
                                                                                                                                          ttcatcGGTAGCGTTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGAAGTGTCCAAGGCGATGATGGCCAAACtctc           >  7:131702/2‑149 (MQ=255)
                                                                                                                                          atcatcGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGAAGTGTCCAAGGCGatgatg                      >  5:311570/1‑138 (MQ=255)
                                                                                                                                          atcatcGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGAAGTGTCCAAGGCGatgatg                      <  6:311570/138‑1 (MQ=255)
                                                                                                                                            catcGGTAGCGTTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGAAGTGTCCAAGGCGATGATGGCCAACCTCTCg          >  7:281647/1‑148 (MQ=255)
                                                                                                                                             atcGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGAAGTGTCTAAGGCGATGATGGCCAACCTCTCGa         >  7:214166/1‑148 (MQ=255)
                                                                                                                                                    aCGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGAAGTGTCCAAGGCGATGATGGCCAACCTCTCGAAAGAGaa  <  2:294109/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CAGTGCTACCGCGTGGTTGGTCTGAAAGAGACCGTTATCTTCGCCGACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGAAGTGTCCAAGGCGATGATGGCCAACCTCTCGAAAGAGAA  >  NC_002947/543579‑543874

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: