Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,105,036 A→T 5.1% intergenic (+125/+28) PP_4488 → / ← phhR hypothetical protein/sigma‑54 dependent transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,105,0360AT 5.1% 663.5 / 15.3 193intergenic (+125/+28)PP_4488/phhRhypothetical protein/sigma‑54 dependent transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (94/89);  new base T (8/2);  total (102/91)
Fisher's exact test for biased strand distribution p-value = 1.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 196 total aligned reads displayed.

GGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCATCAAACTTTCCAGCA  >  NC_002947/5104888‑5105185
                                                                                                                                                    |                                                                                                                                                     
        ctgctCCACCGGCTGATACCTGGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTccc                                                                                                                                              <  2:348812/150‑1 (MQ=255)
          gctCCACCGGCTTATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTTTTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAg                                                                                                                                            <  1:3855788/150‑1 (MQ=255)
          gctCCACCGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACACGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCGTTGCAGTccttt                                                                                                                                            <  2:143371/150‑4 (MQ=38)
            tCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGACCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTAACAGTGGGAGCGGCCATGTGTGGCGAAAGGACTGCAACGCAGTCCCAGGt                                                                                                                                          >  1:5036791/1‑150 (MQ=39)
            tCCAACGGCTGATACCTTGGAATCGGGGGCCCGGCGGGCCCCGATCGCGGCACAAGGCCCCGCCCCCCCGTAATTGGCGCAGCCCCAGCACACCGGGTTTCCTGTGGGGGCGGCCGTGTGTGGCGCAAGGGCTGGCACCCCGTCCCACGt                                                                                                                                          >  2:1746571/1‑150 (MQ=255)
               aagggCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCTACACAAGGCGGATCCCAAAGTTAAGGCGCTCAGCCAGAGCACAGGGTTTTACCTGTGGGAGGGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCt                                                                                                                                       <  2:2307079/147‑1 (MQ=255)
                 cGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACCGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCttt                                                                                                                                     <  2:931477/150‑1 (MQ=255)
                    cTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGAc                                                                                                                                  <  2:6800893/150‑1 (MQ=255)
                    cTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGAc                                                                                                                                  <  2:499171/150‑1 (MQ=255)
                     tGATACCTTGGAATCGGGGGCCCAGCGGGCCCCTATAGCGACACAAGGCCGCTCCCACAGATAATGCGCTCATCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTAGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGAcc                                                                                                                                 >  2:4441579/1‑150 (MQ=255)
                        tACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACTGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACcaac                                                                                                                              <  2:3577018/150‑1 (MQ=255)
                           cTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCaacaac                                                                                                                           <  2:6317186/150‑1 (MQ=255)
                           cTTGGAATCGGGGGCCCGGCGGGCCCCTATCGCGACACACGGCCGCTCCCACAGATAATGCGCTCAGCCAGCGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCaacaac                                                                                                                           <  1:2191757/150‑1 (MQ=255)
                                 aTCGGGGGCCCGGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACTGATAATGCGCTCAGCCAGAGCACAGCTTATTACCTGTGGGAGGGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATcc                                                                                                                     <  1:5493714/150‑1 (MQ=255)
                                    gggggCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATccgcc                                                                                                                  <  2:2235788/150‑1 (MQ=255)
                                     ggggCCCGCAGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGCGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATccgccg                                                                                                                 >  2:6588798/1‑150 (MQ=255)
                                         cccTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATTCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACtt                                                                                                             >  2:1683648/1‑150 (MQ=255)
                                               gggCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTACTCCGCCGACTTGTTAAg                                                                                                       >  2:5152428/1‑150 (MQ=255)
                                                ggCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCACAGCCAGAGCACAGGGTTTTACCTGGGGTAGCGGCCTTGTGACGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGc                                                                                                      <  1:661665/150‑1 (MQ=255)
                                                  ccccTAGCGCGCCCCCCGGCCGCTCCCACAGATAATGCGCTCACCCTTAGAACAGGTTTTTACCTGTGGGCGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGa                                                                                                    <  1:160858/150‑1 (MQ=255)
                                                    ccTATCGCGACACACGGCCGCTCCCACAGATAAGGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGt                                                                                                  <  2:4643190/150‑1 (MQ=255)
                                                       aTCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAAt                                                                                               <  2:6758047/150‑1 (MQ=255)
                                                        tCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATc                                                                                              <  2:5096546/150‑1 (MQ=255)
                                                          gcgACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATccc                                                                                            >  1:4468015/1‑150 (MQ=255)
                                                            gACACAAGGCCGCTCCCACAGATAAGGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGGCCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTc                                                                                          <  2:1572842/150‑1 (MQ=255)
                                                             acacAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCa                                                                                         >  1:2976661/1‑150 (MQ=255)
                                                                 aaGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGTACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACtt                                                                                     >  2:3695649/1‑150 (MQ=255)
                                                                 aaGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACtt                                                                                     >  2:2727527/1‑150 (MQ=255)
                                                                 aaGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACtt                                                                                     >  2:4742237/1‑150 (MQ=255)
                                                                    gCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACCGCGTATTACCTGTGGGAGCGGCCTTGGGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACttgtt                                                                                  <  1:2662916/150‑1 (MQ=255)
                                                                    gCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACttgtt                                                                                  <  2:260990/150‑1 (MQ=255)
                                                                      cgatCCCACAGATAATGCGCTCTGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCTAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCTAGTAATCCCTCAACTTGTTgg                                                                                <  2:5158154/147‑1 (MQ=255)
                                                                      cGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTgg                                                                                >  2:6150678/1‑150 (MQ=255)
                                                                         tCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCaa                                                                             >  2:10933/1‑150 (MQ=255)
                                                                         ccccACAGATATGGCGCACAGCCATAGAACAGGTTTTTACGTGTGGGAGCGGCCTTGTGCCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGCCTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCaa                                                                             <  1:1746589/149‑1 (MQ=255)
                                                                           ccACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAAtt                                                                           >  1:6172857/1‑150 (MQ=255)
                                                                           ccACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAAtt                                                                           >  2:1101085/1‑150 (MQ=255)
                                                                            cacaGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTg                                                                          >  2:4807163/1‑150 (MQ=255)
                                                                             acaGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGt                                                                         >  1:6269024/1‑150 (MQ=255)
                                                                             acaGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGt                                                                         >  1:830992/1‑150 (MQ=255)
                                                                              caGATAATGCGCTCAGCCAGAGCACAGCGTAGTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTc                                                                        <  2:5261753/150‑1 (MQ=255)
                                                                               aGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCg                                                                       >  1:5873768/1‑150 (MQ=255)
                                                                               aGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCg                                                                       >  1:6011041/1‑150 (MQ=255)
                                                                               aGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCg                                                                       >  1:6521323/1‑150 (MQ=255)
                                                                                gATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTaa                                                                                                        >  1:4143531/1‑116 (MQ=255)
                                                                                  tattGCGCTCCGCCCGAGCACTGCTTATTACCTGTGGGCGCGGCATTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCgtgt                                                                    <  1:1841787/147‑1 (MQ=255)
                                                                                  tAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCgtgt                                                                    <  1:1071593/150‑1 (MQ=255)
                                                                                   aaTGCGCTCCGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGGCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGtgtg                                                                   <  1:699754/150‑1 (MQ=255)
                                                                                   aaTGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGtgtg                                                                   <  2:4513215/150‑1 (MQ=255)
                                                                                     tGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCg                                                                 >  2:3412426/1‑150 (MQ=255)
                                                                                          tCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCgt                                                                      >  1:2606592/1‑140 (MQ=255)
                                                                                          tCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCgt                                                                      <  2:2606592/140‑1 (MQ=255)
                                                                                            aGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACAcc                                                          >  2:505730/1‑150 (MQ=255)
                                                                                              ccAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTaa                                                                                                                        >  1:1422836/1‑86 (MQ=255)
                                                                                               ctgagCACAGCTTATTACCGGGGGGAGCGGCCTTGTGCCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAg                                                       <  1:5732469/148‑1 (MQ=255)
                                                                                               cAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAg                                                       >  1:6241402/1‑150 (MQ=255)
                                                                                               cAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAg                                                       >  1:3072224/1‑150 (MQ=255)
                                                                                               cAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAg                                                       >  2:5505831/1‑150 (MQ=255)
                                                                                               cAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAg                                                       >  2:1894583/1‑150 (MQ=255)
                                                                                               cAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTACAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAg                                                       >  1:3129489/1‑150 (MQ=255)
                                                                                                agagCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGt                                                      >  1:1958458/1‑150 (MQ=255)
                                                                                                   gCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCtt                                                   <  1:2912765/150‑1 (MQ=255)
                                                                                                    cacaGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCttt                                                  <  2:6857023/150‑1 (MQ=255)
                                                                                                     accgCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCtttt                                                 <  1:4045896/147‑1 (MQ=255)
                                                                                                     acaGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCtttt                                                 <  2:1757757/150‑1 (MQ=255)
                                                                                                     acaGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCtttt                                                 <  1:4045901/150‑1 (MQ=255)
                                                                                                     acaGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCtttt                                                 <  1:2543133/150‑1 (MQ=255)
                                                                                                           tATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCaa                                           <  1:4459317/150‑1 (MQ=255)
                                                                                                            aTTACCTGTGGGAGCGGCCTTGTGTCGCTAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAg                                          <  2:2331392/150‑1 (MQ=255)
                                                                                                             ttACCTGTGGTCGCGGCATTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAgg                                         <  1:2309921/150‑1 (MQ=255)
                                                                                                             ttACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAgg                                         >  1:3079678/1‑150 (MQ=255)
                                                                                                                 cTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCg                                     >  2:1048256/1‑150 (MQ=255)
                                                                                                                  tgtgGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGa                                    >  2:4171558/1‑150 (MQ=255)
                                                                                                                    tgGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGact                                  >  2:5492342/1‑150 (MQ=255)
                                                                                                                      ggAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGactac                                <  2:2976661/150‑1 (MQ=255)
                                                                                                                        aGCGGCCTTGTGTCGCGAAAGGACAGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTg                              <  1:1036519/150‑1 (MQ=255)
                                                                                                                           ggCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGa                           >  1:3932751/1‑150 (MQ=255)
                                                                                                                             ccTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGaaa                         >  1:5691809/1‑150 (MQ=255)
                                                                                                                                 gtgtCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTcc                     <  1:1682778/150‑1 (MQ=255)
                                                                                                                                 gtgtCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTcc                     <  1:2324822/150‑1 (MQ=255)
                                                                                                                                   gtCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTcccc                   >  2:1125510/1‑150 (MQ=255)
                                                                                                                                       cgAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCATc               <  1:2361915/150‑1 (MQ=255)
                                                                                                                                       cgAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCATc               >  2:5741310/1‑150 (MQ=255)
                                                                                                                                        cctgggACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCATCa              >  1:5525501/5‑150 (MQ=255)
                                                                                                                                         ctgggACTGCGTTGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCATCaa             >  1:4905406/4‑150 (MQ=255)
                                                                                                                                         ctgggACTGCGTTGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCATCaa             >  1:6501012/4‑150 (MQ=255)
                                                                                                                                         aaaGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTATTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCATCaa             >  1:1584192/1‑150 (MQ=255)
                                                                                                                                          tgggACTGCGTTGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGACACTCCCCCATCaaa            >  2:3775077/3‑150 (MQ=255)
                                                                                                                                          tgggACTGCGTTGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCATCaaa            >  2:1418197/3‑150 (MQ=255)
                                                                                                                                          tgggACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGaaa                                                              >  1:3187154/3‑100 (MQ=255)
                                                                                                                                          aaGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCATCaaa            <  1:6422572/150‑1 (MQ=255)
                                                                                                                                          aaGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCATCaaa            <  1:3695649/150‑1 (MQ=255)
                                                                                                                                               cTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCATCAAACTTTc       >  2:3356529/1‑150 (MQ=255)
                                                                                                                                               cTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCATCAAACTTTc       <  1:6303105/150‑1 (MQ=255)
                                                                                                                                                tGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCATCAAACTTTcc      >  2:3125616/1‑150 (MQ=255)
                                                                                                                                                tGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCATCAAACTTTcc      >  2:6007715/1‑150 (MQ=255)
                                                                                                                                                tGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCATCAAACTTTcc      >  1:6310535/1‑150 (MQ=255)
                                                                                                                                                tGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCATCAAACTTTcc      <  2:5036791/150‑1 (MQ=255)
                                                                                                                                                tGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCATCAAACTTTcc      >  2:833816/1‑150 (MQ=255)
                                                                                                                                                   aaCGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCATCAAACTTTCCAGc   >  2:2928345/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCATCAAACTTTCCAGCA  >  NC_002947/5104888‑5105185

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