Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,127,066 A→T 6.1% intergenic (+70/‑112) PP_4511 → / → PP_4512 AraC family transcriptional regulator/DMT superfamily permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,127,0660AT 6.1% 398.0 / 3.7 128intergenic (+70/‑112)PP_4511/PP_4512AraC family transcriptional regulator/DMT superfamily permease
Reads supporting (aligned to +/- strand):  ref base A (42/77);  new base T (6/2);  total (49/79)
Fisher's exact test for biased strand distribution p-value = 5.23e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 131 total aligned reads displayed.

AAAGTGCCTTTTCAGCGGCGATGCTGCGGGCGTTCGGGTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGG  >  NC_002947/5126917‑5127215
                                                                                                                                                     |                                                                                                                                                       
aaaGTGCCTTTTCAGCGGCGATGCTGCGGGCGTTCGGGTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTa                                                                                                                                                         >  2:1382235/1‑150 (MQ=255)
  aGTGCCTTTTCAGCGGCGATGCTGCGGGCGTTCGGGTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTaaa                                                                                                                                                       <  2:2323512/150‑1 (MQ=255)
    tGCCTTTTCAGCGGCGATGCTGCGGGCGTTCGGGTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAAcc                                                                                                                                                     <  2:3191709/150‑1 (MQ=255)
     gCCTTTTCAGCGGCGATGCTGCGGGCGTTCGGGTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAAccc                                                                                                                                                    <  1:2512288/150‑1 (MQ=255)
        ttttCAGCGGCGATGCTGCGGGCGTTCGGGTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCt                                                                                                                                                 <  1:1040163/150‑1 (MQ=255)
        ttttCAGCGGCGATGCTGCGGGCGTTCGGGTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCt                                                                                                                                                 <  1:674608/150‑1 (MQ=255)
         tttCAGCGGCGATGCTGCGGGCGTTCGGGTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACGCTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTc                                                                                                                                                <  2:3167196/150‑1 (MQ=255)
         tttCAGCGGCGATGCTGCGGGCGTTCGGGTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTc                                                                                                                                                <  2:1768013/150‑1 (MQ=255)
          ttCAGCGGCGATGCTGCGGGCGTTCGGGTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTTTACCCGcgaa                                                                                                                                               <  2:1717142/150‑4 (MQ=255)
          ttCAGCGGCGATGCTGCGGGCGTTCGGGTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTTTACCCGcgaa                                                                                                                                               <  2:4481974/150‑4 (MQ=255)
             aGCGGCGATGCTGCGGGCGTTCGGGTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTAc                                                                                                                                            <  2:3611004/150‑1 (MQ=255)
              gcggcgATGCTGCGGGCGTTCGGGTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACa                                                                                                                                           >  1:248661/1‑150 (MQ=255)
              gcggcgATGCTGCGGGCGTTCGGGTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACa                                                                                                                                           <  1:1043921/150‑1 (MQ=255)
               cggcgATGCTGCGGGCGTTCGGGTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGGAAACCCGCGCCTACAg                                                                                                                                          >  2:4240728/1‑150 (MQ=255)
                ggcgATGCTGCGGGCGTTCGGGTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAgg                                                                                                                                         <  1:1829296/150‑1 (MQ=255)
                ggcgATGCTGCGGGCGTTCGGGTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAgg                                                                                                                                         >  2:2539617/1‑150 (MQ=255)
                 gcgATGCTGCGGGCGTTCGGGTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAggt                                                                                                                                        >  2:3013164/1‑149 (MQ=255)
                 gcgATGCTGCGGGCGTTCGGGTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGa                                                                                                                                        <  2:3455603/150‑1 (MQ=255)
                   gATGCTGCGGGCGTTCGGGTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACCGGATg                                                                                                                                      >  2:3512992/1‑150 (MQ=255)
                    aTGCTGCGGGCGTTCGGGTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGc                                                                                                                                    >  2:697181/1‑150 (MQ=255)
                    aTGCTGCGGGCGTTCGGGTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGc                                                                                                                                    >  2:697136/1‑150 (MQ=255)
                    aTGCTGCGGGCGTTCGGGTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAAACCGCGCCGACAGGATGc                                                                                                                                    >  2:697168/1‑150 (MQ=255)
                      gctgcGGGCGTTCGGGTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGccc                                                                                                                                  <  2:727472/150‑1 (MQ=255)
                           gggCGTTCGGGTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGcc                                                                                                                             >  2:2909054/1‑150 (MQ=255)
                            ggCGTTCGGGTGTCCACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCt                                                                                                                            <  1:2751409/150‑1 (MQ=255)
                              cGTTCGGGTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGa                                                                                                                          <  1:4153575/150‑1 (MQ=255)
                               gTTCGGGTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGc                                                                                                                                    >  2:1439446/1‑139 (MQ=255)
                               gTTCGGGTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCCCTCCTACAGGATGCCCCCGCCTGaa                                                                                                                         >  2:1871251/1‑150 (MQ=255)
                                  cGGGTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAAtgt                                                                                                                      <  2:573498/150‑1 (MQ=255)
                                   gggTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAAtgtt                                                                                                                     <  2:248661/150‑1 (MQ=255)
                                     gtgtACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAAtgttgt                                                                                                                   <  1:4204574/150‑1 (MQ=255)
                                     gtgtACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAAtgttgt                                                                                                                   <  1:4497299/150‑1 (MQ=255)
                                     gtgtACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAAtgttgt                                                                                                                   <  1:4204630/150‑1 (MQ=255)
                                      tgtACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTtgtg                                                                                                                  <  1:2909054/150‑1 (MQ=255)
                                      tgtACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTtgtg                                                                                                                  <  2:3039915/150‑1 (MQ=255)
                                      tgtACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTtgtg                                                                                                                  <  1:2824975/150‑1 (MQ=255)
                                         acacCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCaa                                                                                                               >  1:1013828/1‑150 (MQ=255)
                                            ccTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACTGGATGCCCCCGCCTGAATGTTGTGCAAAAc                                                                                                            <  2:3778457/150‑1 (MQ=255)
                                            ccTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAAc                                                                                                            <  2:283045/150‑1 (MQ=255)
                                            ccTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAAc                                                                                                            <  1:4234251/150‑1 (MQ=255)
                                              tttGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCACCTTCGCGGGGAAACCCGCGACAAAAGGAATCCCCCGCCTGAATGTTTTGCAAAAACt                                                                                                          >  2:1284313/1‑150 (MQ=255)
                                                 ggCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTa                                                                                                       <  2:807718/150‑1 (MQ=255)
                                                 ggCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTa                                                                                                       <  1:3767711/150‑1 (MQ=255)
                                                 ggCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTa                                                                                                       <  1:1733627/150‑1 (MQ=255)
                                                 ggCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTa                                                                                                       <  2:1334976/150‑1 (MQ=255)
                                                 ggCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTa                                                                                                       <  1:1913066/150‑1 (MQ=255)
                                                  gCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAAtgttgt                                                                                                                   >  1:718553/1‑137 (MQ=255)
                                                  gCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCACCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAg                                                                                                      >  1:1898049/1‑150 (MQ=255)
                                                    tttGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGa                                                                                                    <  2:3589277/150‑1 (MQ=255)
                                                      tGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGc                                                                                                  <  2:422915/150‑1 (MQ=255)
                                                      tGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGc                                                                                                  <  2:423121/150‑1 (MQ=255)
                                                      tGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCATCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGc                                                                                                  <  2:1225316/150‑1 (MQ=255)
                                                         ggcgAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCggg                                                                                              <  2:1013828/150‑1 (MQ=255)
                                                          gcgAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGt                                                                                             <  2:542785/150‑1 (MQ=255)
                                                          gcgAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTATAGGAGCGGGt                                                                                             <  2:1732930/150‑1 (MQ=255)
                                                           cgAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAAtgtt                                                                                                                     <  1:2233244/126‑1 (MQ=255)
                                                           cgAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCt                                                                                                          >  2:190093/1‑137 (MQ=255)
                                                               agTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCCCTCCTAAAGGATTCCCCCCCCCGGATGTTGTGCAAAACCTGTAGGAGCGGGTTTAcc                                                                                        >  2:1520660/1‑150 (MQ=255)
                                                                gTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTAccc                                                                                       <  1:2570909/150‑1 (MQ=255)
                                                                 tCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGGTAACCCCCGACGACCCGCAGCCCCCCCCCGCAAGTTGTGGAAAACCGGTAGGAGCGGGGTTTcccc                                                                                      >  1:3169434/1‑149 (MQ=255)
                                                                      cGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGaaa                                                                                 <  2:3471049/150‑1 (MQ=255)
                                                                         tAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACCGTATGCCCCCGCCTGAATGTTGTGCAAAACCTGAAGGAGCGGGTTTACCCGCGAAAggg                                                                              <  2:2445338/150‑1 (MQ=255)
                                                                           aCTGACAGCTTTCTTGCGACAGAATGTGTTTGCCTGTGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGGGGGTAAACCCGCTCCTACAGTATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGGGGTTTTCCCCGCGAAAGGGcc                                                                            <  2:1162849/150‑1 (MQ=38)
                                                                           aCTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAAACCGCTCCAACACGATGCCCCCCCCGGAATGTTTTTCAAAAACGGTAGGAGCGGGGTTTACCCCGAAAGGGcc                                                                            >  1:977566/1‑150 (MQ=255)
                                                                            cTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGGTCACCCCCGCCAACCCGCGGCCCCCCCCGGGATGTTGTGGAAAACCAGGTGGAGGGGGTGTTACCGCCGAAGGGCCg                                                                           >  1:3224603/1‑150 (MQ=255)
                                                                             tGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTTAAACCGCGGCTAAAGGATTCCCCCCCCGGAATGTTGTGCACAACCGGTAGGAGCGGGGTTTCCCGCGAAAGGGCCgg                                                                          >  2:2811856/1‑150 (MQ=255)
                                                                             tGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCCCCGGAATGTTGTGCCAAACCTGTAGGAGCGGGGTTACCCGCGCAAGGGCCgg                                                                          >  2:2511459/1‑150 (MQ=38)
                                                                             tGACAGCAGACTTGCTACTTAATGTGTTTGCCTGAGGCTTTACAATTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGATCATACAGCATGCCCCCGCATTAATTTTGTGAAAAACCTTTAGGAGCGGTTTTACCCGCGAAAGGGCCgg                                                                          <  2:179657/150‑1 (MQ=255)
                                                                                 agctgTCTTGCGACAGATTGTGTTTGCCTGCGGCTTTACACTTTGCCTTGGCGGGCCCCTTCGCGGGTAAACCCGCACCTACCGGATGCCCCCGCCTTAATGTTTTGCAAAACCGGTAGGAGGGGGTTTACCCGCGAAAGGGCCGGCACt                                                                      <  2:1166392/150‑1 (MQ=38)
                                                                                 agcagTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCCCGCCTAAAGGATGGCCCCGCCGGAATGTTGTGCAAAACCAGTAGGAGCGGGTTTTCCCCCGACAGGGCCGGCAcg                                                                      >  1:4456477/1‑149 (MQ=255)
                                                                                    cgTCTTGCGACATATTGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCGGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGATTGTTGTGCAAAACCTGTAGGAGGGGGTTTACCCGCGAAAGGGCCGGCACTGcc                                                                   <  2:3093690/149‑1 (MQ=38)
                                                                                    agTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTTAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGcc                                                                   <  2:3093686/150‑1 (MQ=255)
                                                                                      tCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCAACAGGATGCCCCCCCCTGAATGTTGTGCAAAAACGGTATGATCGGGTTTAACCGCGAAAGGGCCGGCACTGccgg                                                                 >  2:285213/1‑148 (MQ=38)
                                                                                        ttGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGGAAAACCGCGCCGAAACGCTGCCCCCCCCGGAATGTTGTGGAAAACCGGTAGGAGCGGGTTTTGACCCGGAAGGGCCCGCCCGGcccggt                                                               >  1:2171134/1‑146 (MQ=255)
                                                                                         tGCTACTGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTag                                                              <  2:2885166/150‑1 (MQ=255)
                                                                                            gACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCCAGTAgaga                                                           >  2:3243509/1‑150 (MQ=255)
                                                                                                       tttGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCGGAATGTTGTGCAAAACCGGGAGGAGCGGGTTTACCCCCGAAAGGGCCGGCCCGGCCGGTAGAGATCTCCATCAAg                                                >  2:365688/1‑150 (MQ=38)
                                                                                                       tttGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAAACCGCTCCAACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCAGGAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCGGTAGAGATCTCCATAAAg                                                >  2:4125458/1‑150 (MQ=38)
                                                                                                              gtggCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGTTAAACCCGATCCTACCGGTTGCCCCCGCCTGATTTTTGTGCAAAACCTGTTGGTGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACccac                                         <  2:621071/148‑1 (MQ=38)
                                                                                                                   tttACACTTTGCGTGTGCCGGCCCGTTCGCGGGTAAACCCGCTCCTCCCGGAGGCCCCGGCTTGATTGTTGTGCAAAACCTGTAGGTGCGGTTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGa                                    <  2:397045/150‑1 (MQ=37)
                                                                                                                   tttACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAAACCGCTCCTACAGGATGCCCCCGCCGGTATTTTGTGCAAAACCAGTAGGAGCGGGTTTAACCGCGAAAGGGCCGGCACGGCCGGTAGAGATATCCATAAAGGACCCACCAAGa                                    >  2:2190309/1‑150 (MQ=255)
                                                                                                                   tttACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAAACCGCACCAAAAGGATGCCCCCGCCGGTATGTTGTTCAAAAGATGTAGGAGCGGGTTTAACCGCGAAAGGGCCGGCACTGGCTGTAGAGATCACCATCAAGGAACCAACAAGa                                    >  2:2190319/1‑150 (MQ=255)
                                                                                                                   tttACACTTTGCCTGTGCCGGCCCCTGCGCGGGTAAACCCGCTCCTACCGGATGCCCCCGCCTGAATGTTGTGCAAAACCTTTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGa                                    <  2:2652627/150‑1 (MQ=255)
                                                                                                                       cacTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCCCTACTACAGGATGCCCCCCGCCGGAATGTTGTGCAAAACCGGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTACATCAAGGAACCACCAAGAacc                                 >  2:2933493/1‑150 (MQ=37)
                                                                                                                          tttgcctttgccGGCCCCTTCGCGGGGAAACCCCCGCCGAAAGGATGCCCCCCCCGGAAAGTTGTGCAAAAACGGTAGGAGCCGGGTTACCCCCCCAAGGGCCCGCCCGGCCCGGAGAGATCGCCAACAACGACCCCCCAAGAaacaccc                             >  2:3425752/1‑149 (MQ=255)
                                                                                                                           ttgcctttgccGGCCCCTTCGCGGGTAAAACCGCTCCTACATGATGCCCCCGCCGGAATGTTTTGGAAAACCTGGTGGAGCGGGTTTACCCGCGGAAGGGCCGGCACTGCCTGTAGAGATCTACATCAAGGACCCACCATGAAAcaacgc                            >  2:737801/1‑150 (MQ=255)
                                                                                                                               ctttgccGGCCCCTTCGCGGGTAAACCCCCGCCTAAAGGATGCCCCCCCCAGAATGTTTTGCCAAAACTGTAGGAGCGGGGTTAACCGCGAAAGGGCCGGCACAGCCGGTAGAGATCTCCATCAAGGACCCACCATGTACcaccgcaccg                        >  2:2695300/1‑150 (MQ=255)
                                                                                                                                       gCCCCTTCGCGGTTAAACCCGCTCCTACAGGATGCCCCCCGCCTTATTGTTGTGCAAAACCTGTAGGAGCGGTTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCgg                 <  1:2933493/150‑1 (MQ=38)
                                                                                                                                         cccTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGGGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGgcgc              <  2:3028250/150‑1 (MQ=255)
                                                                                                                                         cccTTCGCGGGTAAACCCGCGCCTACAGGATGCCCCCCCCGGTATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCGGCCGGTAGTGGTCTCCAGCAAGGGCCCACCAAGAACCACCGCACCGCCCGCGgcgc              >  2:4009315/1‑150 (MQ=37)
                                                                                                                                         cccTTCGCGGGTAAACCCCCGCCTACACGATGCCCCCGCCGGAATGTTGTTCAAAAACAGTAGGAGCGGGGTTAACCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCAACAAGGACCCACCAAGAACCACCGCACCGCCCACGgcgc              >  2:4373702/1‑150 (MQ=255)
                                                                                                                                         cccTTCGCGGGTAAACCCCCACCTAAACGATTCCCCCCCCGGAATGTTGTGCAAAACCAGTAGGAGCGGGTTTACCCGCCCAAGGGCCGGCCCAGCCAGTAGAGAAAACCAACAAAGACCCAACAAGAAACAACCCACCCCCCACGgcgc              >  2:1349076/1‑150 (MQ=255)
                                                                                                                                          ccTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGCGCAAACCCTGTAGGAGCGGGATTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCa             <  1:1382235/150‑1 (MQ=255)
                                                                                                                                            ttCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACt           <  2:4404741/150‑1 (MQ=255)
                                                                                                                                            ttCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCAGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACAGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACt           >  2:1329281/1‑150 (MQ=255)
                                                                                                                                              cgcgGGTAAACCCGCTCCTACAGGAGGCCCCCGCCTGAATGTTGTGCAACACCTGTCGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGc         <  2:3062275/150‑1 (MQ=255)
                                                                                                                                               gcgGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCACGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGca        >  2:1732833/1‑150 (MQ=255)
                                                                                                                                                 gggTAAACCCGCACCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGcaca      >  1:3189343/1‑150 (MQ=255)
                                                                                                                                                  ggTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACAt     >  2:1529066/1‑150 (MQ=255)
                                                                                                                                                    tAAACCCGCTCCTCCAGGATGCCCCCTCCTGTTTGTTGCGCTACCCGTGTCGGAGGGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTg   <  1:4128924/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                       
AAAGTGCCTTTTCAGCGGCGATGCTGCGGGCGTTCGGGTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGG  >  NC_002947/5126917‑5127215

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