Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,484,044 G→C 5.6% I641M (ATC→ATG)  PP_3091 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,484,0440GC 5.6% 252.4 / 4.2 81I641M (ATC→ATG) PP_3091membrane protein
Reads supporting (aligned to +/- strand):  ref base G (42/34);  new base C (2/3);  total (44/37)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGATCACGCCAGCGTCGGCCGGGGTGCCTTCGCCTGCGGCGCTGGCCCAGGCCTGGGGCAGGCTGGCAGCGACGAAGGCCGGGTCTGGATGGGCCTGGGGCTCGGTCAGCAACAGGTACAGCTTGAGGGTGTTGTAGGCCTCGATGATCGAGGCCACCTGCTGCTCGTCCAGGCGCCCGAGCATTTCTTCGGACAGTGACAGCCCACCGGTGGCGGCGGACAGCTCGGCGCCATCGCTGGCAGTGGCCAGGTTGACCTTGGCGGCATAGGCGCCAGCGCGGGTGTTGCCCAGC  >  NC_002947/3483898‑3484191
                                                                                                                                                  |                                                                                                                                                   
gCGATCACGCCAGCGTCGGCCGGGGTGCCTTCGCCTGCGGCGCTGGCCCAGGCCTGGGGCAGGCTGGCAGCGACGAAGGCCGGGTCTGGATGGGCCTGGGGCTCGGTCAGCAACAGGTACAGCTTGAGGGTGTTGTAGGCCTCGATGATc                                                                                                                                                  <  1:667165/150‑1 (MQ=255)
   aTCACGCCAGCGTCGGCCGGGGTGCCTTCGCCTGCGGCGCTGGCCCAGGCCTGGGGCAGGCTGGCAGCGACGAAGGCCGGGTCTGGATGGGCCTGGGGCTCGGTCAGCAACAGGTACAGCTTGAGGGTGTTGTAGGCCTCGATGATCGAg                                                                                                                                               <  2:604628/150‑1 (MQ=255)
     cACGCCAGCGTCGGCCGGGGTGCCTTCGCCTGCGGCGCTGGCCCAGGCCTGGGGCAGGCTGGCAGCGACGAAGGCCGGGTCTGGATGGGCCTGGGGCTCGGTCAGCAACAGGTACAGCTTGAGGGTGTTGTAGGCCTCGATGATCGAGGc                                                                                                                                             <  2:2180852/150‑1 (MQ=255)
      aCGCCAGCGTCGGCCGGGGTGCCTTCGCCTGCGGCGCTGGCCCAGGCCTGGGGCGGGCTGGCAGCGACGAAGGCCGGGTCTGGATGGGCCTGGGGCTCGGTCAGCAACAGGTACAGCTTGAGGGTGTTGTAGGCCTCGATCATCGAGGcc                                                                                                                                            <  1:1893521/150‑1 (MQ=255)
      aCGCCAGCGTCGGCCGGGGTGCCTTCGCCTGCGGCGCTGGCCCAGGCCTGGGGCAGGCTGGCAGCGACGAAGGCCGGGTCTGGATGGGCCTGGGGCTCGGTCAGCAACAGGTACAGCTTGAGGGTGTTGTAGGCCTCGATGATCGAGGcc                                                                                                                                            >  1:2104570/1‑150 (MQ=255)
         ctgccccaggcctgggccagcgccgcaggcgaaggcaccccggccTGGGGCAGGCTGGCAGCGACGAAGGCCGGGTCTGGATGGGCCTGGGGCTCGGTCAGCAACAGGTACAGCTTGAGGGTGTTGTAGGCCTCGATGATCGAGGCCAcc                                                                                                                                         >  2:44719‑M2/46‑150 (MQ=255)
          cAGCGTCGGCCGGGGTGCCTTCGCCTGCGGCGCTGGCCCAGGCCTGGGGCAGGCTGGCAGCGACGAAGGCCGGGTCTGGATGGGCCTGGGGCTCGGTCAGCAACAGGTACAGCTTGAGGGTGTTGTAGGCCTCGATGATCGAGGCCACct                                                                                                                                        >  1:104561/1‑150 (MQ=255)
            gCGTCGGCCGGGGTGCCTTCGCCTGCGGCGCTGGCCCAGGCCTGGGGCAGGCTGGCAGCGACGAAGGCCGGGTCTGGATGGGCCTGGGGCTCGGTCAGCAACAGGTACAGCTTGAGGGTGTTGTAGGCCTCGATGATCGAGGCCACctgc                                                                                                                                      <  2:104561/150‑1 (MQ=255)
                  GCCGGGGTGCCTTCGCCTGCGGCGCTGGCCCAGGCCTGGGGCAGGCTGGCAGCGACGAAGGCCGGGTCTGGATGGGCCTGGGGCTCGGTCAGCAACAGGTACAGCTTGAGGGTGTTGTAGGCCTCGATGATCGAGGCCACCTGCTGCTCg                                                                                                                                <  1:2027936‑M1/150‑1 (MQ=255)
                    tggagTGCCTTCGCCTGCGGCGCTGGCCCAGGCATGGGGCAGGCTGGCAGCGACGAAGGCCGGGTCTGGATGGGCCTGGGGCTCGGTCAGCAACAGGTACAGCTTGAGGGTGTTGTAGGCCTCGATGATCGAGGCCACCTGCTGCTCGTc                                                                                                                              <  1:594151/146‑1 (MQ=255)
                    ctgggccagcgccgcaggcgaaggcaccccggccTGGGGCAGGCTGGCAGCGACGAAGGCCGGGTCTGGATGGGCCTGGGGCTCGGTCAGCAACAGGTACAGCTTGAGGGTGTTGTAGGCCTCGATGATCGAGGCCACctgctg                                                                                                                                    >  1:10054‑M2/35‑144 (MQ=255)
                    ctgggccagcgccgcaggcgaaggcaccccggccTGGGGCAGGCTGGCAGCGACGAAGGCCGGGTCTGGATGGGCCTGGGGCTCGGTCAGCAACAGGTACAGCTTGAGGGTGTTGTAGGCCTCGATGATCGAGGCCACctgctg                                                                                                                                    <  2:10054‑M2/110‑1 (MQ=255)
                       ggTGCCTTCGCCTGCGGCGCTGGCCCAGGCCTGGGGCAGGCTGGCAGCGACGAAGGCCGGGTCTGGATGGGCCTGGGGCTCGGTCAGCAACAGGTACAGCTTGAGGGTGTTGTAGGCCTCgatga                                                                                                                                                    >  2:280932/1‑125 (MQ=255)
                       ggTGCCTTCGCCTGCGGCGCTGGCCCAGGCCTGGGGCAGGCTGGCAGCGACGAAGGCCGGGTCTGGATGGGCCTGGGGCTCGGTCAGCAACAGGTACAGCTTGAGGGTGTTGTAGGCCTCgatga                                                                                                                                                    <  1:280932/125‑1 (MQ=255)
                         tGCCTTCGCCTGCGGCGCTGGCCCAGGCCTGGGGCAGGCTGGCAGCGACGAATGCCGGGTCTGGATGGGCCTGGGGCTCGGTCAGCAACATGTACAGCTTGAGGGTGTTGTAGGCCTCGATGATCGAGGCCACCTGCTGCTCGTCCAGgc                                                                                                                         >  1:25738/1‑150 (MQ=255)
                            cTTCGCCTGCGGCGCTGGCCCAGGCCTGGGGCAGGCTGGCAGCGACGAAGGCCGGGTCTGGATGGGCCTGGGGCTCGGTCAGCAACAGGTACAGCTTGAGGGTGTTGTAGGCCTCGATGATCGAGGCCACCTGCTGCTCGTCCAGGCGcc                                                                                                                      >  1:186458/1‑150 (MQ=255)
                             ttCGCCTGCGGCGCTGGCCCAGGCCTGGGGCAGGCTGGCAGCGACGAAGGCCGGGTCTGGATGGGCCTGGGGCTCGGTCAGCAACAGGTACAGCTTGAGGGTGTTGTAGGCCTCGATGATCGAGGCCACCTGCTGCTCGTCCAGGCGccc                                                                                                                     <  2:186458/150‑1 (MQ=255)
                                 ccTGCGGCGCTGGCCCAGGCCTGGGGCAGGCTGGCAGCGACGAAGGCCGGGTCTGGATGGGCCTGGGGCTCGGTCAGCAACAGGTACAGCTTGAGGGTGTTGTAGGCCTCGATGATCGAGGCCACCTGCTGCTCGTCCAGGCGCCCGAGc                                                                                                                 <  2:2104570/150‑1 (MQ=255)
                                  cTGCGGCGCTGGCCCAGGCCTGGGGCAGGCTGGCAGCGACGAAGGCCGGGTCTGGATGGGCCTGGGGCTCGGTCAGCAACAGGTACAGCTTGAGGGTGTTGTAGGCCTCGATGATCGAGGCCACCTGCTGCTCGTCCAGGCGCCCGAGCa                                                                                                                >  2:157475/1‑150 (MQ=255)
                                         gcTGGCCCAGGCCTGGGGCAGGCTGGCAGCGACGAAGGCCGGGTCTGGATGGGCCTGGGGCTCGGTCAGCAACAGGTACAGCTTGAGGGTGTTGTAGGCCTCGATGATCGAGGCCACCTGCTGCTCGTCCAGGCGCCCGAGCATttcttc                                                                                                         <  1:1079104/150‑1 (MQ=255)
                                          cTGGCCCAGGCCTGGGGCAGGCTGGCAGCGACGAAGGCCGGGTCTGGATGGGCCTGGGGCTCGGTCAGCAACAGGTACAGCTTGAGGGTGTTGTAGGCCTCGATGATCGAGGCCACCTGCTGCTCGTCCAGGCGCCCGAGCATTTCTTCg                                                                                                        <  2:835414/150‑1 (MQ=255)
                                            ggCCCAGGCCTGGGGCAGGCTGGCAGCGACGAAGGCCGGGTCTGGATGGGCCTGGGGCTCGGTCAGCAACAGGTACAGCTTGAGGGTGTTGTAGGCCTCGATGATCGAGGCCACCTGCTGCTCGTCCAGGCGCCCGAGCATTTCTTCGGa                                                                                                      <  1:2201803/150‑1 (MQ=255)
                                               ccAGGCCTGGGGCAGGCTGGCAGCGACGAAGGCCGGGTCTGGATGGGCCTGGGGCTCGGTCAGCAACAGGTACAGCTTGAGGGTGTTGTAGGCCTCGATGATCGAGGCCACCTGCTGCTc                                                                                                                                 <  2:1103248/120‑1 (MQ=255)
                                               ccAGGCCTGGGGCAGGCTGGCAGCGACGAAGGCCGGGTCTGGATGGGCCTGGGGCTCGGTCAGCAACAGGTACAGCTTGAGGGTGTTGTAGGCCTCGATGATCGAGGCCACCTGCTGCTc                                                                                                                                 >  1:1103248/1‑120 (MQ=255)
                                               ccAGGCCTGGGGCAGGCTGGCAGCGACGAAGGCCGGGTCTGGATGGGCCTGGGGCTCGGTCAGCAACAGGTACAGCTTGAGGGTGTTGTAGGCCTCGATGATCGAGGCCACCTGCTGCTCGTCCAGGCGCCCGAGCATTTCTTCGGACAg                                                                                                   >  2:1408146/1‑150 (MQ=255)
                                               ccAGGCCTGGGGCAGGCTGGCAGCGACGAAGGCCGGGTCTGGATGGGCCTGGGGCTCGGTCAGCAACAGGTACAGCTTGAGGGTGTTGTAGGCCTCGATGATCGAGGCCACCTGCTGCTCGTCCAGGCGCCCGAGCATTTCTTCGGACAg                                                                                                   >  2:98841/1‑150 (MQ=255)
                                                cAGGCCTGGGGCAGGCTGGCAGCGACGAAGGCCGGGTCTGGATGGGCCTGGGGCTCGGTCAGCAACAGGTACAGCTTGAGGGTGTTGTAGGCCTCGATGATCGAGGCCACCTGCTGCTCGTCCAGGCGCCCGAGCATTTCTTCGGACAGt                                                                                                  <  1:958768/150‑1 (MQ=255)
                                                cAGGCCTGGGGCAGGCTGGCAGCGACGAAGGCCGGGTCTGGATGGGCCGGGGGATCGGTCAGCAACAGGGACAGCTTGAGGGTGTTGTAGGCCTCGATGATCGAGGCCACCTGCTGCTCGTCCAGGCGCCCGAGCATTTCTTCGGACAGt                                                                                                  >  1:478279/1‑150 (MQ=255)
                                                        gggCAGGCTGGCAGCGACGAAGGCCGGGTCTGGATGGGCCTGGGGCTCGGTCAGCAACAGGTACAGCTTGAGGGTGTTGTAGGCCTCGATGATCGAGGCCACCTGCTGCTCGTCCAGGCGCCCGAGCATttctt                                                                                                          >  1:1963749/1‑134 (MQ=255)
                                                        gggCAGGCTGGCAGCGACGAAGGCCGGGTCTGGATGGGCCTGGGGCTCGGTCAGCAACAGGTACAGCTTGAGGGTGTTGTAGGCCTCGATGATCGAGGCCACCTGCTGCTCGTCCAGGCGCCCGAGCATttctt                                                                                                          <  2:1963749/134‑1 (MQ=255)
                                                              tgggccAGCGCCGAAGGCCGGGTCTGGATGGGCCTGGGGCTCGGTCAGCAACAGGTACAGCTTGAGGGTGTTGTAGGCCTCGATGATCGAGGCCACCTGCTGCTCGTCCAGGCGCCCGAGCATTTCTTCGGACAGTGACAGCCCACCGGt                                                                                    <  1:835604/145‑1 (MQ=255)
                                                              gCTGGCAGCGACGAAGGCCGGGTCTGGATGGGCCTGGGGCTCGGTCAGCAACAGGTACAGCTTGAGGGTGTTGTAGGCCTCGATGATCGAGGCCACCTGCTGCTCGTCCAGGCGCCCGAGCATTTCTTCGGACAGTGACAGCCCACCGGt                                                                                    >  2:1788703/1‑150 (MQ=255)
                                                                    aGCGACGAAGGCCGGGTCTGGATGGGCCTGGGGCTCGGTCAGCAACAGGTACAGCTTGAGGGTGTTGTAGGCCTCGATGATCGAGGCCACCTGCTGCTCGTCCAGGCGCCCGAGCATTTCTTCGGACAGTGACAGCCCACCGGTggcggc                                                                              >  2:2055106/1‑150 (MQ=255)
                                                                       gacgaAGGCCGGGTCTGGATGGGCCTGGGGCTCGGTCAGCAACAGGTACAGCTTGAGGGTGTTGTAGGCCTCGATGATCGAGGCCACCTGCTGCTCGTCCAGGCGCCCGAGCATTTCTTCGGACAGTGACAGCCCACCGGTGGCGGCGGa                                                                           <  1:1408146/150‑1 (MQ=255)
                                                                           aaGGCCGGGTCTGGATGGGCCTGGGGCTCGGTCAGCAACAGGTACAGCTTGAGGGTGTTGTAGGCCTCGATGATCGAGGCCACCTGCTGCTCGTCCAGGCGCCCGAGCATTTCTTCGGACAGTGACAGCCCACCGGTGGCGGCGGACAGc                                                                       >  1:1291366/1‑150 (MQ=255)
                                                                                   gTCTGGATGGGCCTGGGGCTCGGTCAGCAACAGGTACAGCTTGAGGGTGTTGTAGGCCTCGATGATCGAGGCCACCTGCTGCTCGTCCAGGCGCCCGAGCATTTCTTCGGACAGTGACAGCCCACCGGTGGCGGCGGACAGCTCGGCGcc                                                               >  1:1416909/1‑150 (MQ=255)
                                                                                    tCTGGATGGGCCTGGGGCTCGGTCAGCAACAGGTACAGCTTGAGGGTGTTGTAGTCCTCGATGATCGAGGCCACCTGCTGCTCGTCCAGGCGCCCGAGCATTTCTTCGGACAGTGACAGCCCACCGGTGGCGGCGGACAGCTCGGCGCCa                                                              >  1:1266938/1‑150 (MQ=255)
                                                                                       ggATGGGCCTGGGGCTCGGTCAGCAACAGGTACAGCTTGAGGGTGTTGTAGGCCTCGATGATCGAGGCCACCTGCTGCTCGTCCAGGCGCCCGAGCATTTCTTCGGACAGGCACAGCCCACCGGTGGCGGCGGACAGCTCGGCGCCATCg                                                           <  2:1319655/150‑1 (MQ=255)
                                                                                        gATGGGCCTGGGGCTCGGTCAGCAACAGGTACAGCTTGAGGGTGTTGTAGGCCTCGATGATCGAGGCCACCTGCTGCTCGTCCAGGCGCCCGAGCATTTCTTCGGACAGTGACAGCCCACCGGTGGCGGCGGACAGCTCGGCGCCATCGc                                                          >  2:1521502/1‑150 (MQ=255)
                                                                                        gATGGGCCTGGGGCTCGGTCAGCAACAGGTACAGCTTGAGGGTGTTGTAGGCCTCGATGATCGAGGCCACCTGCTGCTCGTCCAGGCGCCCGAGCATTTCTTCGGACAGTGACAGCCCACCGGTGGCGGCGGACAGCTCGGCGCCATCGc                                                          >  2:1968860/1‑150 (MQ=255)
                                                                                        gATGGGCCTGGGGCTCGGTCAGCAACAGGTACAGCTTGAGGGTGTTGTAGGCCTCGATGATCGAGGCCACCTGCTGCTCGTCCAGGCGCCCGAGCATTTCTTCGGACAGTGACAGCCCACCGGTGGCGGCGGACAGCTCGGCGCCATCGc                                                          >  2:1300833/1‑150 (MQ=255)
                                                                                        gATGGGCCTGGAGCTCGGTCAGCAACAGGTACAGCTTGAGGGTGTTGTAGGCCTCGATGATCGAGGCCACCTGCTGCTCGTCCAGGCGCCCGAGCATTTCTTCGGACAGTGACAGCCCACCGGTGGCGGCGGACAGCTCGGCGCCATCGc                                                          >  2:1617175/1‑150 (MQ=255)
                                                                                          tGGGCCTGGGGCTCGGTCAGCAACAGGTACAGCTTGAGGGTGTTGTAGGCCTCGATGATCGAGGCCACCTGCTGCTCGTCCAGGCGCCCGAGCATTTCTTCGGACAGTGACAGCCCACCGGTGGCGGCGGACAGCTCGGCGCCATCGCTg                                                        >  2:2285663/1‑150 (MQ=255)
                                                                                           gggCCTGGGGCTCGGTCAGCAACAGGTACAGCTTGAGGGTGTTGTAGGCCTCGATGATCGAGGCCACCTGCTGCTCGTCCAGGCGCCCGAGCATTTCTTCGGACAGTGACAGCCCACCGGTGGCGGCGGACAGCTCGGCGCCATAGCTgg                                                       <  1:661433/150‑1 (MQ=255)
                                                                                           gggCCTGGGGCTCGGTCAGCAACAGGTACAGCTTGAGGGTGTTGTAGGCCTCGATGACCGAGGCCACCTGCTGCTGGTCCAGGCGCCCGAGCGTTGCTGCGGACAGTGACAGCCCACCGGTGGCGGCGGACAGCTCGGCGCCATCGCTgg                                                       <  1:2102313/150‑1 (MQ=255)
                                                                                             gCCTGGGGCTCGGTCAGCAACAGGTACAGCTTGAGGGTGTTGTAGGCCTCGATGATCGAGGCCACCTGCTGCTCGTCCAGGCGCCCGAGCATTTCTTCGGACAGTGACAGCCCACCGGTGGCGGCGGACAGCTCGGCGCCATCGCTGGCa                                                     >  1:1943291/1‑150 (MQ=255)
                                                                                               cTGGGGCTCGGTCAGCAACAGGTACAGCTTGAGGGTGTTGTAGGCCTCGATGATCGAGGCCACCTGCTGCTCGTCCAGGCGCCCGAGCATTTCTTCGGACAGTGACAGCCCACCGGTGGCGGCGGACAGCTCGGCGCCATCGCTGGCAGt                                                   >  2:1530365/1‑150 (MQ=255)
                                                                                                tGGGGCTCGGTCAGCAACAGGTACAGCTTGAGGGTGTTGTAGGCCTCGATGATCGAGGCCACCTGCTGCTCGTCCAGGCGCCCGAGCATTTCTTCGGACAGTGACAGCCCACCGGTGGCGGCGGACAGCTCg                                                                    <  1:2046913/132‑1 (MQ=255)
                                                                                                tGGGGCTCGGTCAGCAACAGGTACAGCTTGAGGGTGTTGTAGGCCTCGATGATCGAGGCCACCTGCTGCTCGTCCAGGCGCCCGAGCATTTCTTCGGACAGTGACAGCCCACCGGTGGCGGCGGACAGCTCg                                                                    >  2:2046913/1‑132 (MQ=255)
                                                                                                    gCTCGGTCAGCAACAGGTACAGCTTGAGGGTGTTGTAGGCCTCGATGATCGAGGCCACCTGCTGCTCGTCCAGGCGCCCGAGCATTTCTTCGGACAGTGACAGCCCACCGGTGGCGGCGGACAGCTCGGCGCCATCGCTGGCAGTGGCCa                                              >  2:1681472/1‑150 (MQ=255)
                                                                                                    gCTCGGTCAGCAACAGGTACAGCTGGGGGGTGTTGTAGGCCTCGATGATCGAGGCCACCTGCTGCTCGTCCAGGCGCCCGAGCATTTCTTCGGACAGTGACAGCCCACCGGTGGCGGCGGACAGCTCGGCGCCATCGCTGGCAGTGGccc                                              >  1:1655228/1‑149 (MQ=255)
                                                                                                     cTCGGTCAGCAACAGGTACAGCTTGAGGGTGTTGTAGGCCTCGATGATCGAGGCCACCTGCTGCTCGTCCAGGCGCCCGAGCATTTCTTCGGACAGTGACAGCCCACCGGTGGCGGCGGACAGCTCGGCGCCATCGCTGGCAGTGGTCAg                                             <  2:1266938/150‑1 (MQ=255)
                                                                                                      tCGGTCAGCAACGGGTACAGCTTGAGGGTGTTGTAGGCCTCGATGATCGAGGCCACCTGCTGCTCGTCCAGGCGCCCGAGCATTTCTTCGGACAGTGACAGCCCACCGGTGGCGGCGGACAGCTCGGCGCCATCGCTGGCAGTGGCCAgg                                            >  2:472452/1‑150 (MQ=255)
                                                                                                      tCGGTCAGCAACAGGTACAGCTTGAGGGTGTTGTAGGCCTCGATGATCGAGGCCACCTGCTGCTCGTCCAGGCGCCCGAGCATttctt                                                                                                          >  2:211990/1‑88 (MQ=255)
                                                                                                      tCGGTCAGCAACAGGTACAGCTTGAGGGTGTTGTAGGCCTCGATGATCGAGGCCACCTGCTGCTCGTCCAGGCGCCCGAGCATttctt                                                                                                          <  1:211990/88‑1 (MQ=255)
                                                                                                      tCGGTCAGCAACAGGTACAGCTTGAGGGTGTTGTAGGCCTCGATGATCGAGGCCACCTGCTGCTCGTCCAGGCGCCCGAGCATTTCTTCGGACAGTGACAGCCCACCGGTGGCGGCGGACAGCTCGGCGCCATCGCTGGCAGTGGCCAga                                            <  2:228625/150‑2 (MQ=255)
                                                                                                        ggTCAGCAACAGGTACAGCTTGAGGGTGTTGTAGGCCTCGATGATCGAGGCCACCTGCTGCTCGTCCAGGCGCCCGAGCATTTCTTCGGACAGTGACAGCCCACCGGTGGCGGCGGACAGCTCGGCGCCATCGCTGGCAGTGGCCAGGtt                                          <  1:2285663/150‑1 (MQ=255)
                                                                                                               aaCAGGTACAGCTTGAGGGTGTTGTAGGCCTCGATGATCGAGGCCACCTGCTGCTCGTCCAGGCGCCCGAGCATTTCTTCGGACAGTGACAGCCCACCGGTGGCGGCGGACAGCTCGGCGCCATCGCTGGCAGTGGCCAGGTTGACCTTg                                   <  1:1968860/150‑1 (MQ=255)
                                                                                                                aCAGGTACAGCTTGAGGGTGTTGTAGGCCTCGATGATCGAGGCCACCTGCTGCTCGTCCAGGCGCCCGAGCATTTCTTCGGACa                                                                                                    >  2:1469471/1‑84 (MQ=255)
                                                                                                                aCAGGTACAGCTTGAGGGTGTTGTAGGCCTCGATGATCGAGGCCACCTGCTGCTCGTCCAGGCGCCCGAGCATTTCTTCGGACa                                                                                                    <  1:1469471/84‑1 (MQ=255)
                                                                                                                aCAGGTACAGCTTGAGGGTGTTGTAGGCCTCGATGATCGAGGCCACCTGCTGCTCGTCCAGGCGCCCGAGCATTTCTTCGGACAGTGACAGCCCACCGGTGGCGGCGGACAGCTCGGCGCCATCGCTGGCAGTGGCCAGGTTGACCTTgg                                  <  2:1985358/150‑1 (MQ=255)
                                                                                                                aCAGGTACAGCTTGAGGGTGTTGTAGGCCTCGATGATCGAGGCCACCTGCTGCTCGTCCAGGCGCCCGAGCATTTCTTCGGACAGTGACAGCCCACCGGTGGCGGCGGACAGCTCGGCGCCATCGCTGGCAGTGGCCAGGTTGACCTTgg                                  >  2:1642108/1‑150 (MQ=255)
                                                                                                                   ggTACAGCTTGAGGGTGTTGTAGGCCTCGATGATCGAGGCCACCTGCTGCTCGTCCAGGCGCCCGAGCATTTCTTCGGACAGTGACAGCCCACCGGTGGCGGCGGACAGCTCGGCGCCATCGCTGGCAGTGGCCAGGTTGACCTTggcgg                               <  2:1416909/150‑1 (MQ=255)
                                                                                                                      aCAGCTTGAGGGTGTTGTAGGCCTCGATGATCGAGGCCACCTGCTGCTCGTCCAGGCGCCCGAGCATTTCTTCGGACAGTGACAGCCCACCGGTGGCGGCGGACAGCTCGGCGCCATCGCTGGCAGTGGCCAGGTTGACCTTGGCGGCAt                            <  1:1788703/150‑1 (MQ=255)
                                                                                                                      aCAGCTTGAGGGTGTTGTAGGCCTCGATGATCGAGGCCACCTGCTGCTCGTCCAGGCGCCCGAGCATTTCTTCGGACAGTGACAGCCCACCGGTGGCGGCGGACAGCTCGGCGCCATCGCTGGCAGTGACCAGGTTGACCTTGGCGGCAt                            >  2:2206638/1‑150 (MQ=255)
                                                                                                                        ctcTTGAGGGTGTTGTAGGCCTCGATCTTTTCGGCCGTCTGCTGCTCGTCCAGGCGCCCGAGCATTTCTTCGGACAAAGACAGCCCACCGGTGGCGGAGCAAATGACGGTGGTGGGGCTGGCAGTGGCCAGGTTGACCTTGGCGGCATAg                          <  1:2223637/148‑1 (MQ=255)
                                                                                                                          cTTGAGGGTGTTGTAGGCCTCGATGATCGAGGCCACCTGCTGCTCGTCCAGGCGCCCGAGCATTTCTTCGGACAGTGACAGCCCACCGGTGGCGGCGGACAGCTCGGCGCCATCGCTGGCAGTGGCCAGGTTGACCTTGGCGGCATAGgc                        <  1:98841/150‑1 (MQ=255)
                                                                                                                          cTTGAGGGTGTTGTAGGCCTCGATGATCGAGGCCACCTGCTGCTCGTCCAGGCGCCCGAGCATTTCTTCGGACAGTGACAGCCCACCGGTGGCGGCGGACAGCTCGGCGCCATCGCTGGCAGTGGCCAGGTTGACCTTGGCGGCATAGgc                        >  1:1172810/1‑150 (MQ=255)
                                                                                                                           ttGAGGGTGTTGTAGGCCTCGATGATCGAGGCCACCTGCTGCTCGTCCAGGCGCCCGAGCATTTCTTCGGACAGTGACAGCCCACCGGTGGCGGCGGACAGCTCGGCGCCATCGCTGGCAGTGGCCAGGTTGACCTTGGCGGCATAggag                       >  2:2223637/1‑148 (MQ=255)
                                                                                                                                 gtgtTGTAGGCCTCGATGATCGAGGCCACCTGCTGCTCGTCCAGGCGCCCGAGCATTTCTTCGGACAGTGACAGCCCACCGGTGGCGGCGGACAGCTCGGCGCCATCGCTGGCAGTGGCCAGGTTGACCTTGGCGGCATAGGCGCCAgcg                 >  2:2052811/1‑150 (MQ=255)
                                                                                                                                    ttgtAGGCCTCGATGATCGAGGCCACCTGCTGCTCGTCCAGGCGCCCGAGCATTTCTTCGGACAGTGACAGCCCACCGGTGGCGGCGGACAGCTCGGCGCCATCGCTGGCAGTGGCCAGGTTGACCTTGGCGGCATAGGCGCCAGCGCgg              >  2:559354/1‑150 (MQ=255)
                                                                                                                                      gtAGGCCTCGATGATCGAGGCCACCTGCTGCTCGTCCAGGCGCCCGAGCATTTCTTCGGACAGTGACAGCCCACCGGTGGCGGCGGACAGCTCGGCGCCATCGCTg                                                        <  1:1365936/106‑1 (MQ=255)
                                                                                                                                      gtAGGCCTCGATGATCGAGGCCACCTGCTGCTCGTCCAGGCGCCCGAGCATTTCTTCGGACAGTGACAGCCCACCGGTGGCGGCGGACAGCTCGGCGCCATCGCTg                                                        >  2:1365936/1‑106 (MQ=255)
                                                                                                                                      gtAGGCCTCGATGATCGAGGCCACCTGCTGCTCGTCCAGGCGCCCGAGCATTTCTTCGGACAGTGACAGCCCACCGGTGGCGGCGGACAGCTCGGCGCCATCGCTGGCAGTGGCCAGGTTGACCTTGGCGGCATAGGCGACAGcgcgcat            <  2:646147/150‑4 (MQ=255)
                                                                                                                                         ggCCTCGATCATCGAGGCCACCTGCTGCTCGTCCAGGCGCCCGAGCATTTCTTCGGACAGTGACAGCCCACCGGTGGCGGCGGACAGCTCGGCGCCATCGCTGGCAGTGGCCAGGTTGACCTTggcg                                >  1:2215475/1‑127 (MQ=255)
                                                                                                                                         ggCCTCGATCATCGAGGCCACCTGCTGCTCGTCCAGGCGCCCGAGCATTTCTTCGGACAGTGACAGCCCACCGGTGGCGGCGGACAGCTCGGCGCCATCGCTGGCAGTGGCCAGGTTGACCTTggcg                                <  2:2215475/127‑1 (MQ=255)
                                                                                                                                          gCCTCGATGATCGAGGCCACCTGCTGCTCGTCCAGGCGCCCGAGCATTTCTTCGGACAGTGACAGCCCACCGGTGGCGGCGGACAGCTCGGCGCCATCGCTGGCAGTGGCCAGGTTGACCTTGGCGGCATAGGCGCCAGCGCGGGTGTTg        >  2:1605592/1‑150 (MQ=255)
                                                                                                                                          gCCTCGATGATCGAGGCCACCTGCTGCTCGTCCAGGCGCCCGAGCATTTCTTCGGACAGTGACAGCCCACCGGTGGCGGCGGACAGCTCGGCGCCATCGCTGGCAGTGGCCAGGTTGACCTTGGCGGCATAGGCGCCAGCGCGGGTGTTg        <  1:36463/150‑1 (MQ=255)
                                                                                                                                          gCCTCGATCATCGAGGCCACCTGCTGCTCGTCCAGGCGCCCGAGCATTTCTTCAGACAGTTTCAGCCCACCGGTGGCGGCGGACAGCTCGGCGCCATCGCTGGCAGTGGCCAGGTTGACCTTGGCGGCATAGGCGCCAGCGCGGGTGTTg        >  2:824136/1‑150 (MQ=255)
                                                                                                                                           ccTCGATGATCGAGGCCACCTGCTGCTCGTCCAGGCGCCCGAGCATTTCTTCGGACAGTGACAGCCCACCGGTGGCGGCGGACAGCTCGGCGCCATCGCTGGCAGTGGCCAGGTTGACCTTGGCGGCATAGGCGCCAGCGCGGGTGTTGc       >  1:1637825/1‑150 (MQ=255)
                                                                                                                                            cTCGATGATCGAGGCCACCTGCTGCTCGTCCAGGCGCCCGAGCATTTCTTCGGACAGTGACAGCCCACCGGTGGCGGCGGACAGCTCGGCGCCATCGCTGGCAGTGGCCAGGTTGACCTTGGCGGCATAGGCGCCAGCGCGGGTGTTGcc      >  2:2218582/1‑150 (MQ=255)
                                                                                                                                            cTCGATGATCGAGGCCACCTGCTGCTCGTCCAGGCGCCCGAGCATTTCTTCGGACAGTGACAGCCCACCGGTGGCGGCGGACAGCTCGGCGCCATCGCTGGCAGTGGCCAGGTTGACCTTGGCGGCATAGGCGCCAGCGCGGGTGTTGcc      >  1:1827812/1‑150 (MQ=255)
                                                                                                                                             tCGATGATCGAGGCCACCTGCTGCTCGTCCAGGCGCCCGAGCATTTCTTCGGACAGTGACAGCCCACCGGTGGCGGCGGACAGCTCGGCGCCATCGCTGGCAGTGGCCAGGTTGACCTTGGCGGCATAGGCGCCAGCGCGGGTGTTGccc     >  1:1101984/1‑150 (MQ=255)
                                                                                                                                               gatgatCGAGGCCACCTGCTGCTCGTCCAGGCGCCCGAGCATTTCTTCGGACAGTGACAGCCCACCGGTGGCGGCGGACAGCTCGgcg                                                                 >  2:545213/1‑88 (MQ=255)
                                                                                                                                               gatgatCGAGGCCACCTGCTGCTCGTCCAGGCGCCCGAGCATTTCTTCGGACAGTGACAGCCCACCGGTGGCGGCGGACAGCTCGgcg                                                                 <  1:545213/88‑1 (MQ=255)
                                                                                                                                               gatgatCGAGGCCACCTGCTGCTCGTCCAGGCGCCCGAGCATTTCTTCGGACAGTGACAGCCCACCGGTGGCGGCGGACAGCTCGGCGCCATCGCTGGCAGTGGCCAGGTTGACCTTGGCGGCATAGGCGCCAGCGCGGGTGTTGCCCAg   >  2:673279/1‑150 (MQ=255)
                                                                                                                                                atgatCGAGGCCACCTGCTGCTCGTCCAGGCGCCCGAGCATTTCTTCGGACAGTGACAGCCCACCGGTGGCGGCGGACAGCTCGGCGCCATCGCTGGCAGTGGCCAGGTTGACCTTGGCGGCATAGGCGCCAGCGCGGGTGTTGCCCAgc  >  2:1488811/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GCGATCACGCCAGCGTCGGCCGGGGTGCCTTCGCCTGCGGCGCTGGCCCAGGCCTGGGGCAGGCTGGCAGCGACGAAGGCCGGGTCTGGATGGGCCTGGGGCTCGGTCAGCAACAGGTACAGCTTGAGGGTGTTGTAGGCCTCGATGATCGAGGCCACCTGCTGCTCGTCCAGGCGCCCGAGCATTTCTTCGGACAGTGACAGCCCACCGGTGGCGGCGGACAGCTCGGCGCCATCGCTGGCAGTGGCCAGGTTGACCTTGGCGGCATAGGCGCCAGCGCGGGTGTTGCCCAGC  >  NC_002947/3483898‑3484191

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