Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,203,029 C→A 11.1% L430M (CTG→ATG) ‡ PP_4583 → putative Peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,203,0290CA11.1% 318.7 / 23.6 108L430M (CTG→ATG) ‡PP_4583putative Peptidase
Reads supporting (aligned to +/- strand):  ref base C (53/43);  new base A (7/5);  total (60/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.34e-01

Only 100 of 108 total aligned reads displayed.

AAGCGCTCAACCGCCCGGCCCAGGTGCGCCAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTCACATGCGCGCCTGAGCCTGCTGGAGCGCGGCGTGGCCTTT  >  NC_002947/5202883‑5203175
                                                                                                                                                  |                                                                                                                                                  
aaGCGCTCAACCGCCCGGCCCAGGTGCGCCAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTgg                                                                                                                                                 >  1:1841754/1‑150 (MQ=255)
aaGCGCTCAACCGCCCGGCCCAGGTGCGCCAACTGACGTTGCCCACGGGAGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTgg                                                                                                                                                 >  2:2140950/1‑150 (MQ=255)
 aGCGCTCAACCGCCCGGCCCAGGTGCGCCAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGt                                                                                                                                                <  1:801891/150‑1 (MQ=255)
    gcTCAACCGCCCGGCCCAGGTGCGCCAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCACCCGCAGCGGACGGCACCCTGGTGcc                                                                                                                                             <  2:1841754/150‑1 (MQ=255)
     cTCAACCGCCCGGCCCAGGTGCGCCAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCg                                                                                                                                            >  1:2065643/1‑150 (MQ=255)
      tCAACCGCCCGGCCCAGGTGCGCCAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGa                                                                                                                                           >  2:614335/1‑150 (MQ=255)
      tCAACCGCCCGGCCCAGGTGCGCCAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGa                                                                                                                                           <  1:789897/150‑1 (MQ=255)
      tCAACCGCCCGGCCCAGGTGCGCCAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGa                                                                                                                                           >  2:253427/1‑150 (MQ=255)
      tCAACCGCCCGGCCCAGGTGCGCCAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGa                                                                                                                                           >  2:214290/1‑150 (MQ=255)
      tCAACCGCCCGGCCCAGGTGCGCCAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCAGGGTGCCgt                                                                                                                                           <  1:1378776/150‑2 (MQ=255)
             cccGGCCCAGGTGCGCCAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTg                                                                                                                                               <  2:193739/139‑1 (MQ=255)
             cccGGCCCAGGTGCGCCAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTg                                                                                                                                               >  1:193739/1‑139 (MQ=255)
             cccGGCCCAGGTGCGCCAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCt                                                                                                                                    <  2:565492/150‑1 (MQ=255)
             cccGGCCCAGGTGCGCCAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCt                                                                                                                                    <  2:2065643/150‑1 (MQ=255)
             cccGGCCCAGGTGCGCCAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCAGGGTGCCgtc                                                                                                                                          >  1:1303383/1‑142 (MQ=255)
             cccGGCCCAGGTGCGCCAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCAGGGTGCCgtc                                                                                                                                          <  2:1303383/144‑3 (MQ=255)
                  cccAGGTGCGCCAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTAc                                                                                                                               >  1:107834/1‑150 (MQ=255)
                   ccAGGTGCGCCAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTAcg                                                                                                                              >  1:339756/1‑150 (MQ=255)
                      ggTGCGCCAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTAcgccg                                                                                                                           >  2:1659406/1‑150 (MQ=255)
                       gTGCGCCAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCa                                                                                                                                        <  1:2364031/136‑1 (MQ=255)
                       gTGCGCCAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCa                                                                                                                                        >  2:2364031/1‑136 (MQ=255)
                       gTGCGCCAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGACCCGGCTGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTAcgccgc                                                                                                                          <  1:150055/150‑1 (MQ=255)
                              aaCTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAgg                                                                                                                   >  1:1464782/1‑150 (MQ=255)
                              aaCTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAgg                                                                                                                   >  2:2051025/1‑150 (MQ=255)
                                    cGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTgg                                                                                                             <  1:321313/150‑1 (MQ=255)
                                        gCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCaa                                                                                                         <  1:214290/150‑1 (MQ=255)
                                          ccACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGa                                                                                                       >  1:226096/1‑150 (MQ=255)
                                          ccACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGa                                                                                                       >  2:1712077/1‑150 (MQ=255)
                                              gggCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGt                                                                                                   <  1:1659406/150‑1 (MQ=255)
                                                  gcCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCa                                                                                                                                        >  1:709833/1‑109 (MQ=255)
                                                  gcCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCa                                                                                                                                        <  2:709833/109‑1 (MQ=255)
                                                  gcCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCa                                                                                               <  1:485858/150‑1 (MQ=255)
                                                    ccAGGTGGTGCTGCCGCAAACGCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCAct                                                                                             >  2:1192910/1‑150 (MQ=255)
                                                    ccAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCAct                                                                                             >  1:630309/1‑150 (MQ=255)
                                                    ccAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCAct                                                                                             >  1:1709700/1‑150 (MQ=255)
                                                     cAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCActg                                                                                            >  2:716238/1‑150 (MQ=255)
                                                        gtggtgCTCAAGCAACCCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCActgtac                                                                                         <  2:1760218/150‑1 (MQ=255)
                                                              cTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGTCACCCTGGTGCCGCAAACCCTGCTGGGCATTGCCAAGGTTCAGGGCAAGACCGTGCCActgtacctgtac                                                                                   >  1:540554/1‑150 (MQ=255)
                                                              cTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCActgtacctgtac                                                                                   >  1:2344563/1‑150 (MQ=255)
                                                              cTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCActgtacctgtac                                                                                   >  1:1028173/1‑150 (MQ=255)
                                                               tCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGCCAACCGCAGCGGACGGCACCCTGTTACCGATCAGCCTGGTACGCCGCCGCCCGGACGTGGGCAAGACCGTGCCActgtacctgtagg                                                                                  >  2:797179/1‑148 (MQ=255)
                                                                 aaGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGc                                                                                >  1:2144452/1‑150 (MQ=255)
                                                                 aaGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGc                                                                                >  2:2051258/1‑150 (MQ=255)
                                                                   gCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCACTCTCTTTCTCTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTCGTAGGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTa                                                                              <  2:1172779/150‑1 (MQ=255)
                                                                     aaaCCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATg                                                                            <  2:1085889/150‑1 (MQ=255)
                                                                       aCCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGc                                                                                                                                      <  1:853649/90‑1 (MQ=255)
                                                                       aCCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGc                                                                                                                                      >  2:853649/1‑90 (MQ=255)
                                                                         cccGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGc                                                                        >  2:239681/1‑150 (MQ=255)
                                                                           cGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCt                                                                      <  1:253427/150‑1 (MQ=255)
                                                                               gCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATgg                                                                  <  1:2031675/150‑1 (MQ=255)
                                                                                 tCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGa                                                                                                                                           >  2:2199821/1‑75 (MQ=255)
                                                                                        aTTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGcc                                                                                                                                     >  1:457145/1‑74 (MQ=255)
                                                                                        aTTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGcc                                                                                                                                     <  2:457145/74‑1 (MQ=255)
                                                                                         ttCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTc                                                        >  2:1910543/1‑150 (MQ=255)
                                                                                         ttCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTc                                                        <  2:1495072/150‑1 (MQ=255)
                                                                                           cGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGa                                                      <  2:566746/150‑1 (MQ=255)
                                                                                           cGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGa                                                      >  2:1359720/1‑150 (MQ=255)
                                                                                             aTGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTAcgccgcc                                                                                                                         <  1:1391262/81‑1 (MQ=255)
                                                                                             aTGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTAcgccgcc                                                                                                                         >  2:1391262/1‑81 (MQ=255)
                                                                                                  gATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTgg                                               <  1:2172577/150‑1 (MQ=255)
                                                                                                  gATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTgg                                               >  2:2172577/1‑150 (MQ=255)
                                                                                                   aTGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGt                                              <  1:890921/150‑1 (MQ=255)
                                                                                                        tATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTca                                         >  1:720421/1‑150 (MQ=255)
                                                                                                        tATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTca                                         >  1:454597/1‑150 (MQ=255)
                                                                                                         aTGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTcac                                        >  2:439519/1‑150 (MQ=255)
                                                                                                          tGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTgg                                                                                                             >  1:2225598/1‑80 (MQ=255)
                                                                                                          tGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTgg                                                                                                             <  2:2225598/80‑1 (MQ=255)
                                                                                                          tGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTAt                                                                             <  1:802193/112‑1 (MQ=255)
                                                                                                          tGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTAt                                                                             >  2:802193/1‑112 (MQ=255)
                                                                                                            tCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTCACATg                                     >  2:462384/1‑150 (MQ=255)
                                                                                                                tGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTg                                                                                                              >  1:810341/1‑73 (MQ=255)
                                                                                                                tGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTg                                                                                                              <  2:810341/73‑1 (MQ=255)
                                                                                                                   gcgcCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTCACATGCGCGCCt                              <  2:226096/150‑1 (MQ=255)
                                                                                                                     gcCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCGCTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTCACATGCGCGCCTGa                            >  1:1093624/1‑150 (MQ=255)
                                                                                                                     gcCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTCACATGCGCGCCTGa                            <  2:2344563/150‑1 (MQ=255)
                                                                                                                      ccTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCt                                                         <  1:435098/120‑1 (MQ=255)
                                                                                                                      ccTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCt                                                         >  2:435098/1‑120 (MQ=255)
                                                                                                                       cTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGaaa                                                             >  2:1101065/1‑115 (MQ=255)
                                                                                                                       cTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGaaa                                                             <  1:1101065/115‑1 (MQ=255)
                                                                                                                       cTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTCACATGCGCGCCTGAGc                          >  1:2166739/1‑150 (MQ=255)
                                                                                                                             gCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGAccc                                                   >  2:867096/1‑119 (MQ=255)
                                                                                                                             gCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGAccc                                                   <  1:867096/119‑1 (MQ=255)
                                                                                                                             gCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTCACATGCGCGCCTGAGCctgctg                    <  1:1767728/150‑1 (MQ=255)
                                                                                                                                 ccGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCt                                                                      >  2:1312301/1‑96 (MQ=255)
                                                                                                                                 ccGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCt                                                                      <  1:1312301/96‑1 (MQ=255)
                                                                                                                                   gCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTCACATgcgc                                  <  1:1244654/130‑1 (MQ=255)
                                                                                                                                   gCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTCACATgcgc                                  >  2:1244654/1‑130 (MQ=255)
                                                                                                                                    ggctgatcGGCACCAGGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTCACATGCGCGCCTGAGCCTGCTGGAGcgcg             >  2:241337/8‑150 (MQ=255)
                                                                                                                                      ctgatcGGCACCAGGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTCACATGCGCGCCTGAGCCTGCTGGAGCgcggc           >  2:1781630/6‑150 (MQ=255)
                                                                                                                                         atcGGCACCAGGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGaaa                                                             <  2:945928/95‑1 (MQ=255)
                                                                                                                                         atcGGCACCAGGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGaaa                                                             >  1:945928/3‑97 (MQ=255)
                                                                                                                                         atcGGCACCAGGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTtctc                                          <  1:1460206/114‑1 (MQ=255)
                                                                                                                                         atcGGCACCAGGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTtctc                                          >  2:1460206/3‑116 (MQ=255)
                                                                                                                                         atcGGCACCAGGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGAACTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTCACATGCGCGCCTGAGCCTGCTGGAGCGCGGCGTg        >  1:418423/3‑150 (MQ=255)
                                                                                                                                          tcGGCACCAGGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATgg                                                                  <  1:1421820/90‑1 (MQ=255)
                                                                                                                                          tcGGCACCAGGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATgg                                                                  >  2:1421820/2‑91 (MQ=255)
                                                                                                                                            ggCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTCACATGCGCGCCTGAGCCTGCTGGAGCGCGGCGTGGcc     >  2:852860/1‑150 (MQ=255)
                                                                                                                                            ggCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTCACATGCGCGCCTGAGCCTGCTGGAGCGCGGCGTGGcc     >  2:1635286/1‑150 (MQ=255)
                                                                                                                                               aCCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTCACATGCGCGCCTGAGCCTGCTGGAGCGCGGCGTGGCCttt  >  1:1780790/1‑150 (MQ=255)
                                                                                                                                               aCCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTCACATGCGCGCCTGAGCCTGCTGGAGCGCGGCGTGGCCttt  <  2:2355496/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
AAGCGCTCAACCGCCCGGCCCAGGTGCGCCAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTCACATGCGCGCCTGAGCCTGCTGGAGCGCGGCGTGGCCTTT  >  NC_002947/5202883‑5203175

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: