Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,561,243 G→C 7.0% I102M (ATC→ATG)  PP_4048 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,561,2430GC 7.0% 218.4 / 5.2 71I102M (ATC→ATG) PP_4048hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (35/31);  new base C (3/2);  total (38/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.80e-01

ACTGTCTGCAGGATGACCTGACTTTCACGTATCAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCC‑ATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCAGCCAGGTCACGGGCACTGAGCAC  >  NC_002947/4561106‑4561389
                                                                                                                                          |                                                                                                                                                  
aCTGTCTGCAGGATGACCTGACTTTCACGTATCAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCC‑ATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATCATCGgcgctgt                                                                                                                                         <  6:84669/149‑4 (MQ=255)
 cTGTCTGCAGGATGACCTGACTTTCACGTATCAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCC‑ATCGAGCCCACGGGAGTGCCCACAGCCCCCACCACACCGCCGATGATCGGCGCGAt                                                                                                                                         >  7:327477/1‑148 (MQ=255)
    tCTGCAGGATGACCTGACTTTCACGTATCAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCC‑ATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAg                                                                                                                                      <  4:316383/148‑1 (MQ=255)
      tGCAGGATGACCTGACTTTCACGTATCAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCC‑ATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATCATCGgcgctgtg                                                                                                                                        >  4:214612/1‑140 (MQ=255)
      tGCAGGATGACCTGACTTTCACGTATCAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCC‑ATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATCATCGgcgctgtg                                                                                                                                        <  3:214612/144‑5 (MQ=255)
        cAGGATGACCTGACTTTCACGTATCAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCC‑ATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAggg                                                                                                                                    <  3:122275/146‑1 (MQ=255)
          ggATGACCTGACTTTCACGTATCAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCC‑ATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGtcatc                                                                                                                               <  1:177909/149‑1 (MQ=255)
          ggATGACCTGACTTTCACGTATCAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCC‑ATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGtcatc                                                                                                                               <  1:236737/149‑1 (MQ=255)
           gATGACCTGACTTTCACGTATCAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCC‑ATCGAGCCCACGGGAGTGCTCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGtcat                                                                                                                                >  3:219456/1‑147 (MQ=255)
                 cTGACTTTCACGTATCAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCC‑ATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATcaggcagg                                                                                                                        >  1:218172/1‑149 (MQ=255)
                 cTGACTTTCACGTATCAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCC‑ATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATcaggcagg                                                                                                                        >  1:309900/1‑149 (MQ=255)
                  tGACTTTCACGTATCAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCC‑ATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCaggcaggc                                                                                                                       >  3:89944/1‑149 (MQ=255)
                    aCTTTCACGTATCAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGTGCCTGCC‑ATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGcca                                                                                                                     >  7:249940/1‑149 (MQ=255)
                    aCTTTCACGTATCAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCC‑ATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGcc                                                                                                                      <  5:171053/148‑1 (MQ=255)
                          aCGTATCAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCC‑ATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCCATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGc                                                                                                                >  1:137054/1‑148 (MQ=255)
                            gTATCAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCC‑ATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCaa                                                                                                              >  5:220572/1‑148 (MQ=255)
                                cgaCAGCAAGCACAAACCACGCAAGAACGCTGCTGCGTCCTGCCCAACCAGTGCGCCTGCC‑ATCGAGCCCACGGGAGTGCCCCCCGCCCCCACCCCAGCGCCGATGATCGCGGCGATCAGGGTCATCAGGCAGGCCACCAGCAGCAccc                                                                                                         >  8:88242/3‑149 (MQ=255)
                                cAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCC‑ATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACAcc                                                                                                          >  5:137368/1‑148 (MQ=255)
                                 aaCAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCC‑ATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCACCACCCa                                                                                                        >  5:55508/1‑149 (MQ=255)
                                 aaCAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCC‑ATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCa                                                                                                        >  6:34350/1‑149 (MQ=255)
                                    agcaagcaCAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCC‑ATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCActc                                                                                                     >  4:151806/1‑149 (MQ=255)
                                     gcaagcaCAAACCCCGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCC‑ATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCActct                                                                                                    <  3:295268/149‑1 (MQ=255)
                                        agcaCAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCC‑ATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGc                                                                                                 <  6:267454/149‑1 (MQ=255)
                                             aaaCCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCC‑ATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTc                                                                                            <  6:55508/149‑1 (MQ=255)
                                               aCCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCC‑ATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCtt                                                                                                   >  2:180191/1‑140 (MQ=255)
                                               aCCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCC‑ATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCtt                                                                                                   <  1:180191/140‑1 (MQ=255)
                                               aCCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCC‑ATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCcg                                                                                          >  8:190787/1‑149 (MQ=255)
                                               aCCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCC‑ATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCcg                                                                                          <  4:89944/149‑1 (MQ=255)
                                               aCCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCC‑ATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCcg                                                                                          <  3:161805/149‑1 (MQ=255)
                                                ccACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCC‑ATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCcgg                                                                                         <  2:218172/149‑1 (MQ=255)
                                                 cACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCC‑ATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCcggc                                                                                        >  4:300102/1‑149 (MQ=255)
                                                    gCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCC‑ATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCggcggc                                                                                     >  7:29643/1‑149 (MQ=255)
                                                         aaCGCTGCTGCTTCCTGCCCAACCAGTCCGCCTGCC‑ATCGAGCCCACGGGAGTGCCCACCGCCACCACCACAGCGCCGATGATCCGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTgcgg                                                                                >  6:296982/1‑149 (MQ=255)
                                                          aCGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCC‑ATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTgcggc                                                                               <  3:151806/149‑1 (MQ=255)
                                                             ctgctgctTCCTGCCCAACCAGTGCGCCTGCC‑ATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGtt                                                                            <  3:300102/149‑1 (MQ=255)
                                                              tgctgctTCCTGCCCAACCAGTGCGCCTGCCGATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGtt                                                                            >  1:225587/1‑149 (MQ=255)
                                                                   catCCTGCCCAACCAGTGCGCCTGCC‑ATCGAGCCGAGGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCATGGTCATCAGGCAGTCCACCAGCAACACCGACTCTTGCCATTCCGGGGGCTGCGGCGTTGCCTgg                                                                      <  5:296982/147‑1 (MQ=255)
                                                                      ccTGCCCAACCAGTGCGCCTGCC‑ATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCt                                                                    <  8:29643/148‑1 (MQ=255)
                                                                           ccAACCAGTGCGCCTGCC‑ATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTc                                                              <  3:342625/149‑1 (MQ=255)
                                                                            cAACCAGTGTGCCTGCC‑ATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTcc                                                             <  8:249940/149‑1 (MQ=255)
                                                                            cAACCAGTGCGCCTGCC‑ATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTcc                                                             >  1:247629/1‑149 (MQ=255)
                                                                            cAACCAGTGCGCCTGCC‑ATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTcc                                                             <  5:34350/149‑1 (MQ=255)
                                                                             aaCCAGTGCGCCTGCC‑ATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCa                                                                                                                                       <  4:295316/74‑1 (MQ=255)
                                                                             aaCCAGTGCGCCTGCC‑ATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCa                                                                                                                                       >  3:295316/1‑74 (MQ=255)
                                                                              aCCAGTGCGCCTGCC‑ATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTccgc                                                           >  2:346904/1‑149 (MQ=255)
                                                                                  gTGCGCCTGCC‑ATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCATGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCcgcg                                                       >  2:48113/1‑149 (MQ=255)
                                                                                   tGCGCCTGCC‑ATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGCTCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCcgcg                                                       >  6:264977/1‑148 (MQ=255)
                                                                                     cgcCTGCC‑ATCGAGCCCACGGGAGTGCTCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAg                                                    <  4:219456/149‑1 (MQ=255)
                                                                                     cgcCTGCC‑ATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTa                                                     <  8:327477/148‑1 (MQ=255)
                                                                                      gcCTGCC‑ATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAgc                                                   <  7:190787/149‑1 (MQ=255)
                                                                                         tGCC‑ATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCt                                                >  5:347649/1‑149 (MQ=255)
                                                                                           cc‑ATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCtt                                               <  6:137368/148‑1 (MQ=255)
                                                                                                     ccACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGgcgc                                     >  4:217679/1‑149 (MQ=255)
                                                                                                       aCGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGc                                                                          <  5:232201/110‑1 (MQ=255)
                                                                                                       aCGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGc                                                                          >  6:232201/1‑110 (MQ=255)
                                                                                                       aCGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAg                                   <  2:309900/149‑1 (MQ=255)
                                                                                                       aCGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAg                                   >  7:43883/1‑149 (MQ=255)
                                                                                                                cccACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGc                          <  3:217679/149‑1 (MQ=255)
                                                                                                                   aCCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTCGCGCTTGCCCACGCGCAGGCACACCGCAGc                       >  6:77771/1‑149 (MQ=255)
                                                                                                                      gCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCAGCCAg                    <  8:43883/149‑1 (MQ=255)
                                                                                                                      gCCCCCACCACAGCGACGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCAGCCAg                    <  6:220572/149‑1 (MQ=255)
                                                                                                                           caccacAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCAGCCAGGTCAc               >  8:75268/1‑149 (MQ=255)
                                                                                                                              cacaGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCAGCCAGGTCACggg            >  1:147463/1‑149 (MQ=255)
                                                                                                                                tcgcgcCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCAGCCAGGTCACGGGCa          >  6:241315/3‑149 (MQ=255)
                                                                                                                                tcgcgcCGATCATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCAGCCAGGTCACGGGCa          >  6:120213/3‑149 (MQ=255)
                                                                                                                                   cgcCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCAGCCAGGTCACGGGCACt        >  7:210629/1‑148 (MQ=255)
                                                                                                                                   cgcCGATCATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCAGCCAGGTCACGGGCACt        >  8:26383/1‑148 (MQ=255)
                                                                                                                                     ccGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCCCCAGCAACACCCACTCTTGCCATTACGGCGGCTGCGGCGTTGCCTGGCTGACGTCCGCCGCGTAGCGCTTTCCCAGGCGCAGGACCCACGCAGACAGGTCACTGGCACTGAg     >  5:209243/1‑149 (MQ=255)
                                                                                                                                      cGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGGTTGCCCAGGCGCAGGCACACCGCAGCCAGGTCACGGGCACTGAGc    <  6:347649/149‑1 (MQ=255)
                                                                                                                                      cGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCtt                                               <  2:142664/106‑1 (MQ=255)
                                                                                                                                      cGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCtt                                               >  1:142664/1‑106 (MQ=255)
                                                                                                                                       gatgatCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCCCTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCAGCCAGGTCACGGGCACTGAGca   >  7:242126/1‑149 (MQ=255)
                                                                                                                                       gatgatCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCAGCCAGGTCACGGGCACTGAGca   <  1:48113/149‑1 (MQ=255)
                                                                                                                                        atgatCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTcc                                                             <  3:211610/90‑1 (MQ=255)
                                                                                                                                        atgatCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTcc                                                             >  4:211610/1‑90 (MQ=255)
                                                                                                                                         tgatCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCAGCCAGGTCACGGGCACTGAGcac  >  1:174458/1‑148 (MQ=255)
                                                                                                                                          |                                                                                                                                                  
ACTGTCTGCAGGATGACCTGACTTTCACGTATCAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCC‑ATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCAGCCAGGTCACGGGCACTGAGCAC  >  NC_002947/4561106‑4561389

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