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.

GGCGGTCGTTATGACTGGGCAGATCAAGAATCTCTTAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAA  >  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)
                                         cgTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATAttt                                                                                                         <  3:322654/149‑1 (MQ=255)
                                         cgTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATAttt                                                                                                         >  2:140665/1‑149 (MQ=255)
                                           ttGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTc                                                                                                       <  5:26171/149‑1 (MQ=255)
                                           ttGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTc                                                                                                       <  1:131015/149‑1 (MQ=255)
                                                 ggACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCAcct                                                                                                 >  7:168471/1‑148 (MQ=255)
                                                  gacgacCGATAAACGTGATGCCAAACAGGTTCCCTGTCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGGTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGt                                                                                                <  1:296775/149‑1 (MQ=255)
                                                  gacgacCGATAAACGTGATGACAAACAGTTTACCTGGCGTTGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATAGTTTGAACCTTGTTCGCACGTTGGGAAGGTTGGTAATATTTTCGGACCGt                                                                                                <  6:287249/149‑1 (MQ=255)
                                                      acCGATAAACGTGATGACAAACAGTTTACCTTGAGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTAAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTa                                                                                             >  2:276/1‑148 (MQ=255)
                                                        cGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAg                                                                                          >  8:404774/1‑149 (MQ=255)
                                                          aTAAATGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATAtttt                                                                                                        >  6:40011/1‑133 (MQ=255)
                                                          aTAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGt                                                                                        <  7:351631/149‑1 (MQ=255)
                                                             caCGTGATGACAAACAGTTTACCTGGCGTTGTGGTGTTACCTACCTGTTTGATAATGGTGTAACACCTTACTTCGGCAATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAACGTaa                                                                                      <  4:245892/147‑1 (MQ=255)
                                                              aCGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCTCACCGTCTAAAGGTAAg                                                                                     >  4:30381/1‑148 (MQ=255)
                                                              aCGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAg                                                                                     >  7:129526/1‑148 (MQ=255)
                                                              aCGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAg                                                                                     >  7:340833/1‑148 (MQ=255)
                                                              aCGTGATGACAAACAGATTACCTGGCGGGGTGGTGGTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTAAAGCGAATAGTTTGAACCTTCTTCGCACGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAg                                                                                     >  2:30932/1‑148 (MQ=255)
                                                               cGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCa                                                                                   <  1:269350/149‑1 (MQ=255)
                                                                gTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAg                                                                                  >  2:120365/1‑149 (MQ=255)
                                                                      aCAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATg                                                                              >  5:5106/1‑147 (MQ=255)
                                                                       cAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCttcttc                                                                                                                                    >  2:8180/1‑93 (MQ=255)
                                                                       cAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCttcttc                                                                                                                                    <  1:8180/93‑1 (MQ=255)
                                                                        aaaCAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGt                                                                          >  3:327615/1‑149 (MQ=255)
                                                                          aCAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCg                                                                        >  7:116356/1‑149 (MQ=255)
                                                                                tACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGc                                                                                    <  3:64338/131‑1 (MQ=255)
                                                                                tACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGc                                                                                    >  4:64338/1‑131 (MQ=255)
                                                                                  ccTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGaaa                                                                <  5:136210/149‑1 (MQ=255)
                                                                                  ccTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTAAGAAGTCGGCGTGaaa                                                                <  7:382993/149‑1 (MQ=255)
                                                                                  ccTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATAGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGaaa                                                                <  4:151592/149‑1 (MQ=255)
                                                                                  ccTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATAGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGaaa                                                                <  1:58988/149‑1 (MQ=255)
                                                                                    tGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATTGCGAATCTTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAata                                                              >  6:351953/1‑149 (MQ=255)
                                                                                    tGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAata                                                              >  8:227104/1‑149 (MQ=255)
                                                                                    tGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGCATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGAGTGAAata                                                              >  1:392002/1‑149 (MQ=255)
                                                                                      gCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCATCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGTCGTGAAATATg                                                            <  6:145680/149‑1 (MQ=255)
                                                                                      gCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGCCGGCGTGAAata                                                              <  8:303351/147‑1 (MQ=255)
                                                                                      gCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGGTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAatat                                                             <  6:101609/148‑1 (MQ=255)
                                                                                           gtggtgTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTg                                                                                                                            <  6:5544/81‑1 (MQ=255)
                                                                                                gtTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGt                                                                                 <  6:300006/118‑1 (MQ=255)
                                                                                                gtTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGt                                                                                 >  5:300006/1‑118 (MQ=255)
                                                                                                gtTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTGTGAAGTCGGCGTGAAATATGTACCgg                                                      >  2:303253/1‑145 (MQ=255)
                                                                                                 ttAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCgg                                                      >  8:226434/1‑144 (MQ=255)
                                                                                                     cTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTc                                             <  7:371395/149‑1 (MQ=255)
                                                                                                      tACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTcc                                            >  5:327720/1‑149 (MQ=255)
                                                                                                        ccTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagc                                          >  8:192027‑M1/1‑99 (MQ=255)
                                                                                                          tGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagcca                                        >  1:11307‑M1/1‑97 (MQ=255)
                                                                                                            tttGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGGTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGAtt                                        <  1:304432/147‑1 (MQ=255)
                                                                                                               gATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAg                                    <  3:99218/148‑1 (MQ=255)
                                                                                                               gATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGt                                   >  3:264056/1‑149 (MQ=255)
                                                                                                                aTAATGGGGTAACACCTTACTTCAGCTATAGCTAATCGGTTGAACCTTCTTCGCAAGTGGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGt                                   <  4:86805/148‑1 (MQ=255)
                                                                                                                   aTGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTAc                                <  8:116356/148‑1 (MQ=255)
                                                                                                                   aTGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTAc                                <  8:340833/148‑1 (MQ=255)
                                                                                                                    tGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGTATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGc                                                                      >  7:82566/1‑109 (MQ=255)
                                                                                                                    tGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGTATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGc                                                                      <  8:82566/109‑1 (MQ=255)
                                                                                                                    tGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTg                                       >  2:167640/1‑140 (MQ=255)
                                                                                                                     ggTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTgg                             >  5:118662/1‑149 (MQ=255)
                                                                                                                      gtgtAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGt                            <  1:179551/149‑1 (MQ=255)
                                                                                                                      gtgtAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGt                            >  5:370632/1‑149 (MQ=255)
                                                                                                                      gtgtAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGt                            >  8:283730/1‑149 (MQ=255)
                                                                                                                      gggtAACACCTTACTTCAGCAATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAATGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCTATTGTAGTTACTGGt                            <  8:94360/147‑1 (MQ=255)
                                                                                                                       tgtAACACCTTACTTCAGCTATAGCGAATAGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTg                           <  1:348517/149‑1 (MQ=255)
                                                                                                                        gtAACACCTTACTTGAGCTATAGAGAATCGTTTTAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGa                                          <  6:84950/133‑1 (MQ=255)
                                                                                                                        gtAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGa                                          >  5:84950/1‑133 (MQ=255)
                                                                                                                         tAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttca                          >  7:246747‑M1/1‑82 (MQ=255)
                                                                                                                           acacCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTTGTGCCgt                       >  7:40234/1‑149 (MQ=255)
                                                                                                                           acacCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCgt                       >  7:46526/1‑149 (MQ=255)
                                                                                                                             acCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAata                                                              >  3:174454/1‑108 (MQ=255)
                                                                                                                             acCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAata                                                              <  4:174454/108‑1 (MQ=255)
                                                                                                                               cTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTTCGCACCGTCTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatga                    <  5:344389‑M1/149‑73 (MQ=255)
                                                                                                                               cTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGta                    >  7:226533/1‑148 (MQ=255)
                                                                                                                                  aCTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGtata                  <  2:64192/147‑1 (MQ=255)
                                                                                                                                    ttCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGt                                      >  4:33046/1‑125 (MQ=255)
                                                                                                                                      cAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAat                                                               <  6:330716/98‑1 (MQ=255)
                                                                                                                                           ataGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGAt                                                 >  4:3766/1‑107 (MQ=255)
                                                                                                                                           ataGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGAt                                                 <  3:3766/107‑1 (MQ=255)
                                                                                                                                            taGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTaaaa      >  8:307992/1‑149 (MQ=255)
                                                                                                                                            taGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTaaaa      >  6:170977/1‑149 (MQ=255)
                                                                                                                                              gCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAAc     >  3:120727/1‑148 (MQ=255)
                                                                                                                                               cGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAAcc    >  2:182539/1‑148 (MQ=255)
                                                                                                                                                 aaTAGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACcaa  >  4:152292/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GGCGGTCGTTATGACTGGGCAGATCAAGAATCTCTTAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAA  >  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: