Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 543,730 C→T 6.5% A657V (GCT→GTT)  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,7300CT 6.5% 224.5 / 3.5 77A657V (GCT→GTT) rpoCDNA‑directed RNA polymerase subunit beta'
Reads supporting (aligned to +/- strand):  ref base C (33/38);  new base T (3/2);  total (37/40)
Fisher's exact test for biased strand distribution p-value = 6.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTACCGCGTGGTTGGTCTGAAAGAGACCGTTATCTTCGCCGACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGAAGTGTCCAAGGCGATGATGGCCAACCTCTCGAAAGAGAAGG  >  NC_002947/543584‑543876
                                                                                                                                                   |                                                                                                                                                  
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                                                                                                                            <  4:60245/128‑1 (MQ=255)
                                           cAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCga                                                                                                                            >  3:60245/1‑128 (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                                                      <  7:123182/149‑1 (MQ=255)
                                                                                            tCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAg                                                      >  1:294109/1‑149 (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:274633/149‑1 (MQ=255)
                                                                                                        aCTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGaa                                          <  8:197548/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                                                             <  8:90300/119‑1 (MQ=255)
                                                                                                                   ccGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTtgtg                                                             >  7:90300/1‑119 (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)
                                                                                                                                                cGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGAAGTGTCCAAGGCGATGATGGCCAACCTCTCGAAAGAGAAgg  >  4:104385/1‑149 (MQ=255)
                                                                                                                                                 gCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCTGCCTGGTAACCCAGGGCGAGATGTACAACAACGTCAACGACTTGTGGTCTAAGCCCAACGACGAACTGTCCAAGGCGATGATGGCAAACCTCTCGAAAGAGAAgg  >  5:20057/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CTACCGCGTGGTTGGTCTGAAAGAGACCGTTATCTTCGCCGACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGAAGTGTCCAAGGCGATGATGGCCAACCTCTCGAAAGAGAAGG  >  NC_002947/543584‑543876

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: