Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,630,365 T→C 6.8% W315R (TGG→CGG) ‡ PP_3200 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,630,3650TC 6.8% 249.3 / 6.7 89W315R (TGG→CGG) ‡PP_3200hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (41/41);  new base C (3/3);  total (45/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.24e-01

GCTATTTCTCCACAGACGACTTCTTTGGCGATGGCGCCGAGCGCATGTTCGTCGGCGCGCCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCGAACTTCGCCAACCTCGACCCGGTCAGCGGCAAGGCCGGCGAGTTCTACTGGGTGTATCCCGAGGGTATCG  >  NC_002947/3630216‑3630512
                                                                                                                                                     |                                                                                                                                                   
gCTATTTCTCCACAGACGACTTCTTTGGCGATGGCGCCGAGCGCATGTTCGTCGGCGCGCCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCt                                                                                                                                                     <  1:1506480/150‑1 (MQ=255)
 cTATTTCTCCACAGACGACTTCTTTGGCGATGGCGCCGAGCGCATGTTCGTCGGCGCGCCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTg                                                                                                                                                    <  1:451850/150‑1 (MQ=255)
 cTATTTCTCCACAGACGACTTCTTTGGCGATGGCGCCGAGCGCATGTTCGTCGGCGCGCCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTg                                                                                                                                                    <  2:1487499/150‑1 (MQ=255)
  tATTTCTCCACAGACGACTTCTTTGGCGATGGCGCCGAGCGCATGTTCGTCGGCGCGCCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTgg                                                                                                                                                   <  1:1457762/150‑1 (MQ=255)
        tcCATAGACGACTTCTTTGGCGATGGCGCCGAGCGCATGTTCGTCGGCGCGCCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCAcc                                                                                                                                             <  1:1304679/150‑1 (MQ=255)
         ccACAGACGACTTCTTTGGCGATGGCGCCGAGCGCATGTTCGTCGGCGCGCCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCCAGCGCAcct                                                                                                                                            <  1:411601/150‑2 (MQ=255)
          cacaGACGACTTCTTTGGCGATGGCGCCGAGCGCATGTTCGTCGGCGCGCCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGTGCGGCATGCAGGTGCGCTGGCGCAGCga                                                                                                                                           >  1:398732/1‑150 (MQ=255)
          cacaGACGACTTCTTTGGCGATGGCGCCGAGCGCATGTTCGTCGGCGCGCCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCga                                                                                                                                           <  1:972174/150‑1 (MQ=255)
          cacaGACGACTTCTTTGGCGATGGCGCCGAGCGCATGTTCGTCGGCGCGCCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCga                                                                                                                                           <  2:1601948/150‑1 (MQ=255)
          cacaGACGACTTCTTTGGCGATGGCGCCGAGCGCATGTTCGTCGGCGCGCCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCga                                                                                                                                           >  1:174469/1‑150 (MQ=255)
          cacaGACGACTTCTTTGGCGATGGCGCCGAGCGCATGTTCGTCGGCGCGCCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCga                                                                                                                                           >  1:765013/1‑150 (MQ=255)
          cacaGACGACTTCTTTGGCGATGGCGCCGAGCGCATGTTCGTCGGCGCGCCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCCAGCGCAcctg                                                                                                                                           >  1:1714406/1‑148 (MQ=255)
           acaGACGACTTCTTTGGCGATGGCGCCGAGCGCATGTTCGTCGGCGCGCCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCgag                                                                                                                                          >  1:149826/1‑150 (MQ=255)
            caGACGATTTTTTTGGCGATGGCGCCGAGCGCATGTTCGTCGGCGCGCCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCgaga                                                                                                                                         <  2:47912/150‑1 (MQ=255)
             aGACGACTTCTTTGGCGATGGCGCCGAGCGCATGTTCGTCGGCGCGCCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGAc                                                                                                                                        <  1:782905/150‑1 (MQ=255)
             aGACGACTTCTTTGGCGATGGCGCCGAGCGCATGTTCGTCGGCGCGCCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCCAGCGCAcctgcat                                                                                                                                        <  2:1714406/150‑6 (MQ=255)
              gacgacTTCTTTGGCGATGGCGCCGAGCGCATGTTCGTCGGCGCGCCGGTGTTCCCAAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACg                                                                                                                                       <  1:900522/150‑1 (MQ=255)
                 gacTTCTTTGGCGATGGCGCCGAGCGCATGTTCGTCGGCGCGCCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCgaga                                                                                                                                         <  1:1135740/145‑1 (MQ=255)
                 gacTTCTTTGGCGATGGCGCCGAGCGCATGTTCGTCGGCGCGCCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCgaga                                                                                                                                         >  2:1135740/1‑145 (MQ=255)
                 gacTTCTTTGGCGATGGCGCCGAGCGCATGTTCGTCGGCGCGCCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTc                                                                                                                                    >  1:491071/1‑150 (MQ=255)
                 gacTTCTTTGGCGATGGCGCCGAGCGCATGTTCGTCGGCGCGCCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTc                                                                                                                                    >  1:193389/1‑150 (MQ=255)
                 gacTTCTTTGGCGATGGCGCCGAGCGCATGTTCGTCGGCGCGCCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCACTGGCGCACCGAGACGGTc                                                                                                                                    <  1:1681100/150‑1 (MQ=255)
                    ttcttTGGCGATGGCGCCGAGCGCATGTTCGTCGGCGCGCCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGAc                                                                                                                                 >  1:42736/1‑150 (MQ=255)
                    ttcttTGGCGATGGCGCCGAGCGCATGTTCGTCGGCGCGCCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCCAGCGCAcct                                                                                                                                            >  1:106334/1‑138 (MQ=255)
                    ttcttTGGCGATGGCGCCGAGCGCATGTTCGTCGGCGCGCCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCCAGCGCAcct                                                                                                                                            <  2:106334/139‑2 (MQ=255)
                    ttcttTGGCGATGGCGCCGAGCGCATGTTCGTCGGCGCGCCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGCGCGCTGGCGCCCCGAGACGGTCGAc                                                                                                                                 >  2:923792/1‑150 (MQ=255)
                      cttTGGCGATGGCGCCGAGCGCATGTTCGTCGGCGCGCCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTg                                                                                                                               <  2:1510218/150‑1 (MQ=255)
                       tttGGCGATGGCGCCGAGCGCATGTTCGTCGGCGCGCCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTgg                                                                                                                              <  2:157024/150‑1 (MQ=255)
                        ttGGCGATGGCGCCGAGCGCATGTTCGTCGGCGCGCCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTggg                                                                                                                             <  2:1701852/150‑1 (MQ=255)
                          ggCGATGGCGCCGAGCGCATGTTCGTCGGCGCGCCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGc                                                                                                                           <  1:1395341/150‑1 (MQ=255)
                          ggCGATGGCGCCGAGCGCATGTTCGTCGGCGCGCCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGc                                                                                                                           <  2:193389/150‑1 (MQ=255)
                           gCGATGGCGCCGAGCGCATGTTCGTCGGCGCGCCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGcc                                                                                                                          >  2:1333837/1‑150 (MQ=255)
                                   gcCGAGCGCATGTTCGTCGGCGCGCCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGcc                                                                                                                  >  2:198249/1‑150 (MQ=255)
                                   gcCGAGCGCATGTTCGTCGGCGCGCCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGACAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGcc                                                                                                                  >  1:898055/1‑150 (MQ=255)
                                           cATGTTCGTCGGCGCGCCGGTGTTCCCCAGTGGCCAAATACTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGttt                                                                                                          <  1:1040460/150‑1 (MQ=255)
                                            aTGTTCGTCGGCGCGCCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTggg                                                                                                                             >  2:1466583/1‑130 (MQ=255)
                                            aTGTTCGTCGGCGCGCCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTggg                                                                                                                             <  1:1466583/130‑1 (MQ=255)
                                              gTTCGTCGGCGCGCCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCa                                                                                                       <  2:765013/150‑1 (MQ=255)
                                               ttCGTCGGCGCGCCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCAc                                                                                                      <  2:1802180/150‑1 (MQ=255)
                                                 cgtcgGCGCGCCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGa                                                                                                    <  2:42736/150‑1 (MQ=255)
                                                    cgGCGCGCCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACaa                                                                                                 <  2:1836981/150‑1 (MQ=255)
                                                          gcCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGttt                                                                                                          <  2:1456048/135‑1 (MQ=255)
                                                          gcCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGttt                                                                                                          >  1:1456048/1‑135 (MQ=255)
                                                          gcCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTcgcc                                                                                           <  2:1285943/150‑1 (MQ=255)
                                                          gcCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGACGTTTCCTTATGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTcgcc                                                                                           >  1:797639/1‑150 (MQ=255)
                                                                    cccAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCCAGCGAAGGACTCAGGGCAGGGCGGCATGCAGGGGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAAt                                                                                 >  2:1705778/1‑150 (MQ=255)
                                                                    cccAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCCGGGCCGCGCGGCAAGTCGGTTTTCGAGGTGTCTGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAAt                                                                                 >  1:1670613/1‑150 (MQ=255)
                                                                    cccAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTcgccg                                                                                          >  1:1020615/1‑141 (MQ=255)
                                                                    cccAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTcgccg                                                                                          <  2:1020615/141‑1 (MQ=255)
                                                                    cccAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAAt                                                                                 <  2:1211294/150‑1 (MQ=255)
                                                                     ccAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATg                                                                                >  1:1319226/1‑150 (MQ=255)
                                                                        gTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTAc                                                                             <  2:174469/150‑1 (MQ=255)
                                                                           gCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGcc                                                                          >  1:1798674/1‑150 (MQ=255)
                                                                             cAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCg                                                                        >  1:1567444/1‑150 (MQ=255)
                                                                             cAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCg                                                                        >  1:1612554/1‑150 (MQ=255)
                                                                                 ctctGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCGAACt                                                                    >  2:1600555/1‑150 (MQ=255)
                                                                                      gCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACCCGTCGCCGCAACGCAATGTACGCCCGGCCTTCGcc                                                               >  2:605012/1‑150 (MQ=255)
                                                                                        gTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCGAACTTCGCCaa                                                             >  1:587348/1‑150 (MQ=255)
                                                                                         ttcttcCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAAt                                                                                 >  2:23391/1‑129 (MQ=255)
                                                                                         ttcttcCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAAt                                                                                 <  1:23391/129‑1 (MQ=255)
                                                                                         ttcttcCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCGAACTTCGCCAAc                                                            <  2:1612554/150‑1 (MQ=255)
                                                                                         ttcttcCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCGAACTTCGCCAAc                                                            <  2:494083/150‑1 (MQ=255)
                                                                                             tcCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGAGCCTGTACGCCATCCGTTTCCACGACAAGACGACGCAACTGAAGGCACTCGCAAACATCGCCAACCTCg                                                        >  1:1277989/1‑150 (MQ=255)
                                                                                              cgaCGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCGAACTTCGCCAACCTCGa                                                       >  1:367707/3‑150 (MQ=255)
                                                                                                aCGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCGAACTTCGCCAACCTCGAcc                                                     >  1:277454/1‑150 (MQ=255)
                                                                                                        gACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCGAACTTCGCCAACCTCGACCCGGTCAGc                                             >  2:1710067/1‑150 (MQ=255)
                                                                                                        gACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCGAACTTCGCCAACCTCGACCCGGTCAGc                                             <  1:1259386/150‑1 (MQ=255)
                                                                                                             gcaagcCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCGAACTTCGCCAACCTCGACCCGGTCAGCGGCaa                                        <  2:1319226/150‑1 (MQ=255)
                                                                                                             gcaagcCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACCCCCCAACTTCGCCAACCTCGACCCGGTCCGCGGcca                                        >  2:220150/1‑148 (MQ=255)
                                                                                                             gcaagcCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCCCCGTCACTCAATGTACGCCCGAACTTCGCCAACCTCCACCCGGTCAgtcgggg                                        >  2:100166/1‑144 (MQ=255)
                                                                                                                          tCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCGAACTTCGCCAACCTCGACCCGGTCAGCGGCAAGGCCGGCGAGTTc                           >  1:375216/1‑150 (MQ=255)
                                                                                                                            agggcagggcGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCGAACTTCGCCAACCTCGACCCGGTCAGCGGCAAGGCCGGCGAGTTCTa                         >  2:1437894/1‑150 (MQ=255)
                                                                                                                            agggcagggcGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCGAACTTCGCCAACCTCGACCCGGTCAGCGGCAAGGCCGGCGAGTTCTa                         <  1:1333837/150‑1 (MQ=255)
                                                                                                                              ggcagggcGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCGAACTTCGCCAACCTCGACCCGGTCAGCGGCAAGGCCGGCGAGTTCTACt                       >  2:674299/1‑150 (MQ=255)
                                                                                                                               gcagggcGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCGAACTTCGCCAACCTCGACCCGGTCAGCGGCAAGGCCGGCGAGTTCTACTg                      >  1:1575332/1‑150 (MQ=255)
                                                                                                                               gcagggcGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCGAACTTCGCCAACCTCGACCCGGTCAGCGGCAAGGCCGGCGAGTTCTACTg                      <  2:1575332/150‑1 (MQ=255)
                                                                                                                                cagggcGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCGAACTTCGCCAACCTCGACCCGGTCAGCGGCAAGGCCGGCGAGTTCTACTgg                     <  2:1068744/150‑1 (MQ=255)
                                                                                                                                 agggcGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCGAACTTCGCCAACCTCGACCCGGTCAGCGGCAAGGCCGGCGAGTTCTACTggg                    <  1:1705778/150‑1 (MQ=255)
                                                                                                                                  gggcGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCGAACTTCGCCAACCTCGACCCGGTCAGCGGCAAGGCCGGCGAGTTCTACTGGgt                   <  1:196449/150‑1 (MQ=255)
                                                                                                                                  gggcGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCGAACTTCGCCAACCTCGACCCGGTCAGCGGCAAGGCCGGCGAGTTCTACTGGgt                   <  2:1670613/150‑1 (MQ=255)
                                                                                                                                   ggcggcATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCGAACTTCGCCAGCCTCGACCCGGTCAGCGGCAAGGCCGGCGAGTTCTACTGGgtg                  >  2:1656118/1‑150 (MQ=255)
                                                                                                                                   ggcggcATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCGAACTTCGCCAACCTCGACCCGGTCAGCGGCAAGGCCGGCGAGTTCTACTgggat                  >  1:276835/1‑148 (MQ=255)
                                                                                                                                         aTGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCGAACTTCGCCAACCTCGACCCGGTCAGCGGCAAGGCCGGCGAGTTCTACTGGGTGTATccc            >  1:1543249/1‑150 (MQ=255)
                                                                                                                                          tctcggTGCGCCAGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCGAACTTCGCCAACCTCGACCCGGTCAGCGGCAAGGCCGGCGAGTTCTACTGGGTGTATCCCg           >  1:298045/5‑150 (MQ=255)
                                                                                                                                            cAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCGAACTTCGCCAACCTCGACCCGGTCAGCGGCAAGGCCGGCGAGTTCTACTGGGTGTATCCCGAg         <  1:923792/150‑1 (MQ=255)
                                                                                                                                             aGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCTATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCGAACTTCGCCAACCTCGACCCGGTCAGCGGCAAGGCCGGCGAGTTCTACTGGGTGTATCCCGAgg        >  1:238897/1‑150 (MQ=255)
                                                                                                                                             aGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCGAACTTCGCCAACCTCGACCCGGTCAGCGGCAAGGCCGGCGAGTTCTACTGGGTGTATCCCGAgg        <  2:437077/150‑1 (MQ=255)
                                                                                                                                                tGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCTATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCGAACTTCGCCAACCTCGACCCGGTCAGCGGCAAGGCCGGCGAGTTCTACTGGGTGTATCCCGAGGGTa     >  2:1630257/1‑150 (MQ=255)
                                                                                                                                                   gcTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCGAACTTCGCCAACCTCGACCCGGTCAGCGGCAAGGCCGGCGAGTTCTACTGGGTGTATCCCGAGGGTATCg  <  2:492041/150‑1 (MQ=255)
                                                                                                                                                   gcTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCGAACTTCGCCAACCTCGACCCGGTCAGCGGCAAGGCCGGCGAGTTCTACTGGGTGTATCCCGAGGGTATCg  >  2:983969/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GCTATTTCTCCACAGACGACTTCTTTGGCGATGGCGCCGAGCGCATGTTCGTCGGCGCGCCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCGAACTTCGCCAACCTCGACCCGGTCAGCGGCAAGGCCGGCGAGTTCTACTGGGTGTATCCCGAGGGTATCG  >  NC_002947/3630216‑3630512

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