Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,498,638 T→A 5.0% intergenic (+22/+30) PP_3987 → / ← PP_5651 site‑specific recombinase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,498,6380TA 5.0% 699.4 / 17.6 198intergenic (+22/+30)PP_3987/PP_5651site‑specific recombinase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (93/95);  new base A (4/6);  total (97/101)
Fisher's exact test for biased strand distribution p-value = 7.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

Only 100 of 199 total aligned reads displayed.

TACGCCCCAAGGCGGCAAGCGAAATCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGCTCTTTGGACCTTGCCAGAGCTAGT  >  NC_002947/4498489‑4498785
                                                                                                                                                     |                                                                                                                                                   
 aCGCCCCAAGGCGGCAAGCGAAATCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTc                                                                                                                                                    >  1:1396796/1‑150 (MQ=255)
    ccccAAGGCGGCAAGCGAAATCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTcccc                                                                                                                                                 >  1:652493/1‑150 (MQ=255)
    ccccAAGGCGGCAAGCGAAATCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTcccc                                                                                                                                                 >  1:1634174/1‑150 (MQ=255)
     cccAAGGCGGCAAGCGAAATCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCg                                                                                                                                                <  2:1396796/150‑1 (MQ=255)
      ccAAGGCGGCAAGCGAAATCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGc                                                                                                                                               >  2:3388201/1‑150 (MQ=255)
         aGGCGGCAAGCGAAATCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTATGCTCCCCGCttt                                                                                                                                            <  1:1294760/150‑1 (MQ=255)
         aGGCGGCAAGCGAAATCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCAAAGCGGGGAGCGTACCCCGCttc                                                                                                                                            <  1:1116863/150‑2 (MQ=255)
           gcggcAAGCGAAATCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGa                                                                                                                                          <  2:2153471/150‑1 (MQ=255)
           gcggcAAGCGAAATCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGa                                                                                                                                          <  1:3930011/150‑1 (MQ=255)
            cggcAAGCGAAATCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAg                                                                                                                                         <  2:2843481/150‑1 (MQ=255)
               cAAGCGAAATCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGcac                                                                                                                                      >  2:3300926/1‑150 (MQ=255)
                aaGCGAAATCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCAAAGCGGGGAGCGTACCCCGCTTCGaattac                                                                                                                                     <  2:2291406/150‑6 (MQ=255)
                   cGAAATCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAAtt                                                                                                                                  <  1:1099455/150‑1 (MQ=255)
                     aaaTCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTca                                                                                                                                >  1:4254330/1‑150 (MQ=255)
                      aaTCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTcac                                                                                                                               >  1:2778871/1‑150 (MQ=255)
                        tCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACaa                                                                                                                             <  1:659109/150‑1 (MQ=255)
                          ggCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACt                                                                                                                           >  1:821526/1‑150 (MQ=255)
                          ggCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCAAAGCGGGGAGCGTACCCCGCTTCGaat                                                                                                                                        >  1:3405579/1‑135 (MQ=255)
                           gCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTg                                                                                                                          >  1:1909220/1‑150 (MQ=255)
                             gATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGAt                                                                                                                        <  1:3448459/150‑1 (MQ=255)
                                  gggTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTAtt                                                                                                                   <  2:1776651/150‑1 (MQ=255)
                                  gggTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTAtt                                                                                                                   <  2:2178872/150‑1 (MQ=255)
                                   ggTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTg                                                                                                                  <  1:1783517/150‑1 (MQ=255)
                                     tCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAg                                                                                                                >  2:690706/1‑150 (MQ=255)
                                     tCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAg                                                                                                                >  2:4209892/1‑150 (MQ=255)
                                           aaGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATAc                                                                                                          <  1:767538/150‑1 (MQ=255)
                                                   tATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTa                                                                                                  >  1:2371639/1‑150 (MQ=255)
                                                     ttGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATc                                                                                                >  1:3185664/1‑150 (MQ=255)
                                                        ggTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAAt                                                                                             <  2:3727243/150‑1 (MQ=255)
                                                        ggTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAAt                                                                                             >  2:3544190/1‑150 (MQ=255)
                                                         gTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTGCGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATc                                                                                            >  1:3512504/1‑150 (MQ=255)
                                                             cacTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAg                                                                                        >  1:3503542/1‑150 (MQ=255)
                                                              acTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGc                                                                                       >  2:4463141/1‑150 (MQ=255)
                                                              acTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGc                                                                                       >  2:2430822/1‑150 (MQ=255)
                                                              acTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGc                                                                                       >  1:2693020/1‑150 (MQ=255)
                                                              acTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGc                                                                                       >  1:1144144/1‑150 (MQ=255)
                                                              acTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGc                                                                                       >  1:1378167/1‑150 (MQ=255)
                                                                tGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCtt                                                                                     >  2:2809173/1‑150 (MQ=255)
                                                                tGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCtt                                                                                     >  1:3015234/1‑150 (MQ=255)
                                                                 gACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTtg                                                                                    <  1:646/150‑1 (MQ=255)
                                                                 gACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTtg                                                                                    >  1:2502927/1‑150 (MQ=255)
                                                                  aCAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTtgc                                                                                   >  1:3093412/1‑150 (MQ=255)
                                                                   cAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTtgct                                                                                  <  2:1996793/150‑1 (MQ=255)
                                                                    aaGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTtgctg                                                                                 <  1:1189192/150‑1 (MQ=255)
                                                                      gcgcATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTgctgct                                                                               >  1:4330994/1‑150 (MQ=255)
                                                                        gcATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGctct                                                                             <  1:1204040/150‑1 (MQ=255)
                                                                        gcATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGctct                                                                             <  2:976563/150‑1 (MQ=255)
                                                                         cATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTg                                                                            >  2:2635332/1‑150 (MQ=255)
                                                                           tCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACGCGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGcc                                                                          >  2:2094611/1‑150 (MQ=255)
                                                                               acaAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCAccc                                                                      >  1:2052126/1‑150 (MQ=255)
                                                                               acaAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCAccc                                                                      >  2:3425404/1‑150 (MQ=255)
                                                                                 aaaCGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCaa                                                                    >  2:2764567/1‑150 (MQ=255)
                                                                                  aaCGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAat                                                                   >  1:3612428/1‑150 (MQ=255)
                                                                                    cGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAatca                                                                 <  1:402903/150‑1 (MQ=255)
                                                                                    cGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAatca                                                                 <  2:1413103/150‑1 (MQ=255)
                                                                                        atataGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATcatcat                                                             <  2:1154425/150‑1 (MQ=255)
                                                                                         tataGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATa                                                            <  2:1579179/150‑1 (MQ=255)
                                                                                          ataGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAg                                                           >  1:2486009/1‑150 (MQ=255)
                                                                                              aCGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGttt                                                       >  1:980950/1‑150 (MQ=255)
                                                                                              aCGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGttt                                                       <  2:2397660/150‑1 (MQ=255)
                                                                                                gAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCg                                                     >  2:4090314/1‑150 (MQ=255)
                                                                                                 aGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCga                                                    <  2:1261204/150‑1 (MQ=255)
                                                                                                   tAGGTGAAATAGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCTCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCgaga                                                  <  2:151767/150‑1 (MQ=255)
                                                                                                      gTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGaa                                               >  1:2659371/1‑150 (MQ=255)
                                                                                                      gTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGaa                                               <  1:4407955/150‑1 (MQ=255)
                                                                                                         aaaTCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATAc                                            <  1:2282382/150‑1 (MQ=255)
                                                                                                           aTCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCt                                          <  1:4102974/150‑1 (MQ=255)
                                                                                                           aTCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAAAACCt                                          >  1:1757342/1‑150 (MQ=255)
                                                                                                            tCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTc                                         <  2:652323/150‑1 (MQ=255)
                                                                                                              gTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGc                                       >  2:1930140/1‑150 (MQ=255)
                                                                                                                aCAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCAc                                     >  1:1485020/1‑150 (MQ=255)
                                                                                                                    ccAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGcca                                 <  2:3190467/150‑2 (MQ=255)
                                                                                                                     cAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTg                                >  2:1002659/1‑150 (MQ=255)
                                                                                                                      aaCACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGa                               <  1:1683651/150‑1 (MQ=255)
                                                                                                                      aaCACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGa                               <  1:1062517/150‑1 (MQ=255)
                                                                                                                      aaCACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGa                               >  2:1201660/1‑150 (MQ=255)
                                                                                                                       acacGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAg                              <  2:3626053/150‑1 (MQ=255)
                                                                                                                       acacGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAg                              >  1:2824499/1‑150 (MQ=255)
                                                                                                                         acGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTa                            <  1:3069410/150‑1 (MQ=255)
                                                                                                                          cGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAg                           <  1:4464059/150‑1 (MQ=255)
                                                                                                                           gATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGc                          >  2:76911/1‑150 (MQ=255)
                                                                                                                           gATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGc                          >  1:1817853/1‑150 (MQ=255)
                                                                                                                              aGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGctct                       <  1:597655/150‑1 (MQ=255)
                                                                                                                                gTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATTCGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGCTCttt                     <  2:2711953/150‑1 (MQ=255)
                                                                                                                                gTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGCTCttt                     <  1:3883374/150‑1 (MQ=255)
                                                                                                                                gTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGCTCttt                     <  1:3300926/150‑1 (MQ=255)
                                                                                                                                gTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGCTCttt                     >  1:3635219/1‑150 (MQ=255)
                                                                                                                                gTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGCTCttt                     >  2:1800390/1‑150 (MQ=255)
                                                                                                                                gTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGCTCttt                     >  2:3040639/1‑150 (MQ=255)
                                                                                                                                 tAATTCGAAGCGGGGTACGCTCCCCTCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGCTCTTTg                    <  2:3870282/150‑1 (MQ=255)
                                                                                                                                    ttCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGCTCTTTGGAc                 >  1:2241657/1‑150 (MQ=255)
                                                                                                                                    ttCGAACCGGGCTACGCTCCCCGCTTTGAGCACAATACACAACTCATCTATTGAGCTATACGCAATGTATCAATCTTAGCTTGCAGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGCTCTTTGGAc                 <  2:3191628/150‑1 (MQ=255)
                                                                                                                                     tcaaaGCGGGGAGCGTACCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGCTCTTTGGAcc                >  1:3725487/4‑150 (MQ=255)
                                                                                                                                       gAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGCTCTTTGGACCtt              >  2:64390/1‑150 (MQ=255)
                                                                                                                                        aaGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGCTCTTTGGACCTTg             <  2:4435923/150‑1 (MQ=255)
                                                                                                                                        aaGCGGGGAGCGTACCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGCTCTTTGGACCTTg             >  1:3787715/1‑150 (MQ=255)
                                                                                                                                           cGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGCTCTTTGGACCTTGCCa          <  1:1327205/150‑1 (MQ=255)
                                                                                                                                           cGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGCTCTTTGGACCTTGCCa          <  2:3372497/150‑1 (MQ=255)
                                                                                                                                                tACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGCTCTTTGGACCTTGCCAGAGCt     <  1:896194/150‑1 (MQ=255)
                                                                                                                                                   gCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGCTCTTTGGACCTTGCCAGAGCTAGt  <  1:745499/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TACGCCCCAAGGCGGCAAGCGAAATCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGCTCTTTGGACCTTGCCAGAGCTAGT  >  NC_002947/4498489‑4498785

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