Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,155,522 C→G 5.9% P44P (CCC→CCG) ‡ 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,5220CG 5.9% 823.9 / 28.8 247P44P (CCC→CCG) ‡PP_3655cytosine/purine/uracil/thiamine/allantoin permease family protein
Reads supporting (aligned to +/- strand):  ref base C (107/124);  new base G (7/7);  total (115/132)
Fisher's exact test for biased strand distribution p-value = 7.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 262 total aligned reads displayed.

CTTGCCTTGGAGTCAGTCATGACCAGTGCAGCCAATTCGGCACCCCTGATAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCGGCCCAGGGCCCGCAAATGGGCATTCCGCAGATGATCCA  >  NC_002947/4155373‑4155671
                                                                                                                                                     |                                                                                                                                                     
cttgccttgGAGTCCGTCATGACCAGTGCAGCCAATTCGGCACCCCTGATAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGccc                                                                                                                                                       >  1:4067050/1‑150 (MQ=255)
 ttgccttgGAGTCAGTCATGACCAGTGCAGCCAATTCGGCACCCCTGATAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCa                                                                                                                                                      >  1:417031/1‑150 (MQ=255)
    ccttgGAGTCAGTCATGACCAGTGCAGCCAATTCGGCACCCCTGATAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCg                                                                                                                                                   >  2:4604897/1‑150 (MQ=255)
     cttgGAGTCAGTCATGACCAGTGCAGCCAATTCGGCACCCCTGATAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGt                                                                                                                                                  >  2:6519590/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)
              agtcaTGACCAGTGCAGCCAATTCGGCACCCCTGATAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGgcgc                                                                                                                                         >  1:6723843/1‑150 (MQ=255)
              agtcaTGACCAGTGCAGCCAATTCGGCACCCCTGATAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGGTGACGATGGGCAACGgcgc                                                                                                                                         <  1:6281844/150‑1 (MQ=255)
              agtcaTGACCAGTGCAGCCAATTCGGCACCCCTGATAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGGTGACGATGGGCAACGgcgc                                                                                                                                         <  1:991562/150‑1 (MQ=255)
                tcaTGACCAGTGCAGCCAATTCGGCACCCCTGATAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGgcgcgc                                                                                                                                       >  2:2716523/1‑150 (MQ=255)
                tcaTGACCAGTGCAGCCAATTCGGCACCCCTGATAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGgcgcgc                                                                                                                                       >  1:172546/1‑150 (MQ=255)
                tcaTGACCAGTGCAGCCAATTCGGCACCCCTGATAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGgcgcgc                                                                                                                                       >  2:4480311/1‑150 (MQ=255)
                 caTGACCAGTGCAGCCAATTCGGCACCCCTGATAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCt                                                                                                                                      <  1:4480311/150‑1 (MQ=255)
                   tGACCAGTGCAGCCAATTCGGCACCCCTGATAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTgg                                                                                                                                    <  2:4929954/150‑1 (MQ=255)
                      ccAGTGCAGCCAATTCGGCACCCCTGATAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGgtg                                                                                                                                 >  2:373479/1‑150 (MQ=255)
                       cAGTGCAGCCAATTCGGCACCCCTGATAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACAGCGCGCCGGTGACGATGGGCAACGGCGCGATgttgcc                                                                                                                                <  2:1714188/150‑6 (MQ=38)
                          tGCAGCCAATTCGGCACCCCTGATAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCa                                                                                                                             <  2:6314895/150‑1 (MQ=255)
                           gCAGCCAATTCGGCACCCCTGATAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAg                                                                                                                            <  1:5772682/150‑1 (MQ=255)
                             aGCCAATTCGGCACCCCTGATAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCt                                                                                                                          <  1:2369171/150‑1 (MQ=255)
                              gCCAATTCGGCACCCCTGATAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGATCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTg                                                                                                                         >  2:3265977/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:2705430/150‑1 (MQ=255)
                                    tCGGCACCCCTGATAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCAc                                                                                                                   <  1:2261829/150‑1 (MQ=255)
                                    tCGGCACCCCTGATAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCAc                                                                                                                   <  1:2092427/150‑1 (MQ=255)
                                    tCGGCACCCCTGATAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCAc                                                                                                                   <  2:1684229/150‑1 (MQ=255)
                                      ggCACCCCTGATAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCt                                                                                                                 <  2:666767/150‑1 (MQ=255)
                                           cccTGATAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCGGAAcc                                                                                                            >  1:4568643/1‑150 (MQ=255)
                                            ccTGATAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCt                                                                                                           <  1:2564920/150‑1 (MQ=255)
                                             cTGATAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTg                                                                                                          <  2:1317652/150‑1 (MQ=255)
                                             cTGATAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTg                                                                                                          <  1:3047164/150‑1 (MQ=255)
                                               gATAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGt                                                                                                        >  1:4825399/1‑150 (MQ=255)
                                                   gAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATgg                                                                                                    >  1:2225413/1‑150 (MQ=255)
                                                    aaaaaCACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATggg                                                                                                   <  1:2612169/150‑1 (MQ=255)
                                                    aaaaaCACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATggg                                                                                                   <  1:1983717/150‑1 (MQ=255)
                                                     aaaaCACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATgggg                                                                                                  <  2:4098137/150‑1 (MQ=255)
                                                       aaCACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCa                                                                                                <  2:1157568/150‑1 (MQ=255)
                                                            acGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCAtcgtc                                                                                           <  1:1228848/150‑1 (MQ=255)
                                                              gATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGc                                                                                         <  1:3919077/150‑1 (MQ=255)
                                                               aTCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGcc                                                                                        <  1:3169478/150‑1 (MQ=255)
                                                                  ggCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATc                                                                                     <  2:1246250/150‑1 (MQ=255)
                                                                   gCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATcc                                                                                    >  2:5227298/1‑150 (MQ=255)
                                                                   gCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATcc                                                                                    >  2:5986581/1‑150 (MQ=255)
                                                                   gCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATcc                                                                                    >  1:6673589/1‑150 (MQ=255)
                                                                         tGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGGTCGGCGGCAACCTCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCg                                                                              <  2:6480982/150‑1 (MQ=255)
                                                                          gCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCgg                                                                             >  1:6185666/1‑150 (MQ=255)
                                                                             cccGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCa                                                                          <  1:1024482/150‑1 (MQ=255)
                                                                             cccGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCa                                                                          <  2:4836186/150‑1 (MQ=255)
                                                                             cccGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCa                                                                          >  2:1419504/1‑150 (MQ=255)
                                                                             cccGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCa                                                                          <  1:4358508/150‑1 (MQ=255)
                                                                                        gCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTcggcg                                                               <  2:3981617/150‑1 (MQ=255)
                                                                                        gCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCGGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTcggcg                                                               >  2:2139767/1‑150 (MQ=255)
                                                                                          cATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCggg                                                             <  1:1714453/150‑1 (MQ=255)
                                                                                            tGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCCTCGTCGCCAt                                                                                      <  2:739127/123‑1 (MQ=255)
                                                                                                 aGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGa                                                      <  2:2057922/150‑1 (MQ=255)
                                                                                                 aGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGa                                                      <  2:3839303/150‑1 (MQ=255)
                                                                                                 aGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGa                                                      >  2:6132852/1‑150 (MQ=255)
                                                                                                  ggTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGAt                                                     <  1:165023/150‑1 (MQ=255)
                                                                                                   gTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATg                                                    <  1:573077/150‑1 (MQ=255)
                                                                                                   gTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATg                                                    <  1:1337105/150‑1 (MQ=255)
                                                                                                      aGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGcc                                                 <  1:4109385/150‑1 (MQ=255)
                                                                                                      aGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGcc                                                 <  1:4706595/150‑1 (MQ=255)
                                                                                                      aGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGcc                                                 <  1:2877922/150‑1 (MQ=255)
                                                                                                        ggATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCt                                               >  2:1549958/1‑150 (MQ=255)
                                                                                                             tGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACt                                          >  2:6866535/1‑150 (MQ=255)
                                                                                                              gTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTc                                         <  2:5552598/150‑1 (MQ=255)
                                                                                                                tCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCgg                                       <  2:6631331/150‑1 (MQ=255)
                                                                                                                 cacaCTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCGGc                                      >  1:6499032/1‑150 (MQ=255)
                                                                                                                   cacTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCGGccc                                    >  1:1136781/1‑150 (MQ=255)
                                                                                                                    acTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCGGCCCa                                   >  1:1467903/1‑150 (MQ=255)
                                                                                                                    acTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCGGCCCa                                   >  1:5868586/1‑150 (MQ=255)
                                                                                                                    acTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCGGCCCa                                   >  1:3443619/1‑150 (MQ=255)
                                                                                                                     ctctGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCGGCCCAg                                  >  2:482778/1‑150 (MQ=255)
                                                                                                                      tctGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTTCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGTGTGCTTATGGCCCTGCACTCGGACCAgg                                 <  1:5062008/150‑1 (MQ=255)
                                                                                                                      tctGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCGGCCCAgg                                 <  2:5785166/150‑1 (MQ=255)
                                                                                                                        tGGTTCGGCGGCAACATCGCGCCGTTGCCTATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCGGCCCAGGGc                               <  2:23015/150‑1 (MQ=255)
                                                                                                                          gTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCGGCCCAGGGccc                             <  1:1754588/150‑1 (MQ=255)
                                                                                                                                  gcAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCATCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCGGCCCAGGGCCCGCAAATgg                     >  1:3898718/1‑150 (MQ=255)
                                                                                                                                   acccagcgcgCCGGTGACGATGGGCAACGGCGCGATGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCGGCCCAGGGCCCGCAAATggg                    >  2:2327906/7‑150 (MQ=38)
                                                                                                                                      cATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCGGCCCAGGGCCCGCAAATGGGCAt                 >  2:4959009/1‑150 (MQ=255)
                                                                                                                                       agcgcgCCGGTGACGATGGGCAACGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCGGCCCAGGGCCCGCAAATGGGCAtt                >  2:1854352/3‑150 (MQ=255)
                                                                                                                                       aTCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCGGCCCAGGGCCCGCAAATGGGCAtt                >  2:3140740/1‑150 (MQ=255)
                                                                                                                                       aTCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCGGCCCAGGGCCCGCAAATGGGCAtt                >  2:2442321/1‑150 (MQ=255)
                                                                                                                                        tCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATcc                                                                                    >  1:2387732/1‑81 (MQ=255)
                                                                                                                                        tCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATcc                                                                                    <  2:2387732/81‑1 (MQ=255)
                                                                                                                                          gcgcCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGccc                                                >  2:538988/1‑115 (MQ=255)
                                                                                                                                          gcgcCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGccc                                                <  1:538988/115‑1 (MQ=255)
                                                                                                                                          gcgcCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCGGCCCAGGGCCCGCAAATGGGCATTCCg             >  1:4896556/1‑150 (MQ=255)
                                                                                                                                          gcgcCGGTGACGATGGGCAACGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCATTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCGGCCCAGGGCCCGCAAATGGGCATTCCg             >  2:2513873/1‑150 (MQ=255)
                                                                                                                                            gcCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCGGCCCAGGGCCCGCAAATGGGCATTCCGCa           <  1:4304804/150‑1 (MQ=255)
                                                                                                                                            gcCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCGGCCCAGGGCCCGCAAATGGGCATTCCGCa           <  2:417031/150‑1 (MQ=255)
                                                                                                                                              cGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCGGCCCAGGGCCCGCAAATGGGCATTCCGCAga         >  2:6269664/1‑150 (MQ=255)
                                                                                                                                                  gCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCGGCCCAGGGCCCGCAAATGGGCATTCCGCAgatgat     >  1:7444/1‑150 (MQ=255)
                                                                                                                                                   cccATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCGGCCCAGGGCCCGCAAATGGGCATTCCGCAGATGATc    >  2:6165865/1‑150 (MQ=255)
                                                                                                                                                   cccATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCGGCCCAGGGCCCGCAAATGGGCATTCCGCAGATGATc    >  1:799785/1‑150 (MQ=255)
                                                                                                                                                   cccATCGTCACCGGCGCGCTGGGTGGGCAGCTGTTCCACCTGAACCTGGTCTGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCGGCCCCGGGCCCGCAAATGGGCCTTCCGCAGATGATc    >  2:4268/1‑150 (MQ=255)
                                                                                                                                                    ccATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCGGCCCAGGGCCCGCAAATGGGCATTCCGCAGATGATcc   >  1:5011206/1‑150 (MQ=255)
                                                                                                                                                     cATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCGGCCCAGGGCCCGCAAATGGGCATTCCGCAGATGATCCa  >  2:853953/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTTGCCTTGGAGTCAGTCATGACCAGTGCAGCCAATTCGGCACCCCTGATAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGATCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCGGCCCAGGGCCCGCAAATGGGCATTCCGCAGATGATCCA  >  NC_002947/4155373‑4155671

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: