Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 173,566 A→G 5.8% intergenic (+122/+36) 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,5660AG 5.8% 283.0 / 5.0 87intergenic (+122/+36)fhuB/hemLiron(3+)‑hydroxamate import ABC transporter permease/glutamate‑1‑semialdehyde aminotransferase (aminomutase)
Reads supporting (aligned to +/- strand):  ref base A (45/36);  new base G (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 = 9.77e-01

TATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGCCAGGTAAAC  >  NC_000913/173418‑173712
                                                                                                                                                    |                                                                                                                                                  
tATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATa                                                                                                                                                    >  3:728741/1‑149 (MQ=255)
        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:549290/149‑1 (MQ=255)
                         aaTTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACaaa                                                                                                                           <  2:198951/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)
                                                tACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTCGGTCGGATCAg                                                                                                                                                  >  1:650561/1‑103 (MQ=255)
                                                tACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATGCGGCGTGAACGCCTTATCCGAc                                                                                                                              <  1:315939/123‑1 (MQ=38)
                                                tACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATGCGGCGTGAACGCCTTATCCGAc                                                                                                                              >  2:315939/1‑123 (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                                                                                >  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                                     >  3:586969/1‑149 (MQ=255)
                                                                                                               aTCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATa                                     >  2:102953/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                      >  3:275757/1‑149 (MQ=255)
                                                                                                                              ggCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCAAAc                      <  2:349395/149‑1 (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             >  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   >  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)
                                                                                                                                                    |                                                                                                                                                  
TATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGCCAGGTAAAC  >  NC_000913/173418‑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: ‑