Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 169,749 T→A 10.0% intergenic (+22/‑29) fhuA → / → fhuC ferrichrome outer membrane transporter/iron(3+)‑hydroxamate import ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913169,7490TA10.0% 265.5 / 16.9 80intergenic (+22/‑29)fhuA/fhuCferrichrome outer membrane transporter/iron(3+)‑hydroxamate import ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base T (39/33);  new base A (4/4);  total (43/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.35e-01

ACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTG‑CCCATTTCACAAGTTGGCTGTTATGC‑AGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTCCTGCCGGGAAAGTGACC  >  NC_000913/169606‑169897
                                                                                                                                               |                                                                                                                                                      
aCGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGt                                                                                                                                                   <  2:23883/149‑1 (MQ=255)
  gTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGAAATCCGTG‑c                                                                                                                                                <  4:609290/149‑1 (MQ=255)
                   caacCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTG‑CCCATTTCACAAGTTGGc                                                                                                                               >  1:1184570/1‑149 (MQ=255)
                    aacCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTG‑CCCATTTCACAAGTTgg                                                                                                                                >  1:681935/1‑147 (MQ=255)
                    aacCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTG‑CCCATTTCACAAGTTgg                                                                                                                                <  2:681935/147‑1 (MQ=255)
                      ccTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTG‑CCCATTTCACAAGTTGGCTGt                                                                                                                            >  4:194072/1‑149 (MQ=255)
                          ttCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTcc                                                                                                                                                     >  2:1825509/1‑121 (MQ=255)
                          ttCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTcc                                                                                                                                                     <  1:1825509/121‑1 (MQ=255)
                          ttCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTG‑CCCATTTCACAAGTTGGCTGTTATg                                                                                                                        >  1:1237548/1‑149 (MQ=255)
                             gATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTG‑ccc                                                                                                                                              <  4:475831/124‑1 (MQ=255)
                             gATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTG‑ccc                                                                                                                                              >  3:475831/1‑124 (MQ=255)
                              aTCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTG‑CCCATTTCACAAGTTGGCTGtt                                                                                                                           >  3:1437966/1‑142 (MQ=255)
                              aTCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTG‑CCCATTTCACAAGTTGGCTGtt                                                                                                                           <  4:1437966/142‑1 (MQ=255)
                               tCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGAAATCCGTG‑ccc                                                                                                                                              <  2:1596948/122‑1 (MQ=255)
                               tCGTCAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGAAATCCGTG‑ccc                                                                                                                                              >  1:1596948/1‑122 (MQ=255)
                                   gAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTG‑CCCATTTCACAAGTTGGCTGTTATGC‑AGGAATac                                                                                                              <  3:1602629/149‑1 (MQ=255)
                                          tCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTcc                                                                                                                                                     <  4:409623/105‑1 (MQ=255)
                                          tCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTcc                                                                                                                                                     >  3:409623/1‑105 (MQ=255)
                                            gCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTCCCGTTTTTAATTTCTTTTTTGGGGCACGG‑TTTTCGTGCCCCTTTTCAAAGTTTGGCTGTT‑TGCAAGGAAT‑CACGaaataat                                                                                                    >  2:951440/1‑144 (MQ=255)
                                             ccAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTG‑CCCATTTCACAAGTTGGCTGTTATGC‑AGGAATACACGAATCAtt                                                                                                    <  3:194072/149‑1 (MQ=255)
                                               aGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTG‑CCCATTTCACAAGTTGGCTGTTATGC‑AGGAATACACGAATCATTcc                                                                                                  <  2:1184570/149‑1 (MQ=255)
                                                     tttAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTG‑CCCATTTCACAAGTTGGCTGtt                                                                                                                           <  1:1650860/119‑1 (MQ=255)
                                                     tttAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTG‑CCCATTTCACAAGTTGGCTGtt                                                                                                                           >  2:1650860/1‑119 (MQ=255)
                                                       tAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTG‑cc                                                                                                                                               <  4:1846561/97‑1 (MQ=255)
                                                       tAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTG‑cc                                                                                                                                               >  3:1846561/1‑97 (MQ=255)
                                                        aaCACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGt                                                                                                                                                   <  4:1474737/93‑1 (MQ=255)
                                                        aaCACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGt                                                                                                                                                   >  3:1474737/1‑93 (MQ=255)
                                                          cacTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTG‑CCCATTTCACAAGTTGGCTGTTa                                                                                                                          <  2:83184/115‑1 (MQ=255)
                                                          cacTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTG‑CCCATTTCACAAGTTGGCTGTTa                                                                                                                          >  1:83184/1‑115 (MQ=255)
                                                             ttATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTG‑CCCATTTCACAAg                                                                                                                                    <  3:921449/102‑1 (MQ=255)
                                                             ttATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTG‑CCCATTTCACAAg                                                                                                                                    >  4:921449/1‑102 (MQ=255)
                                                               aTGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTG‑CCCATTTCACAAGTTGGCTGTTATGC‑AGGAATACACGAATCATTCCGATACCACTTTTGCAc                                                                                  >  3:1825811/1‑149 (MQ=255)
                                                                   ctgctTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTG‑CCCATTTCACAAGTTGGCTGTTATGC‑AGGAATACACGAATCATTCCGATACCACTTTTGCACTGCg                                                                              <  2:1237548/149‑1 (MQ=255)
                                                                       ttCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGt                                                                                                                                                   <  2:1440510/78‑1 (MQ=255)
                                                                       ttCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGt                                                                                                                                                   >  1:1440510/1‑78 (MQ=255)
                                                                           ggggCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTG‑CCCATTTCACAAGTTGGCTGTTATGC‑AGGAATACACGAATCATTCCGa                                                                                                >  4:1367700/1‑123 (MQ=255)
                                                                           ggggCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTG‑CCCATTTCACAAGTTGGCTGTTATGC‑AGGAATACACGAATCATTCCGa                                                                                                <  3:1367700/123‑1 (MQ=255)
                                                                             ggCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTG‑CCCATTTCACAAGTTgg                                                                                                                                >  4:998188/1‑90 (MQ=255)
                                                                             ggCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTG‑CCCATTTCACAAGTTgg                                                                                                                                <  3:998188/90‑1 (MQ=255)
                                                                                            gTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTG‑CCCATTTCACAAGTTGGCTGTTATGC‑AGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCggg                                                     >  2:1672108/1‑149 (MQ=255)
                                                                                               gTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTG‑CCCATTTCACAAGTTGGCTGTTATGC‑AGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATAtctc                                                                     <  1:420542/130‑1 (MQ=255)
                                                                                               gTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTG‑CCCATTTCACAAGTTGGCTGTTATGC‑AGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATAtctc                                                                     >  2:420542/1‑130 (MQ=255)
                                                                                                    aaacgcaaccTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTG‑CCCATTTCACAAGTTGGCTGTTATGC‑AGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCt                                             <  1:951440/146‑1 (MQ=255)
                                                                                                      ccgcaaccTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTG‑CCCATTTCACAAGTTGGCTGTTATGC‑AGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCtt                                            >  4:1870773/1‑148 (MQ=255)
                                                                                                      ccgcaaccTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTG‑CCCATTTCACAAGTTGGCTGTTATGC‑AGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCtt                                            <  3:1870773/148‑1 (MQ=255)
                                                                                                             cTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTG‑CCCATTTCACAAGTTGGCTGTTATGC‑AGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCtttt                                          >  2:238938/1‑143 (MQ=255)
                                                                                                             cTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTG‑CCCATTTCACAAGTTGGCTGTTATGC‑AGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCtttt                                          <  1:238938/143‑1 (MQ=255)
                                                                                                              ttCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTG‑CCCATTTCACAAGTTGGCTGTTATGC‑AGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCg                                   >  1:1753509/1‑149 (MQ=255)
                                                                                                                   tttCTAATTTCTCTTTTGGGGCACGGATTTCCGTG‑CCCATTTCACAAGTTGGCTGTTATGC‑AGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATAtct                                                                      <  3:1493217/109‑1 (MQ=255)
                                                                                                                   tttCTAATTTCTCTTTTGGGGCACGGATTTCCGTG‑CCCATTTCACAAGTTGGCTGTTATGC‑AGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATAtct                                                                      >  4:1493217/1‑109 (MQ=255)
                                                                                                                   tttCTAATTTCTCTTTTGGGGCACGGATTTCCGTG‑CCCATTTCACAAGTTGGCTGTTATGC‑AGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTc                              <  4:1825811/149‑1 (MQ=255)
                                                                                                                     tCTAATTTCTCTTTTGGGGCACGGATTTCCGTG‑CCCATTTCACAAGTTGGCTGTTATGC‑AGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGgc                                                    >  4:757883/1‑125 (MQ=255)
                                                                                                                     tCTAATTTCTCTTTTGGGGCACGGATTTCCGTG‑CCCATTTCACAAGTTGGCTGTTATGC‑AGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGgc                                                    <  3:757883/125‑1 (MQ=255)
                                                                                                                         aTTTCTCTTTTGGGGCACGGATTTCCGTG‑CCCATTTCACAAGTTGGCTGTTATGC‑AGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCgg                                                      <  1:1415177/119‑1 (MQ=255)
                                                                                                                         aTTTCTCTTTTGGGGCACGGATTTCCGTG‑CCCATTTCACAAGTTGGCTGTTATGC‑AGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCgg                                                      >  2:1415177/1‑119 (MQ=255)
                                                                                                                          tttCTCTTTTGGGGCACGGATTTCCGTG‑CCCATTTCACAAGTTGGCTGTTATGC‑AGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTaa                                                                           <  2:853768/97‑1 (MQ=255)
                                                                                                                          tttCTCTTTTGGGGCACGGATTTCCGTG‑CCCATTTCACAAGTTGGCTGTTATGC‑AGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTaa                                                                           >  1:853768/1‑97 (MQ=255)
                                                                                                                              tctTTTGGGGCACGGATTTCCGTG‑CCCATTTCACAAGTTGGCTGTTATGC‑AGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTc                   <  2:1753509/149‑1 (MQ=255)
                                                                                                                                ttttGGGGCACGGATTTCCGTG‑CCCATTTCACAAGTTGGCTGTTATGC‑AGGAATACACGAATCATTcc                                                                                                  <  3:1523699/68‑1 (MQ=255)
                                                                                                                                ttttGGGGCACGGATTTCCGTG‑CCCATTTCACAAGTTGGCTGTTATGC‑AGGAATACACGAATCATTcc                                                                                                  >  4:1523699/1‑68 (MQ=255)
                                                                                                                                     gggCACGGAAATCCGTG‑CCCATTTCACAAGTTGGCTGTTATGC‑AGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCtttt                                          <  4:347305/119‑1 (MQ=255)
                                                                                                                                     gggCACGGAAATCCGTG‑CCCATTTCACAAGTTGGCTGTTATGC‑AGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCtttt                                          >  3:347305/1‑119 (MQ=255)
                                                                                                                                     gggCACGGAAATCCGTG‑CCCATTTCACAAGTTGGCTGTTATGC‑AGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTaa                         <  2:138153/136‑1 (MQ=255)
                                                                                                                                     gggCACGGAAATCCGTG‑CCCATTTCACAAGTTGGCTGTTATGC‑AGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTaa                         >  1:138153/1‑136 (MQ=255)
                                                                                                                                       gCACGGATTTCCGTG‑CCCATTTCACAAGTTGGCTGTTATGC‑AGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTCCTGCCGGGa          >  1:1724673/1‑149 (MQ=255)
                                                                                                                                       gCACGGATTTCCGTG‑CCCATTTCACAAGTTGGCTGTTATGC‑AGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTCCTGCCGGGa          >  2:756665/1‑149 (MQ=255)
                                                                                                                                       gCACGGAAATCCGTG‑CCCATTTCACAAGTTGGCTGTTATGC‑AGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTCCTGCCGGGa          >  1:1232387/1‑149 (MQ=255)
                                                                                                                                        cACGGATTTCCGTG‑CCCATTTCACAAGTTGGCTGTTATGC‑AGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTCCTGcc              >  4:395355/1‑144 (MQ=255)
                                                                                                                                        cACGGATTTCCGTG‑CCCATTTCACAAGTTGGCTGTTATGC‑AGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTCCTGcc              <  3:395355/144‑1 (MQ=255)
                                                                                                                                        cACGGATTTCCGTG‑CCCATTTCACAAGTTGGCTGTTATGC‑AGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTCCTGCCGGGaa         >  4:581509/1‑149 (MQ=255)
                                                                                                                                         aCGGATTTCCGTG‑CCCATTTCACAAGTTGGCTGTTATGC‑AGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTCCTGCCGGGaaa        >  3:1070081/1‑149 (MQ=255)
                                                                                                                                          cGGATTTCCGTG‑CCCATTTCACAAGTTGGCTGTTATGC‑AGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTCCTGCCGGGAAAg       >  3:185711/1‑149 (MQ=255)
                                                                                                                                           ggATTTCCGTG‑CCCATTTCACAAGTTGGCTGTTATGC‑AGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTCCTGCCGGGaa         >  1:688199/1‑146 (MQ=255)
                                                                                                                                           ggATTTCCGTG‑CCCATTTCACAAGTTGGCTGTTATGC‑AGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTCCTGCCGGGaa         <  2:688199/146‑1 (MQ=255)
                                                                                                                                           ggATTTCCGTG‑CCCATTTCACAAGTTGGCTGTTATGC‑AGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTCCTGCCGGGAAAGt      >  1:1002655/1‑149 (MQ=255)
                                                                                                                                           ggATTTCCGTG‑CCCATTTCACAAGTTGGCTGTTATGC‑AGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTCCTGCCGGGAAAGt      <  2:1002655/149‑1 (MQ=255)
                                                                                                                                           ggATTTCCGTG‑CCCATTTCACAAGTTGGCTGTTATGC‑AGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTCCTGCCGGGAAAGt      >  4:1248059/1‑149 (MQ=255)
                                                                                                                                            gATTTCCGTG‑CCCATTTCACAAGTTGGCTGTTATGC‑AGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGc                                              >  4:853448/1‑108 (MQ=255)
                                                                                                                                            gATTTCCGTG‑CCCATTTCACAAGTTGGCTGTTATGC‑AGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGc                                              <  3:853448/108‑1 (MQ=255)
                                                                                                                                            gATTTCCGTG‑CCCATTTCACAAGTTGGCTGTTATGC‑AGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTCCTGCCGGGAAAGTg     >  3:1739674/1‑149 (MQ=255)
                                                                                                                                              tttCCGTG‑CCCATTTCACAAGTTGGCTGTTATGC‑AGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTCCCGGCAGGAAAggtta   <  2:721827/149‑5 (MQ=255)
                                                                                                                                               ttCCGTG‑CCCATTTCACAAGTTGGCTGTTATGC‑AGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTCCTGCCGGGAAAGTGAcc  >  4:518246/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                      
ACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTG‑CCCATTTCACAAGTTGGCTGTTATGC‑AGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTCCTGCCGGGAAAGTGACC  >  NC_000913/169606‑169897

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑