Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 624,598 A→T 28.6% intergenic (‑88/‑287) 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,5980AT28.6% 130.9 / 23.1 112intergenic (‑88/‑287)fepB/entCferrienterobactin ABC transporter periplasmic binding protein/isochorismate synthase 1
Reads supporting (aligned to +/- strand):  ref base A (36/36);  new base T (2/33);  total (38/74)
Fisher's exact test for biased strand distribution p-value = 3.05e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 123 total aligned reads displayed.

GCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATT  >  NC_000913/624459‑624737
                                                                                                                                           |                                                                                                                                           
gCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGgaga                                                                                                                                             <  2:524335/140‑1 (MQ=255)
gCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGgaga                                                                                                                                             <  1:1329228/140‑1 (MQ=255)
gCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGgaga                                                                                                                                             <  1:1329365/140‑1 (MQ=255)
     aGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGtt                                                                                                                                        <  1:1152321/140‑1 (MQ=255)
        ccTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAgg                                                                                                                                     <  2:1323532/140‑1 (MQ=255)
        ccTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAAATCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAgg                                                                                                                                     <  2:1081784/140‑1 (MQ=255)
            ttAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGTGGGTTAGGGTGa                                                                                                                                 >  2:165450/1‑140 (MQ=255)
            ttAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTgg                                                                                                                                 <  2:565372/140‑2 (MQ=255)
                ttgAAGGGCGTTGGGTTTTAGGGGGGCGAGTCTAACAAATCAGCTTCCTTTTATTAATATGTTTATGGGCTTAATGAAAATTTCGACAACGCGCACAATCCGTCCCCTCGCCCCTTTGGGTTGAGGTTTTGGGTGTgggg                                                                                                                             <  2:310155/138‑1 (MQ=255)
                tAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAgggg                                                                                                                             >  1:1441185/1‑140 (MQ=255)
                  gAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGTTTAAGGGCGTAATGAAAATTTCGAAAAAGCGAACAATCCGTCCCCTCGCCCCTTTGGGGTGTGGGTTAGGGTGAGGGGaa                                                                                                                           <  2:716954/140‑1 (MQ=38)
                     gggCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAAcag                                                                                                                        >  1:997485/1‑140 (MQ=255)
                      ggCTTTGCGGTAGTGCGGGGCGAGTCACACAAAACAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATTACAAATTCGAAAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGTGGTTTAGGGTGAGGGGAAcagc                                                                                                                       <  2:178348/140‑1 (MQ=38)
                      ggCGTTGCGGTTGAGCGGGGCTTGTCTCACAAATCAGGTTCCTGTTATTAATATGGTTAAGGGCTTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGTTGTGGTTTAGGGTGAGGGGAAcagc                                                                                                                       <  2:1454939/140‑1 (MQ=255)
                      ggCGTTGCGGTTGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGAAAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGTGTGGGTTAGGGTGAGGGGAAcagc                                                                                                                       <  2:1338889/140‑1 (MQ=255)
                      ggCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAAcagc                                                                                                                       >  2:1109461/1‑140 (MQ=255)
                        cGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAAcagcca                                                                                                                     >  1:772491/1‑140 (MQ=255)
                                tAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGt                                                                                                             >  1:540909/1‑140 (MQ=255)
                                   tgCGGGGCGATTCTCACAAATCAGCTTCCTTTTATTAATAAGGTAACGGGCGTAATGACAAATTCCACAAAGCCCACAACCCGCCCCCCCGCCCATTTGGGTTGTGGGTTTGGGGGAGGGCAACAGCCAGCACTGGTGCg                                                                                                          <  1:154092/139‑1 (MQ=255)
                                     gggggCGAGTCTAACAATTCAGCTTCCTGTTATTAATATGGTTCTGGGCGTAATGACAAATTCGACAAAGCGCACACTCCGTCCCCTCGCCCCTTTGGGGAGTGGGTTAGGGTGCGGGGAACAGCCAGCACTGGTGCGaa                                                                                                        <  2:34841/139‑1 (MQ=37)
                                     cGGGGCGCGTCTCACAACCCCGCTGCCGTTTCTGACTCTGGTTACGGGGGTACTGACCAATTCGCCAACGCGCACATTGCGTCCCCTCGCGCCTTGGGGGAGCGGGGGAGGGGGCGGGGACCAGCCAGCACTGGTGCGaa                                                                                                        <  2:34269/140‑1 (MQ=255)
                                     cGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGaa                                                                                                        >  1:498878/1‑140 (MQ=255)
                                        ggCGTGTCTCACAAATCCGCTTCCTGTTATTAATCAGGTTAAGGGCGTAATGACAAATACGACAACGCGCACAACCCGTCCCCCCGCCCCTTGGGGGAGCGGGTTAGGGGGCGGGGAACAGCCAGCACTGGTGCGAACAt                                                                                                     <  1:677773/140‑1 (MQ=255)
                                        gcctattctcAAAAAACAGCTTCATTTTATTAATACGTTTACGGGCTTAATCAAAAATCCCAAAACCCCCACAATCCGTCCCCTCGCCCCTTTGGGTATTGGGTTCGGGTGAGGGTAACAGCCAGCACTGGTGCGAACAt                                                                                                     <  1:530181/134‑1 (MQ=16)
                                         gCGTGTCTCACAAATCAGCTTCCTGTTATTAATATGGTTAAGGGCGTAAGGACAATTTCGACAAAGCGCACAATCCGTCCCCCCGCCCCTTTGGGGAGTGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACAtt                                                                                                    <  2:453747/140‑1 (MQ=38)
                                         gCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACAtt                                                                                                    <  2:810813/140‑1 (MQ=255)
                                            aGTCTCACAAATCAGCTTCCTGTTATTAATCAGGTTAAGGGCGTAATGAAAAAGTCGACAAAGCGCACAATCCGTCCCCACGCCCCTTTGGGGAGAGGGTTTGGGTGAGGGGAACAGCCAGCACTGGTTCGAACATTAAc                                                                                                 >  1:564511/1‑140 (MQ=38)
                                            aGTCTCACAAATCAGCTTCCTGTTATTAATATGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCTGCACTGGTGCGACCATTatc                                                                                                 >  1:1293344/1‑138 (MQ=255)
                                            aGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAAc                                                                                                 >  1:994449/1‑140 (MQ=255)
                                            aGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAAc                                                                                                 >  1:26593/1‑140 (MQ=255)
                                               ctcACAAATCAGCTTCCTGTTATTAATAAGGTTATGGGCGTAATGACAAATCCGACAACGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGTGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCt                                                                                              <  2:104180/140‑1 (MQ=255)
                                               ctcACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCATCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACACt                                                                                              <  1:773136/140‑1 (MQ=255)
                                               ataacaAATCAGCTTCCTTTTTTTAATAAGTTTACGGGCTTATTTAAAATTTCTAAAAAGCGCAAAATCCGTCCCCTCGCCCCTTTGGGGTGTGGTTTTGGGTGTGGGGAACAGCCAGCACTGGTGCGAACATTAACCCt                                                                                              <  1:730477/137‑1 (MQ=19)
                                               ataaaaaaTCTGTTTCCTTTTTTTAATTAGTTTTAGGGCTTAATCAAAATTTCCAAAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGTGTGGGTTTGGGTGTGGGGAACAGCCAGCACTGGTGCGAACATTAACCCt                                                                                              <  1:730135/135‑1 (MQ=19)
                                                    aaacccgCTTCCTTTTATTAATATGGTTACGGGGGTAATCAAAATGTCCACAACGCCCACAATCCGCCCCCCCGCCCCTTGGGGGAGCGGGTTGGGGTGTGGGCAACCGCCAGCACTGGTGCGAACATTAACCCTCAccc                                                                                         <  1:1003084/134‑1 (MQ=255)
                                                    aaaTCAGATTCCTTTTATTAATAAGGTTACGGGCGTATTGACAAATTCGACAACGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGTGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCAccc                                                                                         <  1:847670/140‑1 (MQ=37)
                                                    aaaTCACCTCCATTTTATTAATAAGTTTATGGCCTTAATCAAAATTACAACAAAGCGAAAAAACCGTCCCCCCGCCCCTTTGGGGTGTGGTTTTGGGTGCGGGCAACAGCCAGCACTGGTGCGAACATTAACCCTCAccc                                                                                         <  1:1109461/140‑1 (MQ=255)
                                                      aTCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCa                                                                                       >  2:849295/1‑140 (MQ=255)
                                                          gCTTCCTGTTATTAATAAGGTTAATGGCGTGATGACAAATTCGACAAAGCGCACAATACGTCCCCTCGCACCTTTGGGGTGAGGGTTAGGGTGAGGGGTACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCTGccc                                                                                   >  2:755966/1‑140 (MQ=38)
                                                           ttcccTTTTATTAATAAGTTTATGGGCTAAATCAAAAATTCGAAAAAGCGCAAAAACCGTCCCCGCCCCCCTTGGGGTTGTGGTTTTGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCt                                                                                  <  1:1072149/137‑1 (MQ=255)
                                                           gtcccTGTTATTAATATGGTTCTGGGCGTAATGCCAAATGCCCCAAAGCGCAAAATCCGCCCCCTCGCCCCTTGGGGGTGTGGGTTCGGGTGTGGGTAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCt                                                                                  <  1:1141446/137‑1 (MQ=255)
                                                           cTTCCTTTTTTTAAAACGGTTACGGGCTTAATCACAATTTCCACAACGCCCACAATCCGCCCCCCCGCCCCTTGGGGGTGTGGTTTCGGGTGTGGGCAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCt                                                                                  <  1:1095777/140‑1 (MQ=255)
                                                           cTTCCTGTTTTTAATAAGTTTATGGGCTTAATGAAAATTTCTAAAATGCGCACATTCCGTCCCCTCGCCCCTTTGGGGTGTGGGTTTGGGTGTGGGTAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCt                                                                                  <  1:1620690/140‑1 (MQ=21)
                                                           cTTCCTGTTATTAATATGGTTAAGGGCGTAATGACAATTTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGTGCGGTTTTGGGTGTGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCt                                                                                  <  2:1316892/140‑1 (MQ=37)
                                                           cTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCt                                                                                  >  1:313604/1‑140 (MQ=255)
                                                           cTTCCTGTTATGAATACGGTTAAGGGCGTATTGACAAATTCTACAAAGCGCACAATCCGTCCCCACGCCCCTTTGGGGTGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCt                                                                                  <  2:1109700/140‑1 (MQ=37)
                                                             tCCTTTTTTTAATATGTTTACGGGCTTAATGAAAATTTCCAAAAAGCCCAAAAACCGCCCCCCCGCCCCTTTGGGTTGTGGTTTTGGTTGCGGGCAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCa                                                                                <  1:969164/140‑1 (MQ=255)
                                                             tCCTGTAATTAAAAAGGTAAAGGGCTAAATGACAAATTCGACAAAGCCCACAAACCGACCCATCGCCCCTTTGGGGAGAGGGTTAGGGTGCGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCa                                                                                <  2:456722/140‑1 (MQ=25)
                                                             cccgttaattaCTACGTTTACGGGGTTAATCACAAATTCGACAAAGCGCAAAATCCGTCCCCTCCCCCCTTTGGGTAGTGGTTTCGGGTGCGGGCAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCa                                                                                <  1:1417845/133‑1 (MQ=255)
                                                               cTGTTTTTAATACGTTTCTGGGCGTAATGACAACTTCCAAAACGCCCAAAACCGGTCCCCTCCCCCCTTTGGGTTGTGGTTTGGGGTGGGGGCCACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCAcc                                                                              <  1:916950/140‑1 (MQ=255)
                                                                tGTTATTAATAAGGTTAAGGGCGAAATGACAAATTCGACAAAGCGCACAAACCGTCCCCTCGCCCCTTTGCGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTCGAGCGAACATTAACCCACACCCCAGCCCTCAccc                                                                             >  2:1455982/1‑140 (MQ=37)
                                                                  ttattaATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTg                                                                           >  2:152299/1‑140 (MQ=255)
                                                                    attaATCAGTTTATGGGCTTAATCAAAAATTCTAAAAAGCGCACAATCCGTCCCCCCGCCCCTTGGGGTTGTGGTTTTGGGTGTGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGa                                                                         <  1:1029762/140‑1 (MQ=255)
                                                                    attaATAAGGTTCAGGGCGTAATGCCAAATTCGACAAAGCGCACAATCCGTCCCCCCGCCCCTTTGGGGAGAGGGTTAGGGGGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGa                                                                         <  2:1606447/140‑1 (MQ=255)
                                                                    aataataaTGTTAAGGGCGTATAGAAAAAGTCCCCAAAGCGCACACTCCCACCCCTCGCCCCTTATGGGAGAGTGTTAGGTTGAGGGGAACACCCAGCAATTGTGCGAACAATAACACTCACCCCAGCCCTCACCCTGGa                                                                         <  2:105660/138‑1 (MQ=255)
                                                                     ttaATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGaa                                                                        >  2:328092/1‑140 (MQ=255)
                                                                     ttaATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGaa                                                                        <  2:213321/140‑1 (MQ=255)
                                                                        ataaGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAAACCGACCCCTCGCCCCTTAGGGGAGAGGGTTAGGGTGAGGCGAACAGCCAGCACAGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCAGGAAtgg                                                                     >  1:1684826/1‑140 (MQ=37)
                                                                         taaGGTGAAGGGCGTAATGACACATTCGCCAACGCGCACATTCCGTCCCCGCGCCCCTTTGGGGAGAGGGTTAGGGTGCGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGga                                                                    <  2:74504/140‑1 (MQ=25)
                                                                           cgtttATGGGCTTAATCAAAAATTCCAAAATGCGCAAAACCCGTCCCCTCCCCCTTTTGGGTTGTGGTTTTGGGTGTGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGgaga                                                                  <  1:1153969/137‑1 (MQ=255)
                                                                           cggTTACGGGCTTAATCACAAATACCACAAAGCGCACAATCCGCCCCCGCGCCCCTTTGGGGTCTGGGTTAGGGTGTGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGgaga                                                                  <  1:1298282/139‑1 (MQ=255)
                                                                           cggTTAAGGGCGTAACGACAAATACGACAAAGCGCACAACCCGTCCCCACCCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCACACCCCAGCCCTCACCCTGGAAGGgaga                                                                  <  2:686747/139‑1 (MQ=25)
                                                                           agtttACGGGCGTAAAGACAAATTCCCCAAAGCGCACAATCCGTCCCCGCGCCCATTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCCGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGgaga                                                                  <  2:486292/137‑1 (MQ=25)
                                                                           aGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACCCCCACCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGgaga                                                                  <  2:926055/140‑1 (MQ=255)
                                                                           aGGTAAAGGGCGTAATTACAAATTCGACAAAGCGCACAAGCCGTCCCCTCGCCCCTTTGGGGAGAGGTTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGgaga                                                                  <  2:928600/140‑1 (MQ=255)
                                                                            ggTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGTGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGagag                                                                 >  1:1464813/1‑140 (MQ=255)
                                                                            ggTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACCTTAACCCTCACCCCAGCCCTCACCCTGGAAGGGagag                                                                 >  1:1464761/1‑140 (MQ=255)
                                                                            ggTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGagag                                                                 >  2:443109/1‑140 (MQ=255)
                                                                            ggTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGagag                                                                 >  2:443008/1‑140 (MQ=255)
                                                                                aaGGGCGTATTGACAAATACTACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAggggg                                                             <  2:3992/140‑1 (MQ=255)
                                                                                     cGTACTGACAAATTCCACAAAGCGCACAAGCCGTCCCCTCGCCCCTGTGGGGAGAGGGTTAGTCTGAGGGGCACAGCCAGCACCCGTGCGAACATAAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGaa                                                        <  2:1680421/140‑1 (MQ=255)
                                                                                     cGTAATGACAATTTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTTGGGTGTGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGaa                                                        <  2:128459/140‑1 (MQ=255)
                                                                                     cGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTACCCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGaa                                                        <  2:1490244/140‑1 (MQ=255)
                                                                                     cGTAATGACAAATTCGACAAAGCGCACAAGCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGGGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGaa                                                        <  2:1218970/140‑1 (MQ=255)
                                                                                     cGTAATGACAAATTCCAAAACGCGCACACTCCGCCCCCCCGCCCCTTTGGGGAGAGGGTTTGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGaa                                                        <  1:1398491/140‑1 (MQ=25)
                                                                                       tAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACg                                                      >  1:1322805/1‑140 (MQ=255)
                                                                                       tAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACg                                                      >  1:480008/1‑140 (MQ=255)
                                                                                          tGACAACTTCGACAACGCGCACAATCCGTCCCCTCGCCCGTTTGGGGAGTGGGTTAGGGTGAGGGTAACAGCCAGCACTGGTGCGAACTTTAACCCTCACCCCAGCCCTCACCCTGGTAGGGAGAGGGGGCAGAACGgcg                                                   <  1:653463/140‑1 (MQ=255)
                                                                                          tGAAAAATTCTACAAAGCGCACAATCCGTCCCCTCGCCCTTTTGGGTTGTGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGgcg                                                   <  1:36184/140‑1 (MQ=25)
                                                                                            aCAAATACTACAAACCCCAAAAACCGACCCCACCCCCATTTGGGTAGAGGTTTAGGGTGAGGGCAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCa                                                 <  1:1435514/140‑1 (MQ=255)
                                                                                             cAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAg                                                >  2:811853/1‑140 (MQ=255)
                                                                                             cAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCGTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAg                                                >  2:812014/1‑140 (MQ=255)
                                                                                                 ttCGACAACGGGCACAATCCGTCCCCTCGCCCCTTGGGGGAGTGGTTTCGGGTGTGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACa                                            <  1:1461633/140‑1 (MQ=255)
                                                                                                 ttCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGTGTGGGTTTGGGTGAGGGGAACAGCCAGCACTGGTGAGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACa                                            <  1:214968/140‑1 (MQ=255)
                                                                                                 ttCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACa                                            <  2:1441185/140‑1 (MQ=255)
                                                                                                      aaaaGCGCAAAATCCGTCCCCTCGCCCCTTTGGGTTGTGGGTTTGGGTGTGGGTAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATcaca                                       <  1:506279/139‑1 (MQ=255)
                                                                                                             cacaATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTTGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCt                                <  1:1643382/140‑1 (MQ=255)
                                                                                                                     gTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGGACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGaa                        <  2:94365/140‑1 (MQ=255)
                                                                                                                     gTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGaa                        >  2:1317604/1‑140 (MQ=255)
                                                                                                                            cGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCa                 >  1:46790/1‑140 (MQ=255)
                                                                                                                            cGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCa                 >  2:1161756/1‑140 (MQ=255)
                                                                                                                            cGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCa                 >  2:1160734/1‑140 (MQ=255)
                                                                                                                                  tttGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCt                                                                                              <  1:1232381/57‑1 (MQ=255)
                                                                                                                                  tttGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCt                                                                                              >  2:1232381/1‑57 (MQ=255)
                                                                                                                                       ggAGTGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCAAATCCATCAATAATGCTtct      <  2:1317348/140‑1 (MQ=255)
                                                                                                                                         agagGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATTCTTCTCa    >  1:440335/1‑140 (MQ=255)
                                                                                                                                         agagGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCa    >  2:680172/1‑140 (MQ=255)
                                                                                                                                         agagGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCa    >  2:90364/1‑140 (MQ=255)
                                                                                                                                           agGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCATCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCAtt  >  1:840322/1‑140 (MQ=255)
                                                                                                                                           agGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGACGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCAtt  >  1:617334/1‑140 (MQ=255)
                                                                                                                                           |                                                                                                                                           
GCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATT  >  NC_000913/624459‑624737

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