Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,316,022 A→C 10.0% 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,0220AC10.0% 198.7 / 13.0 70A35A (GCA→GCC) PP_3788non‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base A (34/29);  new base C (2/5);  total (36/34)
Fisher's exact test for biased strand distribution p-value = 2.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.33e-01

TCAGGGAATCAGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTGCAACTGCGCTATCCCGT  >  NC_002947/4315876‑4316165
                                                                                                                                                  |                                                                                                                                               
tcaGGGAATCAGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGc                                                                                                                                               >  2:146408/1‑149 (MQ=255)
 caGGGAATCAGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTAGTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGGATCATCCACCGCTACCCATCGAATGCAGCg                                                                                                                                              >  7:279381/1‑149 (MQ=255)
  aGGGAATCAGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACAGCTACCCATCGAATGCAGCgg                                                                                                                                             >  8:30521/1‑149 (MQ=255)
    ggAATCAGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCa                                                                                                                                           >  7:201119/1‑149 (MQ=255)
          aGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCCGCTGCATTcgat                                                                                                                                     <  7:81096/149‑2 (MQ=255)
          aGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCCGCTGCATTcgat                                                                                                                                     <  8:114338/149‑2 (MQ=255)
           gTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCCGCTGCATTcgatg                                                                                                                                    <  5:249773/149‑3 (MQ=255)
           gTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGAACACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCGTCGAATGTCGCTGCATTcgatg                                                                                                                                    <  1:354220/149‑3 (MQ=255)
             ctcGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTcgacgac                                                                                                                                  >  1:309788/1‑148 (MQ=255)
              tcGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAg                                                                                                                                  >  7:176990/1‑148 (MQ=255)
              tcGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAg                                                                                                                                  >  3:199574/1‑148 (MQ=255)
                gTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAg                                                                                                                                  <  4:212625/146‑1 (MQ=255)
                   cGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCATCGGCATTCGACGAGGACCa                                                                                                                             <  6:96845/148‑1 (MQ=255)
                   cGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAg                                                                                                                            >  3:312704/1‑149 (MQ=255)
                    gCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAgc                                                                                                                           <  8:163386/149‑1 (MQ=255)
                      cacacaGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAgcgc                                                                                                                         >  2:287968/1‑149 (MQ=255)
                        cacaGCAAGGAGCGCAGCATGAACCATGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCtt                                                                                                                       <  4:283831/149‑1 (MQ=255)
                        cacaGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCtt                                                                                                                       >  4:315733/1‑149 (MQ=255)
                                  aGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGc                                                                                                             <  7:30521/149‑1 (MQ=255)
                                       agcaTGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAAc                                                                                                         >  3:239281/1‑148 (MQ=255)
                                              aCCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTcgacga                                                                                                                                   <  3:145420/115‑1 (MQ=255)
                                              aCCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTcgacga                                                                                                                                   >  4:145420/1‑115 (MQ=255)
                                                cAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCg                                                                                                >  7:128877/1‑148 (MQ=255)
                                                 aGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGaa                                                                                              <  4:180474/149‑1 (MQ=255)
                                                  ggCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGaa                                                                                              >  8:116025/1‑148 (MQ=255)
                                                  ggCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAAc                                                                                             <  8:201119/149‑1 (MQ=255)
                                                   gCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACt                                                                                            <  8:128877/149‑1 (MQ=255)
                                                     aTTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGc                                                                                          >  2:7543/1‑149 (MQ=255)
                                                         cTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGAGCAGCGCTTGTCCTACAGCGAACTAAGGATCGAACTGCTca                                                                                       >  5:313252/1‑148 (MQ=255)
                                                            aCGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTcaca                                                                                     <  7:116025/147‑1 (MQ=255)
                                                            aCGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCacac                                                                                    <  4:312704/148‑1 (MQ=255)
                                                              gATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTg                                                                                  >  6:92731/1‑148 (MQ=255)
                                                              gATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGc                                                                                 >  5:304553/1‑149 (MQ=255)
                                                                   cGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCt                                                                            <  3:198403/149‑1 (MQ=255)
                                                                    gggAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTc                                                                           >  1:236144/1‑149 (MQ=255)
                                                                          ccctgcctACCCCGTGCACTACTCGCAGCTCGAATACATCGAGCGTATCATCCACCGATACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTCACCTACAGCGAACTAACGATCGAACTGCTCACACTCAACGCTCGGcag                                                                      >  4:279140/2‑146 (MQ=255)
                                                                           cctgcctACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCa                                                                    >  3:210253/1‑149 (MQ=255)
                                                                                ctACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAAc                                                               >  5:79173/1‑149 (MQ=255)
                                                                                 tACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACt                                                              >  1:184787/1‑149 (MQ=255)
                                                                                          cACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTgc                                                      <  3:315733/148‑1 (MQ=255)
                                                                                             ttCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAg                                                  >  1:40255/1‑149 (MQ=255)
                                                                                                tcGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCgg                                               <  4:210253/149‑1 (MQ=255)
                                                                                                  gCTGGTCGAATACATAGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACAGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCggg                                              >  4:288414/1‑148 (MQ=255)
                                                                                                       tCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGAtgtg                                        <  8:176990/149‑1 (MQ=255)
                                                                                                         gAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTgg                                       <  2:309788/148‑1 (MQ=255)
                                                                                                               aTCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGt                                >  6:205328/1‑149 (MQ=255)
                                                                                                               aTCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCACGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGt                                <  1:7543/149‑1 (MQ=255)
                                                                                                                  gatcGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGaa                             <  4:67133/146‑1 (MQ=255)
                                                                                                                  gAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGaa                             >  8:11371/1‑149 (MQ=255)
                                                                                                                   aGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGc                                                                     <  3:173976/108‑1 (MQ=255)
                                                                                                                   aGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGc                                                                     >  4:173976/1‑108 (MQ=255)
                                                                                                                   aGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAAt                            >  8:93643/1‑149 (MQ=255)
                                                                                                                   aGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAAt                            <  4:239281/149‑1 (MQ=255)
                                                                                                                    gCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATg                           <  5:92731/149‑1 (MQ=255)
                                                                                                                       tATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGCGa                        <  1:287968/149‑1 (MQ=255)
                                                                                                                         tcatcCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGCGACa                      <  1:146408/149‑1 (MQ=255)
                                                                                                                            tcCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTg                   >  3:340184/1‑149 (MQ=255)
                                                                                                                               aCCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTGCa                 >  6:370990/1‑148 (MQ=255)
                                                                                                                                      cccATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTGCAACTgcgc          >  6:170397/1‑148 (MQ=255)
                                                                                                                                      cccATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTGCAACTgcg           <  4:199574/147‑1 (MQ=255)
                                                                                                                                      cccATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTGCAACTgcg           <  7:93643/147‑1 (MQ=255)
                                                                                                                                       tcgtCGAATGCCGCTGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTGCAACTGCGCTa        >  5:99636/4‑149 (MQ=255)
                                                                                                                                        cATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTGCAACTGCGCTAt       >  1:141137/1‑149 (MQ=255)
                                                                                                                                         aTCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTGCAACTGCGCTATc      <  2:184787/149‑1 (MQ=255)
                                                                                                                                         aTCGAACGCACCGGGATTCCTCGCATACAAGCGCTTGACCTACAGCGCACTAAGGATCCAACTGCTCACACTGCACGCTCTGCTGCACGAACTCGGCGTGCGCAGCGGGGATGTGGTGGCGGTGAATGCGACAGTGAAACTGCGCTATc      <  8:279381/149‑1 (MQ=255)
                                                                                                                                          tCGAATGCCGCTGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCg                                                         >  4:212683/1‑97 (MQ=255)
                                                                                                                                          tCGAATGCCGCTGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCg                                                         <  3:212683/97‑1 (MQ=255)
                                                                                                                                           cGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCg                                  <  7:204197/119‑1 (MQ=255)
                                                                                                                                            gAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTGCAACTGCGCTATccc    <  5:205328/148‑1 (MQ=255)
                                                                                                                                             aaTGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACGGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTGCAACTGCGCTATCCcgt  >  8:148710/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
TCAGGGAATCAGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTGCAACTGCGCTATCCCGT  >  NC_002947/4315876‑4316165

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