Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 624,584 G→T 5.2% intergenic (‑74/‑301) fepB ← / → entC ferrienterobactin ABC transporter periplasmic binding protein/isochorismate synthase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913624,5840GT 5.2% 756.1 / 22.3 211intergenic (‑74/‑301)fepB/entCferrienterobactin ABC transporter periplasmic binding protein/isochorismate synthase 1
Reads supporting (aligned to +/- strand):  ref base G (101/98);  new base T (4/7);  total (105/106)
Fisher's exact test for biased strand distribution p-value = 5.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.68e-01

Only 100 of 215 total aligned reads displayed.

CCTGAACTGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTC  >  NC_000913/624435‑624732
                                                                                                                                                     |                                                                                                                                                    
ccTGAACTGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCg                                                                                                                                                      >  1:1299852/1‑150 (MQ=255)
 cTGAACTGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATTACCAATTCGACAAAGCGCACAATCCGTCCCCTCGc                                                                                                                                                     >  1:3154975/1‑150 (MQ=255)
 cTGAACTGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGc                                                                                                                                                     <  1:2932843/150‑1 (MQ=255)
 cTGAACTGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGc                                                                                                                                                     <  1:2932684/150‑1 (MQ=255)
 cTGAACTGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGc                                                                                                                                                     <  1:2842431/150‑1 (MQ=255)
 cTGAACTGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGAAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGc                                                                                                                                                     >  2:3234489/1‑150 (MQ=255)
    aaCTGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGTCGAGTCTCACAAATCAGCTTCCTGGTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGcccc                                                                                                                                                  >  2:3904492/1‑150 (MQ=255)
    aaCTGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGcccc                                                                                                                                                  >  1:3100160/1‑150 (MQ=255)
       tGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCttt                                                                                                                                               >  1:1659840/1‑150 (MQ=255)
        gCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTg                                                                                                                                              >  2:1274024/1‑150 (MQ=255)
         cgtcTATTCCTGAAAGCAACAGTCCTGAGAATAGAAGGGCGGTGCGGTAGAGCGTGGCGAGTCTCCCAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTgg                                                                                                                                             <  1:524819/147‑1 (MQ=255)
         cGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTgg                                                                                                                                             <  2:4113235/150‑1 (MQ=255)
          ggCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTggg                                                                                                                                            <  1:407462/150‑1 (MQ=255)
          ggCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTggg                                                                                                                                            <  1:4224858/150‑1 (MQ=255)
             tATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGgag                                                                                                                                         <  1:3528023/150‑1 (MQ=255)
              aTTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGgaga                                                                                                                                        <  1:2244446/150‑1 (MQ=255)
               ttCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGagag                                                                                                                                       <  2:3276653/150‑1 (MQ=255)
                 ccTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCTCCCCAAAGGGGCGAggg                                                                                                                                     <  2:1934268/150‑1 (MQ=255)
                   tGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCTCCCCAAAGGGGCGAgggga                                                                                                                                   <  2:2858290/150‑3 (MQ=255)
                   tGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCTCCCCAAAGGGGCGAgggga                                                                                                                                   <  2:693245/150‑3 (MQ=255)
                   tGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCTCCCCAAAGGGGCGAgggga                                                                                                                                   <  2:668204/150‑3 (MQ=255)
                   tGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGtt                                                                                                                                   >  1:1123595/1‑150 (MQ=255)
                   tGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGtt                                                                                                                                   >  2:2017690/1‑150 (MQ=255)
                      aaGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAgg                                                                                                                                <  1:394186/150‑1 (MQ=255)
                      aaGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAgg                                                                                                                                <  1:1295899/150‑1 (MQ=255)
                         cAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTg                                                                                                                             >  1:902742/1‑150 (MQ=255)
                          aaaaGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGa                                                                                                                            >  1:918621/1‑150 (MQ=255)
                           aaaGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAg                                                                                                                           <  2:3358775/150‑1 (MQ=255)
                                  tGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAAca                                                                                                                    >  2:2780280/1‑150 (MQ=255)
                                    ttAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAAcagc                                                                                                                  >  2:2450063/1‑150 (MQ=255)
                                    ttAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAAcagc                                                                                                                  >  1:1077960/1‑150 (MQ=255)
                                      aaTAGATGGGCGTTGCGGTAGAGGGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGTTTATGGGGTTAATGACAATTTCGAAAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGTTGTGGGTTTGGGTTTGGGGAAcagcca                                                                                                                <  2:3059509/150‑1 (MQ=255)
                                      aaTAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAAcagcca                                                                                                                >  2:3334357/1‑150 (MQ=255)
                                      aaTAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAAcagcca                                                                                                                >  1:26253/1‑150 (MQ=255)
                                       aTAGAGGGGCGTGGCGGTAGAGCGGGGCGATTCTCACAAAGCAGCTTCCTGTTATTAATAAGTTTACGGGGTTATTGACAAATTCGACAAAGCGCAAATTCCGTCCCCTCGCCCGTTTGGGGTGTGGGTTTGGGTGTGGGGAAcagccag                                                                                                               <  2:1384891/150‑1 (MQ=255)
                                          gacggggGTTGCGGTAGTGCGGGGCGTGTCTCACAAATCAGCTTCCTTTTATTAATATGGTTATGGGCTTAATTACAAATTCTAAAACGGGCACAATCCGTCCCCTCGCCCCTTTGGGGTGTGGTTTTGGGTGTGGGGAACAGCCAGCAc                                                                                                            <  2:211532/147‑1 (MQ=255)
                                          gacgggCGGTGCGGTAGAGCGGGGCGAGTCACACAAATCCGCTACCTGTTATTAATACGGTTAAGGGGGTATTGACAAATTCGACAACGCGCACAAACCGTCCCCGCGCCCCTTTGGGGAGCGGGTTTGGGTTAGGGTAACAGCCAGCAc                                                                                                            <  2:219901/147‑1 (MQ=255)
                                          gAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTGTGGGGAGAGGGTGCGGGTGAGGGGAACAGCCAGCAc                                                                                                            >  2:2822744/1‑150 (MQ=255)
                                            aGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTg                                                                                                          >  2:3460152/1‑150 (MQ=255)
                                               gCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTg                                                                                                       >  2:2431047/1‑150 (MQ=255)
                                               gCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTg                                                                                                       >  1:4288578/1‑150 (MQ=255)
                                               gCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTg                                                                                                       >  2:1567807/1‑150 (MQ=255)
                                                   tGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGaa                                                                                                   >  1:1792000/1‑150 (MQ=255)
                                                      ggTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACAt                                                                                                >  2:1345060/1‑150 (MQ=255)
                                                        tAGAGCGGGGCGAGTCCCCCAAAGCAGCTTCCTGTTATTAATACGGTTAAGGGCTTAATGCCAAATTCGACAACGCGCACACTCCGTCCCCCCGCCCCTTTGGGGAGTGGGTTAGGGTGCGGGTAACAGCCAGCACTGGTGCGAACATTa                                                                                              <  2:1239669/150‑1 (MQ=255)
                                                               gggCGAGTCTCACAAATCAGCTTCCTTTTATTAATAATGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCTTCCCCTCGCCCCTTTGGGGAGAGGGTTTGGGTGTGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCa                                                                                       <  2:528417/150‑1 (MQ=38)
                                                                    cgTCACAAAAATCCGCTGCCTTTTATTAATTCTGTTACGGGGGTAATCACAAATTCGACAACGCGCACACGCCGTCCCCTCGCCCCTTTGGGGTGTGGGTTTGGGTGTGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCa                                                                                  <  1:4213242/149‑1 (MQ=255)
                                                                    aGTCTCACAAATCCGCTTCCTGTTATTAATATGGTTAAGGGCGTTATGAGAAATTCGACAAAGCGCACAATCCGTGCCCTCGCGCCTTTGGGGTGTGGGTTAGGGTGTGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCa                                                                                  <  2:2354981/150‑1 (MQ=255)
                                                                       ctcACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGTTGAGGGTTAGGGTGTGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGcc                                                                               <  2:493720/150‑1 (MQ=255)
                                                                        tcACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGccc                                                                              >  1:849878/1‑150 (MQ=255)
                                                                         cacaAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCt                                                                             >  1:3832817/1‑150 (MQ=255)
                                                                         cacaAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCt                                                                             >  1:3132825/1‑150 (MQ=255)
                                                                          acaAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTc                                                                            >  2:4179522/1‑150 (MQ=255)
                                                                           aaaaTCCGCTTCCTTTTTTTAATATGTTTATGGGCTTAATCAAAATTTCTAAAAAGCGAAAAATCCGCCCCCCCGCCCCTTGGGGGTTTGGTTTTGGGTGTGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCa                                                                           <  1:1017774/149‑1 (MQ=255)
                                                                             aaTCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCAcc                                                                         >  1:1709312/1‑150 (MQ=255)
                                                                               tCAGATTCCTTTTATTAATACGGTTATGGGCGTAATGACAAATTCGACAAAGCCCAAAAACCGTCCCCTCGCCCCTTTGGGTAGTGGGTTTGGGTGTGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCt                                                                       <  1:1778012/150‑1 (MQ=255)
                                                                                   cTTCGTGTTATTAATTTGTTTATGGGCTTAATCAAAAATACTACAACGCGCACAATCCGCCCCCCCGCCCGTTGGGGTTGTGGTTTTGGGTGTGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGaa                                                                   <  1:966571/150‑1 (MQ=255)
                                                                                   cTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGaa                                                                   >  1:3285587/1‑150 (MQ=255)
                                                                                   cTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGaa                                                                   >  1:1836859/1‑150 (MQ=255)
                                                                                   cTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGaa                                                                   >  2:664024/1‑150 (MQ=255)
                                                                                     tCCTGTTTTAAATATGTTTATGGGCGTACCCCAAAATCCGACACAGCCCACAACCCGCCCCCCCCCCCGTTTGGGTTGGGGGTTGGGGTGCGGGCAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAgg                                                                 <  1:3331129/150‑1 (MQ=255)
                                                                                     tCCTGTTATTAATATGGTTATGGGCGTAATGACAAATGCGACAACGCGCACAATCCGTCCCCTCGCCCCTTTGGGGTGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAgg                                                                 <  2:3474336/150‑1 (MQ=37)
                                                                                       cTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGga                                                               >  2:2691330/1‑150 (MQ=255)
                                                                                        tGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGgag                                                              >  2:1289689/1‑150 (MQ=255)
                                                                                          ttattaATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGagag                                                            >  2:3715650/1‑150 (MQ=255)
                                                                                                atacGGTTACGGCCTTAATCACAAATCCGAAAAAGCGCAAAATCCGTCCCCGCGCCCCTTTGGGGAGTGGTTTTGGGTGTGGGCAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCa                                                      <  1:2328801/150‑1 (MQ=255)
                                                                                                      ttAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACgg                                                >  1:1933889/1‑150 (MQ=255)
                                                                                                         aGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGACCAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGgcgc                                             >  2:4070736/1‑150 (MQ=255)
                                                                                                         aGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGgcgc                                             >  1:1144261/1‑150 (MQ=255)
                                                                                                          gtgcttAATGACAATTTCGAAAAAGCTCACTTTCCGTCCCCTCGCGCATTTTGGGATTGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCAGCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCTGAACGGCGCa                                            <  2:662046/145‑1 (MQ=21)
                                                                                                          gggCTTAAGGACAAATGCGACAAAGCGCACAATCCGTCCCCGCGCCCCTTTGGGGTGTGGTTTTGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCa                                            <  2:2057302/150‑1 (MQ=255)
                                                                                                           ggCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAg                                           >  2:365341/1‑150 (MQ=255)
                                                                                                           ggCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAg                                           >  2:4327328/1‑150 (MQ=255)
                                                                                                               tAATTACAATTTATACAAAGCGCACAATCCTTCCCCCCGCCCCTTTGGGTATTGGGTTTGGGGGTGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACa                                       <  1:398548/150‑1 (MQ=255)
                                                                                                               tAATTAAAAATTCGACAACGCGCACAAGCCGTCCCCCCGCCCCTTGGGGGAGTGGGTTCGGGTGCGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACa                                       <  1:3110604/150‑1 (MQ=255)
                                                                                                                    aCAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATcaca                                  >  1:4110056/1‑150 (MQ=255)
                                                                                                                    aCAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATcaca                                  >  1:1505791/1‑150 (MQ=255)
                                                                                                                    aCAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCGTTGGGGAGAGGGTTAGGGTGAGGGGACCAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATcaca                                  >  2:2977584/1‑150 (MQ=255)
                                                                                                                     cAAATTCTACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGTAGAGGTTTTGGGTGTGGGTAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACAt                                 <  1:3883770/150‑1 (MQ=255)
                                                                                                                     cAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACAt                                 >  1:2370354/1‑150 (MQ=255)
                                                                                                                      acaTTCGCCACACCGCACAATCCGTCCCCTGGCCCCTTTGGGGCGTGGGTTCGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACAtt                                <  2:941972/148‑1 (MQ=255)
                                                                                                                      aaaTTCGACAAAGCGCACAATCCGACCCCTCGCCCGTTTGGGTTGTGGGTTTGGGTGTGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACAtt                                <  1:3014816/150‑1 (MQ=255)
                                                                                                                        aTTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTgc                              >  1:4330390/1‑150 (MQ=255)
                                                                                                                          tCGACAAAGCGCACAAGCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTTGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTgcgc                            <  2:3799012/150‑1 (MQ=255)
                                                                                                                            gACAACGCGCACAAGCCGTCCCCCCGCCCCTTTGGGGAGAGGGTTTGGGTGCGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCtt                          <  2:902742/150‑1 (MQ=255)
                                                                                                                             aCAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTa                         >  1:3303033/1‑150 (MQ=255)
                                                                                                                                caGCGCACACTCGGGCCCCTCGCCCCGTGGGGGAGTGGGTTAGGGTGCGGGGCACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGc                      <  2:2593653/149‑1 (MQ=255)
                                                                                                                                acgcgcACAACCCGGCCCCGCGCCCCTTGGGGGAGCGGGTTCGGGTGCGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGc                      <  2:1317062/148‑1 (MQ=255)
                                                                                                                                aaGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGc                      <  2:4253561/150‑1 (MQ=255)
                                                                                                                                   cgcACAATCCGTCCCCTCGCCCCTTTGGGGAGTGGGTTTGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGaa                   <  1:450762/150‑1 (MQ=255)
                                                                                                                                   cgaacaAGCCGTCCCCTCGCCCTTTTGGGGTTTGGGTTTGGGTGTGGGCAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGaa                   <  1:3107048/147‑1 (MQ=255)
                                                                                                                                    gcACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAAt                  >  1:563099/1‑150 (MQ=255)
                                                                                                                                     cacaATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATc                 <  2:3591040/150‑1 (MQ=255)
                                                                                                                                     cacaATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATc                 <  2:634866/150‑1 (MQ=255)
                                                                                                                                     cacaATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATc                 >  1:1162230/1‑150 (MQ=255)
                                                                                                                                           cctaacccTCTCCCCAAAGGGGCGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGGAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCaa           >  1:2654062/6‑150 (MQ=37)
                                                                                                                                             gtcgcgtCGCCCTTTGGGGGTGTGGGTTTGGTTTAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGTGAGAGGGGTCAGAACGGCGCAGTACATCACATTGCGCTTATGCGAATCCATCaata         <  1:2375994/144‑1 (MQ=255)
                                                                                                                                             gTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCaata         <  1:2873067/150‑1 (MQ=255)
                                                                                                                                                cccTCGCCCCTTTGGGGAGTGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATg      <  2:3132825/150‑1 (MQ=255)
                                                                                                                                                    ggCCCGTGTGGGGTGTGGGTTTGGGGGCGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTtc  <  2:3699197/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CCTGAACTGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTC  >  NC_000913/624435‑624732

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: