Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,777,050 C→G 8.0% W166C (TGG→TGC)  PP_2429 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,777,0500CG 8.0% 198.5 / 4.6 68W166C (TGG→TGC) PP_2429membrane protein
Reads supporting (aligned to +/- strand):  ref base C (31/30);  new base G (4/2);  total (35/33)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGAGGCGTTGACACGTTGAGCCGAGGCCGCTCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAA‑ACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCATCAACCAGACCTTGGGGT  >  NC_002947/2776905‑2777198
                                                                                                                                                  |                                                                                                                                                     
cgAGGCGTTGACACGTTGAGCCGAGGCCGCTCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAA‑ACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCAcc                                                                                                                                                   >  2:138658/1‑150 (MQ=255)
 gAGGCGTTGACACGTTGAGCCCAGGCCGCTCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAA‑ACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGGCACC‑g                                                                                                                                                 >  1:376524/1‑150 (MQ=255)
      gTTGACACGTTGAGCCGAGGCCGCTCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAA‑ACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑gg                                                                                                                                                >  2:714524/1‑146 (MQ=255)
      gTTGACACGTTGAGCCGAGGCCGCTCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAA‑ACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑gg                                                                                                                                                <  1:714524/146‑1 (MQ=255)
        tGACCCGTTGTGCCTGGCCCGCGAAGGCCCCTTGTCGGTAAAGTAGGGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAA‑ACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTc                                                                                                                                          <  2:376524/150‑1 (MQ=255)
              gTTGAGCCGAGGCCGCTCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAA‑ACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCAcc                                                                                                                                    >  2:218921/1‑150 (MQ=255)
                   gCCGAGGCCGCTCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAA‑ACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCaa                                                                                                                               >  1:526122/1‑150 (MQ=255)
                   gCCGAGGCCGCTCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAA‑ACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCaa                                                                                                                               >  1:4267/1‑150 (MQ=255)
                   gCCGAGGCCCCTCAGGCCCCGTGTCGTTAAAGTAGAGGCTCAGGCCAATACTGCAGCCCATACCGATATCAACAAA‑ACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCaa                                                                                                                               <  2:90049/150‑1 (MQ=255)
                         gCCGCTCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAA‑ACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGc                                                                                                                         >  2:221284/1‑150 (MQ=255)
                           cGCTCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAA‑ACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGccc                                                                                                                       >  2:620973/1‑150 (MQ=255)
                                  gCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAA‑ACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGCCGCTCCCCCCCCCCAACACCGCCCAggcggc                                                                                                                <  2:526122/150‑1 (MQ=255)
                                   ccccTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAA‑ACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGcc                                                                                                               <  1:568421/150‑1 (MQ=255)
                                    cccTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAA‑ACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCa                                                                                                              <  1:138658/150‑1 (MQ=255)
                                       ttGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAA‑ACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGc                                                                                                           >  2:154560/1‑150 (MQ=255)
                                       ttGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAA‑ACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGc                                                                                                           >  1:655551/1‑150 (MQ=255)
                                         gTCGGTCAAGTAGAGGCTGAGGCCAATACTGCAGCCCATGCCGATA‑CAAGAAATAGTCGCTGCAACCACAGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTc                                                                                                         <  2:839790/150‑1 (MQ=255)
                                             gTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAA‑ACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCGCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCa                                                                                                     >  1:93159/1‑150 (MQ=255)
                                              tCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAA‑ACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCag                                                                                                    >  2:531286/1‑150 (MQ=255)
                                                  gTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAA‑ACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCg                                                                                                >  1:354984/1‑150 (MQ=255)
                                                      agGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAA‑ACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCATCACCGCCCAGGCGGCCATGCTCTGCCGAGcatggc                                                                                            <  2:680195/150‑6 (MQ=255)
                                                        gCTTAGGCAAATACTGCAGCCCATGCCCCTATCAACAAA‑ACCCCCGCCAACACTTGGTTGAACACACGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCACAACCCCAACATCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGc                                                                                          <  2:51109/150‑1 (MQ=255)
                                                        gCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAA‑ACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGc                                                                                          >  1:260748/1‑150 (MQ=255)
                                                          tCAGGCCAATACTGCAGCCCATGCCGATATCAACAAA‑ACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGc                                                                                                 <  1:388007/141‑1 (MQ=255)
                                                          tCAGGCCAATACTGCAGCCCATGCCGATATCAACAAA‑ACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGc                                                                                                 >  2:388007/1‑141 (MQ=255)
                                                             ggCCAATACTGCAGCCCATGCCGATATCAACAAA‑ACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCTTTTCCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGaa                                                                                     <  2:422848/150‑1 (MQ=255)
                                                                 aaTACTGCAGGCCATGCCGATATCAACAAA‑ACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCCTGCCCGGCAGAGCGATCAGCAAGAACAcc                                                                                 >  2:491404/1‑150 (MQ=255)
                                                                 aaTACTGCAGCCCATGCCGATATCAACAAA‑ACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCAACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACAcc                                                                                 >  2:667002/1‑150 (MQ=255)
                                                                             cATGCCGATATCAACAAA‑ACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGc                                                                     >  1:268120/1‑150 (MQ=255)
                                                                              aTGCCGATATCAACAAA‑ACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGcc                                                                    >  1:638893/1‑150 (MQ=255)
                                                                                gCCGATATCAACAAA‑ACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATc                                                                                             >  2:297640/1‑123 (MQ=255)
                                                                                gCCGATATCAACAAA‑ACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATc                                                                                             <  1:297640/123‑1 (MQ=255)
                                                                                  cGATATCAACAAA‑ACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCa                                                                                                     <  1:635642/113‑1 (MQ=255)
                                                                                  cGATATCAACAAA‑ACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCa                                                                                                     >  2:635642/1‑113 (MQ=255)
                                                                                  cGATATCAACAAA‑ACCCCCGCCAACACTTGGTTGAACACCCCCACCCGCCCCAGTGCCTGCCGTCACC‑GGGCGCTCCCCACCCCCAACTCCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCacac                                                                <  2:736019/150‑1 (MQ=255)
                                                                                    atatCAACAAA‑ACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGAt                                                                                              >  1:805981/1‑118 (MQ=255)
                                                                                    atatCAACAAA‑ACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGAt                                                                                              <  2:805981/118‑1 (MQ=255)
                                                                                      atCAACAAT‑ACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACTCTGGc                                                            <  2:384041/150‑1 (MQ=255)
                                                                                      atCAACAAA‑ACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCAATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGc                                                            <  1:227288/150‑1 (MQ=255)
                                                                                        caacaaA‑ACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCaa                                                          >  2:126973/1‑150 (MQ=255)
                                                                                         aacaaA‑ACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCa                                                                                            <  2:438537/115‑1 (MQ=255)
                                                                                         aacaaA‑ACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCCGCCATGCTCGGCAGAGCGATCa                                                                                            >  1:438537/1‑115 (MQ=255)
                                                                                          acaaA‑ACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAacg                                                        >  1:439582/1‑150 (MQ=255)
                                                                                            aaa‑aCCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCCGCCACACTGGCAacgac                                                      <  2:439582/150‑1 (MQ=255)
                                                                                                      gCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCa                                                                                >  2:665530/1‑115 (MQ=255)
                                                                                                      gCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCa                                                                                <  1:665530/115‑1 (MQ=255)
                                                                                                             cTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATAccc                                      <  2:268120/150‑1 (MQ=255)
                                                                                                             cTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATAccc                                      <  1:154560/150‑1 (MQ=255)
                                                                                                                   tGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGacca                                >  2:348578/1‑150 (MQ=255)
                                                                                                                     aaCACGCGCAACCGCCCGGGTGCCTGCAGGCACC‑GGGCGCTCCCCGCCCCCAACACCGCCCAGGCGGCCATGCTCGGGAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAATGTTGGCGCGGCTGATACCGCGaccacc                              <  2:506588/150‑1 (MQ=255)
                                                                                                                           cGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCAAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACg                        >  1:625976/1‑150 (MQ=255)
                                                                                                                                 cGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCATc                  <  2:638893/150‑1 (MQ=255)
                                                                                                                                    cccGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCATCAAc               <  1:126973/150‑1 (MQ=255)
                                                                                                                                    cccGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCATCAAc               <  2:260748/150‑1 (MQ=255)
                                                                                                                                     ccGGTGGCTGCAGGCACCGGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCATCAAc               >  1:664572/1‑150 (MQ=255)
                                                                                                                                     ccGGTGGCTGCAGGCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTg                                                              <  1:771969/102‑1 (MQ=255)
                                                                                                                                     ccGGTGGCTGCAGGCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTg                                                              >  2:771969/1‑102 (MQ=255)
                                                                                                                                     ccGGTGGCTGCAGCCACCGGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCAGCAAc               >  2:477128/1‑150 (MQ=255)
                                                                                                                                     ccGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAgc                                                                              >  2:367331/1‑86 (MQ=255)
                                                                                                                                     ccGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAgc                                                                              <  1:367331/86‑1 (MQ=255)
                                                                                                                                      cGGTGGCTGCAGGCACCGGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCATCAAcc              >  1:180782/1‑150 (MQ=255)
                                                                                                                                            cTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCa                                                           >  2:41788/1‑98 (MQ=255)
                                                                                                                                            cTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCa                                                           <  1:41788/98‑1 (MQ=255)
                                                                                                                                            cTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCATCAACCAGACCtt       >  1:720587/1‑150 (MQ=255)
                                                                                                                                             tGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCATCAACCAGACCTTg      <  1:218921/150‑1 (MQ=255)
                                                                                                                                             tGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCATCAACCAGACCTTg      <  1:348578/150‑1 (MQ=255)
                                                                                                                                             tGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCATCAACCAGACCTTg      <  1:620973/150‑1 (MQ=255)
                                                                                                                                              gCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCAAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCATCAACCAGACCTTgg     <  2:625976/150‑1 (MQ=255)
                                                                                                                                                 gccgcc‑GCGCGCGTCCCACCCACACCCCCGCCACGCCGGCCACGTTCGGCAGAGTGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCATCAACCAGACCTTGGGGt  <  2:354984/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
CGAGGCGTTGACACGTTGAGCCGAGGCCGCTCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAA‑ACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCATCAACCAGACCTTGGGGT  >  NC_002947/2776905‑2777198

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