Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,078,925 G→C 9.8% intergenic (‑122/+275) mgtE ← / ← PP_t56 magnesium transporter/tRNA‑Arg

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,078,9250GC 9.8% 215.3 / 11.8 71intergenic (‑122/+275)mgtE/PP_t56magnesium transporter/tRNA‑Arg
Reads supporting (aligned to +/- strand):  ref base G (32/32);  new base C (2/5);  total (34/37)
Fisher's exact test for biased strand distribution p-value = 4.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.85e-01

CTTTTTTAACTTCTACTTCGGTCATAGGCGAACTCCACTCCCAGCAGCGGAGCACGCCGAGGGAATCAATCGGTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAATTTTTGGCCAGAACAATCGTTTGCGCGATAAAGCTGGGACAACGCGTGAAGCTATCAGTGCGACGATTGTTGTGGATGGGAAAGCACATTTG  >  NC_002947/5078779‑5079068
                                                                                                                                                  |                                                                                                                                               
cTTTTTTAACTTCTACTTCGGTCATAGGCGAACTCCACTCCCAGCAGCGGAGCACGCCGAGGGAATCAATCGGTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCcggcgg                                                                                                                                              <  1:625879/150‑1 (MQ=255)
 ttttttAACTTCTACTTCGGTCATAGGCGAACTCCACTCCCAGCAGCGGAGCACGCCGAGGGAATCAATCGGTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCggg                                                                                                                                             <  2:118700/150‑1 (MQ=255)
 ttttttAACTTCTACTTCGGTCATAGGCGAACTCCACTCCCAGCAGCGGAGCACGCCGAGGGAATCAATCGGTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCggg                                                                                                                                             <  2:2877/150‑1 (MQ=255)
  tttttAACTTCTACTTCGGTCATAGGCGAACTCCACTCCCAGCAGCGGAGCACGCCGAGGGAATCAATCGGTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCgggg                                                                                                                                            >  2:758498/1‑150 (MQ=255)
  tttttAACTTCTACTTCGGTCATAGGCGAACTCCACTCCCAGCAGCGGAGCACGCCGAGGGAATCAATCGGTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCgggg                                                                                                                                            <  1:717311/150‑1 (MQ=255)
  tttttAACTTCTACTTCGGTCATAGGCGAACTCCACTCCCAGCAGCGGAGCACGCCGAGGGAATCAATCGGTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCgggg                                                                                                                                            >  2:758461/1‑150 (MQ=255)
       aaCTTCTACTTCGGTCATAGGCGAACTCCACTCCCAGCAGCGGAGCACGCCGAGGGAATCAATCGGTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGCCGGGGCttgt                                                                                                                                       <  2:71467/150‑4 (MQ=255)
       aaCTTCTACTTCGGTCATAGGCGAACTCCACTCCCAGCAGCGGAGCACGCCGAGGGAATCAATCGGTCAGATCTTTCTGTATGAATCGCGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGCCGGGGCttgt                                                                                                                                       <  1:378679/150‑4 (MQ=255)
        aCTTCTACTTCGGTCATAGGCGAACTCCACTCCCAGCAGCGGAGCCCGCCGAGGGAATCAATCGGTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGCCGGGGCttgtg                                                                                                                                      <  1:155622/150‑5 (MQ=255)
        aCTTCTACTTCGGTCATAGGCGAACTCCACTCCCAGCAGCGGAGCACGCCGAGGGAATCAATCGGTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACg                                                                                                                                      >  2:1044190/1‑150 (MQ=255)
         cTTCTACTTCGGTCATAGGCGAACTCCACTCCCAGCAGCGGAGCACGCCGAGGGAATCAATCGGTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGCCGGGGCttgtgg                                                                                                                                     <  1:362796/150‑6 (MQ=255)
          ttCTACTTCGGTCATAGGCGAACTCCACTCCCAGCAGCGGAGCACGCCGAGGGAATCAATCGGTCAGATCTTTCTGTCTGATCGGCATTCCCGTCCAACCACAGGGTAAGGCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCGGACggg                                                                                                                                    >  1:878495/1‑150 (MQ=255)
          ttCTACTTCGGTCATAGGCGAACTCCACTCCCAGCAGCGGAGCACGCCGAGGGAATCAATCGGTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACggg                                                                                                                                    >  2:704239/1‑150 (MQ=255)
           tCTACTTCGGTCATAGGCGAACTCCACTCCCAGCAGCGGAGCACGCCGAGGGAATCAATCGGTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTgccaggc                                                                                                                                   <  2:672288/150‑7 (MQ=255)
           tCTACTTCGGTCATAGGCGAACTCCACTCCCAGCAGCGGAGCACGCCGAGGGAATCAATCGGTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGgc                                                                                                                                   <  1:1305512/150‑1 (MQ=255)
               cTTCGGTCATAGGCGAACTCCACTCCCAGCAGCGGAGCACGCCGAGGGAATCAATCGGTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCaa                                                                                                                               <  1:758461/150‑1 (MQ=255)
               cTTCGGTCATAGGCGAACTCCACTCCCAGCAGCGGAGCACGCCGAGGGAATCAATCGGTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCaa                                                                                                                               <  1:758498/150‑1 (MQ=255)
                      cATAGGCGAACTCCACTCCCAGCAGCGGAGCACGCCGAGGGAATCAATCGGTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATaca                                                                                                                        <  1:704239/150‑1 (MQ=255)
                      cATAGGCGAACTCCACTCCCAGCAGCGGAGCACGCCGAGGGAATCAATCGGTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATaca                                                                                                                        <  1:394402/150‑1 (MQ=255)
                      cATAGGCGAACTCCACTCCCAGCAGCGGAGCACGCCGAGGGAATCAATCGGTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATaca                                                                                                                        <  1:1367701/150‑1 (MQ=255)
                       aTAGGCGAACTCCACTCCCAGCAGCGGAGCACGCCGAGGGAATCAATCGGTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATacac                                                                                                                       >  2:670494/1‑150 (MQ=255)
                        tAGGCGAACTCCACTCCCAGCAGCGGAGCACGCCGAGGGAATCAATCGGTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACt                                                                                                                      >  1:83699/1‑150 (MQ=255)
                           gCGAACTCCACTCCCAGCAGCGGAGCACGCCGAGGGAATCAATCGGTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATaca                                                                                                                        >  2:1407111/1‑145 (MQ=255)
                           gCGAACTCCACTCCCAGCAGCGGAGCACGCCGAGGGAATCAATCGGTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATaca                                                                                                                        <  1:1407111/145‑1 (MQ=255)
                           gCGAACTCCACTCCCAGCAGCGGAGCACGCCGAGGGAATCAATCGGTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCt                                                                                                                   >  1:1048082/1‑150 (MQ=255)
                                 tccactccCAGCAGCGGAGCACGCCGAGGGAATCAATCGGTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGCCGGGGCttgtg                                                                                                                                      >  2:626962/1‑121 (MQ=255)
                                 tccactccCAGCAGCGGAGCACGCCGAGGGAATCAATCGGTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGCCGGGGCttgtg                                                                                                                                      <  1:626962/125‑5 (MQ=255)
                                   cactccCAGCAGCGGAGCACGCCGAGGGCCAAAATCGGTCAGATGTGTCGGTAGGAATCGTGGGAGCGAGGAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGt                                                                                                           >  2:375018/1‑150 (MQ=255)
                                   cactccCAGCAGCGGAGCACGCCGAGGGAATCAATCGGTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGt                                                                                                           <  2:581255/150‑1 (MQ=255)
                                        ccAGCAGCGGAGCACGCCGAGGGAATCAATCGGTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGAt                                                                                                      >  2:1152145/1‑150 (MQ=255)
                                            cagcGGAGCACGCCGAGGGAATCAATCGGTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCt                                                                                                  <  2:781371/150‑1 (MQ=255)
                                                 gAGCACGCCGAGGGAATCAATCGGTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAAt                                                                                                                              <  2:1282451/117‑1 (MQ=255)
                                                 gAGCACGCCGAGGGAATCAATCGGTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAAt                                                                                                                              >  1:1282451/1‑117 (MQ=255)
                                                     aCGCCGAGGGAATCAATCGGTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAAttttt                                                                                         <  1:670494/150‑1 (MQ=255)
                                                     aCGCCGAGGGAATCAATCGGTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAAttttt                                                                                         <  2:83699/150‑1 (MQ=255)
                                                      cGCCGAGGGAATCAATCGGTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAATTTTTg                                                                                        <  1:1152145/150‑1 (MQ=255)
                                                       gCCGAGGGAATCAATCGGTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAATTTTTgg                                                                                       >  1:338822/1‑150 (MQ=255)
                                                          gAGGGAATCAATCGGTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGTCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAATTTTTGGCCa                                                                                    >  2:529041/1‑150 (MQ=255)
                                                          gAGGGAATCAATCGGTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGCAGATCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAATTTTTGGCCa                                                                                    <  1:1106044/150‑1 (MQ=255)
                                                                 tcaatcGGTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAATTTTTGGCCATAACAAt                                                                             >  1:374247/1‑150 (MQ=255)
                                                                      cGGTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAATTTTTGGCCAGAACAATCGttt                                                                        <  1:1044190/150‑1 (MQ=255)
                                                                        gTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCggg                                                                                                                                             <  1:946682/79‑1 (MQ=255)
                                                                        gTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCggg                                                                                                                                             >  2:946682/1‑79 (MQ=255)
                                                                        gTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGTCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAATTTTTGGCCAGAACAATCGTTTgc                                                                      <  1:529041/150‑1 (MQ=255)
                                                                           aGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGAt                                                                                                      <  1:293885/115‑1 (MQ=255)
                                                                           aGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGAt                                                                                                      >  2:293885/1‑115 (MQ=255)
                                                                                    tGTATGAATCGTGTTAGCGAGTAACTACTGGGTCAGTCCCTGGTGGAATCCCACAACGCCCGGCGGGGGGGAGGGCCGGAAGCATACACTGGTGAAGCTGTCAGTTCGCTTAAATTTTTGGCCAGCACAACCGTTTACGCGATAAAGCTg                                                          <  2:518674/150‑1 (MQ=255)
                                                                                    tGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAATTTTTGGCCAGAACAATCGTTTGCGCGATAAAGCTg                                                          <  2:878495/150‑1 (MQ=255)
                                                                                           agctttttAGCGCGTAACTCATGGATCAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAATTTTTGGCCAGAACAATCGTTTGCGCGATAAAGCTGGGACAAc                                                   <  2:374247/144‑1 (MQ=255)
                                                                                                   aGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAATTTTTGGCCAGAACAATCGTTTGCGCGATAAAGCTGGGACAACGCGTGAAg                                           >  2:1166348/1‑150 (MQ=255)
                                                                                                           actactGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAATTTTTGGCCAGAACAATCGTTTGCGCGATAAAGCTGGGACAACGCGTGAAGCTATCAGt                                   <  1:375018/150‑1 (MQ=255)
                                                                                                                  ggTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAATTTTTGGCCAGAACAATCGTTTGCGCGATAAAGCTGGGACAACGCGTGaa                                            >  1:1070400/1‑134 (MQ=255)
                                                                                                                  ggTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAATTTTTGGCCAGAACAATCGTTTGCGCGATAAAGCTGGGACAACGCGTGaa                                            <  2:1070400/134‑1 (MQ=255)
                                                                                                                       gTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAATTTTTGGCCAGAACAATCGTTTGCGCGATAAAGCTGGGACAACGCGTGAAGCTATCAGTGCGACGAttgtt                       >  2:344693/1‑150 (MQ=255)
                                                                                                                        tCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAAtttt                                                                                          <  2:525665/82‑1 (MQ=255)
                                                                                                                        tCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAAtttt                                                                                          >  1:525665/1‑82 (MQ=255)
                                                                                                                        tCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAATTTTTGGCCAGAACAATCGTTTGCGCGATAAAGCTGGGACAACGCGTGAAGCTATCAGTGCGACGAttgttg                      >  2:117258/1‑150 (MQ=255)
                                                                                                                              gtggtTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAATTTTTGGCCAGAACAATCGTTTGCGCGATAAAGCTGGGACAACGCGTGAAg                                           <  1:464070/123‑1 (MQ=255)
                                                                                                                              gtggtTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAATTTTTGGCCAGAACAATCGTTTGCGCGATAAAGCTGGGACAACGCGTGAAg                                           >  2:464070/1‑123 (MQ=255)
                                                                                                                               tggtTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAATTTTTg                                                                                        >  2:996445/1‑77 (MQ=255)
                                                                                                                               tggtTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAATTTTTg                                                                                        <  1:996445/77‑1 (MQ=255)
                                                                                                                                      ccgtcaGCCCCGCCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAATTTTTGGCCAGAACAATCGTTTGCGCGATAAAGCTGGGACAACGCGTGAAGCTATCAGTGCGACGATTGTTGTGGATGGGAAAGca        >  2:313932/6‑150 (MQ=255)
                                                                                                                                      ccACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAATTTTTGGCCAGAACAATCGTTTGCGCGATAAAGCTGGGACAACGCGTGAAGCTATCAGTGCGACGATTGTTGTGGATGGGAAAGca        >  2:376647/1‑150 (MQ=255)
                                                                                                                                        acaAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAATTTTTGGCCAGAACAATCGTTTGCGCGATAAAGCTGGGACAACGCGTGAAGCTATCAGTGCGACGATTGTTGTGGATGGGAAAGcaca      >  2:594581/1‑150 (MQ=255)
                                                                                                                                        acaAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAATTTTTGGCCAGAACAATCGTTTGCGCGATAAAGCTGGGACAACGCGTGAAGCTATCAGTGCGACGATTGTTGTGGATGGGAAAGcaca      >  1:588602/1‑150 (MQ=255)
                                                                                                                                         caAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAATTTTTGGCCAGAACAATCGTTTGCGCGATAAAGCTGGGACAACGCGTGAAGCTATCAGTGCGACGATTGTTGTGGATGGGAAAGCACAt     >  1:1212257/1‑150 (MQ=255)
                                                                                                                                          aaGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAATTTTTGGCCAGAACAATCGTTTGCGCGATAAAGCTGGGACAACGCGTGAAGCTATCAGTGCGACGATTGTTGTGGATGGGAAAGCACAtt    >  1:643285/1‑150 (MQ=255)
                                                                                                                                            gCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAATTTTTGGCCAGAACAATCGTTTGCGCGATAAAGCTGGGACAACGCGTGAAGCTATCAGTGCGACGATTGTTGTGGATGGGAAAGCACATTTg  >  2:1234485/1‑150 (MQ=255)
                                                                                                                                             ccccGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAATTTTTGGCCAGAACAATCGTTTGCGCGATAAAGCTGGGACAACGCGTGAAGCTATCAGTGCGACGAttgtt                       <  1:1214307/128‑1 (MQ=255)
                                                                                                                                             ccccGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAATTTTTGGCCAGAACAATCGTTTGCGCGATAAAGCTGGGACAACGCGTGAAGCTATCAGTGCGACGAttgtt                       >  2:1214307/1‑128 (MQ=255)
                                                                                                                                              cccGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAATTTTTGGCCAGAACAATCGTTTGCGCGATAAAGCTGGGACAACGCGTGAAGCTATCAGTGCGACGATTGTTGTGGATGGGaaa           <  2:873819/139‑1 (MQ=255)
                                                                                                                                              cccGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAATTTTTGGCCAGAACAATCGTTTGCGCGATAAAGCTGGGACAACGCGTGAAGCTATCAGTGCGACGATTGTTGTGGATGGGaaa           >  1:873819/1‑139 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
CTTTTTTAACTTCTACTTCGGTCATAGGCGAACTCCACTCCCAGCAGCGGAGCACGCCGAGGGAATCAATCGGTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAATTTTTGGCCAGAACAATCGTTTGCGCGATAAAGCTGGGACAACGCGTGAAGCTATCAGTGCGACGATTGTTGTGGATGGGAAAGCACATTTG  >  NC_002947/5078779‑5079068

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