Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,630,366 G→A 7.7% W315* (TGG→TAG) ‡ PP_3200 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,630,3660GA 7.7% 224.2 / 7.3 86W315* (TGG→TAG) ‡PP_3200hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (41/38);  new base A (4/3);  total (45/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.79e-01

CTATTTCTCCACAGACGACTTCTTTGGCGATGGCGCCGAGCGCATGTTCGTCGGCGCGCCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCGAACTTCGCCAACCTCGACCCGGTCAGCGGCAAGGCCGGCGAGTTCTACTGGGTGTATCCCGAGGGTATCGAGG  >  NC_002947/3630217‑3630515
                                                                                                                                                     |                                                                                                                                                     
cTATTTCTCCACAGACGACTTCTTTGGCGATGGCGCCGAGCGCATGTTCGTCGGCGCGCCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTg                                                                                                                                                       <  2:1487499/150‑1 (MQ=255)
cTATTTCTCCACAGACGACTTCTTTGGCGATGGCGCCGAGCGCATGTTCGTCGGCGCGCCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTg                                                                                                                                                       <  1:451850/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:174469/1‑150 (MQ=255)
         cacaGACGACTTCTTTGGCGATGGCGCCGAGCGCATGTTCGTCGGCGCGCCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCga                                                                                                                                              <  2:1601948/150‑1 (MQ=255)
         cacaGACGACTTCTTTGGCGATGGCGCCGAGCGCATGTTCGTCGGCGCGCCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCga                                                                                                                                              >  1:765013/1‑150 (MQ=255)
         cacaGACGACTTCTTTGGCGATGGCGCCGAGCGCATGTTCGTCGGCGCGCCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCga                                                                                                                                              <  1:972174/150‑1 (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                                                                                                                                            >  2:1135740/1‑145 (MQ=255)
                gacTTCTTTGGCGATGGCGCCGAGCGCATGTTCGTCGGCGCGCCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCgaga                                                                                                                                            <  1:1135740/145‑1 (MQ=255)
                gacTTCTTTGGCGATGGCGCCGAGCGCATGTTCGTCGGCGCGCCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTc                                                                                                                                       >  1:193389/1‑150 (MQ=255)
                gacTTCTTTGGCGATGGCGCCGAGCGCATGTTCGTCGGCGCGCCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTc                                                                                                                                       >  1:491071/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                                                                                                                              <  2:193389/150‑1 (MQ=255)
                         ggCGATGGCGCCGAGCGCATGTTCGTCGGCGCGCCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGc                                                                                                                              <  1:1395341/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:1612554/1‑150 (MQ=255)
                                                                            cAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCg                                                                           >  1:1567444/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                                                <  1:1259386/150‑1 (MQ=255)
                                                                                                       gACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCGAACTTCGCCAACCTCGACCCGGTCAGc                                                >  2:1710067/1‑150 (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                            <  1:1333837/150‑1 (MQ=255)
                                                                                                                           agggcagggcGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCGAACTTCGCCAACCTCGACCCGGTCAGCGGCAAGGCCGGCGAGTTCTa                            >  2:1437894/1‑150 (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                      <  2:1670613/150‑1 (MQ=255)
                                                                                                                                 gggcGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCGAACTTCGCCAACCTCGACCCGGTCAGCGGCAAGGCCGGCGAGTTCTACTGGgt                      <  1:196449/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)
                                                                                                                                                     ggCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCGAACTTCGCCAACCTCGACCCGGTCAGCGGCAAGGCCGGCGAGTTCTACTGGGTGTATCCCGAGGGTATCGAgg  <  1:198249/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTATTTCTCCACAGACGACTTCTTTGGCGATGGCGCCGAGCGCATGTTCGTCGGCGCGCCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCGAACTTCGCCAACCTCGACCCGGTCAGCGGCAAGGCCGGCGAGTTCTACTGGGTGTATCCCGAGGGTATCGAGG  >  NC_002947/3630217‑3630515

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: ‑