Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,316,022 A→C 9.8% A35A (GCA→GCC)  PP_3788 → non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,316,0220AC 9.8% 202.5 / 10.8 69A35A (GCA→GCC) PP_3788non‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base A (28/34);  new base C (3/4);  total (31/38)
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

AGGGAATCAGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTGCAACTGCGCTATC  >  NC_002947/4315878‑4316161
                                                                                                                                                |                                                                                                                                           
atggAATCAGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCATGCATTACTGACGATCACGGGAGGGCCTTCCTACCCCGTTCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCgg                                                                                                                                         <  5:365805/147‑1 (MQ=255)
aGGGAATCAGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCgg                                                                                                                                         >  8:149495/1‑149 (MQ=255)
     aTCAGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGACTACCCCGTGCACTTCTCGGTGTGCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAAGTCCGCGGGATTc                                                                                                                                    >  6:20473/1‑149 (MQ=255)
        aGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCCGCTGCATTcgat                                                                                                                                 <  6:265071/149‑2 (MQ=255)
        aGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCCGCTGCATTcgat                                                                                                                                 <  2:122162/149‑2 (MQ=255)
        aGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCCGCTGCATTcgat                                                                                                                                 <  7:198989/149‑2 (MQ=255)
        aGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTcga                                                                                                                                  <  4:229265/148‑1 (MQ=255)
        aGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATAGAGCGTATCATCCACCGCTACCCATCGAATGCCGCTGCATTcga                                                                                                                                  <  6:351212/148‑1 (MQ=255)
            tcGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAg                                                                                                                              >  4:78773/1‑148 (MQ=255)
              gTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAgg                                                                                                                             <  5:116551/147‑1 (MQ=255)
              gTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAg                                                                                                                              <  3:78773/146‑1 (MQ=255)
              gTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGAc                                                                                                                           >  5:63177/1‑149 (MQ=255)
              gTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGAc                                                                                                                           <  6:63177/149‑1 (MQ=255)
                  gCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATGCATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAgc                                                                                                                       >  3:146741/1‑149 (MQ=255)
                   ccACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAgcg                                                                                                                      <  8:183191/149‑1 (MQ=255)
                    cacacaGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCGTCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAgcgc                                                                                                                     >  2:12725/1‑149 (MQ=255)
                    cacacaGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAgcgc                                                                                                                     <  3:117956/149‑1 (MQ=255)
                      cacaGCAAGGAGCGCAGCATGAACCAGGCATTACTGACTATCACGGGACGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCTAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCtt                                                                                                                   <  7:173107/149‑1 (MQ=255)
                       acaGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTg                                                                                                                  <  3:59272/149‑1 (MQ=255)
                         aGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATGCATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGAc                                                                                                                <  4:146741/149‑1 (MQ=255)
                         aGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGAc                                                                                                                <  6:72138/149‑1 (MQ=255)
                           cAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGAcc                                                                                                               >  1:88180/1‑148 (MQ=255)
                             aGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTa                                                                                                             >  5:277783/1‑148 (MQ=255)
                              ggATCGCAGCATGAACCAGGGATTACTGACGATCACGGGAGGGCCTGCCTACGGCGTGCTCTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGATTGCAGCGGCATTCGACGAGGACCAGCGCTTGTCCTACa                                                                                                           <  4:117453/149‑1 (MQ=255)
                                    cagcaTGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAAc                                                                                                     <  6:277783/149‑1 (MQ=255)
                                     agcaTGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAAc                                                                                                     >  7:85940/1‑148 (MQ=255)
                                             ccAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATc                                                                                             >  3:240734/1‑148 (MQ=255)
                                                 gCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACt                                                                                        >  3:194906/1‑149 (MQ=255)
                                                   aTTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTTGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGc                                                                                      >  4:239380/1‑149 (MQ=255)
                                                    ttACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCt                                                                                     >  7:124728/1‑149 (MQ=255)
                                                       cTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTca                                                                                   >  4:46740/1‑148 (MQ=255)
                                                        tGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTcaca                                                                                 >  2:199653/1‑149 (MQ=255)
                                                            gATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTg                                                                              >  7:92117/1‑148 (MQ=255)
                                                                aCGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGc                                                                         <  8:293543/149‑1 (MQ=255)
                                                                  gggAGGGCCTGACTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTc                                                                       >  2:213395/1‑149 (MQ=255)
                                                                    gAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCg                                                                      >  4:69423/1‑148 (MQ=255)
                                                                          ctgcctACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCAt                                                               >  3:348121/1‑149 (MQ=255)
                                                                          ctgcctACCCCGCGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCAt                                                               <  6:331171/149‑1 (MQ=255)
                                                                             cctACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCGTCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGaa                                                            <  1:12725/149‑1 (MQ=255)
                                                                                  cccGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCgg                                                       >  5:5308/1‑149 (MQ=255)
                                                                                   ccGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGc                                                      <  3:46740/149‑1 (MQ=255)
                                                                                     gTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGt                                                    <  1:199653/149‑1 (MQ=255)
                                                                                           ttCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAg                                              >  1:362074/1‑149 (MQ=255)
                                                                                           ttCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAg                                              >  6:290709/1‑149 (MQ=255)
                                                                                                tcTGGTCGAATAGATCGAGCGTATCATCCTCCGCTACCCATCGAATTCAGCGGCATTCGACGAGGACCCGCGCTTGACCTACAGCGAACTATGGATCGAACTGCTCACACTGCACGCTCGGCTGCATTAACTCGGCGGGGGCAGCgg                                           <  8:240423/146‑1 (MQ=255)
                                                                                                    gTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGAtgt                                     <  8:124728/149‑1 (MQ=255)
                                                                                                    gTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGAtg                                      <  4:348121/148‑1 (MQ=255)
                                                                                                    gTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGAt                                       <  1:213395/147‑1 (MQ=255)
                                                                                                     tCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGGGCGCAGCGGGGAtgtg                                    <  3:239380/149‑1 (MQ=255)
                                                                                                         aTACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTg                                <  8:85940/149‑1 (MQ=255)
                                                                                                            cATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCgg                             <  2:88180/149‑1 (MQ=255)
                                                                                                            cATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCg                              <  4:194906/148‑1 (MQ=255)
                                                                                                             aTCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGAt                                       >  8:231582/1‑138 (MQ=255)
                                                                                                             aTCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGAt                                       <  7:231582/138‑1 (MQ=255)
                                                                                                             aTCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGt                            >  6:116828/1‑149 (MQ=255)
                                                                                                             aTCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGt                            >  1:164948/1‑149 (MQ=255)
                                                                                                             aTCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGt                            >  1:19210/1‑149 (MQ=255)
                                                                                                              tCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTGGCGGt                            >  8:101646/1‑148 (MQ=255)
                                                                                                                   cGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGc                      <  6:5308/149‑1 (MQ=255)
                                                                                                                   cGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGc                      <  3:69423/149‑1 (MQ=255)
                                                                                                                       tcatcCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGCGACa                  <  5:290709/149‑1 (MQ=255)
                                                                                                                          tcCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGCGCCAGTg               >  7:306707/1‑149 (MQ=255)
                                                                                                                               cGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTGCAACt          <  5:116828/149‑1 (MQ=255)
                                                                                                                                     tcgtCGAATGCCGCTGCATTCGATGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTGCAACTGCGCTa    >  1:205161/4‑149 (MQ=255)
                                                                                                                                      cATCGACTGCAGCGGCAGTCGACGAGGACCAGCTCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTGCAACTGCGCTAt   <  5:20473/149‑1 (MQ=255)
                                                                                                                                       gtCGAATGCCGCTGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTGCAACTGCGCTATc  >  7:296754/2‑149 (MQ=255)
                                                                                                                                       aTCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGCGa                    <  7:298351/131‑1 (MQ=255)
                                                                                                                                       aTCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGCGa                    >  8:298351/1‑131 (MQ=255)
                                                                                                                                       aTCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTGCAACTGCGCTATc  <  8:92117/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                           
AGGGAATCAGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTGCAACTGCGCTATC  >  NC_002947/4315878‑4316161

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: ‑