Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 503,835 A→T 9.0% intergenic (+35/‑116) PP_0414 → / → rpe polyamine ABC transporter permease/ribulose‑5‑phosphate 3‑epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947503,8350AT 9.0% 286.0 / 11.5 101intergenic (+35/‑116)PP_0414/rpepolyamine ABC transporter permease/ribulose‑5‑phosphate 3‑epimerase
Reads supporting (aligned to +/- strand):  ref base A (49/42);  new base T (4/5);  total (53/48)
Fisher's exact test for biased strand distribution p-value = 7.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

Only 100 of 101 total aligned reads displayed.

CGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTATTGCCCCCTCCATTC  >  NC_002947/503688‑503981
                                                                                                                                                   |                                                                                                                                                   
cGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAg                                                                                                                                                   <  2:1984573/150‑1 (MQ=255)
cGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGAATCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAg                                                                                                                                                   >  1:1984573/1‑150 (MQ=255)
 gACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGc                                                                                                                                                  <  1:689872/150‑1 (MQ=255)
        gCCGCGGCGGCGCTCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAaga                                                                                                                                          <  2:160731/150‑1 (MQ=255)
        gCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAaga                                                                                                                                          <  1:704301/150‑1 (MQ=255)
         ccGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAagag                                                                                                                                         >  2:653405/1‑150 (MQ=255)
         ccGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCCGGCTTTAGCCGCGAagag                                                                                                                                         <  1:1797174/150‑1 (MQ=255)
           gcggcggcgACCTTGCTGATTGCGTTCTCGGTGTGGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCTGCCTCTTCGCGGCTAAAGCCGGAAGAGGc                                                                                                                                       <  1:1494700/150‑1 (MQ=255)
            cggcggcgACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGc                                                                                                                                       <  2:1620401/149‑1 (MQ=255)
            cggcggcgACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGc                                                                                                                                       >  1:1620401/1‑149 (MQ=255)
             ggcggcgACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCg                                                                                                                                     <  1:653405/150‑1 (MQ=255)
             ggcggcgACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCGGTACCGGCCTCTTCCGGCTTTAGCCGCGAAGAGGCCg                                                                                                                                     <  2:1044025/150‑1 (MQ=255)
                       ttGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTGCCGGCCTCTTCCGGCTTTAGCCgcgaa                                                                                                                                            >  2:2216432/1‑129 (MQ=38)
                       ttGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTGCCGGCCTCTTCCGGCTTTAGCCgcgaa                                                                                                                                            <  1:2216432/133‑5 (MQ=38)
                        tGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTACCTGTGCCGGCCTCTTCCGGCTTTAGCCGCGAAGAGGCCGGTACAggcaag                                                                                                                          <  1:2106693/150‑5 (MQ=25)
                          cTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGt                                                                                                                        <  1:1201460/150‑1 (MQ=255)
                           tGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCTCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGtt                                                                                                                       >  1:981108/1‑150 (MQ=255)
                            gATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTg                                                                                                                      <  1:1886205/150‑1 (MQ=255)
                                 cGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGa                                                                                                                 <  2:131469/150‑1 (MQ=255)
                                  gTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGac                                                                                                                >  2:1335930/1‑150 (MQ=255)
                                    tctcGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGacac                                                                                                              <  2:1987445/150‑1 (MQ=255)
                                      tcGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGcc                                                                                                                                                 >  2:1466346/1‑114 (MQ=255)
                                      tcGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGcc                                                                                                                                                 <  1:1466346/114‑1 (MQ=255)
                                           gttgttGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGcc                                                                                                                                      <  1:615179/119‑1 (MQ=255)
                                           gttgttGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGcc                                                                                                                                      >  2:615179/1‑119 (MQ=255)
                                           gttgttGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGa                                                                                                       <  2:381032/150‑1 (MQ=255)
                                            ttgttgCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGGTGACTTGCCGTTTCCGGCCTCTTTGCTGCTAAGGCCGGTGGAGGCCTGGACGTGTTTCTTTTGGCTACACCAATTGAt                                                                                                      >  2:2343132/1‑150 (MQ=255)
                                            ttgttgCTGCTGACCCTGGAGGGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGAt                                                                                                      <  1:1160147/150‑1 (MQ=255)
                                             tgttgCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATa                                                                                                     <  2:2027046/150‑1 (MQ=255)
                                              gttgCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTc                                                                                                                   >  2:2340025/1‑135 (MQ=255)
                                              gttgCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTc                                                                                                                   <  1:2340025/135‑1 (MQ=255)
                                              gttgCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGAATGAGGCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATAc                                                                                                    <  1:2275402/150‑1 (MQ=255)
                                               ttgCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGcc                                                                                                                                                 <  2:1932375/105‑1 (MQ=255)
                                               ttgCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGcc                                                                                                                                                 >  1:1932375/1‑105 (MQ=255)
                                                tgctgcTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCa                                                                                                  >  1:1646524/1‑150 (MQ=255)
                                                  ctgctgACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGcc                                                                                                                                      <  1:1072112/112‑1 (MQ=255)
                                                  ctgctgACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGcc                                                                                                                                      >  2:1072112/1‑112 (MQ=255)
                                                  ctgctgACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAAt                                                                                                <  2:2206581/150‑1 (MQ=255)
                                                          ccTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTc                                                                                        >  1:2184676/1‑150 (MQ=255)
                                                          ccTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGGTGCGTTGGACGGCACGCATTGATACCAATAAGCCTTc                                                                                        >  2:2397582/1‑150 (MQ=255)
                                                            tGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCAt                                                                                      <  2:889158/150‑1 (MQ=255)
                                                               aGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGc                                                                                   <  1:1042503/150‑1 (MQ=255)
                                                                gTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAg                                                                                                                               <  2:1414390/105‑1 (MQ=255)
                                                                gTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAg                                                                                                                               >  1:1414390/1‑105 (MQ=255)
                                                                        gTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGAGGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGa                                                                          <  1:2397582/150‑1 (MQ=255)
                                                                         tGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTAAg                                                                         >  1:2280074/1‑150 (MQ=255)
                                                                            gCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTACCTGTGCCGGCCTCTTCCGGCTTTAGCCgcgaa                                                                                                                                            <  1:2242482/80‑5 (MQ=18)
                                                                               tAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACAcc                                                                                                             >  2:1936566/1‑108 (MQ=255)
                                                                               tAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACAcc                                                                                                             <  1:1936566/108‑1 (MQ=255)
                                                                                    agaAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATgcgc                                                              <  2:2405981/150‑1 (MQ=255)
                                                                                     gaAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATgcgcg                                                             >  1:2164574/1‑150 (MQ=255)
                                                                                        gCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCa                                                          >  1:1940549/1‑150 (MQ=255)
                                                                                            aGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGt                                                      <  2:2005527/150‑1 (MQ=255)
                                                                                            aGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGt                                                      >  1:2005527/1‑150 (MQ=255)
                                                                                            aGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTAAGTTACAATGCGCGCCACCGt                                                      <  2:2280074/150‑1 (MQ=255)
                                                                                                cTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCACCCTCCTCTCGGCTAAAGCTGCCAGCTGCCGGCACAGGTTGCATTAGTCAACACCAATTCATACTAATCACCCTTCCTCTCCTTCCCTGAATTACAATGCGCCCCCCCGtagct                                                  >  2:217441/1‑146 (MQ=255)
                                                                                                 tCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGGTGCGTTGGTAGACACCAATTGATACCAATCAGCCTTCATCGCCTTCCCTGCGTTACAATGCGCGCCACCGTGAAtt                                                 >  1:857420/1‑150 (MQ=255)
                                                                                                   aaCAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGc                                                                                            >  2:1580051/1‑105 (MQ=255)
                                                                                                   aaCAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGc                                                                                            <  1:1580051/105‑1 (MQ=255)
                                                                                                      aGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCa                                            >  1:353329/1‑150 (MQ=255)
                                                                                                      aGCCGGCTTGATGATGGCTTGCCTGTAACGGCCCCTTCGCGGCTAAAGCCGGAACAGGCCGGCAGCGGTTGCGGTGGTCGACACCAATTCATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCCCCGTGCATTTGCCa                                            >  2:281491/1‑150 (MQ=255)
                                                                                                             ttGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATgcgcgc                                                            <  1:1875000/127‑1 (MQ=255)
                                                                                                             ttGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATgcgcgc                                                            >  2:1875000/1‑127 (MQ=255)
                                                                                                               gatgatGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTgcg                                   >  1:2313230/1‑150 (MQ=255)
                                                                                                                  gatgGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTa                                                                      <  1:1668155/112‑1 (MQ=255)
                                                                                                                  gatgGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTa                                                                      >  2:1668155/1‑112 (MQ=255)
                                                                                                                  gatgGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCa                                >  2:1469797/1‑150 (MQ=255)
                                                                                                                     ggCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGc                             >  2:397661/1‑150 (MQ=255)
                                                                                                                      gCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCa                            <  2:2313230/150‑1 (MQ=255)
                                                                                                                        ttGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAAtt                                                                                                         >  2:521333/1‑71 (MQ=255)
                                                                                                                        ttGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAAtt                                                                                                         <  1:521333/71‑1 (MQ=255)
                                                                                                                        ttGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCaa                                           >  1:1056223/1‑133 (MQ=255)
                                                                                                                        ttGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCaa                                           <  2:1056223/133‑1 (MQ=255)
                                                                                                                         tGCCCGTACCGGCCTCTCCGCCGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTAGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCCCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGcc                         <  1:463786/150‑1 (MQ=255)
                                                                                                                          gCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGccc                        >  1:1464283/1‑150 (MQ=255)
                                                                                                                           ccTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGAGCGCCATGCAGCCCt                       >  2:1704688/1‑150 (MQ=255)
                                                                                                                            cTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCCTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAACTGCGCCCTGCAGCCCTa                      >  2:906038/1‑150 (MQ=255)
                                                                                                                                 ccGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAg                                                                         <  2:2369064/94‑1 (MQ=255)
                                                                                                                                 ccGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAg                                                                         >  1:2369064/1‑94 (MQ=255)
                                                                                                                                 ccGGCCTCTTCCGGCTTTAGCCGCGAAGGGGCCGGTACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTa                 >  1:1066975/1‑150 (MQ=38)
                                                                                                                                  cGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATGAGCCTTCATCGCTTTCCCTGAGGTACAATGCGCGCCACCGTTAATTTGCCAACAGGGGCGCCATTCAGCCCAACGCTAt                >  2:237798/1‑150 (MQ=255)
                                                                                                                                  cGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTAt                >  1:2081169/1‑150 (MQ=255)
                                                                                                                                   ggCCTCTTCCGGCTTTAGCCGCGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTAtt               >  2:749468/1‑150 (MQ=255)
                                                                                                                                    gCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAAtt                                                                                                         >  2:330476/1‑59 (MQ=255)
                                                                                                                                    gCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAAtt                                                                                                         <  1:330476/59‑1 (MQ=255)
                                                                                                                                      ctctTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCAATGCAGccctaccctaatacc            >  2:734476/1‑145 (MQ=255)
                                                                                                                                      ctctTCGCGGCTAAAGCCGGAAGAGACCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATCCGCGCCACCGTGAATCTGCCAACAGGTGCGCCATAAATCCCTACACTATTGcc            >  2:731556/1‑150 (MQ=255)
                                                                                                                                       tctTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTg                                                                                                                      <  1:76258/43‑1 (MQ=255)
                                                                                                                                       tctTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTg                                                                                                                      >  2:76258/1‑43 (MQ=255)
                                                                                                                                       tctTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTATTGccc           >  1:1433550/1‑150 (MQ=255)
                                                                                                                                       tctTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTATTGccc           >  1:1377711/1‑150 (MQ=255)
                                                                                                                                        ctgcgcgGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTATTGcccc          <  1:493004/147‑1 (MQ=255)
                                                                                                                                        ctTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTATTGcccc          <  1:734476/150‑1 (MQ=255)
                                                                                                                                           cgcgGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATAAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTAc                     >  1:973873/1‑136 (MQ=255)
                                                                                                                                           cgcgGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATAAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTAc                     <  2:973873/136‑1 (MQ=255)
                                                                                                                                            gcgGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTATTGCCCCCTcc      >  1:1447851/1‑150 (MQ=255)
                                                                                                                                            ccgGCTTTAGCCGCGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCGGCCCTACGCTATTGCCCCCTc       >  2:588398/2‑150 (MQ=255)
                                                                                                                                             cgGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTATTGCCCCCTCCa     >  1:1776191/1‑150 (MQ=255)
                                                                                                                                              ggCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTATTGCCCCCTCCAt    >  1:1248810/1‑150 (MQ=255)
                                                                                                                                                cTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTATTGCCCCCTCCCTtc  >  2:2175194/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTATTGCCCCCTCCATTC  >  NC_002947/503688‑503981

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