Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,638,137 T→G 5.4% V67G (GTG→GGG)  PP_2308 → acyl‑CoA thioesterase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,638,1370TG 5.4% 270.1 / 6.2 93V67G (GTG→GGG) PP_2308acyl‑CoA thioesterase II
Reads supporting (aligned to +/- strand):  ref base T (42/46);  new base G (3/2);  total (45/48)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.54e-01

CAGCATTCCGCCCAGCTGGGCCCAGGGGCGCGCTGCCTTTGGCGGGCTGATGGCGGCCATGGTCTATGAAACCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCTCATTGGCCATCAGCTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGCCAGGTGGTGACTTTGGTGCAGGGCAATTTCGGTGCGGGCCGCCCTTCGGTGGTCGAAGTGC  >  NC_002947/2637988‑2638277
                                                                                                                                                       |                                                                                                                                            
cAGCATTCCGCCCAGCTGGGCCCAGGGGCGCGCTGCCTTTGGCGGGCTGATGGCGGCCATGGTCTATGAAACCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCTCATTGGCCATCAGCTTCGTGGCACCCGCGGCGGCGGAtgt                                                                                                                                              <  1:1380436/150‑1 (MQ=255)
 aGCATTCCGCCCAGCTGGGCCCAGGGGCGCGCTGCCTTTGGCGGGCTGATGGCGGCCATGGTCTATGAAACCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCTCATTGGCCATCAGCTTCGTGGCACCCGCGGCGGCGGAtgtg                                                                                                                                             <  2:872871/150‑1 (MQ=255)
     ttCCGCCCAGCTGGGCCCAGGGGCGCGCTGCCTTTGGCGGGCTGATGGCGGCCATGGTCTATGAAACCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCTCATTGGCCATCAGCTTCGTGGCACCCGCGGCGGCGGATGTGCCCa                                                                                                                                         >  2:1621024/1‑150 (MQ=255)
       tctcccAGCTGGGCCCAGGGGCGCGCTGCCTTTGGCGGGCTGATGGCGGCCATGGTCTATGAAACCATGCGCCTCAAGATCAGCGCCAACCGCCCGGTACGCTCATTTTCCATCATCTTCGTGGCACCCGCGGCGGCGGATGTGCCCATc                                                                                                                                       <  1:647888/147‑1 (MQ=255)
        cGCCCAGCTGGGCCCAGGGGCGCGCTGCCTTTGGCGGGCTGATGGCGGCCATGGTCTATGAAACCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCTCATTGGCCATCAGCTTCGTGGCACCCGCGGCGGCGGATGTGCCCATcc                                                                                                                                      >  2:879707/1‑150 (MQ=255)
         gCCCAGCTGGGCCCAGGGGCGCGCTGCCTTTGGCGGGCTGATGGCGGCCATGGTCTATGAAACCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCTCATTGGCCATCAGCTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCg                                                                                                                                     >  2:529280/1‑150 (MQ=255)
           ccAGCTGGGCCCAGGGGCGCGCTGCCTTTGGCGGGCTGATGGCGGCCATGGTCTATGAAACCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCTCATTGGCCATCAGCTTCGTGGCACCCGCGGCGGCGGATGTGCCCATccaga                                                                                                                                   >  2:647888/1‑147 (MQ=255)
             aGCTGGGCCCAGGGGCGCGCTGCCTTTGGCGGGCTGATGGCGGCCATGGTCTATGAAACCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCTCATTGGCCATCAGCTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTc                                                                                                                                 >  2:314291/1‑150 (MQ=255)
                tGGGCCCAGGGGCGCGCTGCCTTTGGCGGGCTGATGGCGGCCATGGTCTATGAAACCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCTCATTGGCCATCAGCTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAg                                                                                                                              <  2:1699522/150‑1 (MQ=255)
                  ggCCCAGGGGCGCGCTGCCTTTGGCGGGCTGATGGCGGCCATGGTCTATGAAACCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCTCATTGGCCATCAGCTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGt                                                                                                                            <  2:554273/150‑1 (MQ=255)
                   gCCCAGGGGCGCGCTGCCTTTGGCGGGCTGATGGCGGCCATGGTCTATGAAACCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCTCATTGGCCATCAGCTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTg                                                                                                                           >  2:1087169/1‑150 (MQ=255)
                     ccAGGGGCGCGCTGCCTTTGGCGGGATCATGGCGGACACGGCCGATGAAACCATGCGCCTCAAGATCAGCGAGTGCCGCCCGGTACGCTCATTGGCCATCAGCTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGa                                                                                                                         <  1:1738818/150‑1 (MQ=255)
                      cAGGGGCGCGCTGCCTTTGGCGGGCTGATGGCGGCCATGGTCTATGAAACCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCTCATTGGCCATCAGCTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAg                                                                                                                        <  1:758381/150‑1 (MQ=255)
                           gcgcgcTGCCTTTGGCGGGCTGATGGCGGCCATGGTCTATGAAACCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCTCATTGGCCATCAGCTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGtttt                                                                                                                   >  2:2403229/1‑150 (MQ=255)
                           gcgcgcTGCCTTTGGCGGGCTGATGGCGGCCATGGTCTATGAAACCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCTCATTGGCCATCAGCTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGtttt                                                                                                                   <  1:1621024/150‑1 (MQ=255)
                           gcgcgcTGCCTTTGGCGGGCTGATGGCGGCCATGGTCTATGAAACCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCTCATTGGCCATCAGCTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGtttt                                                                                                                   <  2:1757311/150‑1 (MQ=255)
                              cgcTGCCTTTGGCGGGCTGATGGCGGCCATGGTCTATGAAACCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCTCATTGGCCATCAGCTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTgcg                                                                                                                >  1:1015817/1‑150 (MQ=255)
                                 tGCCTTTGGCGGGCTGATGGCGGCCATGGTCTATGAAACCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCTCATTGGCCATCAGCTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAgg                                                                                                                             >  1:2222013/1‑134 (MQ=255)
                                 tGCCTTTGGCGGGCTGATGGCGGCCATGGTCTATGAAACCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCTCATTGGCCATCAGCTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAgg                                                                                                                             <  2:2222013/134‑1 (MQ=255)
                                 tGCCTTTGGCGGGCTGATGGCGGCCATGGTCTATGAAACCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCTCATTGGCCATCAGCTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGa                                                                                                             >  2:1660970/1‑150 (MQ=255)
                                  gCCTTTGGCGGGCTGATGGCGGCCATGGTCTATGAAACCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCTCATTGGCCATCAGCTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGaa                                                                                                            <  1:650244/150‑1 (MQ=255)
                                  gCCTTTGGCGGGCTGATGGCGGCCATGGTCTATGAAACCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCTCATTGGCCATCAGCTTCGTGGCAACCGTGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGaa                                                                                                            <  1:551208/150‑1 (MQ=255)
                                     tttGGCGGGCTGATGGCGGCCATGGTCTATGAAACCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCTCATTGGCCATCAGCTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGaaggc                                                                                                         >  2:1602953/1‑150 (MQ=255)
                                     tttGGCGGGCTGATGGCGGCCATGGTCTATGAAACCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCTCATTGGCCATCAGCTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGaaggc                                                                                                         >  2:2203354/1‑150 (MQ=255)
                                             gCTGATGGCGGCCATGGTCTATGAAACCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCTCATTGGCCATCAGCTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGt                                                                                                 <  1:954143/150‑1 (MQ=255)
                                                gATGGCGGCCATGGTCTATGAAACCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCTCATTGGCCATCAGCTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAg                                                                                              <  1:1471163/150‑1 (MQ=255)
                                                  tGGCGGCCATGGTCTATGAAACCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCTCATTGGCCATCAGCTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGaaggcaag                                                                                                      <  2:1089467/140‑1 (MQ=255)
                                                  tGGCGGCCATGGTCTATGAAACCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCTCATTGGCCATCAGCTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGaaggcaag                                                                                                      >  1:1089467/1‑140 (MQ=255)
                                                    gcggcCATGGTCTATGAAACCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCTCATTGGAATTAGTGTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACg                                                                                          <  1:1333400/150‑1 (MQ=255)
                                                      ggcCATGGTCTATGAAACCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCTCATTGGCCATCAGCTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGtt                                                                                        >  2:1673200/1‑150 (MQ=255)
                                                        ccATGGTCTATGAAACCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCTCATTGGCCATCAGCTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGaaggcaa                                                                                                       >  1:182592/1‑133 (MQ=255)
                                                        ccATGGTCTATGAAACCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCTCATTGGCCATCAGCTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGaaggcaa                                                                                                       <  2:182592/133‑1 (MQ=255)
                                                          aTGGTCTATGAAACCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCTCATTGGCCATCAGCTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTg                                                                                    >  2:2219509/1‑150 (MQ=255)
                                                                tATGAAACCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCTCATTGGCCATCAGCTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCcgc                                                                              >  2:724694/1‑150 (MQ=255)
                                                                tATGAAACCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCTCATTGGCCATCAGCTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCcgc                                                                              <  2:1837553/150‑1 (MQ=255)
                                                                 aTGAAACCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCTCATTGGCCATCAGCTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCcgcg                                                                             >  2:407636/1‑150 (MQ=255)
                                                                 aTGAAACCATGCGCCTCAAGATCAGCCCCAACCCCCAGGTAAAATGCGACTTTCATAACCTCGTGGCACTCTCGGCGGCGGATGTGCCCATCCGCTTCGAGGGGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCcgcg                                                                             >  1:1863959/1‑150 (MQ=255)
                                                                    aaaCCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCTCATTGGCCATCAGCTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGaa                                                                                                            >  1:526257/1‑116 (MQ=255)
                                                                    aaaCCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCTCATTGGCCATCAGCTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGaa                                                                                                            <  2:526257/116‑1 (MQ=255)
                                                                      aCCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCTCATTGGCCATCAGCTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGtt                                                                        >  1:929427/1‑150 (MQ=255)
                                                                       ccATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCTCATTGGCCATCAGCTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTc                                                                       <  1:529280/150‑1 (MQ=255)
                                                                        cATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCTCATTGGCCATCAGCTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTgcgc                                                                                                               <  1:1335662/109‑1 (MQ=255)
                                                                        cATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCTCATTGGCCATCAGCTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTgcgc                                                                                                               >  2:1335662/1‑109 (MQ=255)
                                                                        cATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCTCATTGGCCATCAGCTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCa                                                                      <  1:1246843/150‑1 (MQ=255)
                                                                        cATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCTCATTGGCCATCAGCTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCa                                                                      <  1:1884671/150‑1 (MQ=255)
                                                                         aTGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCTCATTGGCCATCAGCTTCGTGGCACCCGCGGCGGCGGATGTGACCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAg                                                                     >  1:2179614/1‑150 (MQ=255)
                                                                            cgcCTCAAGATCAGCGACAACCGCCCGGTACGCTCATTGGCCATCAGCTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGAt                                                                  <  1:1087169/150‑1 (MQ=255)
                                                                              ccTCAAGATCAGCGACAACCGCCCGGTACGCTCATTGGCCATCAGCTTCGTGGCACCAGCGGCGGCGGATGTGGCTTGTCGCTTCTAGGTGGGGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATgg                                                                >  1:289948/1‑150 (MQ=255)
                                                                               cTCAAGATCAGCGACAACCGCCCGGTACGCTCATTGGCCATCAGCTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGc                                                               <  1:1653628/150‑1 (MQ=255)
                                                                               cTCAAGATCAGCGACAACCGCCCGGTACGCTCATTGGCCATCAGCTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGc                                                               >  2:1148860/1‑150 (MQ=255)
                                                                                tCAAGATCAGCGACAACCGCCCGGTACGCTAATTGGCCATCAGCTTCGTGGAACTCGCGGCGGCGGATGTCCCCATCCGGTTCGAAGAGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCGGGGCGGCGCTGTTCAGGATGGcc                                                              >  1:714746/1‑150 (MQ=255)
                                                                                  aaGATCAGCGACAACCGCCCGGTACGCTCATTGGCCATCAGCTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGtt                                                                                        >  2:1652872/1‑122 (MQ=255)
                                                                                  aaGATCAGCGACAACCGCCCGGTACGCTCATTGGCCATCAGCTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGtt                                                                                        <  1:1652872/122‑1 (MQ=255)
                                                                                  aaGATCAGCGACAACCGCCCGGTACGCTCATTGGCCATCAGCTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGCCAg                                                            <  2:1863959/150‑1 (MQ=255)
                                                                                  aaGATCAGCGACAACCGCCCGGTACGCTCATTGGCCATCAGCTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGCCAg                                                            <  1:407636/150‑1 (MQ=255)
                                                                                  aaGATCAGCGACAACCGCCCGGTACGCTCATTGGCCATCAGCTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGCCAg                                                            >  1:1191587/1‑150 (MQ=255)
                                                                                       cAGCGACAACCGCCCGGTACGCTCATTGGCCATCAGCTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGCCAggtggt                                                       <  1:954410/150‑1 (MQ=255)
                                                                                         gCGACAACCGCCCGGTACGCTCATTGGCCATCAGCTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCgcgc                                                                            >  2:936388/1‑127 (MQ=255)
                                                                                         gCGACAACCGCCCGGTACGCTCATTGGCCATCAGCTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCgcgc                                                                            <  1:936388/127‑1 (MQ=255)
                                                                                         gCGACAACCGCCCGGTACGCTCATTGGCCATCAGCTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGCCAGGTGGTGa                                                     <  1:1054704/150‑1 (MQ=255)
                                                                                               aCCGCCCGGTACGCTCATTGGCCATCAGCTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCg                                                                                                   <  2:1607430/98‑1 (MQ=255)
                                                                                               aCCGCCCGGTACGCTCATTGGCCATCAGCTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCg                                                                                                   >  1:1607430/1‑98 (MQ=255)
                                                                                                 cGCCCGGTACGCTCATTGGCCATCAGCTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGCCAGGTGGTGACTTTGGTg                                             >  1:590968/1‑150 (MQ=255)
                                                                                                 cGCCCGGTACGCTCATTGGCCATCAGCTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGCCAGGTGGTGACTTTGGTg                                             >  2:200536/1‑150 (MQ=255)
                                                                                                    ccGGTACGCTCATTGGCCATCAGCTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGCCAGGTGGTGACTTTGGTGCAg                                          <  1:314291/150‑1 (MQ=255)
                                                                                                      ggTACGCTCATTGGCCATCAGCTTCGTGGCACCCGCGGCGGCGGAGGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGCCAGGTGGTGACTTTGGTGCAggg                                        <  2:444190/150‑1 (MQ=255)
                                                                                                          cGCTCATTGGCCATCAGCTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGCCAGGTGGTGACTTTGGTGCAGGGCAAt                                    <  2:613609/150‑1 (MQ=255)
                                                                                                           gCTCATTGGCCATCAGCTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGCCAGGTGGTGACTTTGGTGCAGGGCAAtt                                   >  2:2122173/1‑150 (MQ=255)
                                                                                                             tCATTGGCCATCAGCTTCGTGGCACCGGAATCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGCCAGGTGGTGACTTTGGTGCAGGGCAATTTc                                 >  1:665661/1‑150 (MQ=255)
                                                                                                                aTTGGCCATCAGCTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCGAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGCCAGGTGGTGACTTTGGTGCAgg                                         <  2:213927/140‑1 (MQ=255)
                                                                                                                aTTGGCCATCAGCTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCGAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGCCAGGTGGTGACTTTGGTGCAgg                                         >  1:213927/1‑140 (MQ=255)
                                                                                                                aTTGGCCATCAGCTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGCCAGGTGGTGACTTTGGTGCAgg                                         <  1:670427/140‑1 (MQ=255)
                                                                                                                aTTGGCCATCAGCTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGCCAGGTGGTGACTTTGGTGCAgg                                         >  2:670427/1‑140 (MQ=255)
                                                                                                                aTTGGCCATCAGCTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGCCAGGTGGTGACTTTGGTGCAGGGCAATTTCgg                               >  1:1607453/1‑150 (MQ=255)
                                                                                                                 ttGGCCATCAGCTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGCCAGGTGGTGACTTTGGTGCAGGGCAATTTCGGt                              >  1:86196/1‑150 (MQ=255)
                                                                                                                 ttGGCCATCAGCTTCGTCAGACGCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGCCAGGTGGTGACTTTGGTGCAGGGCAATTTCGGt                              <  1:1306335/150‑1 (MQ=255)
                                                                                                                    gCCATCAGCTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTCGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGCCAggtggt                                                       >  2:2278154/1‑122 (MQ=255)
                                                                                                                    gCCATCAGCTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTCGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGCCAggtggt                                                       <  1:2278154/122‑1 (MQ=255)
                                                                                                                     ccATCAGCTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGCCAGGTGGTGACTTTGGTGCAGGGCAATTTCGGTGCgg                          <  1:2203354/150‑1 (MQ=255)
                                                                                                                           gCTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGCCAGGTGGTGACTTTGGTGCAGGGCAATTTCGGTGCGGgccgcc                    >  2:1645039/1‑150 (MQ=255)
                                                                                                                              ttCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGCCAGGTGGTGACTTTGGTGCAGGGCAATTTCGGTGCGGGCCGCCCt                  >  2:1637719/1‑150 (MQ=255)
                                                                                                                                cGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGCCAGGTGGTGACTTTGGTGCAGGGCAATTTCGGTGCGGGCCGCCCTTc                <  2:929427/150‑1 (MQ=255)
                                                                                                                                   ggCCCCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGCCAGGTGGTGACTTTCGTGCAGGGCAATTTCGGTGCGCGCCGTCCTTCggt             <  2:690637/150‑1 (MQ=255)
                                                                                                                                   ggCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGCGGCAGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGCCAGGTGGTGACTTTGGTGCAGGGCAATTTCGGTGCGGGCCGCCCTTCggt             <  1:724694/150‑1 (MQ=255)
                                                                                                                                   ggCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGCCAGGTGGTGACTTTGGTGCAGGGCAATTTCGGTGCGGGCCGCCCTTCggt             <  1:1602953/150‑1 (MQ=255)
                                                                                                                                   ggCACCCGCGGCGGCGGATGTGCCCATCAGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGCCAGGTGGTGACTTTGGTGCAGGGCAATTTCGGTGCGGGCCGCCCTTCggt             <  1:182875/150‑1 (MQ=255)
                                                                                                                                     cACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGCCAGGTGGTGACTTTGGTGCAGGGCAATTTCGGTGCGGGCCGCCCTTCggtgg           >  2:1050866/1‑150 (MQ=255)
                                                                                                                                          gcggcggcgGATGTGCCTATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGCCAGGTGGTGACTTTGGTGCAGGGCAATTTCGGTGCGGGCCGCCCTTCGGTGGTCGaa      <  2:289948/150‑1 (MQ=255)
                                                                                                                                          gcggcggcgGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGCCAGGTGGTGACTTTGGTGCAGGGCAATTTCGGTGCGGGCCGCCCTTCGGTGGTCGaa      <  1:671507/150‑1 (MQ=255)
                                                                                                                                          gcggcggcgGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGCCAGGTGGTGACTTTGGTGCAGGGCAATTTCGGTGCGGGCCGCCCTTCGGTGGTCGaa      <  2:1015817/150‑1 (MQ=255)
                                                                                                                                          ctcgaagcggATGGGCACATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGCCAGGTGGTGACTTTGGTGCAGGGCAATTTCGGTGCGGGCCGCCCTTCGGTGGTCGaa      >  2:2148621/7‑150 (MQ=255)
                                                                                                                                           tcgaagcggATGGGCACATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGCCAGGTGGTGACTTTGGTGCAGGGCAATTTCGGTg                             >  2:1929533/6‑126 (MQ=255)
                                                                                                                                           tcgaagcggATGGGCACATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGCCAGGTGGTGACTTTGGTGCAGGGCAATTTCGGTg                             <  1:1929533/121‑1 (MQ=255)
                                                                                                                                           tcgaagcggATGGGCACATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGCCAGGTGGTGACTTTGGTGCAGGGCAATTTCGGTGCGGGCCGCCCTTCGGTGGTCGAAg     >  2:5164/6‑150 (MQ=255)
                                                                                                                                              aagcggATGGGCACATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGCCAGGTGGTGACTTTGGTGCAGGGCAATTTCGGTGCGGGCCGCCCTTCGGTGGTCGAAGTGc  <  1:2148621/148‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                            
CAGCATTCCGCCCAGCTGGGCCCAGGGGCGCGCTGCCTTTGGCGGGCTGATGGCGGCCATGGTCTATGAAACCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCTCATTGGCCATCAGCTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGCCAGGTGGTGACTTTGGTGCAGGGCAATTTCGGTGCGGGCCGCCCTTCGGTGGTCGAAGTGC  >  NC_002947/2637988‑2638277

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