Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 173,453 T→A 9.3% intergenic (+9/+149) fhuB → / ← hemL iron(3+)‑hydroxamate import ABC transporter permease/glutamate‑1‑semialdehyde aminotransferase (aminomutase)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913173,4530TA 9.3% 664.6 / 6.8 253intergenic (+9/+149)fhuB/hemLiron(3+)‑hydroxamate import ABC transporter permease/glutamate‑1‑semialdehyde aminotransferase (aminomutase)
Reads supporting (aligned to +/- strand):  ref base T (124/95);  new base A (30/2);  total (155/98)
Fisher's exact test for biased strand distribution p-value = 2.22e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 275 total aligned reads displayed.

TGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGAT  >  NC_000913/173305‑173602
                                                                                                                                                     |                                                                                                                                                     
tggtggGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCata                                                                                                                                                      >  2:1673619/1‑150 (MQ=255)
tggtggGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCata                                                                                                                                                      >  1:214451/1‑150 (MQ=255)
 ggtgggtgGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCatat                                                                                                                                                     >  1:3027914/1‑150 (MQ=255)
 ggtgggtgGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCatat                                                                                                                                                     >  1:3322136/1‑150 (MQ=255)
    gggtggTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAAc                                                                                                                                                  >  2:1626137/1‑150 (MQ=255)
    gggtggTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTACATATAAc                                                                                                                                                  >  2:2869643/1‑150 (MQ=255)
     ggtggtTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTGATTTTTGCATATaaag                                                                                                                                                 >  1:2531516/1‑148 (MQ=255)
     ggtggtTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTAGTTGTTGAGAAAGCAGAGCCGTTAAATTTTAAATATaaag                                                                                                                                                 >  1:483966/1‑148 (MQ=255)
       tggtTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGc                                                                                                                                               >  2:229650/1‑150 (MQ=255)
       tggtTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGc                                                                                                                                               >  2:1278465/1‑150 (MQ=255)
        ggtTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCTGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCGATATTTTATGTTTTTGTTGAGAATGCAGAGCCGTTAATTTTTCCATATTACGGcc                                                                                                                                              >  2:3145215/1‑150 (MQ=255)
        ggtTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCAAATAACGGcc                                                                                                                                              >  2:1890825/1‑150 (MQ=255)
           ttACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTg                                                                                                                                           >  1:173645/1‑150 (MQ=255)
               tgctgGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACg                                                                                                                                       >  1:3758119/1‑150 (MQ=255)
               tgctgGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACg                                                                                                                                       >  2:1374495/1‑150 (MQ=255)
               tgctgGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCAGTACg                                                                                                                                       >  1:281366/1‑150 (MQ=255)
                gctgGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGc                                                                                                                                      <  2:3801184/150‑1 (MQ=255)
                  tgGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCt                                                                                                                                    >  1:4194908/1‑150 (MQ=255)
                    gtgtTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTACATATAAAGGCCCAGTACGCCAgg                                                                                                                                  >  1:797607/1‑150 (MQ=255)
                          gCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTACATATAACGGCCCTGTACGCCTGGACGGAt                                                                                                                            >  2:2946435/1‑150 (MQ=255)
                            tGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATaa                                                                                                                          >  2:3733623/1‑150 (MQ=255)
                               cTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGc                                                                                                                       <  2:1881942/150‑1 (MQ=255)
                                tGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCAAATAACGGCCCAGTACGCCTGGACGGATAAAGCg                                                                                                                      >  1:2508348/1‑150 (MQ=255)
                                tGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTACATATAACGGCCCTGTACGCCTGGACGGATAAAGCg                                                                                                                      >  2:3981903/1‑150 (MQ=255)
                                  gtgtGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGtt                                                                                                                    <  1:996185/150‑1 (MQ=255)
                                      gggCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTGCATATAACGGGCCTGTACGCCTGGACGGATAAACCGTTTaag                                                                                                                >  1:2492737/1‑148 (MQ=255)
                                       ggCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCCTATAACGGCCCGGTACGCCGGGACGGATAAAGCGGTTAcgc                                                                                                               >  1:2042550/1‑150 (MQ=255)
                                             tGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTGGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAAc                                                                                                         <  1:3804390/150‑1 (MQ=255)
                                               gTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACgg                                                                                                       <  2:1942708/150‑1 (MQ=255)
                                                tGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATAAAACGGCCCAGTACGCCAGGGCGGATAAAGCGTTTACGCGACAACGGa                                                                                                      >  1:2680608/1‑150 (MQ=255)
                                                tGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCAAATAAAGGCCCAGTACGCCAGGACGGCTAAAGCGGTTAAGCGACAACCGa                                                                                                      >  1:933010/1‑150 (MQ=255)
                                                tGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTACATATAACGGCCCTGTACGCCAGGACGGATAAAGCGTTTACGCGACCACGGa                                                                                                      >  1:699298/1‑150 (MQ=255)
                                                tGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTACAAATAAAGGCCCTGTACGCCCCGACGGATAAAGCGTTTACGCGACAACGGa                                                                                                      >  1:2190576/1‑150 (MQ=255)
                                                tGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTAAATATAAAGGCCCAGTACGCCAGGACGGCTAAAGCGTTTAAGCGACAACGGa                                                                                                      >  1:8342/1‑150 (MQ=255)
                                                 gCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGGTTACGCGACAACGGAt                                                                                                     >  2:2775546/1‑150 (MQ=255)
                                                   tGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATCACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGc                                                                                                   <  1:1802068/150‑1 (MQ=255)
                                                   tGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGc                                                                                                   <  1:874639/150‑1 (MQ=255)
                                                     tttCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTTATTTTTCCATATAAACGCCCAGTACGCCTGGACGGATAAAGCGTTTACGCGACCAAGGATTccc                                                                                                 >  1:1147781/1‑150 (MQ=255)
                                                     tttCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGccc                                                                                                 >  2:2044901/1‑150 (MQ=255)
                                                     tttCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTAAAAATAAACGCCCAGTAAGCCTCGAAGGAAAAAAAGTTTACACGAAAAAGGGTGccc                                                                                                 >  1:2263828/1‑150 (MQ=255)
                                                       tccattccaGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTTCAAAAAAAGGACCAGTTGGCCGTGACGGATAAAACGTGTAAACGACAACCGATTCCCGa                                                                                               >  1:1817404/1‑150 (MQ=255)
                                                       tccattccaGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCGGGACGGATAAATCGTTTACGCGACAACGGATGCCCGa                                                                                               >  1:2253865/1‑150 (MQ=255)
                                                       tccattccaGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTACATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGa                                                                                               >  2:2561300/1‑150 (MQ=255)
                                                       tccattccaGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTACATATAAAGGCCCAGTACGCCGGGACGGATAAAACGTTTAAGCGACAACCGATGcccca                                                                                               >  1:2540039/1‑148 (MQ=255)
                                                       tccattccaGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATAGTTTATCGATTTGTTGAGAAAGCAGAGCCGTTAATTTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCg                                                                                                >  2:3804030/1‑150 (MQ=255)
                                                        ccattccaGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGAt                                                                                              <  2:3655315/150‑1 (MQ=255)
                                                          attccaGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGGCGGATAAAGCGTTTACGCGACAACGGATGCCCGATGc                                                                                            >  2:198089/1‑150 (MQ=255)
                                                           ttccaGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTGCATATAAAGGCCCAGTACGCCTGGACGGATAAAGCGTTTACGCGACAACCGATGCCCGATGCg                                                                                           >  1:3048865/1‑150 (MQ=255)
                                                           ttccaGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAAAGGCCCGGTACGCCTGGACGGATAAAGCGTTTACGCCCCCACGGATGCCCGAGGCg                                                                                           >  1:88725/1‑150 (MQ=255)
                                                           ttccaGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTACATATAACGGCCCGGTACGCCGGGACGGGTAAAGCGTTTACGCGACAACGGATTCCCCCTGCg                                                                                           >  1:1577989/1‑150 (MQ=255)
                                                            tccaGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGa                                                                                          <  2:3422958/150‑1 (MQ=255)
                                                             ccaGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGAc                                                                                         >  2:640670/1‑150 (MQ=255)
                                                             ccaGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAAAGGCCCTGTACACCTGGACGGATAAATCGTTTAAGCGACAACGGATGCCCGATGCGAc                                                                                         >  1:2386290/1‑150 (MQ=255)
                                                                     cggcggGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCAAATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGgcgc                                                                                 >  2:762910/1‑150 (MQ=255)
                                                                     cggcggGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTACATATAACGGCCCATTACGCCAGGACGGAAAAAGAGTTTACGCGACAACGGATGCCCGATTCGACGCAGgcgc                                                                                 >  1:271345/1‑150 (MQ=255)
                                                                      ggcggGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGcgcg                                                                                >  2:3159569/1‑150 (MQ=255)
                                                                      ggcggGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAAAGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGcgcg                                                                                >  2:554860/1‑150 (MQ=255)
                                                                      ggcggGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTACAAATAAAGGCCCAGTAAGCCAGGACGGATAAAACGTTTAAGCGACAACCGATGCCCCATGCGACGCAGGcgcg                                                                                >  2:1815982/1‑150 (MQ=255)
                                                                                    cAACCTTGATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTa                                                                  <  1:3426607/150‑1 (MQ=255)
                                                                                       ccTTTCTCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACaa                                                               <  2:3703786/150‑1 (MQ=255)
                                                                                        cTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTACATAAAAAGGGCCAGTAAGACAGGTCGGATAAAAACTTTTAGAGACAACCGATGACCCATGCGAAGCTGGCGCGTCTTATAGGGCCTACaaa                                                              >  1:3370988/1‑150 (MQ=255)
                                                                                            aTCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGc                                                          <  2:3826842/150‑1 (MQ=255)
                                                                                            aTCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGc                                                          <  2:2630599/150‑1 (MQ=255)
                                                                                               ggCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGcccc                                                       <  2:10855/150‑1 (MQ=255)
                                                                                                    gcCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCAAATAACGGCCCAGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCAACAAACGGCCCCGAATc                                                  >  2:1354752/1‑150 (MQ=255)
                                                                                                    gcCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTACATATAACGGCCCTGTACGCCAGGACGGATAAAGCGTTTACGCGACAAAGGAAGCCCGATGCGACGCATGCGCCTCTTATAGGGCCAACAAACCGCCCCGaaac                                                  >  2:4050875/1‑148 (MQ=255)
                                                                                                      cATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCAGTAAGCCAGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGt                                                >  2:918139/1‑150 (MQ=255)
                                                                                                      cATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAAATTTTACAAAAAACGGCCCAGTACGCCAGGACCGATAAATCGTTTACGCGACAACGGATGCCCCATTCGACCCTGTCGCGTCATATAGGGCCTACAAACAGCCCCCAAACGt                                                >  2:469706/1‑150 (MQ=255)
                                                                                                          ttttATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGt                                            <  2:3953943/150‑1 (MQ=255)
                                                                                                               tCTATTTGTTGAGAAAGCAGAGCCTTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGAt                                       <  1:3686086/150‑1 (MQ=255)
                                                                                                               tCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTTCGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGAt                                       <  2:692960/150‑1 (MQ=255)
                                                                                                                 tATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCAGTACGCCGGGCCGGATAAAGCCGTTAAGCGACAACGGATGCCCGATGCGCCGCTGGCGCGTCTTATCGGGCCAACAAAAGGCCCCCAAACGTAGGTCGGATaa                                     >  2:2504192/1‑150 (MQ=255)
                                                                                                                  aTTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAg                                    <  2:2416653/150‑1 (MQ=255)
                                                                                                                  aTTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAg                                    <  1:1045098/150‑1 (MQ=255)
                                                                                                                  aTTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAg                                    >  2:2332280/1‑150 (MQ=255)
                                                                                                                  aTTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAAAGGCCCAGTAAGCCTGGACGGATAAAGCGTTTAAGCGACAACCGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAAAGGCCCCGAATCGTAGGTCGGATAAg                                    >  2:3224861/1‑150 (MQ=39)
                                                                                                                  aTTTGTTGAGAAAGCAGAGCCGTTAAATTTTACAAATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCAGGCGCGTCTTATAGGGCCGACAAACGGCCCCGAAACGTAGGTCCGATAAg                                    >  2:896714/1‑150 (MQ=255)
                                                                                                                      gttgAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGt                                <  1:404914/150‑1 (MQ=255)
                                                                                                                         gagaAAGCAGAGCCGTTAATTTTTCCATATCACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCa                             <  2:1953404/150‑1 (MQ=255)
                                                                                                                          agaAAGCAGAGCCGTTAATTTTTACATATAACCGCCCTGTACCCCTGGCCGGATAAAGCGTTTACGCGACCACGGATGCCCGATGCCACGCTGGCGCGTCGTATCGGGCCTACCAACGGCCCCGAAACGTAGGTCGGATAAGGCCTTCAc                            >  2:3602856/1‑150 (MQ=255)
                                                                                                                             aaccAGAGCCTTTAATTTTACAAAAAAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCTACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCAcgcc                         <  1:1871954/147‑1 (MQ=38)
                                                                                                                               gCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATGCGGCGTGAAcgcctt                       <  1:2044901/150‑3 (MQ=255)
                                                                                                                               ccAGACCCGTTAATTTATCCAAATAACGGCCCTGTACGCCTGGACGCATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATGCGGCGTGAAcgcctt                       <  2:2263828/149‑3 (MQ=38)
                                                                                                                                  gagCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATc                    <  1:2917712/150‑1 (MQ=255)
                                                                                                                                   agCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATcc                   <  1:1325872/150‑1 (MQ=255)
                                                                                                                                     ccGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGa                 <  2:3040771/150‑1 (MQ=255)
                                                                                                                                     ccGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGa                 <  2:464899/150‑1 (MQ=255)
                                                                                                                                      cGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGAc                <  2:3869960/150‑1 (MQ=255)
                                                                                                                                       gTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACACCGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACa               <  2:4138610/150‑1 (MQ=255)
                                                                                                                                        ttAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACaa              >  2:1453760/1‑150 (MQ=255)
                                                                                                                                        ttAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACaa              >  2:1307857/1‑150 (MQ=255)
                                                                                                                                        ttAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACaa              >  1:257691/1‑150 (MQ=255)
                                                                                                                                         tAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACaaa             <  2:2716194/150‑1 (MQ=255)
                                                                                                                                          aaTTTTTCCATATAACGGCCCTGTACGCCTGGACGTATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCTTCTTATCGGGCCTAAAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAAc            <  1:3625495/150‑1 (MQ=255)
                                                                                                                                                    atatAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAAGGGTTGCCCGTTGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGa   <  1:554860/150‑1 (MQ=255)
                                                                                                                                                    atatAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGa   >  1:3900272/1‑150 (MQ=255)
                                                                                                                                                    atatAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGa   >  1:314975/1‑150 (MQ=255)
                                                                                                                                                     tataACGGCCCTGTACGCCTGGTCGGATAAAGCGTTTACGCGACCACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGAt  <  2:2466128/150‑1 (MQ=255)
                                                                                                                                                     tataACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGAt  >  2:211289/1‑150 (MQ=255)
                                                                                                                                                     tataACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGAt  >  2:296835/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGAT  >  NC_000913/173305‑173602

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: