Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 182,035 A→T 10.3% intergenic (+24/+41) PP_5SB → / ← PP_0160 5S ribosomal RNA/ferrioxamine receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947182,0350AT10.3% 545.8 / 40.6 174intergenic (+24/+41)PP_5SB/PP_01605S ribosomal RNA/ferrioxamine receptor
Reads supporting (aligned to +/- strand):  ref base A (64/92);  new base T (10/8);  total (74/100)
Fisher's exact test for biased strand distribution p-value = 3.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 177 total aligned reads displayed.

CAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAGCTTTCATTGGTAAGGTTGTTGGCG  >  NC_002947/181886‑182184
                                                                                                                                                     |                                                                                                                                                     
caGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGa                                                                                                                                                       >  1:4177530/1‑150 (MQ=255)
caGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGa                                                                                                                                                       <  2:3477419/150‑1 (MQ=255)
  gAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATg                                                                                                                                                     <  2:3293933/150‑1 (MQ=255)
   aaTTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGc                                                                                                                                                    >  1:2574484/1‑150 (MQ=255)
     tttCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAg                                                                                                                                                  >  2:3851567/1‑150 (MQ=255)
       tCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGAc                                                                                                                                                >  1:2475851/1‑150 (MQ=255)
         ttGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAg                                                                                                                                              >  1:1878123/1‑150 (MQ=255)
          tGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCATCGCAGACAgg                                                                                                                                             <  2:1760491/150‑1 (MQ=255)
           gacgagCATAGCGCATTGGACCCACCTGTTCCCTTCCGGACCTCAGCAGTGTCACGATGCCACGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCATCGCAGACAggg                                                                                                                                            <  1:895834/150‑1 (MQ=17)
            acgacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAgggg                                                                                                                                           <  1:2797484/150‑1 (MQ=255)
                ccATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCttt                                                                                                                                       <  1:2029692/150‑1 (MQ=255)
                ccATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCttt                                                                                                                                       <  1:1989507/150‑1 (MQ=255)
                ccAGAGCGCCTTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCAGCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCttt                                                                                                                                       <  2:919869/150‑1 (MQ=255)
                           tGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCgg                                                                                                                            <  1:1814086/150‑1 (MQ=255)
                             gAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGt                                                                                                                          <  2:2196266/150‑1 (MQ=255)
                             gAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGt                                                                                                                          <  2:1110392/150‑1 (MQ=255)
                                ccaccTGATCCCATCCCTAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTc                                                                                                                       <  1:242688/150‑1 (MQ=255)
                                ccaccTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTc                                                                                                                       <  2:4094360/150‑1 (MQ=255)
                                 caccTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCa                                                                                                                      <  1:883259/150‑1 (MQ=255)
                                 caccTGATCACATCCCGAAATCAGCAGTGAAACGATGAATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAATCCCAATGACCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCa                                                                                                                      >  2:2501830/1‑150 (MQ=18)
                                  accTGATCCCATCCCGAACTGAGCAGTGAACCGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCaa                                                                                                                     <  1:2333275/150‑1 (MQ=255)
                                  accTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGAt                                                                                                                                                      >  1:420788/1‑117 (MQ=255)
                                   ccTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAat                                                                                                                    <  2:465012/150‑1 (MQ=255)
                                   ccTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCGTTGTCTTTTCCGGGTTTCAat                                                                                                                    >  1:2037264/1‑150 (MQ=255)
                                       atcccatcccGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCg                                                                                                                <  1:478285/150‑1 (MQ=255)
                                       atcccatcccGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCATTGTCTTTTCCGGGTTTCAATATCg                                                                                                                <  2:2489700/150‑1 (MQ=255)
                                         cccatcccGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGt                                                                                                              <  2:1603089/150‑1 (MQ=255)
                                         cacatccaGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTACCCATGTGAGAGTAGGTCATTGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTACGGGTTTCAATATAGGt                                                                                                              >  2:3256047/3‑150 (MQ=21)
                                          ccatcccGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGtt                                                                                                             >  1:710974/1‑150 (MQ=255)
                                          ccatcccGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGtt                                                                                                             >  1:4078685/1‑150 (MQ=255)
                                          ccatcccGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGtt                                                                                                             >  2:599712/1‑150 (MQ=255)
                                           catcccGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTa                                                                                                            <  1:1752533/150‑1 (MQ=255)
                                           catcccGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTa                                                                                                            <  2:2982059/150‑1 (MQ=255)
                                           catcccGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCAGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTa                                                                                                            <  2:4044575/150‑1 (MQ=255)
                                               ccGAACTCAGCAGTGAAACGGTGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAAc                                                                                                        >  1:1645819/1‑150 (MQ=255)
                                               ccGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAAc                                                                                                        <  1:3563733/150‑1 (MQ=255)
                                                  aaCTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGa                                                                                                     <  2:1318467/150‑1 (MQ=255)
                                                  aaCTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGa                                                                                                     >  2:3288552/1‑150 (MQ=255)
                                                    cTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATa                                                                                                   <  1:3179250/150‑1 (MQ=255)
                                                     tCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCGTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAg                                                                                                  >  1:599395/1‑150 (MQ=255)
                                                         cagTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGa                                                                                              <  1:966253/150‑1 (MQ=255)
                                                         cagTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGa                                                                                              >  1:1689712/1‑150 (MQ=255)
                                                                         cGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACAtt                                                                              <  2:433167/150‑1 (MQ=255)
                                                                         cGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACAtt                                                                              <  2:3739354/150‑1 (MQ=255)
                                                                           ccGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTgc                                                                            <  1:3985045/150‑1 (MQ=255)
                                                                           ccGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTgc                                                                            <  1:4284645/150‑1 (MQ=255)
                                                                             gATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCAGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGGTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTgcgc                                                                          <  1:1092712/150‑1 (MQ=255)
                                                                               tGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGctc                                                                        <  2:2254586/150‑1 (MQ=255)
                                                                                 ggaGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCtt                                                                      <  2:2498468/148‑1 (MQ=255)
                                                                                 gTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCtt                                                                      <  2:521880/150‑1 (MQ=255)
                                                                                  tAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTc                                                                     >  1:1322354/1‑150 (MQ=255)
                                                                                   aGTGTGGGGTGTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCa                                                                    <  1:619251/150‑1 (MQ=255)
                                                                                      gtgGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCAccc                                                                 <  1:1279410/150‑1 (MQ=255)
                                                                                       tgGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCa                                                                >  2:1964386/1‑150 (MQ=255)
                                                                                       tgGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCa                                                                >  2:1996324/1‑150 (MQ=255)
                                                                                        ggggTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAATCCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCAt                                                               >  2:2302067/1‑150 (MQ=255)
                                                                                        ggggTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCAt                                                               >  1:1351600/1‑150 (MQ=255)
                                                                                         gggTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATg                                                              <  1:3376643/150‑1 (MQ=255)
                                                                                         gggTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATg                                                              <  1:3376639/150‑1 (MQ=255)
                                                                                             ttCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGc                                                          <  2:2754073/150‑1 (MQ=255)
                                                                                                 ccATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCagta                                                      <  1:2600700/150‑1 (MQ=255)
                                                                                                  cATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCACCGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCagtag                                                     <  2:3980767/150‑1 (MQ=255)
                                                                                                  cATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCagtag                                                     <  1:2742120/150‑1 (MQ=255)
                                                                                                     gtgAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGATGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTc                                                  >  2:4452129/1‑150 (MQ=255)
                                                                                                      tgAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCa                                                 <  1:3764507/150‑1 (MQ=255)
                                                                                                      tgAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCa                                                 <  2:2072604/150‑1 (MQ=255)
                                                                                                        agagTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAgg                                               >  1:3873756/1‑150 (MQ=255)
                                                                                                        agagTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAgg                                               <  1:2474725/150‑1 (MQ=255)
                                                                                                          agTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCt                                             >  2:530186/1‑150 (MQ=255)
                                                                                                             aGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGc                                          <  2:3571856/150‑1 (MQ=255)
                                                                                                                   tCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGa                                    <  1:3726765/150‑1 (MQ=255)
                                                                                                                      tCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGt                                 >  1:1107596/1‑150 (MQ=255)
                                                                                                                        aaGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGa                               <  2:177876/150‑1 (MQ=255)
                                                                                                                        aaGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGa                               <  1:3289958/150‑1 (MQ=255)
                                                                                                                        aaGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGa                               <  2:2818768/150‑1 (MQ=255)
                                                                                                                             tCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAg                          >  2:2271930/1‑150 (MQ=255)
                                                                                                                                aaCCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAGCtt                       <  1:2004285/150‑1 (MQ=255)
                                                                                                                                 aCCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAGCttt                      <  2:2319935/150‑1 (MQ=255)
                                                                                                                                   ccAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAGCTTTCa                    <  2:4177530/150‑1 (MQ=255)
                                                                                                                                    cAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTa                                                            <  1:1037492/109‑1 (MQ=255)
                                                                                                                                     aaaGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAAc                                                                                                        >  2:3900960/1‑64 (MQ=255)
                                                                                                                                     aaaGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAGCTTTCAtt                  >  1:1726042/1‑150 (MQ=255)
                                                                                                                                     aaaGCCCCTGTCTGCATCGCAGACAGGGGCTTTGTCTTTTCCGGGtt                                                                                                                         <  1:3459831/47‑1 (MQ=14)
                                                                                                                                     aaaGCCCCTGTCTGCATCGCAGACAGGGGCTTTGTCTTTTCCGGGtt                                                                                                                         >  2:3459831/1‑47 (MQ=14)
                                                                                                                                     aaaGCCCCTGTCTGCATCGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAGCTTTCAtt                  >  2:3673428/1‑150 (MQ=255)
                                                                                                                                     aaaGCCCCTGTCTGCATCGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAGCTTTCAtt                  >  2:1536493/1‑150 (MQ=255)
                                                                                                                                     aaaGCCACTGTCTGCATCGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGTTAGCTGACCGTGGCGCTCACATTGAGCTCTTCTCCCTTGTAGCAGTAGTTCAGGCTTGCACATGAGGTGATGTAGCTTTCAtt                  >  2:4307054/1‑150 (MQ=255)
                                                                                                                                      aaGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAGCTTTCATTg                 <  2:784191/150‑1 (MQ=255)
                                                                                                                                       aGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAGCTTTCATTgg                <  1:1918597/150‑1 (MQ=255)
                                                                                                                                       aGCCCCTGTCTGCATCGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAGCTTTCATTgg                >  2:1892399/1‑150 (MQ=255)
                                                                                                                                         ccccTGTCTGCATCGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAGCTTTCATTGGTa              >  1:2076/1‑150 (MQ=255)
                                                                                                                                         ccccTGTCTGCATCGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAGCTTTCATTGGTa              >  2:1965465/1‑150 (MQ=255)
                                                                                                                                          cccTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAGCTTTCATTGGTaa             <  1:4340728/150‑1 (MQ=255)
                                                                                                                                           ccTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAGCTTTCATTGGTAAg            <  1:2512983/150‑1 (MQ=255)
                                                                                                                                            cTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAGCTTTCATTGGTAAgg           <  2:2574484/150‑1 (MQ=255)
                                                                                                                                            cTGTCTGCATCGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAGCTTTCATTGGTAAgg           >  2:1265603/1‑150 (MQ=255)
                                                                                                                                             tGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAGCTTTCATTGGTAAGgt          >  1:1417920/1‑150 (MQ=255)
                                                                                                                                              gTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAGCTTTCATTGGTAAGgtt         <  1:90993/150‑1 (MQ=255)
                                                                                                                                                  gCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAGCTTTCATTGGTAAGGttgttg     <  2:1140794/150‑1 (MQ=255)
                                                                                                                                                    gATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAGCTTTCATTGGTAAGGTTGTTGGc   >  2:3127062/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAGCTTTCATTGGTAAGGTTGTTGGCG  >  NC_002947/181886‑182184

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