Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,155,525 C→G 5.3% I45M (ATC→ATG PP_3655 → cytosine/purine/uracil/thiamine/allantoin permease family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,155,5250CG 5.3% 899.4 / 28.7 265I45M (ATC→ATGPP_3655cytosine/purine/uracil/thiamine/allantoin permease family protein
Reads supporting (aligned to +/- strand):  ref base C (122/129);  new base G (7/7);  total (129/136)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

Only 100 of 266 total aligned reads displayed.

GCCTTGGAGTCAGTCATGACCAGTGCAGCCAATTCGGCACCCCTGATAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCGGCCCAGGGCCCGCAAATGGGCATTCCGCAGATGATCCAGA  >  NC_002947/4155376‑4155673
                                                                                                                                                     |                                                                                                                                                    
 ccttgGAGTCAGTCATGACCAGTGCAGCCAATTCGGCACCCCTGATAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCg                                                                                                                                                     >  2:4604897/1‑150 (MQ=255)
       agtcagtcATGACCAGTGCAGCCAATTCGGCACCCCTGATAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCg                                                                                                                                               <  2:4984469/150‑1 (MQ=255)
         tcagtcaTGACCAGTGCAGCCAATTCGGCACCCCTGATAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGgc                                                                                                                                             <  2:2876190/150‑1 (MQ=255)
           agtcaTGTCCTGTGCAGCCATGTCGGGAACCCTTTTTGAAAAACACACGATCGGCTTCGTGCCCCAGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACAGTCTGGTTCGGCGGCATCATCGCGCCGGTGACGATGGGCAACGgcgc                                                                                                                                           <  1:6379431/150‑1 (MQ=255)
           agtcaTGACCAGTGCAGCCAATTCGGCACCCCTGATAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGGTGACGATGGGCAACGgcgc                                                                                                                                           <  1:6281844/150‑1 (MQ=255)
             tcaTGACCAGTGCAGCCAATTCGGCACCCCTGATAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGgcgcgc                                                                                                                                         >  2:2716523/1‑150 (MQ=255)
                tGACCAGTGCAGCCAATTCGGCACCCCTGATAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTgg                                                                                                                                      <  1:1778417/150‑1 (MQ=255)
                tGACCAGTGCAGCCAATTCGGCACCCCTGATAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTgg                                                                                                                                      <  2:4929954/150‑1 (MQ=255)
                 gACCAGTGCAGCCAATTCGGCACCCCTGATAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTggg                                                                                                                                     >  1:1656853/1‑150 (MQ=255)
                   ccAGTGCAGCCAATTCGGCACCCCTGATAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGgtg                                                                                                                                   >  2:373479/1‑150 (MQ=255)
                    cAGTGCAGCCAATTCGGCACCCCTGATAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGgtgt                                                                                                                                  >  1:903287/1‑150 (MQ=255)
                      gTGCAGCCAATTCGGCACCCCTGATAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGc                                                                                                                                <  1:1429254/150‑1 (MQ=255)
                       tGCAGCCAATTCGGCACCCCTGATAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCa                                                                                                                               <  2:6314895/150‑1 (MQ=255)
                        gCAGCCAATTCGGCACCCCTGATAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAg                                                                                                                              <  1:5772682/150‑1 (MQ=255)
                           gCCAATTCGGCACCCCTGATAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGATCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTg                                                                                                                           >  2:3265977/1‑150 (MQ=255)
                            ccAATTCGGCACCCCTGATAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGt                                                                                                                          >  2:4366254/1‑150 (MQ=255)
                            ccAATTCGGCACCCCTGATAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGt                                                                                                                          <  1:6230610/150‑1 (MQ=255)
                            ccAATTCGGCACCCCTGATAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGGTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGt                                                                                                                          >  2:6020563/1‑150 (MQ=255)
                                ttCGGCACCCCTGATAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCa                                                                                                                      <  1:855927/150‑1 (MQ=255)
                                 tCGGCACCCCTGATAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCAc                                                                                                                     <  2:6269733/150‑1 (MQ=255)
                                 tCGGCACCCCTGATAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCAc                                                                                                                     <  2:2705430/150‑1 (MQ=255)
                                 tCGGCACCCCTGATAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCAc                                                                                                                     <  2:1684229/150‑1 (MQ=255)
                                 tCGGCACCCCTGATAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCAc                                                                                                                     <  1:2261709/150‑1 (MQ=255)
                                 tCGGCACACATGATAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGACATGGAAATGTAAGGGATCTGTTCACTCTCTTGTTCGGCGGCAACATCGAGCCGTTGCCCATCGTCACCGGAGCGCTGGGTGTGCAGCTGTTCCAc                                                                                                                     <  1:3045358/150‑1 (MQ=255)
                                  cGGCACCCCTGATAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCAcc                                                                                                                    <  1:2520197/150‑1 (MQ=255)
                                   ggCACCCCTGATAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCt                                                                                                                   <  2:4085315/150‑1 (MQ=255)
                                        cccTGATAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCGGAAcc                                                                                                              >  1:4568643/1‑150 (MQ=255)
                                          cTGATAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTg                                                                                                            <  2:1317652/150‑1 (MQ=255)
                                          cTGAGAGAAACACACACGAGCGGCGACGTGCCCCCGCAAGACCGCCATGGAACGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTg                                                                                                            <  2:4941979/150‑1 (MQ=255)
                                                gAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATgg                                                                                                      >  1:2225413/1‑150 (MQ=255)
                                                gAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATgg                                                                                                      <  2:2077534/150‑1 (MQ=255)
                                                 aaaaaCACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATggg                                                                                                     <  1:2612169/150‑1 (MQ=255)
                                                 aaaaaCACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATggg                                                                                                     >  2:6648251/1‑150 (MQ=255)
                                                         acGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCAtcgtc                                                                                             <  1:1228848/150‑1 (MQ=255)
                                                           gATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGc                                                                                           >  2:5364757/1‑150 (MQ=255)
                                                           gATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGc                                                                                           <  1:3919077/150‑1 (MQ=255)
                                                            aTCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGcc                                                                                          <  1:3169478/150‑1 (MQ=255)
                                                             tCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCa                                                                                         >  2:174813/1‑150 (MQ=255)
                                                                gCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATcc                                                                                      >  2:5986581/1‑150 (MQ=255)
                                                                gCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATcc                                                                                      >  1:6673589/1‑150 (MQ=255)
                                                                gCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATcc                                                                                      >  2:5227298/1‑150 (MQ=255)
                                                                     gTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTc                                                                                 >  1:2308322/1‑150 (MQ=255)
                                                                       gCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCgg                                                                               >  1:6185666/1‑150 (MQ=255)
                                                                          cccGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCa                                                                            <  1:1024482/150‑1 (MQ=255)
                                                                          cccGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCa                                                                            <  1:4358508/150‑1 (MQ=255)
                                                                          cccGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCa                                                                            <  2:4836186/150‑1 (MQ=255)
                                                                                     gCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTcggcg                                                                 <  2:3981617/150‑1 (MQ=255)
                                                                                      ccATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTcggcgg                                                                >  1:3732477/1‑150 (MQ=255)
                                                                                       cATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCggg                                                               <  1:1714899/150‑1 (MQ=255)
                                                                                        aTGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCgggg                                                              <  2:2352579/150‑1 (MQ=255)
                                                                                         tGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCAt                                                                                        >  1:739127/1‑123 (MQ=255)
                                                                                          ggAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTg                                                            <  2:675575/150‑1 (MQ=255)
                                                                                              aGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGa                                                        >  2:6132852/1‑150 (MQ=255)
                                                                                              aGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGa                                                        <  2:3839303/150‑1 (MQ=255)
                                                                                                gTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATg                                                      <  1:573077/150‑1 (MQ=255)
                                                                                                gTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATg                                                      <  1:1337105/150‑1 (MQ=255)
                                                                                                gTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATg                                                      <  2:1459251/150‑1 (MQ=255)
                                                                                                   aGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGcc                                                   <  1:2877922/150‑1 (MQ=255)
                                                                                                   aGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGcc                                                   <  1:4706595/150‑1 (MQ=255)
                                                                                                     ggATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCt                                                 <  2:648252/150‑1 (MQ=255)
                                                                                                     ggATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCt                                                 >  2:1549958/1‑150 (MQ=255)
                                                                                                        tCTGTTCACACTCTGGTTCGGCGGCAACTTCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCa                                              <  1:6824627/150‑1 (MQ=255)
                                                                                                          tGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACt                                            >  2:6866678/1‑150 (MQ=255)
                                                                                                           gTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTc                                           <  2:5552598/150‑1 (MQ=255)
                                                                                                            ttCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCg                                          >  2:6566461/1‑150 (MQ=255)
                                                                                                              cacaCTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCGGc                                        >  1:6499032/1‑150 (MQ=255)
                                                                                                              cacaCTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCGGc                                        >  1:2072471/1‑150 (MQ=255)
                                                                                                                cgttctGGGTCGGCGGCAAGATCGCTCCGTTGCCCATCTTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGTTCCTGGTTTGGGGCATCGTCGCCATCCTTGTCGGCCTCCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCGGccc                                      <  1:1648508/147‑1 (MQ=255)
                                                                                                                cacTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCGGccc                                      >  1:1136781/1‑150 (MQ=255)
                                                                                                                 acTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCCTCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGACGGCGGGGTGCTGATGGCCCTGCACTCGGCCCa                                     >  1:4105580/1‑150 (MQ=255)
                                                                                                                 acTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCGGCCCa                                     >  1:3443619/1‑150 (MQ=255)
                                                                                                                 acTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCGGCCCa                                     >  2:4267560/1‑150 (MQ=255)
                                                                                                                 acTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCGGCCCa                                     >  1:5868586/1‑150 (MQ=255)
                                                                                                                 acTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCGGCCCa                                     >  1:506067/1‑150 (MQ=255)
                                                                                                                   tctGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCGGCCCAgg                                   <  2:5785166/150‑1 (MQ=255)
                                                                                                                     tGGTTCGGCGGCAACATCGCGCCGTTGCCTATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCGGCCCAGGGc                                 <  2:23015/150‑1 (MQ=255)
                                                                                                                       gTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCGGCCCAGGGccc                               <  1:1754588/150‑1 (MQ=255)
                                                                                                                       gTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCGGCCCAGGGccc                               >  1:445516/1‑150 (MQ=255)
                                                                                                                        ttCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCGGCCCAGGGCCCg                              >  1:4262281/1‑150 (MQ=255)
                                                                                                                        ttCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCGGCCCAGGGCCCg                              >  2:852003/1‑150 (MQ=255)
                                                                                                                               gcAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCGGCCCCGGGCCCGCAAATgg                       >  1:764806/1‑150 (MQ=255)
                                                                                                                                cAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCGGCCCAGGGCCCGCAAATggg                      <  1:4155913/150‑1 (MQ=255)
                                                                                                                                  ccagcgcgCCGGTGACGATGGGCAACGGCGCGATGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCGGCCCAGGGCCCGCAAATGGGCa                    >  2:1115146/5‑150 (MQ=255)
                                                                                                                                    agcgcgCCGGTGACGATGGGCAACGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCGGCCCAGGGCCCGCAAATGGGCAtt                  >  2:1854352/3‑150 (MQ=255)
                                                                                                                                    aTCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCGGCCCAGGGCCCGCAAATGGGCAtt                  >  2:3140740/1‑150 (MQ=255)
                                                                                                                                     tCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATcc                                                                                      >  1:2387732/1‑81 (MQ=255)
                                                                                                                                     tCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATcc                                                                                      <  2:2387767/81‑1 (MQ=255)
                                                                                                                                       gcgcCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGccc                                                  <  1:538988/115‑1 (MQ=255)
                                                                                                                                       gcgcCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCGGCCCAGGGCCCGCAAATGGGCATTCCg               >  1:4896556/1‑150 (MQ=255)
                                                                                                                                       gcgcCGGTGACGATGGGCAACGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCATTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCGGCCCAGGGCCCGCAAATGGGCATTCCg               >  2:2513873/1‑150 (MQ=255)
                                                                                                                                         gcCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCGGCCCAGGGCCCGCAAATGGGCATTCCGCa             <  2:417031/150‑1 (MQ=255)
                                                                                                                                                cccATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCGGCCCAGGGCCCGCAAATGGGCATTCCGCAGATGATc      >  2:6165865/1‑150 (MQ=255)
                                                                                                                                                 ccATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCGGCCCAGGGCCCGCAAATGGGCATTCCGCAGATGATcc     >  2:1270196/1‑150 (MQ=255)
                                                                                                                                                 ccATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCGGCCCAGGGCCCGCAAATGGGCATTCCGCAGATGATcc     >  2:4787794/1‑150 (MQ=255)
                                                                                                                                                 ccATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCGGCCCAGGGCCCGCAAATGGGCATTCCGCAGATGATcc     >  2:2155304/1‑150 (MQ=255)
                                                                                                                                                 ccATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCGGCCCAGGGCCCGCAAATGGGCATTCCGCAGATGATcc     >  1:5011206/1‑150 (MQ=255)
                                                                                                                                                 ccATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCGGCCCAGGGCCCGCAAATGGGCATTCCGCAGATGATcc     >  1:4793594/1‑150 (MQ=255)
                                                                                                                                                   aTCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCGGCCCAGGGCCCGCAAATGGGCATTCCGCAGATGATCCag   >  2:2428192/1‑150 (MQ=255)
                                                                                                                                                   aTCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCGGCCCAGGGCCCGCAAATGGGCATTCCGCAGATGATCCag   <  2:903287/150‑1 (MQ=255)
                                                                                                                                                    tCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCGGCCCAGGGCCCGCAAATGGGCATTCCGCAGATGATCCaga  <  2:6215225/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCCTTGGAGTCAGTCATGACCAGTGCAGCCAATTCGGCACCCCTGATAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCGGCCCAGGGCCCGCAAATGGGCATTCCGCAGATGATCCAGA  >  NC_002947/4155376‑4155673

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