Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,141,212 T→A 6.1% L44Q (CTG→CAG)  PP_3645 → L‑proline amide hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,141,2120TA 6.1% 353.1 / 12.5 114L44Q (CTG→CAG) PP_3645L‑proline amide hydrolase
Reads supporting (aligned to +/- strand):  ref base T (54/53);  new base A (4/3);  total (58/56)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.74e-01

Only 100 of 127 total aligned reads displayed.

CCTTCAAGGAGCGCCCTCCATGATTCAGCCTGTCGCTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGACCACCTGCAAATTGCACAAGCCACCGTGATCGGC  >  NC_002947/4141063‑4141360
                                                                                                                                                      |                                                                                                                                                    
ccTTCAAGGAGCGCCCTCCATGATTCAGCCTGTCGCTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTtggcct                                                                                                                                                      <  1:1585455/150‑1 (MQ=255)
 cTTCAAGGAGCGCCCTCCATGATTCAGCCTGTCGCTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTtggcctg                                                                                                                                                     >  2:1827895/1‑150 (MQ=255)
 cTTCAAGGAGCGCCCTCCATGATTCAGCCTGTCGCTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTtggcctg                                                                                                                                                     <  2:723877/150‑1 (MQ=255)
  ttCAAGGAGCGCCCTCCATGATTCAGCCTGTCGCTGAACGCACACCTTCCGGCACCAGCTCACTGGATGTCCGAGAGGCCGTTACGTTAACGCTTATTCACCGCGTGGGCCTGAACAAGGAAATGTGGGGCGGACAGTTCGTtggcctgg                                                                                                                                                    >  1:2994705/1‑150 (MQ=255)
  ttCAAGGAGCGCCCTCCATGATTCAGCCTGTCGCTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTtggcctgg                                                                                                                                                    >  1:830941/1‑150 (MQ=255)
  ttCAAGGAGCGCCCTCCATGATTCAGCCTGTCGCTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTtggcctgg                                                                                                                                                    >  1:1113021/1‑150 (MQ=255)
  ttCAAGGAGCGCCCTCCATGATTCAGCCTGTCGCTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTtggcctgg                                                                                                                                                    >  1:764570/1‑150 (MQ=255)
  ttCAAGGAGCGCCCTCCATGATTCAGCCTGTCGCTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTtggcctgg                                                                                                                                                    >  2:868178/1‑150 (MQ=255)
       ggAGCGCCCTCCATGATTCAGCCTGTCGCTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAAc                                                                                                                                               <  2:234292/150‑1 (MQ=255)
             cccTCCATGATTCAGCCTGTCGCTGAACGCGCACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGCAATGTGGGGCGGGCAGTTTGTTGGCCTGGCCAACGACTAc                                                                                                                                         >  2:1363950/1‑150 (MQ=255)
              ccTCCATGATTCAGCCTGTCGCTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTtggcct                                                                                                                                                      <  2:97296/136‑1 (MQ=255)
                  cATGATTCAGCCTGTCGCTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGt                                                                                                                                    >  2:2298266/1‑150 (MQ=255)
                   aTGATTCAGCCTGTCGCTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTc                                                                                                                                   <  2:1184960/150‑1 (MQ=255)
                    tGATTCAGCCTGTCGCTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCa                                                                                                                                  >  1:620913/1‑150 (MQ=255)
                    tGATTCAGCCTGTCGCTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCa                                                                                                                                  >  1:900243/1‑150 (MQ=255)
                      aTTCAGCCTGTCGCTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATc                                                                                                                                <  2:710378/150‑1 (MQ=255)
                       ttCAGCCTGTCGCTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAAc                                                                                                                                               >  2:158213/1‑134 (MQ=255)
                       ttCAGCCTGTCGCTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTAAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAAc                                                                                                                                               <  1:158213/134‑1 (MQ=255)
                          aGCCTTGCGCGGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGTCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCt                                                                                                                            >  2:766993/1‑150 (MQ=255)
                           cccTGTCGGTGAACGCACCCCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACTGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTa                                                                                                                           <  2:459216/149‑1 (MQ=255)
                            ccTGTCGCTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTAc                                                                                                                          <  1:868178/150‑1 (MQ=255)
                            ccTGTCGCTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTAc                                                                                                                          <  1:2101206/150‑1 (MQ=255)
                              tGTCGCTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGa                                                                                                                        <  2:764570/150‑1 (MQ=255)
                              tGTCGCTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGa                                                                                                                        <  1:1827895/150‑1 (MQ=255)
                              tGTCGCTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGa                                                                                                                        <  1:1764616/150‑1 (MQ=255)
                                    tGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCt                                                                                                                  <  2:880440/150‑1 (MQ=255)
                                        cGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCAGGCCAAc                                                                                                                                               >  2:1759438/1‑117 (MQ=255)
                                        cGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCAGGCCAAc                                                                                                                                               <  1:1759438/117‑1 (MQ=255)
                                           acacCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCAc                                                                                                           >  2:945069/1‑150 (MQ=255)
                                             acCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAAc                                                                                                                                               >  2:399997/1‑112 (MQ=255)
                                             acCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACgg                                                                                                         >  2:488259/1‑150 (MQ=255)
                                             acCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACgg                                                                                                         >  1:386452/1‑150 (MQ=255)
                                               cTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTc                                                                                                       >  2:2444844/1‑150 (MQ=255)
                                                    ggCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCgg                                                                                                               <  2:680900/137‑1 (MQ=255)
                                                    ggCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCgg                                                                                                               >  1:680900/1‑137 (MQ=255)
                                                     gCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGcc                                                                                                 >  2:777657/1‑150 (MQ=255)
                                                       aCCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCTTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCcgc                                                                                               <  1:1032155/150‑1 (MQ=255)
                                                       aCCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCcgc                                                                                               >  2:647034/1‑150 (MQ=255)
                                                       aCCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCcgc                                                                                               <  1:647034/150‑1 (MQ=255)
                                                       aCCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCcgc                                                                                               <  1:1445738/150‑1 (MQ=255)
                                                       aCCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCcgc                                                                                               >  2:435470/1‑150 (MQ=255)
                                                        ccAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCcgcg                                                                                              <  1:41490/150‑1 (MQ=255)
                                                         cAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGt                                                                                             >  1:1088909/1‑150 (MQ=255)
                                                                 tGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTAACGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGgccg                                                                                     >  1:2469781/1‑150 (MQ=255)
                                                                    aTGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGccgccg                                                                                  <  2:1420340/150‑1 (MQ=255)
                                                                       tCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGaca                                                                               <  1:845187/150‑1 (MQ=255)
                                                                         ggCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGgccgc                                                                                    <  2:1696847/143‑1 (MQ=255)
                                                                         ggCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGgccgc                                                                                    >  1:1696847/1‑143 (MQ=255)
                                                                          gCCAGGGCCAACTCGTGTTGCTAACCAAGGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGacacac                                                                            <  2:386452/150‑1 (MQ=255)
                                                                          gCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGacacac                                                                            >  1:1476213/1‑150 (MQ=255)
                                                                           ccAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACAcc                                                                           <  1:435470/150‑1 (MQ=255)
                                                                            cAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCg                                                                          >  1:987383/1‑150 (MQ=255)
                                                                               ggCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAg                                                                                                   >  1:1499162/1‑122 (MQ=255)
                                                                                gCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTgg                                                                      >  2:2655019/1‑150 (MQ=255)
                                                                                    aCCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAgg                                                                  <  1:920385/150‑1 (MQ=255)
                                                                                    aCCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAgg                                                                  <  2:457088/150‑1 (MQ=255)
                                                                                       cGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAAc                                                                                                                                               <  1:2750666/70‑1 (MQ=255)
                                                                                       cGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAAc                                                                                                                                               >  2:2750666/1‑70 (MQ=255)
                                                                                        gtggtgCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCaaa                                                                                                    >  1:2196271/1‑112 (MQ=255)
                                                                                         tggtgCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATg                                                             <  1:2655019/150‑1 (MQ=255)
                                                                                          ggtgCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGc                                                            <  2:2019196/150‑1 (MQ=255)
                                                                                            gtgCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTAc                                                                                                                          <  1:1808713/87‑1 (MQ=255)
                                                                                            gtgCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTAc                                                                                                                          >  2:1808713/1‑87 (MQ=255)
                                                                                              gCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGa                                                         <  1:945069/150‑1 (MQ=255)
                                                                                                  aTCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAg                                                     >  1:1203986/1‑150 (MQ=255)
                                                                                                  aTCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAg                                                     >  1:3018885/1‑150 (MQ=255)
                                                                                                      aCGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGc                                                                         <  1:1255272/126‑1 (MQ=255)
                                                                                                       cGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGc                                                >  1:235412/1‑150 (MQ=255)
                                                                                                        ggCGTGGGCATGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGcc                                               <  1:1660536/150‑1 (MQ=255)
                                                                                                           gTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGcc                                                                                                                             >  1:2494459/1‑69 (MQ=255)
                                                                                                           gTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGcc                                                                                                                             <  2:2494459/69‑1 (MQ=255)
                                                                                                               gCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGc                                        <  2:2680366/150‑1 (MQ=255)
                                                                                                               gCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGc                                        <  2:1152223/150‑1 (MQ=255)
                                                                                                                ccTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCt                                       <  2:1626633/150‑1 (MQ=255)
                                                                                                                ccTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCGCGCCGAACTGCt                                       >  1:1397853/1‑150 (MQ=255)
                                                                                                                   gAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCcgc                                                                                               >  2:18422/1‑91 (MQ=255)
                                                                                                                   gAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCcgc                                                                                               <  1:18422/91‑1 (MQ=255)
                                                                                                                      cAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGaccc                                 >  2:312264/1‑149 (MQ=255)
                                                                                                                      cAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGacca                                 >  1:2896247/1‑150 (MQ=255)
                                                                                                                       aaGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACCCACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGaccac                                >  1:47398/1‑150 (MQ=255)
                                                                                                                         ggAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGACCACCt                              <  2:2896247/150‑1 (MQ=255)
                                                                                                                          gAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGACCACCTg                             >  1:1973942/1‑150 (MQ=255)
                                                                                                                           aaaTGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACCCACCGCTGGAAGGCTATGCCGACCAGCTCCCCGAACTGCTTGACCACCTGc                            >  2:1137875/1‑150 (MQ=255)
                                                                                                                                tgGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGACCACCTGCAAAtt                       <  1:777657/150‑1 (MQ=255)
                                                                                                                                   ggCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGACCACCTGCAAATTGca                    <  2:1088909/150‑1 (MQ=255)
                                                                                                                                    gCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGACCACCTGCAAATTGcac                   <  2:830941/150‑1 (MQ=255)
                                                                                                                                    gCGGTCAGTTCGTTGGCCTGCCCAACGACTACCGCGTCATCGCCTACGTCATGCTCGGGCAGGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGACCCCCTGCAAATTGcac                   <  2:1949486/150‑1 (MQ=255)
                                                                                                                                        tCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGa                                             >  2:1489611/1‑120 (MQ=255)
                                                                                                                                        tCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGa                                             <  1:1489611/120‑1 (MQ=255)
                                                                                                                                          tagtCGTTGGCCAGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCt                                                                        >  2:159971/4‑91 (MQ=255)
                                                                                                                                          tagtCGTTGGCCAGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCt                                                                        <  1:159971/88‑1 (MQ=255)
                                                                                                                                          tagtCGTTGGCCAGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGACCACCTGCAAATTGCACAAGCCa             >  1:111995/4‑150 (MQ=255)
                                                                                                                                           gTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGACCACCTGCAAATTGCACAAGCCAc            >  2:856193/1‑150 (MQ=255)
                                                                                                                                            ttCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCcgc                                                                                               <  1:1199972/66‑1 (MQ=255)
                                                                                                                                            ttCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCcgc                                                                                               >  2:1199972/1‑66 (MQ=255)
                                                                                                                                             tCGTTGGCCAGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGAc                                                        <  2:1430716/104‑1 (MQ=255)
                                                                                                                                               gTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGaa                                            >  2:976838/1‑114 (MQ=255)
                                                                                                                                               gTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGaa                                            <  1:976838/114‑1 (MQ=255)
                                                                                                                                                     ctggccAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGACCACCTGCAAATTGCACAAGCCACCGTGATCGGc  <  1:856193/150‑1 (MQ=255)
                                                                                                                                                     ctggccAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGACCACCTGCAAATTGCACAAGCCACCGTGATCGGc  <  2:2181200/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
CCTTCAAGGAGCGCCCTCCATGATTCAGCCTGTCGCTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGACCACCTGCAAATTGCACAAGCCACCGTGATCGGC  >  NC_002947/4141063‑4141360

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: