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.

TGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAA‑TTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGAT  >  NC_000913/173305‑173602
                                                                                                                                                     |                                                                                                                                                     
tggtggGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAA‑TTTTTCCata                                                                                                                                                      >  2:1673619/1‑150 (MQ=255)
tggtggGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAA‑TTTTTCCata                                                                                                                                                      >  1:214451/1‑150 (MQ=255)
 ggtgggtgGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAA‑TTTTTCCatat                                                                                                                                                     >  1:3027914/1‑150 (MQ=255)
 ggtgggtgGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAA‑TTTTTCCatat                                                                                                                                                     >  1:3322136/1‑150 (MQ=255)
    gggtggTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAA‑TTTTTCCATATAAc                                                                                                                                                  >  2:1626137/1‑150 (MQ=255)
    gggtggTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAA‑TTTTTACATATAAc                                                                                                                                                  >  2:2869643/1‑150 (MQ=255)
     ggtggtTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTGA‑TTTTTGCATATaaag                                                                                                                                                 >  1:2531516/1‑148 (MQ=255)
     ggtggtTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTAGTTGTTGAGAAAGCAGAGCCGTTAA‑ATTTTAAATATaaag                                                                                                                                                 >  1:483966/1‑148 (MQ=255)
       tggtTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAA‑TTTTTCCATATAACGGc                                                                                                                                               >  2:229650/1‑150 (MQ=255)
       tggtTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAA‑TTTTTCCATATAACGGc                                                                                                                                               >  2:1278465/1‑150 (MQ=255)
        ggtTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCTGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCGATATTTTATGTTTTTGTTGAGAATGCAGAGCCGTTAA‑TTTTTCCATATTACGGcc                                                                                                                                              >  2:3145215/1‑150 (MQ=255)
        ggtTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAA‑TTTTTCCAAATAACGGcc                                                                                                                                              >  2:1890825/1‑150 (MQ=255)
           ttACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAA‑TTTTTCCATATAACGGCCCTg                                                                                                                                           >  1:173645/1‑150 (MQ=255)
               tgctgGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAA‑TTTTTCCATATAACGGCCCTGTACg                                                                                                                                       >  1:3758119/1‑150 (MQ=255)
               tgctgGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAA‑TTTTTCCATATAACGGCCCTGTACg                                                                                                                                       >  2:1374495/1‑150 (MQ=255)
               tgctgGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAA‑TTTTTCCATATAACGGCCCAGTACg                                                                                                                                       >  1:281366/1‑150 (MQ=255)
                gctgGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAA‑TTTTTCCATATAACGGCCCTGTACGc                                                                                                                                      <  2:3801184/150‑1 (MQ=255)
                  tgGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAA‑TTTTTCCATATAACGGCCCTGTACGCCt                                                                                                                                    >  1:4194908/1‑150 (MQ=255)
                    gtgtTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAA‑TTTTTACATATAAAGGCCCAGTACGCCAgg                                                                                                                                  >  1:797607/1‑150 (MQ=255)
                          gCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAA‑TTTTTACATATAACGGCCCTGTACGCCTGGACGGAt                                                                                                                            >  2:2946435/1‑150 (MQ=255)
                            tGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAA‑TTTTTCCATATAACGGCCCTGTACGCCTGGACGGATaa                                                                                                                          >  2:3733623/1‑150 (MQ=255)
                               cTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAA‑TTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGc                                                                                                                       <  2:1881942/150‑1 (MQ=255)
                                tGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAA‑TTTTTCCAAATAACGGCCCAGTACGCCTGGACGGATAAAGCg                                                                                                                      >  1:2508348/1‑150 (MQ=255)
                                tGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAA‑TTTTTACATATAACGGCCCTGTACGCCTGGACGGATAAAGCg                                                                                                                      >  2:3981903/1‑150 (MQ=255)
                                  gtgtGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAA‑TTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGtt                                                                                                                    <  1:996185/150‑1 (MQ=255)
                                      gggCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAA‑TTTTTGCATATAACGGGCCTGTACGCCTGGACGGATAAACCGTTTaag                                                                                                                >  1:2492737/1‑148 (MQ=255)
                                       ggCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAA‑TTTTTCCCTATAACGGCCCGGTACGCCGGGACGGATAAAGCGGTTAcgc                                                                                                               >  1:2042550/1‑150 (MQ=255)
                                             tGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTGGAGAAAGCAGAGCCGTTAA‑TTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAAc                                                                                                         <  1:3804390/150‑1 (MQ=255)
                                               gTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAA‑TTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACgg                                                                                                       <  2:1942708/150‑1 (MQ=255)
                                                tGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAA‑TTTTTCCATAAAACGGCCCAGTACGCCAGGGCGGATAAAGCGTTTACGCGACAACGGa                                                                                                      >  1:2680608/1‑150 (MQ=255)
                                                tGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAA‑TTTTTCCAAATAAAGGCCCAGTACGCCAGGACGGCTAAAGCGGTTAAGCGACAACCGa                                                                                                      >  1:933010/1‑150 (MQ=255)
                                                tGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAA‑TTTTTACATATAACGGCCCTGTACGCCAGGACGGATAAAGCGTTTACGCGACCACGGa                                                                                                      >  1:699298/1‑150 (MQ=255)
                                                tGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAA‑TTTTTACAAATAAAGGCCCTGTACGCCCCGACGGATAAAGCGTTTACGCGACAACGGa                                                                                                      >  1:2190576/1‑150 (MQ=255)
                                                tGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAA‑TTTTTAAATATAAAGGCCCAGTACGCCAGGACGGCTAAAGCGTTTAAGCGACAACGGa                                                                                                      >  1:8342/1‑150 (MQ=255)
                                                 gCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAA‑TTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGGTTACGCGACAACGGAt                                                                                                     >  2:2775546/1‑150 (MQ=255)
                                                   tGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAA‑TTTTTCCATATCACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGc                                                                                                   <  1:1802068/150‑1 (MQ=255)
                                                   tGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAA‑TTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGc                                                                                                   <  1:874639/150‑1 (MQ=255)
                                                     tttCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTTA‑TTTTTCCATATAAACGCCCAGTACGCCTGGACGGATAAAGCGTTTACGCGACCAAGGATTccc                                                                                                 >  1:1147781/1‑150 (MQ=255)
                                                     tttCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAA‑TTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGccc                                                                                                 >  2:2044901/1‑150 (MQ=255)
                                                     tttCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAA‑TTTTTAAAAATAAACGCCCAGTAAGCCTCGAAGGAAAAAAAGTTTACACGAAAAAGGGTGccc                                                                                                 >  1:2263828/1‑150 (MQ=255)
                                                       tccattccaGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAA‑TTTTTTCAAAAAAAGGACCAGTTGGCCGTGACGGATAAAACGTGTAAACGACAACCGATTCCCGa                                                                                               >  1:1817404/1‑150 (MQ=255)
                                                       tccattccaGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAA‑TTTTTCCATATAACGGCCCTGTACGCCGGGACGGATAAATCGTTTACGCGACAACGGATGCCCGa                                                                                               >  1:2253865/1‑150 (MQ=255)
                                                       tccattccaGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAA‑TTTTTACATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGa                                                                                               >  2:2561300/1‑150 (MQ=255)
                                                       tccattccaGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAA‑TTTTTACATATAAAGGCCCAGTACGCCGGGACGGATAAAACGTTTAAGCGACAACCGATGcccca                                                                                               >  1:2540039/1‑148 (MQ=255)
                                                       tccattccaGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATAGTTTATCGATTTGTTGAGAAAGCAGAGCCGTTAATTTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCg                                                                                                >  2:3804030/1‑150 (MQ=255)
                                                        ccattccaGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAA‑TTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGAt                                                                                              <  2:3655315/150‑1 (MQ=255)
                                                          attccaGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAA‑TTTTTCCATATAACGGCCCTGTACGCCTGGGCGGATAAAGCGTTTACGCGACAACGGATGCCCGATGc                                                                                            >  2:198089/1‑150 (MQ=255)
                                                           ttccaGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAA‑TTTTTGCATATAAAGGCCCAGTACGCCTGGACGGATAAAGCGTTTACGCGACAACCGATGCCCGATGCg                                                                                           >  1:3048865/1‑150 (MQ=255)
                                                           ttccaGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAA‑TTTTTCCATATAAAGGCCCGGTACGCCTGGACGGATAAAGCGTTTACGCCCCCACGGATGCCCGAGGCg                                                                                           >  1:88725/1‑150 (MQ=255)
                                                           ttccaGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAA‑TTTTTACATATAACGGCCCGGTACGCCGGGACGGGTAAAGCGTTTACGCGACAACGGATTCCCCCTGCg                                                                                           >  1:1577989/1‑150 (MQ=255)
                                                            tccaGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAA‑TTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGa                                                                                          <  2:3422958/150‑1 (MQ=255)
                                                             ccaGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAA‑TTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGAc                                                                                         >  2:640670/1‑150 (MQ=255)
                                                             ccaGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAA‑TTTTTCCATATAAAGGCCCTGTACACCTGGACGGATAAATCGTTTAAGCGACAACGGATGCCCGATGCGAc                                                                                         >  1:2386290/1‑150 (MQ=255)
                                                                     cggcggGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAA‑TTTTTCCAAATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGgcgc                                                                                 >  2:762910/1‑150 (MQ=255)
                                                                     cggcggGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAA‑TTTTTACATATAACGGCCCATTACGCCAGGACGGAAAAAGAGTTTACGCGACAACGGATGCCCGATTCGACGCAGgcgc                                                                                 >  1:271345/1‑150 (MQ=255)
                                                                      ggcggGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAA‑TTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGcgcg                                                                                >  2:3159569/1‑150 (MQ=255)
                                                                      ggcggGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAA‑TTTTTCCATATAAAGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGcgcg                                                                                >  2:554860/1‑150 (MQ=255)
                                                                      ggcggGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAA‑TTTTTACAAATAAAGGCCCAGTAAGCCAGGACGGATAAAACGTTTAAGCGACAACCGATGCCCCATGCGACGCAGGcgcg                                                                                >  2:1815982/1‑150 (MQ=255)
                                                                                    cAACCTTGATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAA‑TTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTa                                                                  <  1:3426607/150‑1 (MQ=255)
                                                                                       ccTTTCTCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAA‑TTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACaa                                                               <  2:3703786/150‑1 (MQ=255)
                                                                                        cTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAA‑TTTTTACATAAAAAGGGCCAGTAAGACAGGTCGGATAAAAACTTTTAGAGACAACCGATGACCCATGCGAAGCTGGCGCGTCTTATAGGGCCTACaaa                                                              >  1:3370988/1‑150 (MQ=255)
                                                                                            aTCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAA‑TTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGc                                                          <  2:3826842/150‑1 (MQ=255)
                                                                                            aTCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAA‑TTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGc                                                          <  2:2630599/150‑1 (MQ=255)
                                                                                               ggCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAA‑TTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGcccc                                                       <  2:10855/150‑1 (MQ=255)
                                                                                                    gcCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAA‑TTTTTCCAAATAACGGCCCAGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCAACAAACGGCCCCGAATc                                                  >  2:1354752/1‑150 (MQ=255)
                                                                                                    gcCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAA‑TTTTTACATATAACGGCCCTGTACGCCAGGACGGATAAAGCGTTTACGCGACAAAGGAAGCCCGATGCGACGCATGCGCCTCTTATAGGGCCAACAAACCGCCCCGaaac                                                  >  2:4050875/1‑148 (MQ=255)
                                                                                                      cATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAA‑TTTTTCCATATAACGGCCCAGTAAGCCAGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGt                                                >  2:918139/1‑150 (MQ=255)
                                                                                                      cATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAA‑ATTTTACAAAAAACGGCCCAGTACGCCAGGACCGATAAATCGTTTACGCGACAACGGATGCCCCATTCGACCCTGTCGCGTCATATAGGGCCTACAAACAGCCCCCAAACGt                                                >  2:469706/1‑150 (MQ=255)
                                                                                                          ttttATCTATTTGTTGAGAAAGCAGAGCCGTTAA‑TTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGt                                            <  2:3953943/150‑1 (MQ=255)
                                                                                                               tCTATTTGTTGAGAAAGCAGAGCCTTTAA‑TTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGAt                                       <  1:3686086/150‑1 (MQ=255)
                                                                                                               tCTATTTGTTGAGAAAGCAGAGCCGTTAA‑TTTTTCCATATAACGGCCCTGTTCGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGAt                                       <  2:692960/150‑1 (MQ=255)
                                                                                                                 tATTTGTTGAGAAAGCAGAGCCGTTAA‑TTTTTCCATATAACGGCCCAGTACGCCGGGCCGGATAAAGCCGTTAAGCGACAACGGATGCCCGATGCGCCGCTGGCGCGTCTTATCGGGCCAACAAAAGGCCCCCAAACGTAGGTCGGATaa                                     >  2:2504192/1‑150 (MQ=255)
                                                                                                                  aTTTGTTGAGAAAGCAGAGCCGTTAA‑TTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAg                                    <  2:2416653/150‑1 (MQ=255)
                                                                                                                  aTTTGTTGAGAAAGCAGAGCCGTTAA‑TTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAg                                    <  1:1045098/150‑1 (MQ=255)
                                                                                                                  aTTTGTTGAGAAAGCAGAGCCGTTAA‑TTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAg                                    >  2:2332280/1‑150 (MQ=255)
                                                                                                                  aTTTGTTGAGAAAGCAGAGCCGTTAA‑TTTTTCCATATAAAGGCCCAGTAAGCCTGGACGGATAAAGCGTTTAAGCGACAACCGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAAAGGCCCCGAATCGTAGGTCGGATAAg                                    >  2:3224861/1‑150 (MQ=39)
                                                                                                                  aTTTGTTGAGAAAGCAGAGCCGTTAA‑ATTTTACAAATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCAGGCGCGTCTTATAGGGCCGACAAACGGCCCCGAAACGTAGGTCCGATAAg                                    >  2:896714/1‑150 (MQ=255)
                                                                                                                      gttgAGAAAGCAGAGCCGTTAA‑TTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGt                                <  1:404914/150‑1 (MQ=255)
                                                                                                                         gagaAAGCAGAGCCGTTAA‑TTTTTCCATATCACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCa                             <  2:1953404/150‑1 (MQ=255)
                                                                                                                          agaAAGCAGAGCCGTTAA‑TTTTTACATATAACCGCCCTGTACCCCTGGCCGGATAAAGCGTTTACGCGACCACGGATGCCCGATGCCACGCTGGCGCGTCGTATCGGGCCTACCAACGGCCCCGAAACGTAGGTCGGATAAGGCCTTCAc                            >  2:3602856/1‑150 (MQ=255)
                                                                                                                             aaccAGAGCCTTTAA‑TTTTACAAAAAAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCTACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCAcgcc                         <  1:1871954/147‑1 (MQ=38)
                                                                                                                               gCAGAGCCGTTAA‑TTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATGCGGCGTGAAcgcctt                       <  1:2044901/150‑3 (MQ=255)
                                                                                                                               ccAGACCCGTTAA‑TTTATCCAAATAACGGCCCTGTACGCCTGGACGCATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATGCGGCGTGAAcgcctt                       <  2:2263828/149‑3 (MQ=38)
                                                                                                                                  gagCCGTTAA‑TTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATc                    <  1:2917712/150‑1 (MQ=255)
                                                                                                                                   agCCGTTAA‑TTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATcc                   <  1:1325872/150‑1 (MQ=255)
                                                                                                                                     ccGTTAA‑TTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGa                 <  2:3040771/150‑1 (MQ=255)
                                                                                                                                     ccGTTAA‑TTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGa                 <  2:464899/150‑1 (MQ=255)
                                                                                                                                      cGTTAA‑TTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGAc                <  2:3869960/150‑1 (MQ=255)
                                                                                                                                       gTTAA‑TTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACACCGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACa               <  2:4138610/150‑1 (MQ=255)
                                                                                                                                        ttAA‑TTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACaa              >  2:1453760/1‑150 (MQ=255)
                                                                                                                                        ttAA‑TTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACaa              >  2:1307857/1‑150 (MQ=255)
                                                                                                                                        ttAA‑TTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACaa              >  1:257691/1‑150 (MQ=255)
                                                                                                                                         tAA‑TTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACaaa             <  2:2716194/150‑1 (MQ=255)
                                                                                                                                          aa‑TTTTTCCATATAACGGCCCTGTACGCCTGGACGTATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCTTCTTATCGGGCCTAAAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAAc            <  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)
                                                                                                                                                     |                                                                                                                                                     
TGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAA‑TTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGAT  >  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: ‑