Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 173,573 T→G 5.5% intergenic (+129/+29) fhuB → / ← hemL iron(3+)‑hydroxamate import ABC transporter permease/glutamate‑1‑semialdehyde aminotransferase (aminomutase)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913173,5730TG 5.5% 301.5 / 7.2 86intergenic (+129/+29)fhuB/hemLiron(3+)‑hydroxamate import ABC transporter permease/glutamate‑1‑semialdehyde aminotransferase (aminomutase)
Reads supporting (aligned to +/- strand):  ref base T (45/36);  new base G (3/2);  total (48/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.90e-01

GAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGCCAGGTAAACACCTTC  >  NC_000913/173426‑173718
                                                                                                                                                   |                                                                                                                                                 
cagaAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACACCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTc                                                                                                                                                  <  3:69247/148‑1 (MQ=255)
          agCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATGCGGCGTGAAcgccttatc                                                                                                                                        <  4:666314/149‑6 (MQ=39)
            ccGTTAATTTTTCCATATAACGGCCCTGTACGCCTGAACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCg                                                                                                                                      >  4:582599/1‑149 (MQ=255)
               ttAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATGCGGCGTGAAcgccttat                                                                                                                                         >  4:241844/1‑139 (MQ=39)
               ttAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACa                                                                                                                                   >  2:527833/1‑149 (MQ=255)
               ttAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACa                                                                                                                                   >  3:361454/1‑149 (MQ=255)
                 aaTTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACaaa                                                                                                                                 <  2:198951/149‑1 (MQ=255)
                 aaTTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACaaa                                                                                                                                 <  2:549290/149‑1 (MQ=255)
                  aTTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAAc                                                                                                                                >  2:658502/1‑149 (MQ=255)
                  aTTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAAc                                                                                                                                >  2:730921/1‑149 (MQ=255)
                          atatAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTgg                                                                                                                        >  3:543364/1‑149 (MQ=255)
                          atatAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGAGGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTgg                                                                                                                        <  3:582599/149‑1 (MQ=255)
                           tataACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGa                                                                                                                       <  1:527833/149‑1 (MQ=255)
                            ataACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGAt                                                                                                                      <  1:291903/149‑1 (MQ=255)
                             taACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATc                                                                                                                     <  2:443774/149‑1 (MQ=255)
                               aCGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATcac                                                                                                                   >  2:90077/1‑149 (MQ=255)
                                 ggCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACaa                                                                                                                 >  4:771566/1‑149 (MQ=255)
                                    ccTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGGCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACtt                                                                                                              >  2:559411/1‑149 (MQ=255)
                                    ccTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCCGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACtt                                                                                                              >  4:535188/1‑149 (MQ=255)
                                       gTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAGCTTCGc                                                                                                           <  1:730921/149‑1 (MQ=255)
                                       gTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGc                                                                                                           <  1:658502/149‑1 (MQ=255)
                                        tACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATGCGGCGTGAACGCCTTATCCGAc                                                                                                                                    >  2:315939/1‑123 (MQ=38)
                                        tACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATGCGGCGTGAACGCCTTATCCGAc                                                                                                                                    <  1:315939/123‑1 (MQ=38)
                                             cTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAacac                                                                                                     >  3:183173/1‑149 (MQ=255)
                                                      tAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTgc                                                                                            <  1:90077/149‑1 (MQ=255)
                                                         aGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTgcagc                                                                                         >  1:609887/1‑149 (MQ=255)
                                                         aGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTgcagc                                                                                         >  4:713913/1‑149 (MQ=255)
                                                          gCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTgcagca                                                                                        >  3:126048/1‑149 (MQ=255)
                                                                  gcgACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGgt                                                                                <  4:89133/149‑1 (MQ=255)
                                                                    gACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATc                                                                                      >  1:814082/1‑141 (MQ=255)
                                                                    gACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATc                                                                                      <  2:814082/141‑1 (MQ=255)
                                                                    gACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGgtgt                                                                              >  1:405934/1‑149 (MQ=255)
                                                                       aaCGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTAt                                                                           >  1:303005/1‑149 (MQ=255)
                                                                        aCGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTAtt                                                                          <  4:126048/149‑1 (MQ=255)
                                                                           gATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGat                                                                       >  2:671607/1‑149 (MQ=255)
                                                                                ccGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTAtt                                                                          >  2:277510/1‑141 (MQ=255)
                                                                                ccGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTAtt                                                                          <  1:277510/141‑1 (MQ=255)
                                                                                ccGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATc                                                                    <  4:133986/147‑1 (MQ=255)
                                                                                ccGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATc                                                                    >  3:133986/1‑147 (MQ=255)
                                                                                   aTGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCa                                                               <  2:405934/149‑1 (MQ=255)
                                                                                          gCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCtgtg                                                        >  1:349395/1‑149 (MQ=255)
                                                                                             ggCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTgcgc                                                     >  3:571928/1‑149 (MQ=255)
                                                                                              gcgcGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGtt                                                                             <  1:479644/124‑1 (MQ=255)
                                                                                              gcgcGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGtt                                                                             >  2:479644/1‑124 (MQ=255)
                                                                                                   tCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGa                                               >  3:519016/1‑149 (MQ=255)
                                                                                                       aTCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATa                                           >  2:102953/1‑149 (MQ=255)
                                                                                                       aTCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATa                                           >  3:586969/1‑149 (MQ=255)
                                                                                                        tCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTcc                                                                >  1:598683/1‑127 (MQ=255)
                                                                                                        tCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTcc                                                                <  2:598683/127‑1 (MQ=255)
                                                                                                            gtcTACAAACGGCCCCGAATCGTAGTTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGATAtt                                                                          >  2:253451/3‑113 (MQ=255)
                                                                                                            gtcTACAAACGGCCCCGAATCGTAGGTCGGATAAGGAGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTAtt                                                                          <  1:253451/111‑1 (MQ=255)
                                                                                                               tACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGc                                   >  4:208528/1‑149 (MQ=255)
                                                                                                               tACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGCGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGc                                   >  3:629198/1‑149 (MQ=255)
                                                                                                                aCAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCt                                  >  3:374552/1‑149 (MQ=255)
                                                                                                                  aaaCGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTc                                >  2:434608/1‑149 (MQ=255)
                                                                                                                     cGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCaaa                             <  3:771566/149‑1 (MQ=255)
                                                                                                                      ggCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCAAAc                            <  2:349395/149‑1 (MQ=255)
                                                                                                                      ggCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCAAAc                            >  3:275757/1‑149 (MQ=255)
                                                                                                                         cccGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCt                         <  4:728741/149‑1 (MQ=255)
                                                                                                                          ccGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGTAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTg                        >  4:463232/1‑149 (MQ=255)
                                                                                                                           cGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTGa                       <  3:535188/149‑1 (MQ=255)
                                                                                                                           cGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTGCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTGa                       <  1:559411/149‑1 (MQ=255)
                                                                                                                            gAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCtt                                                                  >  1:201371/1‑105 (MQ=255)
                                                                                                                            gAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCtt                                                                  <  2:201371/105‑1 (MQ=255)
                                                                                                                               tCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGt                   >  1:259627/1‑149 (MQ=255)
                                                                                                                                 gTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGc                 <  4:361454/149‑1 (MQ=255)
                                                                                                                                  tttgTCGGATGCGGCGTGAACGCCTTATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGcc                >  3:389067/4‑149 (MQ=255)
                                                                                                                                  tAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGcc                <  4:183173/149‑1 (MQ=255)
                                                                                                                                  tAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGcc                <  4:586969/149‑1 (MQ=255)
                                                                                                                                  tAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGCGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGcc                <  4:629198/149‑1 (MQ=255)
                                                                                                                                    ggTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCTATGGTGTTATTGATCTCTTCCATGCTGTACGCCACGTACATAAAGCCCGCTTCAAACGCTGACGGTGCCAg              >  2:131780/1‑149 (MQ=255)
                                                                                                                                    ggTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGCCAg              <  2:303005/149‑1 (MQ=255)
                                                                                                                                     gTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGt                   <  3:337205/143‑1 (MQ=255)
                                                                                                                                     gTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGt                   >  4:337205/1‑143 (MQ=255)
                                                                                                                                     gTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGCCAgg             >  3:526014/1‑149 (MQ=255)
                                                                                                                                        ggATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGCCAGGTaa          >  1:283199/1‑149 (MQ=255)
                                                                                                                                         gATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGCCa               <  3:579320/143‑1 (MQ=255)
                                                                                                                                         gATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGCCa               >  4:579320/1‑143 (MQ=255)
                                                                                                                                         gATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGCCAGGTaaa         <  4:374552/149‑1 (MQ=255)
                                                                                                                                         gATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGCCAGGTaaa         >  1:533580/1‑149 (MQ=255)
                                                                                                                                         gATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGCCAGGTaaa         >  4:162891/1‑149 (MQ=255)
                                                                                                                                          aTAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGCCAGGTAAac        >  3:475284/1‑149 (MQ=255)
                                                                                                                                          aTAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGCCAGGTAAac        <  1:434608/149‑1 (MQ=255)
                                                                                                                                               gCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGCCAGGTAAACACCtt   <  4:571928/149‑1 (MQ=255)
                                                                                                                                               gCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGCCAGGTAAACACCtt   >  4:709770/1‑149 (MQ=255)
                                                                                                                                                cGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGCCAGGTAAACACCTTc  <  1:671607/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGCCAGGTAAACACCTTC  >  NC_000913/173426‑173718

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: