Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 169,078 C→A 5.0% S532* (TCG→TAG)  fhuA → ferrichrome outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913169,0780CA 5.0% 642.3 / 14.3 188S532* (TCG→TAG) fhuAferrichrome outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base C (99/79);  new base A (3/7);  total (102/86)
Fisher's exact test for biased strand distribution p-value = 1.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.83e-01

Only 100 of 189 total aligned reads displayed.

GGCGGTCGTTATGACTGGGCAGATCAAGAATCTCTTAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATA‑TTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAA  >  NC_000913/168930‑169222
                                                                                                                                                    |                                                                                                                                                 
ggCGGTCGTTATGACTGGGCAGATCAAGAATCTCTTAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATc                                                                                                                                                   <  5:380930/149‑1 (MQ=255)
         tATGACTGGGCAGATCAAGAATCTCTTAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAAcc                                                                                                                                          >  8:322374/1‑149 (MQ=255)
                gggCAGATCAAGAATCTCTTAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCg                                                                                                                                   >  4:92679/1‑149 (MQ=255)
                   cAGATCAAGAATCTCTTAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATAGTTTGAACCTTCTTCGCaa                                                                                                                                <  7:31532/149‑1 (MQ=255)
                      aTCAAGAATCTCTTAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGtt                                                                                                                             >  2:348158/1‑149 (MQ=255)
                      aTCAAGAATCTCTTAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGtt                                                                                                                             <  6:313353/149‑1 (MQ=255)
                      aTCAAGAATCTCTTAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGtt                                                                                                                             >  5:59534/1‑149 (MQ=255)
                       tCAAGAATCTCTTAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGATg                                                                                                                            >  2:304432/1‑149 (MQ=255)
                        cAAGAATCTCTTAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTgg                                                                                                                           >  7:303351/1‑149 (MQ=255)
                         aaGAATCTCTTAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTgg                                                                                                                           >  1:242149/1‑148 (MQ=255)
                         aaGAATCTCTTAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTgg                                                                                                                           <  1:155771/148‑1 (MQ=255)
                           gAATCTCTTAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGttt                                                                                                                                               >  2:33470/1‑126 (MQ=255)
                           gAATCTCTTAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGttt                                                                                                                                               <  1:33470/126‑1 (MQ=255)
                            aaTCTCTTAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGGGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATAGTTTGAACCTTCTTCGCAAGTTGGGAAg                                                                                                                       <  2:376316/149‑1 (MQ=255)
                                  ttAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGAt                                                                                                                    <  2:350242/146‑1 (MQ=255)
                                  ttAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATg                                                                                                                   <  4:221382/147‑1 (MQ=255)
                                  ttAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATAGTTTGAACCTTCTTCGCAAGTTGGGAAGGATg                                                                                                                   <  8:251488/147‑1 (MQ=255)
                                  ttAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACATCAGCTATAGCGAATCGTTTGAACCTTCATCGCAAGTTGGGAAGGAt                                                                                                                    <  7:350437/146‑1 (MQ=255)
                                  ctcaCCGCGTTGCCGGGACGCCCGATAAACCTTATCACAATCAGTTTACCTGGCGTGGGGGTGTTAACGACCTGTTTGATATTGGTGTAACACCTTACTTCAGCTATAGCGATTCGTTTGAACCTTCTTCGCAAGTTGGGCAGGATGGt                                                                                                                 <  6:299204/146‑1 (MQ=255)
                                    aaCCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATgg                                                                                                                  <  7:322374/146‑1 (MQ=255)
                                    aaCCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACAACTTACTTCAGCTATAGCGAATCGTTTTAACCTTCTTCGCAAGTTGGGAATGATGGTaa                                                                                                               >  3:230724/1‑149 (MQ=255)
                                     aCCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGAGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTaa                                                                                                               <  3:120086/148‑1 (MQ=255)
                                         cgTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATA‑ttt                                                                                                         <  3:322654/149‑1 (MQ=255)
                                         cgTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATA‑ttt                                                                                                         >  2:140665/1‑149 (MQ=255)
                                           ttGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATA‑TTTTc                                                                                                       <  5:26171/149‑1 (MQ=255)
                                           ttGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATA‑TTTTc                                                                                                       <  1:131015/149‑1 (MQ=255)
                                                 ggACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATA‑TTTTCGCAcct                                                                                                 >  7:168471/1‑148 (MQ=255)
                                                  gacgacCGATAAACGTGATGCCAAACAGGTTCCCTGTCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGGTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATA‑TTTTCGCACCGt                                                                                                <  1:296775/149‑1 (MQ=255)
                                                  gacgacCGATAAACGTGATGACAAACAGTTTACCTGGCGTTGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATAGTTTGAACCTTGTTCGCACGTTGGGAAGGTTGGTAATA‑TTTTCGGACCGt                                                                                                <  6:287249/149‑1 (MQ=255)
                                                      acCGATAAACGTGATGACAAACAGTTTACCTTGAGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTAAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATA‑TTTTCGCACCGTCTa                                                                                             >  2:276/1‑148 (MQ=255)
                                                        cGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATA‑TTTTCGCACCGTCTAAAg                                                                                          >  8:404774/1‑149 (MQ=255)
                                                          aTAAATGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATA‑tttt                                                                                                        >  6:40011/1‑133 (MQ=255)
                                                          aTAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATA‑TTTTCGCACCGTCTAAAGGt                                                                                        <  7:351631/149‑1 (MQ=255)
                                                             caCGTGATGACAAACAGTTTACCTGGCGTTGTGGTGTTACCTACCTGTTTGATAATGGTGTAACACCTTACTTCGGCAATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATA‑TTTTCGCACCGTCTAAACGTaa                                                                                      <  4:245892/147‑1 (MQ=255)
                                                              aCGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATA‑TTTTCTCACCGTCTAAAGGTAAg                                                                                     >  4:30381/1‑148 (MQ=255)
                                                              aCGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATA‑TTTTCGCACCGTCTAAAGGTAAg                                                                                     >  7:129526/1‑148 (MQ=255)
                                                              aCGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATA‑TTTTCGCACCGTCTAAAGGTAAg                                                                                     >  7:340833/1‑148 (MQ=255)
                                                              aCGTGATGACAAACAGATTACCTGGCGGGGTGGTGGTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTAAAGCGAATAGTTTGAACCTTCTTCGCACGTTGGGAAGGATGGTAATA‑TTTTCGCACCGTCTAAAGGTAAg                                                                                     >  2:30932/1‑148 (MQ=255)
                                                               cGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATA‑TTTTCGCACCGTCTAAAGGTAAGCa                                                                                   <  1:269350/149‑1 (MQ=255)
                                                                gTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATA‑TTTTCGCACCGTCTAAAGGTAAGCAg                                                                                  >  2:120365/1‑149 (MQ=255)
                                                                      aCAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATA‑TTTTCGCACCGTCTAAAGGTAAGCAGTATg                                                                              >  5:5106/1‑147 (MQ=255)
                                                                       cAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCttcttc                                                                                                                                    >  2:8180/1‑93 (MQ=255)
                                                                       cAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCttcttc                                                                                                                                    <  1:8180/93‑1 (MQ=255)
                                                                        aaaCAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATA‑TTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGt                                                                          >  3:327615/1‑149 (MQ=255)
                                                                          aCAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATA‑TTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCg                                                                        >  7:116356/1‑149 (MQ=255)
                                                                                tACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATA‑TTTTCGCACCGTCTAAAGGTAAGc                                                                                    <  3:64338/131‑1 (MQ=255)
                                                                                tACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATA‑TTTTCGCACCGTCTAAAGGTAAGc                                                                                    >  4:64338/1‑131 (MQ=255)
                                                                                  ccTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATA‑TTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGaaa                                                                <  5:136210/149‑1 (MQ=255)
                                                                                  ccTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATA‑TTTTCGCACCGTCTAAAGGTAAGCAGTAAGAAGTCGGCGTGaaa                                                                <  7:382993/149‑1 (MQ=255)
                                                                                  ccTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATAGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATA‑TTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGaaa                                                                <  4:151592/149‑1 (MQ=255)
                                                                                  ccTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATAGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATA‑TTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGaaa                                                                <  1:58988/149‑1 (MQ=255)
                                                                                    tGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATTGCGAATCTTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATA‑TTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAata                                                              >  6:351953/1‑149 (MQ=255)
                                                                                    tGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATA‑TTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAata                                                              >  8:227104/1‑149 (MQ=255)
                                                                                    tGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGCATGGTAATA‑TTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGAGTGAAata                                                              >  1:392002/1‑149 (MQ=255)
                                                                                      gCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCATCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATA‑TTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGTCGTGAAATATg                                                            <  6:145680/149‑1 (MQ=255)
                                                                                      gCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATA‑TTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGCCGGCGTGAAata                                                              <  8:303351/147‑1 (MQ=255)
                                                                                      gCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGGTGGGAAGGATGGTAATA‑TTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAatat                                                             <  6:101609/148‑1 (MQ=255)
                                                                                           gtggtgTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTg                                                                                                                            <  6:5544/81‑1 (MQ=255)
                                                                                                gtTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATA‑TTTTCGCACCGTCTAAAGGTAAGCAGt                                                                                 <  6:300006/118‑1 (MQ=255)
                                                                                                gtTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATA‑TTTTCGCACCGTCTAAAGGTAAGCAGt                                                                                 >  5:300006/1‑118 (MQ=255)
                                                                                                gtTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATA‑TTTTCGCACCGTCTAAAGGTAAGCAGTGTGAAGTCGGCGTGAAATATGTACCgg                                                      >  2:303253/1‑145 (MQ=255)
                                                                                                 ttAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATA‑TTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCgg                                                      >  8:226434/1‑144 (MQ=255)
                                                                                                     cTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATA‑TTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTc                                             <  7:371395/149‑1 (MQ=255)
                                                                                                      tACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATA‑TTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTcc                                            >  5:327720/1‑149 (MQ=255)
                                                                                                        ccTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATA‑TTTTCGCACCGTCTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagc                                          >  8:192027‑M1/1‑99 (MQ=255)
                                                                                                          tGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATA‑TTTTCGCACCGTCTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagcca                                        >  1:11307‑M1/1‑97 (MQ=255)
                                                                                                            tttGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGGTGGGAAGGATGGTAATA‑TTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGAtt                                        <  1:304432/147‑1 (MQ=255)
                                                                                                               gATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATA‑TTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAg                                    <  3:99218/148‑1 (MQ=255)
                                                                                                               gATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATA‑TTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGt                                   >  3:264056/1‑149 (MQ=255)
                                                                                                                aTAATGGGGTAACACCTTACTTCAGCTATAGCTAATCGGTTGAACCTTCTTCGCAAGTGGGGAAGGATGGTAATA‑TTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGt                                   <  4:86805/148‑1 (MQ=255)
                                                                                                                   aTGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATA‑TTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTAc                                <  8:116356/148‑1 (MQ=255)
                                                                                                                   aTGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATA‑TTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTAc                                <  8:340833/148‑1 (MQ=255)
                                                                                                                    tGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGTATGGTAATA‑TTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGc                                                                      >  7:82566/1‑109 (MQ=255)
                                                                                                                    tGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGTATGGTAATA‑TTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGc                                                                      <  8:82566/109‑1 (MQ=255)
                                                                                                                    tGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATA‑TTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTg                                       >  2:167640/1‑140 (MQ=255)
                                                                                                                     ggTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATA‑TTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTgg                             >  5:118662/1‑149 (MQ=255)
                                                                                                                      gtgtAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATA‑TTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGt                            <  1:179551/149‑1 (MQ=255)
                                                                                                                      gtgtAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATA‑TTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGt                            >  5:370632/1‑149 (MQ=255)
                                                                                                                      gtgtAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATA‑TTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGt                            >  8:283730/1‑149 (MQ=255)
                                                                                                                      gggtAACACCTTACTTCAGCAATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAATGATGGTAATA‑TTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCTATTGTAGTTACTGGt                            <  8:94360/147‑1 (MQ=255)
                                                                                                                       tgtAACACCTTACTTCAGCTATAGCGAATAGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATA‑TTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTg                           <  1:348517/149‑1 (MQ=255)
                                                                                                                        gtAACACCTTACTTGAGCTATAGAGAATCGTTTTAACCTTCTTCGCAAGTTGGGAAGGATGGTAATA‑TTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGa                                          <  6:84950/133‑1 (MQ=255)
                                                                                                                        gtAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATA‑TTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGa                                          >  5:84950/1‑133 (MQ=255)
                                                                                                                         tAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATA‑TTTTCGCACCGTCTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttca                          >  7:246747‑M1/1‑82 (MQ=255)
                                                                                                                           acacCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATA‑TTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTTGTGCCgt                       >  7:40234/1‑149 (MQ=255)
                                                                                                                           acacCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATA‑TTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCgt                       >  7:46526/1‑149 (MQ=255)
                                                                                                                             acCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATA‑TTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAata                                                              >  3:174454/1‑108 (MQ=255)
                                                                                                                             acCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATA‑TTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAata                                                              <  4:174454/108‑1 (MQ=255)
                                                                                                                               cTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTTCGCACCGTCTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatga                    <  5:344389‑M1/149‑73 (MQ=255)
                                                                                                                               cTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATA‑TTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGta                    >  7:226533/1‑148 (MQ=255)
                                                                                                                                  aCTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATA‑TTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGtata                  <  2:64192/147‑1 (MQ=255)
                                                                                                                                    ttCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATA‑TTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGt                                      >  4:33046/1‑125 (MQ=255)
                                                                                                                                      cAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATA‑TTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAat                                                               <  6:330716/98‑1 (MQ=255)
                                                                                                                                           ataGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATA‑TTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGAt                                                 >  4:3766/1‑107 (MQ=255)
                                                                                                                                           ataGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATA‑TTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGAt                                                 <  3:3766/107‑1 (MQ=255)
                                                                                                                                            taGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATA‑TTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTaaaa      >  8:307992/1‑149 (MQ=255)
                                                                                                                                            taGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATA‑TTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTaaaa      >  6:170977/1‑149 (MQ=255)
                                                                                                                                              gCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATA‑TTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAAc     >  3:120727/1‑148 (MQ=255)
                                                                                                                                               cGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATA‑TTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAAcc    >  2:182539/1‑148 (MQ=255)
                                                                                                                                                 aaTAGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATA‑TTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACcaa  >  4:152292/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GGCGGTCGTTATGACTGGGCAGATCAAGAATCTCTTAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATA‑TTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAA  >  NC_000913/168930‑169222

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑