Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,638,140 C→A 5.4% P68H (CCC→CAC)  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,1400CA 5.4% 295.8 / 4.6 93P68H (CCC→CAC) PP_2308acyl‑CoA thioesterase II
Reads supporting (aligned to +/- strand):  ref base C (42/46);  new base A (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 = 7.02e-01

TTCCGCCCAGCTGGGCCCAGGGGCGCGCTGCCTTTGGCGGGCTGATGGCGGCCATGGTCTATGAAACCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCT‑CATTGGCCATCAG‑CTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGCCAGGTGGTGACTTTGGTGCAGGGCAATTTCGGTGCGGGCCGCCCTTCGGTGGTCGAAGTGCCGGCGTTGCCGG  >  NC_002947/2637993‑2638289
                                                                                                                                                     |                                                                                                                                                     
ttCCGCCCAGCTGGGCCCAGGGGCGCGCTGCCTTTGGCGGGCTGATGGCGGCCATGGTCTATGAAACCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCT‑CATTGGCCATCAG‑CTTCGTGGCACCCGCGGCGGCGGATGTGCCCa                                                                                                                                                     >  2:1621024/1‑150 (MQ=255)
  tctcccAGCTGGGCCCAGGGGCGCGCTGCCTTTGGCGGGCTGATGGCGGCCATGGTCTATGAAACCATGCGCCTCAAGATCAGCGCCAACCGCCCGGTACGCT‑CATTTTCCATCAT‑CTTCGTGGCACCCGCGGCGGCGGATGTGCCCATc                                                                                                                                                   <  1:647888/147‑1 (MQ=255)
   cGCCCAGCTGGGCCCAGGGGCGCGCTGCCTTTGGCGGGCTGATGGCGGCCATGGTCTATGAAACCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCT‑CATTGGCCATCAG‑CTTCGTGGCACCCGCGGCGGCGGATGTGCCCATcc                                                                                                                                                  >  2:879707/1‑150 (MQ=255)
    gCCCAGCTGGGCCCAGGGGCGCGCTGCCTTTGGCGGGCTGATGGCGGCCATGGTCTATGAAACCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCT‑CATTGGCCATCAG‑CTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCg                                                                                                                                                 >  2:529280/1‑150 (MQ=255)
      ccAGCTGGGCCCAGGGGCGCGCTGCCTTTGGCGGGCTGATGGCGGCCATGGTCTATGAAACCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCT‑CATTGGCCATCAG‑CTTCGTGGCACCCGCGGCGGCGGATGTGCCCATccaga                                                                                                                                               >  2:647888/1‑147 (MQ=255)
        aGCTGGGCCCAGGGGCGCGCTGCCTTTGGCGGGCTGATGGCGGCCATGGTCTATGAAACCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCT‑CATTGGCCATCAG‑CTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTc                                                                                                                                             >  2:314291/1‑150 (MQ=255)
           tGGGCCCAGGGGCGCGCTGCCTTTGGCGGGCTGATGGCGGCCATGGTCTATGAAACCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCT‑CATTGGCCATCAG‑CTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAg                                                                                                                                          <  2:1699522/150‑1 (MQ=255)
             ggCCCAGGGGCGCGCTGCCTTTGGCGGGCTGATGGCGGCCATGGTCTATGAAACCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCT‑CATTGGCCATCAG‑CTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGt                                                                                                                                        <  2:554273/150‑1 (MQ=255)
              gCCCAGGGGCGCGCTGCCTTTGGCGGGCTGATGGCGGCCATGGTCTATGAAACCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCT‑CATTGGCCATCAG‑CTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTg                                                                                                                                       >  2:1087169/1‑150 (MQ=255)
                ccAGGGGCGCGCTGCCTTTGGCGGGATCATGGCGGACACGGCCGATGAAACCATGCGCCTCAAGATCAGCGAGTGCCGCCCGGTACGCT‑CATTGGCCATCAG‑CTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGa                                                                                                                                     <  1:1738818/150‑1 (MQ=255)
                 cAGGGGCGCGCTGCCTTTGGCGGGCTGATGGCGGCCATGGTCTATGAAACCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCT‑CATTGGCCATCAG‑CTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAg                                                                                                                                    <  1:758381/150‑1 (MQ=255)
                      gcgcgcTGCCTTTGGCGGGCTGATGGCGGCCATGGTCTATGAAACCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCT‑CATTGGCCATCAG‑CTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGtttt                                                                                                                               >  2:2403229/1‑150 (MQ=255)
                      gcgcgcTGCCTTTGGCGGGCTGATGGCGGCCATGGTCTATGAAACCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCT‑CATTGGCCATCAG‑CTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGtttt                                                                                                                               <  2:1757311/150‑1 (MQ=255)
                      gcgcgcTGCCTTTGGCGGGCTGATGGCGGCCATGGTCTATGAAACCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCT‑CATTGGCCATCAG‑CTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGtttt                                                                                                                               <  1:1621024/150‑1 (MQ=255)
                         cgcTGCCTTTGGCGGGCTGATGGCGGCCATGGTCTATGAAACCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCT‑CATTGGCCATCAG‑CTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTgcg                                                                                                                            >  1:1015817/1‑150 (MQ=255)
                            tGCCTTTGGCGGGCTGATGGCGGCCATGGTCTATGAAACCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCT‑CATTGGCCATCAG‑CTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAgg                                                                                                                                         <  2:2222013/134‑1 (MQ=255)
                            tGCCTTTGGCGGGCTGATGGCGGCCATGGTCTATGAAACCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCT‑CATTGGCCATCAG‑CTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAgg                                                                                                                                         >  1:2222013/1‑134 (MQ=255)
                            tGCCTTTGGCGGGCTGATGGCGGCCATGGTCTATGAAACCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCT‑CATTGGCCATCAG‑CTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGa                                                                                                                         >  2:1660970/1‑150 (MQ=255)
                             gCCTTTGGCGGGCTGATGGCGGCCATGGTCTATGAAACCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCT‑CATTGGCCATCAG‑CTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGaa                                                                                                                        <  1:650244/150‑1 (MQ=255)
                             gCCTTTGGCGGGCTGATGGCGGCCATGGTCTATGAAACCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCT‑CATTGGCCATCAG‑CTTCGTGGCAACCGTGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGaa                                                                                                                        <  1:551208/150‑1 (MQ=255)
                                tttGGCGGGCTGATGGCGGCCATGGTCTATGAAACCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCT‑CATTGGCCATCAG‑CTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGaaggc                                                                                                                     >  2:1602953/1‑150 (MQ=255)
                                tttGGCGGGCTGATGGCGGCCATGGTCTATGAAACCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCT‑CATTGGCCATCAG‑CTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGaaggc                                                                                                                     >  2:2203354/1‑150 (MQ=255)
                                        gCTGATGGCGGCCATGGTCTATGAAACCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCT‑CATTGGCCATCAG‑CTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGt                                                                                                             <  1:954143/150‑1 (MQ=255)
                                           gATGGCGGCCATGGTCTATGAAACCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCT‑CATTGGCCATCAG‑CTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAg                                                                                                          <  1:1471163/150‑1 (MQ=255)
                                             tGGCGGCCATGGTCTATGAAACCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCT‑CATTGGCCATCAG‑CTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGaaggcaag                                                                                                                  >  1:1089467/1‑140 (MQ=255)
                                             tGGCGGCCATGGTCTATGAAACCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCT‑CATTGGCCATCAG‑CTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGaaggcaag                                                                                                                  <  2:1089467/140‑1 (MQ=255)
                                               gcggcCATGGTCTATGAAACCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCT‑CATTGG‑AATTAGTGTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACg                                                                                                      <  1:1333400/150‑1 (MQ=255)
                                                 ggcCATGGTCTATGAAACCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCT‑CATTGGCCATCAG‑CTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGtt                                                                                                    >  2:1673200/1‑150 (MQ=255)
                                                   ccATGGTCTATGAAACCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCT‑CATTGGCCATCAG‑CTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGaaggcaa                                                                                                                   >  1:182592/1‑133 (MQ=255)
                                                   ccATGGTCTATGAAACCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCT‑CATTGGCCATCAG‑CTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGaaggcaa                                                                                                                   <  2:182592/133‑1 (MQ=255)
                                                     aTGGTCTATGAAACCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCT‑CATTGGCCATCAG‑CTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTg                                                                                                >  2:2219509/1‑150 (MQ=255)
                                                           tATGAAACCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCT‑CATTGGCCATCAG‑CTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCcgc                                                                                          >  2:724694/1‑150 (MQ=255)
                                                           tATGAAACCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCT‑CATTGGCCATCAG‑CTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCcgc                                                                                          <  2:1837553/150‑1 (MQ=255)
                                                            aTGAAACCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCT‑CATTGGCCATCAG‑CTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCcgcg                                                                                         >  2:407636/1‑150 (MQ=255)
                                                            aTGAAACCATGCGCCTCAAGATCAGCCCCAACCCCCAGGTAAAATGCGACTTTCAT‑AA‑CCTCGTGGCACTCTCGGCGGCGGATGTGCCCATCCGCTTCGAGGGGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCcgcg                                                                                         >  1:1863959/1‑150 (MQ=255)
                                                               aaaCCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCT‑CATTGGCCATCAG‑CTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGaa                                                                                                                        <  2:526257/116‑1 (MQ=255)
                                                               aaaCCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCT‑CATTGGCCATCAG‑CTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGaa                                                                                                                        >  1:526257/1‑116 (MQ=255)
                                                                 aCCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCT‑CATTGGCCATCAG‑CTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGtt                                                                                    >  1:929427/1‑150 (MQ=255)
                                                                  ccATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCT‑CATTGGCCATCAG‑CTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTc                                                                                   <  1:529280/150‑1 (MQ=255)
                                                                   cATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCT‑CATTGGCCATCAG‑CTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTgcgc                                                                                                                           >  2:1335662/1‑109 (MQ=255)
                                                                   cATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCT‑CATTGGCCATCAG‑CTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTgcgc                                                                                                                           <  1:1335662/109‑1 (MQ=255)
                                                                   cATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCT‑CATTGGCCATCAG‑CTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCa                                                                                  <  1:1246843/150‑1 (MQ=255)
                                                                   cATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCT‑CATTGGCCATCAG‑CTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCa                                                                                  <  1:1884671/150‑1 (MQ=255)
                                                                    aTGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCT‑CATTGGCCATCAG‑CTTCGTGGCACCCGCGGCGGCGGATGTGACCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAg                                                                                 >  1:2179614/1‑150 (MQ=255)
                                                                       cgcCTCAAGATCAGCGACAACCGCCCGGTACGCT‑CATTGGCCATCAG‑CTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGAt                                                                              <  1:1087169/150‑1 (MQ=255)
                                                                         ccTCAAGATCAGCGACAACCGCCCGGTACGCT‑CATTGGCCATCAG‑CTTCGTGGCACCAGCGGCGGCGGATGTGGCTTGTCGCTTCTAGGTGGGGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATgg                                                                            >  1:289948/1‑150 (MQ=255)
                                                                          cTCAAGATCAGCGACAACCGCCCGGTACGCT‑CATTGGCCATCAG‑CTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGc                                                                           >  2:1148860/1‑150 (MQ=255)
                                                                          cTCAAGATCAGCGACAACCGCCCGGTACGCT‑CATTGGCCATCAG‑CTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGc                                                                           <  1:1653628/150‑1 (MQ=255)
                                                                           tCAAGATCAGCGACAACCGCCCGGTACGCT‑AATTGGCCATCAG‑CTTCGTGGAACTCGCGGCGGCGGATGTCCCCATCCGGTTCGAAGAGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCGGGGCGGCGCTGTTCAGGATGGcc                                                                          >  1:714746/1‑150 (MQ=255)
                                                                             aaGATCAGCGACAACCGCCCGGTACGCT‑CATTGGCCATCAG‑CTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGtt                                                                                                    >  2:1652872/1‑122 (MQ=255)
                                                                             aaGATCAGCGACAACCGCCCGGTACGCT‑CATTGGCCATCAG‑CTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGtt                                                                                                    <  1:1652872/122‑1 (MQ=255)
                                                                             aaGATCAGCGACAACCGCCCGGTACGCT‑CATTGGCCATCAG‑CTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGCCAg                                                                        >  1:1191587/1‑150 (MQ=255)
                                                                             aaGATCAGCGACAACCGCCCGGTACGCT‑CATTGGCCATCAG‑CTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGCCAg                                                                        <  1:407636/150‑1 (MQ=255)
                                                                             aaGATCAGCGACAACCGCCCGGTACGCT‑CATTGGCCATCAG‑CTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGCCAg                                                                        <  2:1863959/150‑1 (MQ=255)
                                                                                  cAGCGACAACCGCCCGGTACGCT‑CATTGGCCATCAG‑CTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGCCAggtggt                                                                   <  1:954410/150‑1 (MQ=255)
                                                                                    gCGACAACCGCCCGGTACGCT‑CATTGGCCATCAG‑CTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCgcgc                                                                                        >  2:936388/1‑127 (MQ=255)
                                                                                    gCGACAACCGCCCGGTACGCT‑CATTGGCCATCAG‑CTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCgcgc                                                                                        <  1:936388/127‑1 (MQ=255)
                                                                                    gCGACAACCGCCCGGTACGCT‑CATTGGCCATCAG‑CTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGCCAGGTGGTGa                                                                 <  1:1054704/150‑1 (MQ=255)
                                                                                          aCCGCCCGGTACGCT‑CATTGGCCATCAG‑CTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCg                                                                                                               <  2:1607430/98‑1 (MQ=255)
                                                                                          aCCGCCCGGTACGCT‑CATTGGCCATCAG‑CTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCg                                                                                                               >  1:1607430/1‑98 (MQ=255)
                                                                                            cGCCCGGTACGCT‑CATTGGCCATCAG‑CTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGCCAGGTGGTGACTTTGGTg                                                         >  2:200536/1‑150 (MQ=255)
                                                                                            cGCCCGGTACGCT‑CATTGGCCATCAG‑CTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGCCAGGTGGTGACTTTGGTg                                                         >  1:590968/1‑150 (MQ=255)
                                                                                               ccGGTACGCT‑CATTGGCCATCAG‑CTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGCCAGGTGGTGACTTTGGTGCAg                                                      <  1:314291/150‑1 (MQ=255)
                                                                                                 ggTACGCT‑CATTGGCCATCAG‑CTTCGTGGCACCCGCGGCGGCGGAGGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGCCAGGTGGTGACTTTGGTGCAggg                                                    <  2:444190/150‑1 (MQ=255)
                                                                                                     cGCT‑CATTGGCCATCAG‑CTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGCCAGGTGGTGACTTTGGTGCAGGGCAAt                                                <  2:613609/150‑1 (MQ=255)
                                                                                                      gCT‑CATTGGCCATCAG‑CTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGCCAGGTGGTGACTTTGGTGCAGGGCAAtt                                               >  2:2122173/1‑150 (MQ=255)
                                                                                                        t‑CATTGGCCATCAG‑CTTCGTGGCACCGGAATCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGCCAGGTGGTGACTTTGGTGCAGGGCAATTTc                                             >  1:665661/1‑150 (MQ=255)
                                                                                                           aTTGGCCATCAG‑CTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCGAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGCCAGGTGGTGACTTTGGTGCAgg                                                     >  1:213927/1‑140 (MQ=255)
                                                                                                           aTTGGCCATCAG‑CTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCGAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGCCAGGTGGTGACTTTGGTGCAgg                                                     <  2:213927/140‑1 (MQ=255)
                                                                                                           aTTGGCCATCAG‑CTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGCCAGGTGGTGACTTTGGTGCAgg                                                     <  1:670427/140‑1 (MQ=255)
                                                                                                           aTTGGCCATCAG‑CTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGCCAGGTGGTGACTTTGGTGCAgg                                                     >  2:670427/1‑140 (MQ=255)
                                                                                                           aTTGGCCATCAG‑CTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGCCAGGTGGTGACTTTGGTGCAGGGCAATTTCgg                                           >  1:1607453/1‑150 (MQ=255)
                                                                                                            ttGGCCATCAG‑CTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGCCAGGTGGTGACTTTGGTGCAGGGCAATTTCGGt                                          >  1:86196/1‑150 (MQ=255)
                                                                                                            ttGGCCATCAG‑CTTCGTCAGACGCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGCCAGGTGGTGACTTTGGTGCAGGGCAATTTCGGt                                          <  1:1306335/150‑1 (MQ=255)
                                                                                                               gCCATCAG‑CTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTCGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGCCAggtggt                                                                   >  2:2278154/1‑122 (MQ=255)
                                                                                                               gCCATCAG‑CTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTCGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGCCAggtggt                                                                   <  1:2278154/122‑1 (MQ=255)
                                                                                                                ccATCAG‑CTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGCCAGGTGGTGACTTTGGTGCAGGGCAATTTCGGTGCgg                                      <  1:2203354/150‑1 (MQ=255)
                                                                                                                      g‑CTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGCCAGGTGGTGACTTTGGTGCAGGGCAATTTCGGTGCGGgccgcc                                >  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                  <  2:1015817/150‑1 (MQ=255)
                                                                                                                                     gcggcggcgGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGCCAGGTGGTGACTTTGGTGCAGGGCAATTTCGGTGCGGGCCGCCCTTCGGTGGTCGaa                  <  1:671507/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)
                                                                                                                                                    cccATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGc                                                                           <  2:29054/78‑1 (MQ=255)
                                                                                                                                                    cccATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGc                                                                           >  1:29054/1‑78 (MQ=255)
                                                                                                                                                    cccATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGCCAGGTGGTGACTTTGGTGCAGGGCAATTTCGGTGCggg                                     <  1:1675935/116‑1 (MQ=255)
                                                                                                                                                    cccATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGCCAGGTGGTGACTTTGGTGCAGGGCAATTTCGGTGCggg                                     >  2:1675935/1‑116 (MQ=255)
                                                                                                                                                    cccATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGCCAGGTGGTGACTTTGGTGCAGGGCAATTTCGGTGCGGGCGGCCCTTCGTTGGTCGAACTGCCGGCGTTGCCg   >  2:301609/1‑150 (MQ=255)
                                                                                                                                                     ccATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGCCAGGTGGTGACTTTGGTGCAGGGCAATTTCGGTGCGGGCCGCCCTTCGGTGGTCGAAGTGCCGGCGTTGCCgg  >  1:1274054/1‑150 (MQ=255)
                                                                                                                                                     ccANCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGCCAGGTGGTGACTTTGGTTCAGGGCAATTTCGGTGCGGGCCGCCCTTCGGTGGTCGAAGTACCGGCGTTGCCgg  >  2:1519886/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTCCGCCCAGCTGGGCCCAGGGGCGCGCTGCCTTTGGCGGGCTGATGGCGGCCATGGTCTATGAAACCATGCGCCTCAAGATCAGCGACAACCGCCCGGTACGCT‑CATTGGCCATCAG‑CTTCGTGGCACCCGCGGCGGCGGATGTGCCCATCCGCTTCGAGGTGGAGGTTTTGCGCGAAGGCAAGGCGGTTAGCACGTTGCTGGGCCGCGCTGTTCAGGATGGCCAGGTGGTGACTTTGGTGCAGGGCAATTTCGGTGCGGGCCGCCCTTCGGTGGTCGAAGTGCCGGCGTTGCCGG  >  NC_002947/2637993‑2638289

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