Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 173,567 A→C 5.8% intergenic (+123/+35) 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,5670AC 5.8% 303.7 / 10.5 87intergenic (+123/+35)fhuB/hemLiron(3+)‑hydroxamate import ABC transporter permease/glutamate‑1‑semialdehyde aminotransferase (aminomutase)
Reads supporting (aligned to +/- strand):  ref base A (46/36);  new base C (3/2);  total (49/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 = 4.46e-01

GAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGCCAGGTAAAC  >  NC_000913/173426‑173712
                                                                                                                                             |                                                                                                                                                 
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)
               ttAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAg                                                                                                                                                  >  3:349424/1‑128 (MQ=255)
               ttAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAg                                                                                                                                                  <  4:349424/128‑1 (MQ=255)
               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:730921/1‑149 (MQ=255)
                  aTTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAAc                                                                                                                          >  2:658502/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)
                                        tACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTCGGTCGGATCAg                                                                                                                                                  >  1:650561/1‑103 (MQ=255)
                                        tACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATGCGGCGTGAACGCCTTATCCGAc                                                                                                                              >  2:315939/1‑123 (MQ=38)
                                        tACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATGCGGCGTGAACGCCTTATCCGAc                                                                                                                              <  1:315939/123‑1 (MQ=38)
                                        tACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAg                                                                                                                                                  <  2:650561/103‑1 (MQ=255)
                                             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                                                                                <  2:814082/141‑1 (MQ=255)
                                                                    gACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATc                                                                                >  1:814082/1‑141 (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                                                                    <  1:277510/141‑1 (MQ=255)
                                                                                ccGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTAtt                                                                    >  2:277510/1‑141 (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                                                                       >  2:479644/1‑124 (MQ=255)
                                                                                              gcgcGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGtt                                                                       <  1:479644/124‑1 (MQ=255)
                                                                                                   tCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGa                                         >  3:519016/1‑149 (MQ=255)
                                                                                                       aTCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATa                                     >  2:102953/1‑149 (MQ=255)
                                                                                                       aTCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATa                                     >  3:586969/1‑149 (MQ=255)
                                                                                                        tCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTcc                                                          <  2:598683/127‑1 (MQ=255)
                                                                                                        tCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTcc                                                          >  1:598683/1‑127 (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:586969/149‑1 (MQ=255)
                                                                                                                                  tAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGcc          <  4:183173/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             >  4:337205/1‑143 (MQ=255)
                                                                                                                                     gTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGt             <  3:337205/143‑1 (MQ=255)
                                                                                                                                     gTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGCCAgg       >  3:526014/1‑149 (MQ=255)
                                                                                                                                        ggATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGCCAGGTaa    >  1:283199/1‑149 (MQ=255)
                                                                                                                                         gATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGCCa         >  4:579320/1‑143 (MQ=255)
                                                                                                                                         gATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGCCa         <  3:579320/143‑1 (MQ=255)
                                                                                                                                         gATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGCCAGGTaaa   <  4:374552/149‑1 (MQ=255)
                                                                                                                                         gATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGCCAGGTaaa   >  4:162891/1‑149 (MQ=255)
                                                                                                                                         gATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGCCAGGTaaa   >  1:533580/1‑149 (MQ=255)
                                                                                                                                          aTAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGCCAGGTAAac  <  1:434608/149‑1 (MQ=255)
                                                                                                                                          aTAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGCCAGGTAAac  >  3:475284/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                 
GAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGCCAGGTAAAC  >  NC_000913/173426‑173712

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