Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 987,561 G→C 15.3% intergenic (‑579/+24) ompF ← / ← asnS outer membrane porin 1a (Ia;b;F)/asparaginyl tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913987,5610GC15.3% 192.7 / 25.6 66intergenic (‑579/+24)ompF/asnSouter membrane porin 1a (Ia;b;F)/asparaginyl tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base G (27/28);  new base C (6/4);  total (33/33)
Fisher's exact test for biased strand distribution p-value = 7.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.18e-01

AATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCGGAACAGTACCGTAGCGACGCAGATCGCGATACCACCAGTAATCTTCTTTATTCAGGCCCATTTCCAGCATACGCTCGTCCAGCACGTCCAGACGTTCTTCACGCTGGGAGCCACCA  >  NC_000913/987287‑987813
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                            
aaTTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACCAACAGATAACTTACAAGAAAAGGTAAAAAAAACGCCAACggg                                                                                                                                                                                                                                                              >  2:205753/1‑275 (MQ=255)
                                                                      tCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAAATTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTTACGACGCCAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatcacatcc                                                                                                                                                                                         >  1:141366/1‑273 (MQ=255)
                                                                                                                                     aaGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGc                                                                                                                                                                                                                                                       <  6:416628/149‑1 (MQ=255)
                                                                                                                                      aGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGct                                                                                                                                                                                                                                                      >  6:202658/1‑149 (MQ=255)
                                                                                                                                              gATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGTATGATCGTTATGGTGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGa                                                                                                                                                                                                                                              <  7:329739/149‑1 (MQ=255)
                                                                                                                                                ttAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACg                                                                                                                                                                                                                                            >  8:265869/1‑149 (MQ=255)
                                                                                                                                                  aGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCt                                                                                                                                                                                                                                          <  5:202658/149‑1 (MQ=255)
                                                                                                                                                             gTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGCCCGCTGGCTCttt                                                                                                                                                                                                                                                  >  5:145691/1‑129 (MQ=255)
                                                                                                                                                             gTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGCCCGCTGGCTCttt                                                                                                                                                                                                                                                  <  6:145691/130‑2 (MQ=255)
                                                                                                                                                              tAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGATTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCGGAACAGTACCGTAGCGACGCa                                                                                                 <  2:265256/274‑1 (MQ=255)
                                                                                                                                                                        aaTGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCt                                                                                                                                                                                                                                                          >  8:349336/1‑111 (MQ=255)
                                                                                                                                                                        aaTGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCt                                                                                                                                                                                                                                                          <  7:349336/111‑1 (MQ=255)
                                                                                                                                                                         aTGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAAAAGCTGGCGTTAcgcg                                                                                                                                                                                                                   >  6:176031/1‑149 (MQ=255)
                                                                                                                                                                          tGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCgg                                                                                                                                                                                                                  >  7:132635/1‑149 (MQ=255)
                                                                                                                                                                          tGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGCCCGCTGGctct                                                                                                                                                                                                                                                    >  7:188068/1‑115 (MQ=255)
                                                                                                                                                                          tGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGCCCGCTGGctct                                                                                                                                                                                                                                                    <  8:188068/115‑1 (MQ=255)
                                                                                                                                                                                  tCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGAAAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTGGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATGACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCGGAACAGTACCGTAGCGACGCAGATCGCGATACCACCAGTAAt                                                                            <  2:330053/275‑1 (MQ=255)
                                                                                                                                                                                     gTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGc                                                                                                                                                                                                                           >  7:113484/1‑129 (MQ=255)
                                                                                                                                                                                     gTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGc                                                                                                                                                                                                                           <  8:113484/129‑1 (MQ=255)
                                                                                                                                                                                      tAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCCCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTGACGTAAGCAATCAGACGTTCAAAACCAATACCGAAACCTGAATGCGGAACAGTACCGTAGCGACGCAAATCGCGATACCGCCAGGAAccttt                                                                        >  2:473479/1‑274 (MQ=255)
                                                                                                                                                                                          aGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTg                                                                                                                                                                              <  2:406270/169‑1 (MQ=255)
                                                                                                                                                                                           ggATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGaa                                                                                                                                                                                                 >  8:232387/1‑149 (MQ=255)
                                                                                                                                                                                             atgatCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatc                                                                                                                                                                                               >  7:334329/1‑149 (MQ=255)
                                                                                                                                                                                              tgatCGTTATGGCGGGGTTGCAAATAGATTGCTTGCCTTAACAAACCGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatcacatccagcacgt                                                                                                                                                                                  >  1:131957/1‑153 (MQ=255)
                                                                                                                                                                                               gatCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCt                                                                                                                                                                                                                              >  8:140098/1‑116 (MQ=255)
                                                                                                                                                                                               gatCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCt                                                                                                                                                                                                                              <  7:140098/116‑1 (MQ=255)
                                                                                                                                                                                                atCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatcaca                                                                                                                                                                                            <  7:265869/149‑1 (MQ=255)
                                                                                                                                                                                                  cGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACGTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCGGAACAGTACCGTAGCGACGCAGATCGCGATACCACCAGTAATCTTCTTTATTCAGGcc                                                            <  2:163452/275‑1 (MQ=255)
                                                                                                                                                                                                  cGTTAGGGCGGGGTTGAAAATANATTGCTTGACTTAAAAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGTCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTg                                                                                                                                                                              <  3:131957/161‑1 (MQ=255)
                                                                                                                                                                                                      aTGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGGAGACGCTAAATTAGAAGGTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCacgtacga                                                                                                                                                                                  >  2:502193/1‑152 (MQ=255)
                                                                                                                                                                                                                tGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTAcgc                                                                                                                                                                                                                    >  2:27832/1‑109 (MQ=255)
                                                                                                                                                                                                                    aaTAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAg                                                                                                                                                                                                                <  4:27832/109‑1 (MQ=255)
                                                                                                                                                                                                                    aaTAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatcac                                                                                                                                                                                             <  8:212929/128‑1 (MQ=255)
                                                                                                                                                                                                                    aaTAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatcac                                                                                                                                                                                             >  7:212929/1‑128 (MQ=255)
                                                                                                                                                                                                                         aTTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACg                                                                                                                                                                                                            >  4:302187/1‑108 (MQ=255)
                                                                                                                                                                                                                         aTTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCGGAACAGTACCGTAGCGACGCAGATCGCGATACCACCAGTAATCTTCTTTATTCAGGCCCATTTCCAGCATACGCTCGTCCa                                     >  2:237252/1‑275 (MQ=255)
                                                                                                                                                                                                                             cTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTggg                                                                                                                                                                                                        <  2:302191/108‑1 (MQ=255)
                                                                                                                                                                                                                             cTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCGGAACAGTACCGTAGCGACGCa                                                                                                 <  1:11152/211‑1 (MQ=255)
                                                                                                                                                                                                                                 aCTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGa                                                                                                                                                                                                                                              >  4:401590/1‑66 (MQ=255)
                                                                                                                                                                                                                                     aacaaacaGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCt                                                                                                                                                                                                                                          <  2:401595/66‑1 (MQ=255)
                                                                                                                                                                                                                                     aacaaacaGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCGGAACAGTACCGTAGCGACGCAGATCGCGATACCACCAGTAATCTTCTTTATTCAGGCCCATTTCCAGCATACGCTCGTCCAGCACGTCCacac                         >  2:166363/1‑275 (MQ=255)
                                                                                                                                                                                                                                        aaacaGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAg                                                                                                                                                             <  7:216603/140‑1 (MQ=255)
                                                                                                                                                                                                                                        aaacaGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAg                                                                                                                                                             >  8:216603/1‑140 (MQ=255)
                                                                                                                                                                                                                                                aaCTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGc                                                                                                                                                                         >  6:221616/1‑120 (MQ=255)
                                                                                                                                                                                                                                                aaCTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGc                                                                                                                                                                         <  5:221616/120‑1 (MQ=255)
                                                                                                                                                                                                                                                aaCTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCaaaa                                                                                                                                            <  8:132635/149‑1 (MQ=255)
                                                                                                                                                                                                                                                    tGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCaa                                                                                                                                        <  5:176031/149‑1 (MQ=255)
                                                                                                                                                                                                                                                      aaaannaaGATAAAAAAAGAGCCCGCGGGCGCTGCATCGTGAAGACGCTAAATTNGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATAACGTACGTTTAGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGGCCGAAACCTGAATGCGGAACAGTACCGTAGCGACGCAGATCGCGATACCACCAGTAAt                                                                            <  4:512160/201‑1 (MQ=255)
                                                                                                                                                                                                                                                         gAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCg                                                                                                                                   >  6:449860/1‑149 (MQ=255)
                                                                                                                                                                                                                                                          aaaaGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCGGAACAGTACCGTAGCGACGCAGATcgcg                                                                                          >  3:1685/1‑189 (MQ=255)
                                                                                                                                                                                                                                                          aaaaGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCGGAACAGTACCGTAGCGACGCAGATCGCGATaccacc                                                                                  >  2:329312/1‑197 (MQ=255)
                                                                                                                                                                                                                                                            aaGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCaa                                                                                                                                                          >  4:520823/1‑123 (MQ=255)
                                                                                                                                                                                                                                                            aaGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCAGAGTACGTGGGAACGGAATCACATCCCGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGCCGTTCAAAACCAAGACCGAAACCTGAATTCGGAACGGTACCGTAGGGAAGCAGATCGCGATCCCAACAGTAATCTTCTTTATGCAGGCCAATTGCGCGCATACGATCGTCAGTGACGTCCCAGCGTTCGTAACGGTGGGAgaccaca  >  4:26103/1‑269 (MQ=255)
                                                                                                                                                                                                                                                             agagAAAAAAAGGGCCAGCGGGCGCTGGCTCTTGAAGACGAGAAATTAGTAGCTGGCGTTACGAGGAGTACGGGGGAACGGAATCAACTCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCGGAACAGTACCGTAGCGACGCAGATCGCGATACCACCAGTAATCTTCTTTATTCAGGCCCATTTCCAGCATACGCTCGTCCAGCACGTCCAGACGTTCTTCACGCTGGGAGccacc   <  4:473474/273‑1 (MQ=255)
                                                                                                                                                                                                                                                             aGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCGGAACAGTACCGTAGCGACGCAGATCGCGATACCACCAGTAATCTTTTTTATTCAGGCCCATTTCCAGCATACGCTCGTCCAGCACGTCCAGACGTTCTTCACGCTGGGAGccacc   <  2:92361/273‑1 (MQ=255)
                                                                                                                                                                                                                                                              gATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCGGAACAGTACCGTAGCGACGCAGATCGCGATac                                                                                      <  1:1685/189‑1 (MQ=255)
                                                                                                                                                                                                                                                              gATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCGGAACAGTACCGTAGCGACGCAGATCGCGATACCACCAGTa                                                                              <  4:329308/197‑1 (MQ=255)
                                                                                                                                                                                                                                                               aTAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCt                                                                                                                             <  8:334329/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                tAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAg                                                                                                                                                      <  2:520828/123‑1 (MQ=255)
                                                                                                                                                                                                                                                                      aaGAGCCAGCGCCCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGa                                                                                                                                      <  5:445748/133‑1 (MQ=255)
                                                                                                                                                                                                                                                                      aaGAGCCAGCGCCCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGa                                                                                                                                      >  6:445748/1‑133 (MQ=255)
                                                                                                                                                                                                                                                                       agagCCAGCGCCCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCaa                                                                                                                                        >  7:274418/1‑130 (MQ=255)
                                                                                                                                                                                                                                                                       agagCCAGCGCCCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCaa                                                                                                                                        <  8:274418/130‑1 (MQ=255)
                                                                                                                                                                                                                                                                       agagCCAGCGCCCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCgg                                                                                                                     >  5:124539/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                         agCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCGGaa                                                                                                                   <  7:232387/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                         agCCAGCGCCCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCGGaa                                                                                                                   >  7:331941/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                              gCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCGGAACAGTa                                                                                                              >  6:293915/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                            
AATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCGGAACAGTACCGTAGCGACGCAGATCGCGATACCACCAGTAATCTTCTTTATTCAGGCCCATTTCCAGCATACGCTCGTCCAGCACGTCCAGACGTTCTTCACGCTGGGAGCCACCA  >  NC_000913/987287‑987813

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