Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 169,748 T→A 10.0% intergenic (+21/‑30) 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,7480TA10.0% 262.2 / 17.6 80intergenic (+21/‑30)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 = 6.25e-01

ACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTG‑CCCATTTCACAAGTTGGCTGTTATGC‑AGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTCCTGCCGGGAAAGTGAC  >  NC_000913/169606‑169896
                                                                                                                                              |                                                                                                                                                      
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                                                                                                                                             >  3:475831/1‑124 (MQ=255)
                             gATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTG‑ccc                                                                                                                                             <  4:475831/124‑1 (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                                                                                                                         >  1:83184/1‑115 (MQ=255)
                                                          cacTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTG‑CCCATTTCACAAGTTGGCTGTTa                                                                                                                         <  2:83184/115‑1 (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                                                                                                                               <  3:998188/90‑1 (MQ=255)
                                                                             ggCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTG‑CCCATTTCACAAGTTgg                                                                                                                               >  4:998188/1‑90 (MQ=255)
                                                                                            gTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTG‑CCCATTTCACAAGTTGGCTGTTATGC‑AGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCggg                                                    >  2:1672108/1‑149 (MQ=255)
                                                                                               gTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTG‑CCCATTTCACAAGTTGGCTGTTATGC‑AGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATAtctc                                                                    >  2:420542/1‑130 (MQ=255)
                                                                                               gTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTG‑CCCATTTCACAAGTTGGCTGTTATGC‑AGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATAtctc                                                                    <  1:420542/130‑1 (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                                         <  1:238938/143‑1 (MQ=255)
                                                                                                             cTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTG‑CCCATTTCACAAGTTGGCTGTTATGC‑AGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCtttt                                         >  2:238938/1‑143 (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                                                     >  2:1415177/1‑119 (MQ=255)
                                                                                                                         aTTTCTCTTTTGGGGCACGGATTTCCGTG‑CCCATTTCACAAGTTGGCTGTTATGC‑AGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCgg                                                     <  1:1415177/119‑1 (MQ=255)
                                                                                                                          tttCTCTTTTGGGGCACGGATTTCCGTG‑CCCATTTCACAAGTTGGCTGTTATGC‑AGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTaa                                                                          >  1:853768/1‑97 (MQ=255)
                                                                                                                          tttCTCTTTTGGGGCACGGATTTCCGTG‑CCCATTTCACAAGTTGGCTGTTATGC‑AGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTaa                                                                          <  2:853768/97‑1 (MQ=255)
                                                                                                                              tctTTTGGGGCACGGATTTCCGTG‑CCCATTTCACAAGTTGGCTGTTATGC‑AGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTc                  <  2:1753509/149‑1 (MQ=255)
                                                                                                                                ttttGGGGCACGGATTTCCGTG‑CCCATTTCACAAGTTGGCTGTTATGC‑AGGAATACACGAATCATTcc                                                                                                 >  4:1523699/1‑68 (MQ=255)
                                                                                                                                ttttGGGGCACGGATTTCCGTG‑CCCATTTCACAAGTTGGCTGTTATGC‑AGGAATACACGAATCATTcc                                                                                                 <  3:1523699/68‑1 (MQ=255)
                                                                                                                                     gggCACGGAAATCCGTG‑CCCATTTCACAAGTTGGCTGTTATGC‑AGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCtttt                                         >  3:347305/1‑119 (MQ=255)
                                                                                                                                     gggCACGGAAATCCGTG‑CCCATTTCACAAGTTGGCTGTTATGC‑AGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCtttt                                         <  4:347305/119‑1 (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             <  3:395355/144‑1 (MQ=255)
                                                                                                                                        cACGGATTTCCGTG‑CCCATTTCACAAGTTGGCTGTTATGC‑AGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTCCTGcc             >  4:395355/1‑144 (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     >  4:1248059/1‑149 (MQ=255)
                                                                                                                                           ggATTTCCGTG‑CCCATTTCACAAGTTGGCTGTTATGC‑AGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTCCTGCCGGGAAAGt     <  2:1002655/149‑1 (MQ=255)
                                                                                                                                           ggATTTCCGTG‑CCCATTTCACAAGTTGGCTGTTATGC‑AGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTCCTGCCGGGAAAGt     >  1:1002655/1‑149 (MQ=255)
                                                                                                                                            gATTTCCGTG‑CCCATTTCACAAGTTGGCTGTTATGC‑AGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGc                                             <  3:853448/108‑1 (MQ=255)
                                                                                                                                            gATTTCCGTG‑CCCATTTCACAAGTTGGCTGTTATGC‑AGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGc                                             >  4:853448/1‑108 (MQ=255)
                                                                                                                                            gATTTCCGTG‑CCCATTTCACAAGTTGGCTGTTATGC‑AGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTCCTGCCGGGAAAGTg    >  3:1739674/1‑149 (MQ=255)
                                                                                                                                              tttCCGTG‑CCCATTTCACAAGTTGGCTGTTATGC‑AGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTCCCGGCAGGAAAggtta  <  2:721827/149‑5 (MQ=255)
                                                                                                                                              |                                                                                                                                                      
ACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTG‑CCCATTTCACAAGTTGGCTGTTATGC‑AGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTCCTGCCGGGAAAGTGAC  >  NC_000913/169606‑169896

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: ‑