Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 794,958 T→A 5.7% L366Q (CTG→CAG)  PP_0680 → ATP‑dependent protease‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947794,9580TA 5.7% 333.3 / 8.9 107L366Q (CTG→CAG) PP_0680ATP‑dependent protease‑like protein
Reads supporting (aligned to +/- strand):  ref base T (51/49);  new base A (2/4);  total (54/53)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.86e-01

Only 100 of 109 total aligned reads displayed.

GCGCCGGTGGTGTTCGAGCCGCACCCGACCTATGACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATCGGTGAGCTGGGCCGTGTTGCCACTGTCAGCCTGCAGCCGCAGGTGATTCCGCTGAAC  >  NC_002947/794813‑795106
                                                                                                                                                 |                                                                                                                                                    
 cgcCGGTGGTGTTCGAGCCGCACCCGACCTATGACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAg                                                                                                                                                 <  2:1864673/150‑1 (MQ=255)
  gcCGGTGGTGTTCGAGCCGCACCCGACCTATGACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCAGAATCAGGAAGc                                                                                                                                                <  1:817021/150‑1 (MQ=255)
  gcCGGTGGTGGTCGAGCCGCACCCGACCTATGACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGc                                                                                                                                                <  1:98310/150‑1 (MQ=255)
   ccGGTGGTGTTCGAGCCGCACCCGACCTATGACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCAGAATCAGGAAGcc                                                                                                                                               <  2:2357257/150‑1 (MQ=255)
     ggtggtGTTCGAGCCGCACCCGACCTATGACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCga                                                                                                                                             >  1:1996501/1‑150 (MQ=255)
     ggtggtGTTCGAGCCGCACCCAGCCTATGACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCga                                                                                                                                             <  1:1118472/150‑1 (MQ=255)
       tggtgTTCGAGCCGCACCCGACCTATGACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCgaga                                                                                                                                           >  2:2059156/1‑150 (MQ=255)
         gtgtTCGAGCCGCACCCGACCTATGACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGaaa                                                                                                                                         >  1:2255583/1‑150 (MQ=255)
           gtTCGAGCCGCACCCGACCTATGACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAAt                                                                                                                                       >  1:2387456/1‑150 (MQ=255)
                 gCCGCACCCGACCTATGACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTg                                                                                                                                                     <  1:321754/130‑1 (MQ=255)
                 gCCGCACCCGACCTATGACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTggg                                                                                                                                 >  2:1098784/1‑150 (MQ=255)
                 gCCGCACCCGACCTATGACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTggg                                                                                                                                 >  1:1168346/1‑150 (MQ=255)
                  ccGCACCCGACCTATGACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGt                                                                                                                                <  2:1764060/150‑1 (MQ=255)
                     cACCCGACCTATGACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGaa                                                                                                                                                  <  2:1141110/129‑1 (MQ=255)
                     cACCCGACCTATGACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAAAGGCGGCTTCCTGATTCTGGaa                                                                                                                                                  >  1:1141110/1‑129 (MQ=255)
                             cTATGACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGt                                                                                                                     >  2:526457/1‑150 (MQ=255)
                              tATGACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTa                                                                                                                    <  1:147166/150‑1 (MQ=255)
                              tATGACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTa                                                                                                                    >  2:147166/1‑150 (MQ=255)
                               aTGACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTAt                                                                                                                   <  1:1385237/150‑1 (MQ=255)
                                tGACAATTTGTTTGGCCTTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCCAATCCCCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAGGCCGAGAAAATGCTGGGTGAGCCGTTCGTATg                                                                                                                  <  1:1355058/150‑1 (MQ=255)
                                 gACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTTGGTGAGCCGTTCGTATgg                                                                                                                 >  2:2033265/1‑150 (MQ=255)
                                 gACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATgg                                                                                                                 >  1:2165369/1‑150 (MQ=255)
                                 gACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATgg                                                                                                                 <  1:2239017/150‑1 (MQ=255)
                                    aaTTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATgg                                                                                                                 >  1:15025/1‑147 (MQ=255)
                                    aaTTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATgg                                                                                                                 <  2:15025/147‑1 (MQ=255)
                                    aaTTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGAc                                                                                                              >  2:1087443/1‑150 (MQ=255)
                                      tttgtttgGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTTGTAGCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGAcgc                                                                                                            >  1:1978239/1‑150 (MQ=255)
                                      tttctctgGCCGCATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGAcgc                                                                                                            <  2:1978239/144‑1 (MQ=255)
                                        tgtttgGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACgcgc                                                                                                          <  1:1098784/150‑1 (MQ=255)
                                         gtttgGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGa                                                                                                                              >  2:1128773/1‑129 (MQ=255)
                                         gtttgGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGa                                                                                                                              <  1:1128773/129‑1 (MQ=255)
                                            tgGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGAcgc                                                                                                            <  1:2215727/144‑1 (MQ=255)
                                            tgGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGAcgc                                                                                                            >  2:2215727/1‑144 (MQ=255)
                                                 gTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGTGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAgcgtg                                                                                                 >  1:263796/1‑150 (MQ=255)
                                                 gTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCt                                                                                                                                    <  2:1598278/115‑1 (MQ=255)
                                                 gTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCt                                                                                                                                    >  1:1598278/1‑115 (MQ=255)
                                                  tATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAgcgtg                                                                                                 <  1:1735059/149‑1 (MQ=255)
                                                  tATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAgcgtg                                                                                                 >  2:1735059/1‑149 (MQ=255)
                                                    tCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAgcgtgcgt                                                                                              >  1:2385098/1‑150 (MQ=255)
                                                     ggAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGtt                                                                                             <  2:2385098/149‑1 (MQ=255)
                                                       aaTACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGaaaa                                                                                                                                        <  2:792275/105‑1 (MQ=255)
                                                       aaTACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGaaaa                                                                                                                                        >  1:792275/1‑105 (MQ=255)
                                                         tACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAg                                                                                         >  2:735495/1‑150 (MQ=255)
                                                           cAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGc                                                                                           >  1:1355523/1‑146 (MQ=255)
                                                           cAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGc                                                                                           <  2:1355523/146‑1 (MQ=255)
                                                             gCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTcgc                                                                                     <  2:2148404/150‑1 (MQ=255)
                                                                  gACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTa                                                                                                                    <  1:1945440/114‑1 (MQ=255)
                                                                  gACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTa                                                                                                                    >  2:1945440/1‑114 (MQ=255)
                                                                  gACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAg                                                                                <  2:1536696/150‑1 (MQ=255)
                                                                     cAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTg                                                                                                                                   <  2:1342922/96‑1 (MQ=255)
                                                                     cAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTg                                                                                                                                   >  1:1342922/1‑96 (MQ=255)
                                                                        ggTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTg                                                                                                                                   <  1:105756/93‑1 (MQ=255)
                                                                        ggTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTg                                                                                                                                   >  2:105756/1‑93 (MQ=255)
                                                                         gTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTTGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGa                                                                         <  1:2033265/150‑1 (MQ=255)
                                                                          tGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGAt                                                                        <  1:1087443/150‑1 (MQ=255)
                                                                           gcgcTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCt                                                                                                                                                      <  2:308544/71‑1 (MQ=255)
                                                                           gcgcTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCt                                                                                                                                                      >  1:308544/1‑71 (MQ=255)
                                                                           gcgcTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTg                                                                                            >  1:915300/1‑129 (MQ=255)
                                                                           gcgcTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTg                                                                                            <  2:915300/129‑1 (MQ=255)
                                                                           gcgcTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATg                                                                       <  2:2255583/150‑1 (MQ=255)
                                                                             gcTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGa                                                                     <  2:2165369/150‑1 (MQ=255)
                                                                              cTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATgg                                                                                                                 >  1:960584/1‑105 (MQ=255)
                                                                              cTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATgg                                                                                                                 <  2:960584/105‑1 (MQ=255)
                                                                              cTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAg                                                                    >  2:981237/1‑150 (MQ=255)
                                                                                       tcctATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATc                                                           <  2:1996501/150‑1 (MQ=255)
                                                                                       tcctATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATc                                                           >  1:2249866/1‑150 (MQ=255)
                                                                                            tCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTg                                                                                            >  1:1614876/1‑112 (MQ=255)
                                                                                            tCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTg                                                                                            <  2:1614876/112‑1 (MQ=255)
                                                                                            tCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATCGGTGa                                                      >  2:1385797/1‑150 (MQ=255)
                                                                                                 aaCTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATCGGTGAGCTgg                                                 >  2:217223/1‑150 (MQ=255)
                                                                                                    tGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATCGGt                                                        <  2:33093/140‑1 (MQ=255)
                                                                                                    tGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATCGGt                                                        >  1:33093/1‑140 (MQ=255)
                                                                                                     gcgcCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTcgc                                                                                     >  2:1457180/1‑110 (MQ=255)
                                                                                                     gcgcCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTcgc                                                                                     <  1:1457180/110‑1 (MQ=255)
                                                                                                      ggcCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATCGGTGAGCTGGGCCgt                                            <  1:217223/149‑1 (MQ=255)
                                                                                                      cgcCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATCGGTGAGCTGGGCCgt                                            <  2:2249866/150‑1 (MQ=255)
                                                                                                       gcCCAGGAGCGATGCACCGTTCCAACGGCGGCTTCCTGCTTGCGTGTGGTTAAAGAAGGCTGGGTGAGCCGTTCGGATGTGACGCGCTCAAGCGTGCGTATCAGTCGCGCAAGCTGAAGTTGGAGTCGCCGATCGGTGAGCgggtcagtt                                           >  2:339614/1‑144 (MQ=255)
                                                                                                          cAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATCGGTGAGCTGGGCCGTGTTg                                        >  2:549295/1‑150 (MQ=255)
                                                                                                           aGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATCGGTGa                                                      <  2:292828/135‑1 (MQ=255)
                                                                                                           aGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATCGGTGa                                                      >  1:292828/1‑135 (MQ=255)
                                                                                                                 gcTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGGGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATCGGTGAGCTGGGCCGTGTTGCCACTGt                                 <  1:2140359/150‑1 (MQ=255)
                                                                                                                   tGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAg                                                                    >  2:1356988/1‑113 (MQ=255)
                                                                                                                   tGCACCGGGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTAGGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAg                                                                    <  1:1356988/113‑1 (MQ=255)
                                                                                                                     cACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATCGGTGAGCTGGGCCGTGTTGCCACTGTCAGc                             <  1:526457/150‑1 (MQ=255)
                                                                                                                             cAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATCGGTGAGCTGCAGCGTGGTGCCACTGTCAGCCTGCAGcc                     >  2:1161503/1‑150 (MQ=255)
                                                                                                                              aaCGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATCGGTGAGCTGGGCCGTGTTGCCACTGTCAGCCTGCAGCCg                    >  1:653687/1‑150 (MQ=255)
                                                                                                                                  tcggcTTCCAGAATCAGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATCGGTGAGCTGGGCCGTGTTGCCACTGTCAGCCTGCAGCCGCAgg                >  1:2286121/2‑150 (MQ=255)
                                                                                                                                    ggcTTCCAGAATCAGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTCGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATCGGTGAGCTGGGCCGTGTTGCCACTGTCAGCCTGCAGCCGCAGGTg              >  2:1534094/1‑150 (MQ=255)
                                                                                                                                       ttCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGAAGTCGCGCAAGCTGAAGATGGAGTCGCCGATCGTTGAGCTGGGCCGTGTTGCCACTGTGAGGCTGCAGCCGCAGGTGAtt           >  1:1929576/1‑150 (MQ=255)
                                                                                                                                         ccTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATCGGTGAGCTGGGCCGTGTTGCCACTGTCAGCCTGCAGACGCAGGTGATTcc         >  2:136503/1‑150 (MQ=255)
                                                                                                                                          cTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATg                                                                       <  1:1698549/87‑1 (MQ=255)
                                                                                                                                          cTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATg                                                                       >  2:1698549/1‑87 (MQ=255)
                                                                                                                                          cTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATCGGTGAGCTgg                                                 <  1:623802/109‑1 (MQ=255)
                                                                                                                                          cTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATCGGTGAGCTgg                                                 >  2:623802/1‑109 (MQ=255)
                                                                                                                                          cTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATCGGTGAGCTGGGCCGTGTTGCCACTGTCAGCCTACAGc                      >  1:785053/1‑136 (MQ=255)
                                                                                                                                          cTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATCGGTGAGCTGGGCCGTGTTGCCACTGTCAGCCTACAGc                      <  2:785053/136‑1 (MQ=255)
                                                                                                                                          cTGATTCTGGAAGCCGAGAAAATGCTGGGTAAGCAGTTAGTCTCCGCCACCCACAAGCTGGCGTTGCAGTTGCGCAAGCTGAAGATGGAGTCGCCGATCGGTGAGCTGGGCCGTGTTGCCACTGTCAGCCTGCAGCCGCAGGTGATTCCg        >  1:1023282/1‑150 (MQ=255)
                                                                                                                                           tGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATCGGTGAGCTGGGCCGTGTTGCCACTGTCAGCCTGCAGCCGCAGGTGATTCCGc       <  1:1594813/150‑1 (MQ=255)
                                                                                                                                               tCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATCGGTGAGCTGGGCCGTGTTGCCACTGTCAGCCTGCAGCCGCAGGTGATTCCGCTGaa   <  2:2387456/150‑1 (MQ=255)
                                                                                                                                                cTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATCGGTGAGCTGGGCCGTGTTGCCACTGTCAGCCTGCAGCCGCAGGTGATTCCGCTGAAc  <  1:446357/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GCGCCGGTGGTGTTCGAGCCGCACCCGACCTATGACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATCGGTGAGCTGGGCCGTGTTGCCACTGTCAGCCTGCAGCCGCAGGTGATTCCGCTGAAC  >  NC_002947/794813‑795106

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: