Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,159,126 A→G 6.1% intergenic (+30/‑179) PP_3661 → / → PP_3662 membrane protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,159,1260AG 6.1% 286.7 / 6.1 98intergenic (+30/‑179)PP_3661/PP_3662membrane protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (44/48);  new base G (4/2);  total (48/50)
Fisher's exact test for biased strand distribution p-value = 4.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.11e-01

GCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAA  >  NC_002947/4158978‑4159274
                                                                                                                                                    |                                                                                                                                                    
gCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTaa                                                                                                                                                     >  2:1619333/1‑150 (MQ=255)
  ttCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGTCTCTTTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGGAAAc                                                                                                                                                   >  1:1646319/1‑150 (MQ=255)
        cAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTc                                                                                                                                             <  1:690531/150‑1 (MQ=255)
            ttCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTAc                                                                                                                                         >  2:163952/1‑150 (MQ=255)
              cGCCAGGCTGGGCCACGGGCGCTGCCCCCTTTTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAg                                                                                                                                       <  1:479672/150‑1 (MQ=39)
                 cAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAggtg                                                                                                                                    >  1:1497010/1‑150 (MQ=255)
                  aGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGGTGGGTGAAGTCTGATTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAggtgg                                                                                                                                   <  1:331010/150‑1 (MQ=255)
                  aGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAggtgg                                                                                                                                   <  1:163952/150‑1 (MQ=255)
                  aGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAggtgg                                                                                                                                   <  2:1823751/150‑1 (MQ=255)
                  aGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAggtgg                                                                                                                                   <  2:30240/150‑1 (MQ=255)
                   ggCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAggtggg                                                                                                                                  <  1:1880763/150‑1 (MQ=255)
                    gCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGTGTGCAGAAGCTGTAGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAggtgggt                                                                                                                                 <  1:1181300/150‑1 (MQ=255)
                         gCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGtgctgc                                                                                                                            >  1:1396478/1‑150 (MQ=255)
                         aggcatctaagtctttcctgtgcgggtaaaacatccttctttgacaccaccataccgagaaatgtggacggcaacgaggggggctgtgaaaaaccttgagaccgagctgcCTGCTTCGCGGGTGAACCCGCTCCTACAGGTGGGtgctgc                                                                                                                            <  1:1118167‑M1/40‑1 (MQ=255)
                           cACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGccc                                                                                                                          >  1:262375/1‑150 (MQ=255)
                           cACGGGCGCTGGCAACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGccc                                                                                                                          <  1:1169115/150‑1 (MQ=255)
                            catctaagtctttcctgtgcgggtaaaacatccttctttgacaccaccataccgagaaatgtggacggcaacgaggggggctgtgaaaaaccttgagaccgagctgcCTGCTTCGCGGGTGAACCCGCTCCTACAGGTGGGTGCTGCCCt                                                                                                                         <  1:857508‑M1/43‑1 (MQ=255)
                            aCGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCt                                                                                                                         <  2:799935/150‑1 (MQ=255)
                              tctaagtctttcctgtgcgggtaaaacatccttctttgacaccaccataccgagaaatgtggacggcaacgaggggggctgtgaaaaaccttgagaccgagtcggCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGa                                                                                                                       <  2:1119828‑M1/45‑1 (MQ=255)
                                taagtctttcctgtgcgggtaaaacatccttctttgacaccaccataccgagaaatgtggacggcaacgaggggggctgtgaaaaaccttgagaccgagctgcCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAg                                                                                                                     <  2:1714550‑M1/47‑1 (MQ=255)
                                gcgcTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAg                                                                                                                     >  2:1463096/1‑150 (MQ=255)
                                   gtctttcctgtgcgggtaaaacatccttctttgacaccaccataccgagaaatgtggacggcaacgaggggggctgtgaaaaaccttgagaccgagctgcCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACg                                                                                                                  <  1:1550871‑M1/50‑1 (MQ=255)
                                       ccACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCaa                                                                                                              <  2:1396478/150‑1 (MQ=255)
                                        cACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAg                                                                                                             >  1:1487180/1‑150 (MQ=255)
                                          ccGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGc                                                                                                           <  2:170770/150‑1 (MQ=255)
                                            ggTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGccc                                                                                                         <  1:1494085/150‑1 (MQ=255)
                                                 ggtaaaacatccttctttgacaccaccataccgagaaatgtggacggcaacgaggggggctgtgaaaaaccttgagaccgagctgcCTGCTTCGCGGGTGAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTa                                                                                                    >  2:857508‑M1/87‑150 (MQ=255)
                                                 ggtaaaacatccttctttgacaccaccataccgagaaatgtggacggcaacgaggggggctgtgaaaaaccttgagaccgagctgcCTGCTTCGCGGGTGAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTa                                                                                                    >  2:827922‑M1/87‑150 (MQ=255)
                                                 ggtaaaacatccttctttgacaccaccataccgagaaatgtggacggcaacgaggggggctgtgaaaaaccttgagaccgagctgcCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTa                                                                                                    >  2:1550871‑M1/87‑150 (MQ=255)
                                                  gtaaaacatccttctttgacaccaccataccgagaaatgtggacggcaacgaggggggctgtgaaaaaccttgagaccgagctgcCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAg                                                                                                   <  1:672782‑M1/65‑1 (MQ=255)
                                                  aTCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAg                                                                                                   <  2:1213083/150‑1 (MQ=255)
                                                   tCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAgg                                                                                                  <  1:1619333/150‑1 (MQ=255)
                                                    ccTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGa                                                                                                 <  1:829455/150‑1 (MQ=255)
                                                      tCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGc                                                                                               <  2:140483/150‑1 (MQ=255)
                                                      aacatccttctttgacaccaccataccgagaaatgtggacggcaacgaggggggctgtgaaaaaccttgagaccgagctgcCTGCTTCGCGGGTGAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGc                                                                                               >  2:1118167‑M1/82‑150 (MQ=255)
                                                        catccttctttgacaccaccataccgagaaatgtggacggcaacgaggggggctgtgaaaaaccttgagaccgagctgcCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCgg                                                                                             >  1:1714550‑M1/80‑150 (MQ=255)
                                                        aaCCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCgg                                                                                             <  1:688043/150‑1 (MQ=255)
                                                         atccttctttgacaccaccataccgagaaatgtggacggcaacgaggggggctgtgaaaaaccttgagcccgagctgcCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCggg                                                                                            >  2:1697550‑M1/79‑150 (MQ=255)
                                                         aCCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCggg                                                                                            <  2:1598297/150‑1 (MQ=255)
                                                         aCCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCggg                                                                                            <  2:1497010/150‑1 (MQ=255)
                                                           ccttctttgacaccaccataccgagaaatgtggacggcaacgaggggggctgtgaaaaacctggagaccgagtcggCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGtt                                                                                          >  1:770391‑M1/77‑150 (MQ=255)
                                                              tgGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCAc                                                                                       <  1:1463096/150‑1 (MQ=255)
                                                                gCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCAccc                                                                                     >  1:648002/1‑150 (MQ=255)
                                                                  tgacaccaccataccgagaaatgtggacggcaacgaggggggctgtgaaaaaccttgagaccgagctgcCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCcgc                                                                                   >  2:672782‑M1/70‑150 (MQ=255)
                                                                  tGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTACTGCCCTGAAGGCGGCAAGGCCCCTGTAGGAGCGGGTTCACCcgc                                                                                   >  1:8378/1‑150 (MQ=255)
                                                                     ttGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGaa                                                                                <  1:1907287/150‑1 (MQ=255)
                                                                       ggCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAgc                                                                              <  1:675276/150‑1 (MQ=255)
                                                                       ccaccataccgagaaatgtggacggcaacgaggggggctgtgaaaaaccttgagaccgagctgcCTGCTTCGCGGGTGAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAgc                                                                              >  1:334674‑M1/65‑150 (MQ=255)
                                                                           gggggCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAgcaggc                                                                          >  1:1905992/1‑150 (MQ=255)
                                                                                 tGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCCGCTCg                                                                    >  1:24443/1‑150 (MQ=255)
                                                                                  gcggATTACCCTGGGGGTGCAGAGGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGGTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGccc                                                                                                         >  2:949047/4‑112 (MQ=255)
                                                                                  gcggATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGccc                                                                                                         <  1:949047/109‑1 (MQ=255)
                                                                                   gCGATTACTCTGGGGGTGCAGACGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGt                                                                  <  2:1041677/150‑1 (MQ=255)
                                                                                    cGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGtc                                                                 <  2:648002/150‑1 (MQ=255)
                                                                                    aaatgtggacggcaacgaggggggctgtgaaaaaccttgagaccgagtcggCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGtc                                                                 >  1:1119828‑M1/52‑150 (MQ=255)
                                                                                     gATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGtct                                                                >  2:1275710/1‑150 (MQ=255)
                                                                                          ggacggcaacgaggggggctgtgaaaaacctggagaccgagtcggCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAgg                                                           >  2:362077‑M1/46‑150 (MQ=255)
                                                                                          cccTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAg                                                                                                             <  2:121559/100‑1 (MQ=255)
                                                                                          cccTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAg                                                                                                             >  1:121559/1‑100 (MQ=255)
                                                                                           ccTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACGGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCAccc                                                                                     >  2:1055645/1‑123 (MQ=255)
                                                                                           ccTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCggg                                                                                            <  2:396706/116‑1 (MQ=255)
                                                                                           ccTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCggg                                                                                            >  1:396706/1‑116 (MQ=255)
                                                                                           ccTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCAccc                                                                                     <  1:1055645/123‑1 (MQ=255)
                                                                                           ccTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGt                                                          <  1:1275710/150‑1 (MQ=255)
                                                                                           ccTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGt                                                          >  1:1541003/1‑150 (MQ=255)
                                                                                           ccTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGt                                                          <  2:1541003/150‑1 (MQ=255)
                                                                                           ccTGGGGGTGCAGAAGCTGTTGGACTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGt                                                          >  2:946767/1‑150 (MQ=255)
                                                                                                      aGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAGGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGcc                                               >  1:1633081/1‑150 (MQ=255)
                                                                                                            tGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTc                                         >  1:485292/1‑150 (MQ=255)
                                                                                                              ttGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTGACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGt                                       >  2:1130730/1‑150 (MQ=255)
                                                                                                               tGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGtt                                      <  2:1685230/150‑1 (MQ=255)
                                                                                                                 ggCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGc                                    >  1:934066/1‑150 (MQ=255)
                                                                                                                  gCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGcc                                   >  2:1460250/1‑150 (MQ=255)
                                                                                                                       cTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCca                              <  1:1667212/150‑1 (MQ=255)
                                                                                                                        tGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTa                                                                                                    <  2:1549665/79‑1 (MQ=255)
                                                                                                                        tGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTa                                                                                                    >  1:1549665/1‑79 (MQ=255)
                                                                                                                        tGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCcac                             >  1:227221/1‑150 (MQ=255)
                                                                                                                         ggAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCcaca                            <  2:1487180/150‑1 (MQ=255)
                                                                                                                          gagaCCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACAt                           >  2:1878878/1‑150 (MQ=255)
                                                                                                                            gaCCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACAttt                         <  2:262375/150‑1 (MQ=255)
                                                                                                                               cGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTtctc                      >  1:469791/1‑150 (MQ=255)
                                                                                                                                gAGTCGTCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCGGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACTGCCCCCCTCGTTGCCGTCCACATTTCTCg                     <  2:764560/150‑1 (MQ=255)
                                                                                                                                gAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGcc                                   >  1:1248961/1‑136 (MQ=255)
                                                                                                                                gAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGcc                                   <  2:1248961/136‑1 (MQ=255)
                                                                                                                                  gTCGGCTGGTTCGCGGGTAAGCCCGCTCCTGCAGGTGGGTGCTGCCCTGAAGGGGGCAAGGCCCCTGTAGGGGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGt                   <  2:228213/150‑1 (MQ=255)
                                                                                                                                      gctgctTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGGAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCg                     <  1:346426/144‑1 (MQ=255)
                                                                                                                                      gctgctTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGGAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCg                     >  2:346953/1‑144 (MQ=255)
                                                                                                                                      gctgctTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGGAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCg                     >  2:346426/1‑144 (MQ=255)
                                                                                                                                      gctgctTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGGAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCg                     <  1:346953/144‑1 (MQ=255)
                                                                                                                                      gctgctTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTAtgg               >  1:477018/1‑150 (MQ=255)
                                                                                                                                      gctgctTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTAtgg               <  2:1659216/150‑1 (MQ=255)
                                                                                                                                      gctgctTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTAtgg               >  1:373782/1‑150 (MQ=255)
                                                                                                                                        TGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTAtggtg             >  1:946767‑M1/1‑150 (MQ=255)
                                                                                                                                            TCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGACGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGgtgt         >  2:1633081‑M1/1‑150 (MQ=255)
                                                                                                                                                GGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGtcaca     <  1:214835‑M1/150‑3 (MQ=255)
                                                                                                                                                GGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCaaa     >  2:469791‑M1/1‑150 (MQ=255)
                                                                                                                                                  GTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGt                                                          >  2:351470‑M1/1‑95 (MQ=255)
                                                                                                                                                  GTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGt                                                          <  1:351470‑M1/95‑1 (MQ=255)
                                                                                                                                                   TAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTGACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAaagaa  >  1:1130730‑M1/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAA  >  NC_002947/4158978‑4159274

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