Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,768,103 G→C 8.8% P252A (CCT→GCT) ‡ PP_4218 ← lipase/esterase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,768,1030GC 8.8% 348.6 / 21.6 113P252A (CCT→GCT) ‡PP_4218lipase/esterase family protein
Reads supporting (aligned to +/- strand):  ref base G (52/51);  new base C (4/6);  total (56/57)
Fisher's exact test for biased strand distribution p-value = 7.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.03e-01

Only 100 of 116 total aligned reads displayed.

CTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAGCGCTGCACCGAAGGCCG  >  NC_002947/4767954‑4768250
                                                                                                                                                     |                                                                                                                                                   
ctcGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAgg                                                                                                                                                     >  1:2161407/1‑150 (MQ=255)
cgcGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAgg                                                                                                                                                     >  2:1543961/3‑150 (MQ=255)
 tcGATGATCGAACCCATACGCATGAAGTCGTGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAggg                                                                                                                                                    >  1:601621/1‑149 (MQ=255)
           aaCCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCAACAGTGCACCTGCACTGATGc                                                                                                                                          <  2:548064/150‑1 (MQ=255)
           aaCCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCACAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGc                                                                                                                                          <  2:2280253/150‑1 (MQ=255)
           aaCCCATACGCATGAAGTCATGGGTCACTCCCGGGACCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGc                                                                                                                                          <  1:2343198/150‑1 (MQ=255)
            aCCCATACGCATGAAGTCATGGGTCACTCCCGGGCTCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCt                                                                                                                                         <  2:2364667/150‑1 (MQ=255)
              ccATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCACACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTgg                                                                                                                                       >  1:799346/1‑150 (MQ=255)
                aTACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCAACAGTGCACCTGCACTGATGCTgatt                                                                                                                                     <  1:1259618/150‑4 (MQ=255)
                     cATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCt                                                                                                                                >  2:1420927/1‑150 (MQ=255)
                       tGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGc                                                                                                                              >  2:2135414/1‑150 (MQ=255)
                        gAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTg                                                                                                                                                   <  2:2215220/128‑1 (MQ=255)
                           gTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAg                                                                                                                          >  2:825822/1‑150 (MQ=255)
                               tGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAAc                                                                                                                      <  2:881577/150‑1 (MQ=255)
                               tGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAAc                                                                                                                      <  1:1055866/150‑1 (MQ=255)
                                         ccGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAAc                                                                                                            >  2:1293408/1‑150 (MQ=255)
                                         ccGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAAc                                                                                                            >  1:792034/1‑150 (MQ=255)
                                         ccGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAAc                                                                                                            >  1:1626508/1‑150 (MQ=255)
                                           gggATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCAATATGGACACCCTGATCATGCAGCGTGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAAcgc                                                                                                          <  2:1775209/150‑1 (MQ=255)
                                             gATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCgg                                                                                                        <  2:1626508/150‑1 (MQ=255)
                                                cACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATc                                                                                                     <  2:725904/150‑1 (MQ=255)
                                                 aCCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTTTTGCTGGCGACGCTGCCCAGCAACGGGGAAAAACGCggctag                                                                                                    >  2:2343538/1‑146 (MQ=255)
                                                 aCCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGt                                                                                                                   >  1:657089/1‑135 (MQ=255)
                                                 aCCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGt                                                                                                                   <  2:657089/135‑1 (MQ=255)
                                                   cGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTg                                                                                                  >  2:2161776/1‑150 (MQ=255)
                                                   cGCAAGCTGCACCGCCACCCCGCCTCGCTCCGAGTGCCGTGCATAAGAGACACCCTGGTCATGCAGCGGGGGGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGACGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTg                                                                                                  >  2:33213/1‑150 (MQ=255)
                                                    gCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGa                                                                                                 <  2:804621/150‑1 (MQ=255)
                                                    gCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGa                                                                                                 <  2:2440597/150‑1 (MQ=255)
                                                          tGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCAACAGTGCACCTGCACTgatg                                                                                                                                           >  1:1260605/1‑99 (MQ=255)
                                                          tGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCAACAGTGCACCTGCACTGATg                                                                                                                                           <  2:1260605/102‑1 (MQ=255)
                                                           gCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCt                                                                                                                                         >  2:638918/1‑103 (MQ=255)
                                                           gCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCt                                                                                                                                         <  1:638918/103‑1 (MQ=255)
                                                           gCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAAcgcg                                                                                                         >  2:1683245/1‑135 (MQ=255)
                                                            cACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCCAATTAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGc                                                                                         >  1:1307298/1‑150 (MQ=255)
                                                             aCCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCTACGCCCTGGTCATGCAGCGGGTCGCACTCGACTCGCATCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGcc                                                                                        <  1:756581/150‑1 (MQ=255)
                                                             aCCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGcc                                                                                        <  1:1293408/150‑1 (MQ=255)
                                                             aCCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCAACAGTGCACCTGCACTGATGc                                                                                                                                          <  2:1000431/100‑1 (MQ=38)
                                                             aCCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCAACAGTGCACCTGCACTGATGc                                                                                                                                          >  1:1000431/1‑100 (MQ=38)
                                                                     cccGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGc                                                                                         <  1:1688561/141‑1 (MQ=255)
                                                                     cccGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGAGAAAAACGCGGATCGTGACGGTCTGc                                                                                         >  2:1688561/1‑141 (MQ=255)
                                                                      ccGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTg                                                                               <  2:480385/150‑1 (MQ=255)
                                                                        ggCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGc                                                                             <  2:1113452/150‑1 (MQ=255)
                                                                          cTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCACCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGt                                                                           >  2:515491/1‑150 (MQ=255)
                                                                             gCTCCAAGCGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTAtt                                                                        <  1:18628/150‑1 (MQ=255)
                                                                               tCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGcc                                                                                        <  1:748377/132‑1 (MQ=255)
                                                                               tCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGcc                                                                                        >  2:748377/1‑132 (MQ=255)
                                                                               tCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCCCCGTGGCGTATTGc                                                                      >  1:973782/1‑150 (MQ=255)
                                                                               tCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGc                                                                             >  1:952291/1‑143 (MQ=255)
                                                                               tCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGGGACGGTGTGCCGGCACCGTGGCGTATTGc                                                                      >  2:29881/1‑150 (MQ=255)
                                                                                ccAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCt                                                                     <  1:1635014/150‑1 (MQ=255)
                                                                                 cAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTTGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTg                                                                    <  2:601621/150‑1 (MQ=255)
                                                                                  aaGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTgg                                                                   <  1:2111390/150‑1 (MQ=255)
                                                                                  aaGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTgg                                                                   <  2:2124532/150‑1 (MQ=255)
                                                                                       ccGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCACACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGaa                                                              <  2:799346/150‑1 (MQ=255)
                                                                                        cGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAAc                                                             <  1:1385218/150‑1 (MQ=255)
                                                                                          tGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCa                                                           >  2:661561/1‑150 (MQ=255)
                                                                                           gCATAAGCGACACCCTGGTGATGCAGGGGGTCGCACTCGGCAATCAGGATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAGAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCTTATTGCTGGTAGAACCAc                                                          <  2:826269/150‑1 (MQ=255)
                                                                                             aTAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGatcggaag                                                        >  1:1717338/1‑143 (MQ=255)
                                                                                              tAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTcc                                                       >  1:987359/1‑150 (MQ=255)
                                                                                                 gCGACACCCTGGTCATGCAGCGGGTCGCATTCGGCGCCCGGCTTCAGTGTATGTGCACTGTTGCTGGCGCAGTATTCTAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAgg                                                    <  1:532504/150‑1 (MQ=255)
                                                                                                 gCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAgg                                                    >  2:1125/1‑150 (MQ=255)
                                                                                                   gACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGt                                                  >  1:1031466/1‑150 (MQ=255)
                                                                                                    acacCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGccc                                                                                                                            <  1:1797858/75‑1 (MQ=255)
                                                                                                    acacCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTc                                                 <  1:1543961/150‑1 (MQ=255)
                                                                                                    acacCCTGGTCATGACGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGccc                                                                                                                            >  2:1797858/1‑75 (MQ=255)
                                                                                                      acCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTg                                                                       <  2:902311/126‑1 (MQ=255)
                                                                                                      acCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTg                                                                       >  1:902311/1‑126 (MQ=255)
                                                                                                        ccTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGtt                                                                                                                                            <  2:2251268/55‑1 (MQ=255)
                                                                                                        ccTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGtt                                                                                                                                            >  1:2251268/1‑55 (MQ=255)
                                                                                                        ccTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTgcg                                             <  2:752887/150‑1 (MQ=255)
                                                                                                        ccTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTgcg                                             >  2:385193/1‑150 (MQ=255)
                                                                                                         cTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTgg                                                                   >  2:831866/1‑127 (MQ=255)
                                                                                                         cTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTgg                                                                   <  1:831866/127‑1 (MQ=255)
                                                                                                            gTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGGAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTc                                         <  1:2343538/150‑1 (MQ=255)
                                                                                                            gTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGaa                                                              >  1:2109134/1‑129 (MQ=255)
                                                                                                               aTGCAGCGGGACCCACCCGGCATTCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAg                                      >  2:2112276/1‑150 (MQ=255)
                                                                                                                  cAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAg                                   <  1:2161776/150‑1 (MQ=255)
                                                                                                                           gCACTCGGCAATCAGCATCAGTGCAGGTGCGCTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCt                          >  2:77970/1‑150 (MQ=255)
                                                                                                                               tCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCa                      >  2:890322/1‑150 (MQ=255)
                                                                                                                               tCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGGGAAAAACGCGGAGCGTGACGGTCAGGAGGCACCTTGGCTTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAGCCGCTGCCa                      >  1:1750394/1‑150 (MQ=255)
                                                                                                                                 ggCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATa                    >  1:1754016/1‑150 (MQ=255)
                                                                                                                                  gCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAg                   >  2:520878/1‑150 (MQ=255)
                                                                                                                                      tCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAGCGCt               >  2:1758436/1‑150 (MQ=255)
                                                                                                                                      ccAGCAACAGTGCACCTGCACTGATGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGc                                     <  1:2396069/127‑1 (MQ=255)
                                                                                                                                      ccAGCAACAGTGCACCTGCACTGATGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGc                                     >  2:2396069/2‑128 (MQ=255)
                                                                                                                                        aGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATc                                                                                                     >  2:1974471/1‑62 (MQ=255)
                                                                                                                                        aGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATc                                                                                                     <  1:1974471/62‑1 (MQ=255)
                                                                                                                                        aGCAACAGTGCACCTGCACTGATGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGAGGCGTATTGCTGGGAGAACCACTCCAGGGTCTGCGCTTCaa                                       >  2:1304754/1‑124 (MQ=255)
                                                                                                                                         gCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCATGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAGCGCTGCa            >  1:466750/1‑150 (MQ=255)
                                                                                                                                          cATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAGCGCTGCAc           <  1:2135414/150‑1 (MQ=255)
                                                                                                                                           aTCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAGCGCTGCAcc          >  2:2233327/1‑150 (MQ=255)
                                                                                                                                             cAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAACAACACAGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAGCGCTGCACCGa        >  1:2274721/1‑150 (MQ=255)
                                                                                                                                             cAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAAc                                                             <  2:517270/97‑1 (MQ=255)
                                                                                                                                             cAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAAc                                                             >  1:517270/1‑97 (MQ=255)
                                                                                                                                               gTGCAGGTGCGCTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAGCGCTGCACCGAAg      <  1:77970/150‑1 (MQ=255)
                                                                                                                                               gTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAGCGCTGCACCGAAg      <  2:2161407/150‑1 (MQ=255)
                                                                                                                                               gTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAGCGCTGCACCGAAg      <  1:1336140/150‑1 (MQ=255)
                                                                                                                                                tGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAGCGCTGCACCGAAgg     >  1:2345346/1‑150 (MQ=255)
                                                                                                                                                tGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCGGCGCTTCAAGCAGGTAACCGCTGCCATAGCGCTGCACCGAAgg     <  1:825822/150‑1 (MQ=255)
                                                                                                                                                   aGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAGCGCTGCACCGAAGgccg  >  1:109096/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAGCGCTGCACCGAAGGCCG  >  NC_002947/4767954‑4768250

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