Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 998,954 T→C 6.8% A155A (GCT→GCC elfD → putative fimbrial chaperone ElfD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913998,9540TC 6.8% 226.2 / 2.3 75A155A (GCT→GCCelfDputative fimbrial chaperone ElfD
Reads supporting (aligned to +/- strand):  ref base T (34/35);  new base C (4/2);  total (38/37)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

CCCCAGATAATAAAAAATCAGCAGACTTTATTATTACACCACCGCTATTTGTGCTGAATCCGGCAAATGAGAATCTGTTACGCATTATGTACATTGGAGCGCCGTTGGCGAAAGACAGAGAAACCCTTTTCTTCACTAGCGTACGGGCAGTCCCTTCAACAACGAAGCGGAAAGAGGGAAATACCCTGAAGATTGCCACACAAAGCGTCATCAAACTTTTCTGGCGACCAAAAGGTTTAGCGTATCCCTTAGGCGAGGCTCCGGCGAAACTGCGTTGCACTTCGTCAGCTGACATGGTTACGGTCAGTAACCCAACACCTTATTTCATTACCCTGACAGACCTGAAAATAGGTGGAAAAGTAGTTAAAAATCAAATGATTTCCCCCTTTGATAAATACCAATTTTCTCTGCCAAAGGGGGCCAAAAATAGCAGCGTAACGTATCGAACCATCAATGACTACGGGGCGGAAACGCCGCAACTCAACTGTAAATCGTAAGCCGTCTTCAGTTAAGAGAGCGAGATGTATAGA  >  NC_000913/998695‑999224
                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                              
ccccAGATAATAAAAAATCAGCAGACTTTATTATTACACCACCGCTATTTGTGCTGAATCCGGCAAATGAGAATCTGTTACGCATTATGTACATTGGAGCGCCGTTGGCGAAAGACAGAGAAACCCTTTTCTTCACTAGCGTACGGGCAGTCCCTTCAACAACGAAGCGGAAAGAGGGAAATACCCTGAAGATTGCCACACAAAGCGTCATCAAACTTTTCTGGCCACCAAAAGGTTTAGCGTAACCCTTAGGCGAGGCCCCGGAGaaaag                                                                                                                                                                                                                                                                     >  2:94074/1‑269 (MQ=255)
              aaaTCAGCAGACTTTATTATTACACCACCGCTATTTGTGCTGAGTCCGGCAAATGAGAATCTGTTACGCATTATGTACATTGGAGCGCCGTTGGCGAAAGACAGAGAAACCCTTTTCTTCACTAGCGTACGGGCAGTCCCTTCAACAACGAAGCGGAAAGAGGGAAATACCCTGAAGATTGCCACACAAAGCTGCATCAAACTTTTCTGGCGACCAAAAGGTTTAGCATATCCCTTAGTCGAGGCTCCGGCGAAACTGCGTTCCACTTCGt                                                                                                                                                                                                                                                       >  1:802490/1‑271 (MQ=255)
                                       caccGCTATTTGTGCTGAATCCGGCAAATGAGAATCTGTTACGCATTATGTACATTGGAGCGCCGTTGGCGAAAGACAGAGAAACCCTTTTCTTCACTAGCGTACGGGCAGTCCCTTCAACAACGAAGCGGAAAGAGGGAAATACCCTGAAGATTGCCACACAAAGCGTCATCAAACTTTTATGGCGACCAAAAGGTTTAGCGTATCCCTTAGGCGAGGCTCCGGCGAAACTGCGTTGCACTTCGTCAGCTGACATGGTTACGGTCAGCAAc                                                                                                                                                                                                                             >  1:964553/1‑272 (MQ=255)
                                        accGCTATTTGTGCTGAATCCGGCAAATGAGAATCTGTTACGCATTATGTACATTGGAGCGCCGTTGGCGAAAGACAGAGAAACCCTTTTCTTCACTAGCGTACGGGCAGTCCCTTCAACAACGAAGCGGAAAGAGGGAAATACCCTGAAGATTGCCACACAAAGCGTCATCAAACTTTTCTGGCGACCAAAAGGTTTAGCGTATCCCTTAGGCGAGGCTCCGGCGAAACTGCGTTGCACTTCGTCAGCTGACATGGTTACGGTCAGTAAcc                                                                                                                                                                                                                            <  1:94074/272‑1 (MQ=255)
                                               ttttTGTTGAATCCGGCAATGGGATGTTGGTTACGCTTTGTTTAAATTGGGGGGCCGTTGGCGAAAGACAGAGAAACCCTTTGCTTCACTAGCGTACGGGCAGTCCCTTCAACAACGAAGCGGAAAGAGGGAAATACCCTGAAGATTGCCACACAAAGCGTCATCAAACTTTTCTGGCGACCAAAAGGTTTAGCGTATCCCTTAGGCGAGGCTCCGGCGAAACTGCGTTGCACTTCGTCAGCTGACATGGTTACGGTCAGTAACCCAACAcc                                                                                                                                                                                                                     <  2:964552/272‑1 (MQ=255)
                                                    gCTGAATCCGGCAAATGATAATCTGGTACGCATTATGTACATTGGGGCGCCGTTGGCGAAAGACAGAGAAACCCTTTTCTTCACGAGCGTACGGGCAGTCCCTTCAACAACGAAGCGGAAAGAGGGAAATACCCTGAAGATTGCCACACAAAGCGTCATCAAACTTTTCTGGCGACCAAAAGGTTTAGCGTATCCCTTAGGCGAGGCTCCGGCGAAACTGCGTTGCACTTCGTCAGCTGACATGGTTACGGTCAGTAACCCAACACCTTAtt                                                                                                                                                                                                                <  1:593441/272‑1 (MQ=255)
                                                      tGAATCCGGCAAATGAGAATCTGTTACGCATTATGTACATTGGAGCGCCGTTGGCGAAAGACAGAGAAACCCTTTTCTTCACTAGCGTACGGGCAGTCCCTTCAACAACGAAGCGGAAAGAGGGAAATACCCTGAAGATTGCCACACAAAGCGTCATCAAACTTTTCTGGCGACCAAAAGGTTTAGCGTATCCCTTAGGCGAGGCCCCGGCGAAACTGCGTTGCACTTCGTCAGCTGACATGGTTACGGTCAGTAACCCAACACCTTAt                                                                                                                                                                                                                 >  1:593304/1‑269 (MQ=255)
                                                          tCCGGCAAATGAGAATCTGTTACGCATTATGTACATTGGAGCGCCGTTGGCGAAAGACAGAGAAACCCTTTTCTTCACTAGCGTACGGGCAGTCCCTTCAACAACGAAGCGGAAAGAGGGAAATACCCTGAAGATTGCCACACAAAGCGTCATCAAACTTTTCTGGCGACCAAAAGGTTTAGCGTATCCCTTAGGCGAGGCt                                                                                                                                                                                                                                                                                >  1:1013160/1‑202 (MQ=255)
                                                          tCCGGCAAATGAGAATCTGTTACGCATTATGTACATTGGAGCGCCGTTGGCGAAAGACAGAGAAACCCTTTTCTTCACTAGCGTACGGGCAGTCCCTTCAACAACGAAGCGGAAAGAGGGAAATACCCTGAAGATTGCCACACAAAGCGTCATCAAACTTTTCTGGCGACCAAAAGGTTTAGCGTATCCCTTAGGCGAGGCTCCGGCGAAACTGCGTTGCACTTCGt                                                                                                                                                                                                                                                       >  1:501014/1‑227 (MQ=255)
                                                          tCCGGCAAATGAGAATCTGTTACGCATTATGTACATTGGAGCGCCGTTGGCGAAAGACAGAGAAACCCTTTTCTTCACTAGCGTACGGGCAGTCCCTTCAACAACGAAGCGGAAAGAGGGAAATACCCTGAAGATTGCCACACAAAGCGTCATCAAACTTTTCTGGCGACCAAAAGGTTTAGCGTATCCCTTAGGCGAGGCTCCGGCGAAACTGCGTTGCACTTCGt                                                                                                                                                                                                                                                       <  2:501013/227‑1 (MQ=255)
                                                           ccGGCAAATGGGAATCTGTTACGCATTATGTACATTGGAGCGCCGTTGGCGAAAGACAGAGAAACCCTTTTCTTCACTAGCGTACGGGCAGTCCCTTCAACAACGAAGCGGAAAGAGGGAAATACCCTGAAGATTGCCACACAAAGCGTCATCAAACTTTTCTGGCGACCAAAAGGTTTAGCGTATCCCTTAGGCGAGGCTc                                                                                                                                                                                                                                                                               <  2:1013159/202‑1 (MQ=255)
                                                                                         gTACATTGGAGCGCCGTTGGCGAAAGACAGAGAAACCCTTTTCTTCACTAGCGTACGGGCAGTCCCTTCAACAACGAAGCGGAAAGAGGGAAATACCCTGAAGATTGCCACACAAAGCGTCATCAAACTTTTCTGGCGACCAAAAGGTTTAGCGTATCCCTTAGGCGAGGCTCCGGCGAAACTGCGTTGCACTTCGTCAGCt                                                                                                                                                                                                                                                  >  1:491137/1‑202 (MQ=255)
                                                                                          tACATTGGAGCGCCGTTGGCGAAAGACAGAGAAACCCTTTTCTTCACTAGCGTACGGGCAGTCCCTTCAACAACGAAGCGGAAAGAGGGAAATACCCTGAAGATTGCCACACAAAGCGTCATCAAACTTTTCTGGCGACCAAAAGGTTTAGCGTATCCCTTAGGCGAGGCTCCGGCGAAACTGCGTTGCACTTCGTCAGCTg                                                                                                                                                                                                                                                 <  2:491136/202‑1 (MQ=255)
                                                                                                ggAGCGCCGTTGGCGAAAGACAGAGAAACCCTTTTCTTCACTAGCGTACGGGCAGTCCCTTCAACAACGAAGCGGAAAGAGGGAAATACCCTGAAGATTGCCACACAAAGCGTCATCAAACTTTTCTGGCGACCAAAAGGTTTAGCGTATCCCTTAGGCGAGGCCCCGGCGAAACTGCGTTGCACTTCGTCAGCTGACATGGTTCCGGTCAGTAACCCAACCACTTATTTCATTACCCCGACAAACCTAAAAAAAGGCGGGAAAGAGTTaaa                                                                                                                                                                    >  2:994835/1‑272 (MQ=255)
                                                                                                       cGTTGGCGAAAGACAGAGAAACCCTTTTCTTCACTAGCGTACGGGCAGTCCCTTCAACAACGAAGCGGAAAGAGGGAAATACCCTGAAGATTGCCACACAAAGCGTCATCAAACTTTTCTGGCGACCAAAAGGTTTAGCGTATCCCTTAGGCGAGGCTCCGGc                                                                                                                                                                                                                                                                           >  2:40133/1‑163 (MQ=255)
                                                                                                        gTTGGCGAAAGACAGAGAAACCCTTTTCTTCACTAGCGTACGGGCAGTCCCTTCAACAACGAAGCGGAAAGAGGGAAATACCCTGAAGATTGCCACACAAAGCGTCATCAAACTTTTCTGGCGACCAAAAGGTTTAGCGTATCCCTTAGGCGAGGCTCCGGCg                                                                                                                                                                                                                                                                          <  1:40133/163‑1 (MQ=255)
                                                                                                            gCGAAAGACAGAGAAACCCTTTTCTTCACTAGCGTACGGGCAGTCCCTTCAACAACGAAGCGGAAAGAGGGAAATACCCTGAAGATTGCCACACAAAGCGTCATCAAACTTTTCTGGCGACCAAAAGGTTTAGCGTATCCCTTAGGCGAGGCTCCGGCGAAACTGCGTTGCACTTCGTCAGCTGACATGGTTACGGTCAGTAACCCAacac                                                                                                                                                                                                                      >  2:964916/1‑211 (MQ=255)
                                                                                                             cGAAAGACAGAGAAACCCTTTTCTTCACTAGCGTACGGGCAGTCCCTTCAACAACGAAGCGGAAAGAGGGAAATACCCTGAAGATTGCCACACAAAGCGTCATCAAACTTTTCTGGCGACCAAAAGGTTTAGCGTATCCCTTAGGCGAGGCTCCGGCGAAACTGCGTTGCACTTCGTCAGCTGACATGGTTACGGTCAGTAACCCAACAcc                                                                                                                                                                                                                     <  1:964917/211‑1 (MQ=255)
                                                                                                                 aaaaagggaAACCCCTTTCTTACAGAGCGTACGGGAAGTCCCTTCAACAACGAAGCAGAAAGAGGGAAATACCCTGAAGATTGCCACACAAAGCGTCATCAAACCGTTCTGGCGACCAAAAGGTTTAGCGTATCCCTTAGGCGAGGCTCCGGCGAAACTGCGTTGCACTTCGTCAGCTGACATGGTTACGGTCAGTAACCCAACACCTTATTTCATTACCCTGACAGACCTGAAAATAGGTGGAAAAGTAGTTAAAAATCAAATGATTTccc                                                                                                                                                    <  2:802489/265‑1 (MQ=255)
                                                                                                                  agcAGAGAAACCCTTTTCTGTACTAGCGTACGGGCAGTCCCATAAACAACGAAGCGGAAAGAGGGAAATACCCTGAAGATTGCCACACAAAGCGTCATCAAACTTTTCTGTCGACCAAAAGGTTTAGCGTATCCCTTAGGCGAGGCCCCGGCGAAACTGCGTTGCACTTCGTCAGCTGACATGGTTACGGTCAGTAACCCAACACCTTATTTCATTACCCTGACAGACCTGAAAATAGGTGGAAAAGTAGTTAAAAATCAAATGATTTccc                                                                                                                                                    <  1:994836/269‑1 (MQ=255)
                                                                                                                                    ttcACTAGCGTACGGGCAGTCCCTTCAACAACGAAGCGGAAAGAGGGAAATACCCTGAAGATTGCCACACAAAGCGTCATCAAACTTTTCTGGCGACCAAAAGGTTTAGCGTATCCCTTAGGCGAGGCTCCGGCGAAACTGCGTTGCACTTCGTCAGCTGACATGGTTACGGTCAGTAACCCAACACCTTATTTCATTACCCTGACAGACCTGAAAATAGGTGGAAAAGAAGTTAAAAATCAAATGATTTCCCCCTTTGATAAATACCAAtt                                                                                                                                 >  1:1099058/1‑272 (MQ=255)
                                                                                                                                     tcACTAGCGTACGGGAAGTCCCTACAACAACGAAGCGGAAAGAGGGAAATACCCTAAAGATTGCCAACAAAAGCTTCATCAAACTTTTCTGGCGACCAAAAGGTTTAGCGTATCCCTTAGGCGAGGCTCCGGCGAAACTGCGTTGCACTTCGTCAGCTGACATGGTTACGGTCAGTAACCCAACCCCTTATTTCATTACCCTGACAGACCTGAAAATAGGTGGAAAAGTAGTTAAAAATCAAATGATTTCCCCCTTTGATAAATACCAAttt                                                                                                                                <  1:669184/272‑1 (MQ=255)
                                                                                                                                            cGTACGGGCAGTCCCTTCAACAACGAAGCGGAAAGAGGGAAATACCCTGAAGATTGCCACACAAAGCGTCATCAAACTTTTCTGGCGACCAAAAGGTTTAGCGTATCCCTTAGGCGAGGCTCCGGCGAAACTGCGTTGCACTTCGTCAGCTGACATGGTTACGGTCAGTAAc                                                                                                                                                                                                                             >  2:908971/1‑172 (MQ=255)
                                                                                                                                            cGTACGGGCAGTCCCTTCAACAACGAAGCGGAAAGAGGGAAATACCCTGAAGATTGCCACACAAAGCGTCATCAAACTTTTCTGGCGACCAAAAGGTTTAGCGTATCCCTTAGGCGAGGCTCCGGCGAAACTGCGTTGCACTTCGTCAGCTGACATGGTTACGGTCAGTAACCCAACACCTTATTTCATTACCCTGACAGACCTGAAAATAGGTGGAAAAGTAGTTAAAAATCAAATGATTTCACCCTTTTAAAAATACCaaaatttcctgt                                                                                                                         >  2:917802/1‑262 (MQ=255)
                                                                                                                                             gTACGGGCAGTCCCTTCAACAACGAAGCGGAAAGAGGGAAATACCCTGAAGATTGCCACACAAAGCGTCATCAAACTTTTCTGGCGACCAAAAGGTTTAGCGTATCCCTTAGGCGAGGCTCCGGCGAAACTGCGTTGCACTTCGTCAGCTGACATGGTTACGGTCAGTAAcc                                                                                                                                                                                                                            <  1:908972/172‑1 (MQ=255)
                                                                                                                                                         ccTTCAACAACGAAGCGGAAAGAGGGAAATACCCTGAAGATTGCCACACAAAGCGTCATCAAACTTTTCTGGCGACCAAAAGGTTTAGCGTATCCCTTAGGCGAGGCTCCGGCGAAACTGCGTTGCACTTCGTCAGCt                                                                                                                                                                                                                                                  >  2:232639/1‑138 (MQ=255)
                                                                                                                                                          cTTCAACAACGAAGCGGAAAGAGGGAAATACCCTGAAGATTGCCACACAAAGCGTCATCAAACTTTTCTGGCGACCAAAAGGTTTAGCGTATCCCTTAGGCGAGGCTCCGGCGAAACTGCGTTGCACTTCGTCAGCTg                                                                                                                                                                                                                                                 <  1:232639/138‑1 (MQ=255)
                                                                                                                                                          cTTCAACAACGAAGCGGAAAGAGGGAAATACCCTGAAGATTGCCACACAAAGCGTCATCAAACTTTTCTGGCGACCAAAAGGTTTAGCGTATCCCTTAGGCGAGGCTCCGGCGAAACTGCGTTGCACTTCGTCAGCTGACATGGTTACgg                                                                                                                                                                                                                                     >  1:956683/1‑150 (MQ=255)
                                                                                                                                                          cTTCAACAACGAAGCGGAAAGAGGGAAATACCCTGAAGATTGCCACACAAAGCGTCATCAAACTTTTCTGGCGACCAAAAGGTTTAGCGTATCCCTTAGGCGAGGCTCCGGCGAAACTGCGTTGCACTTCGTCAGCTGACATGGTTACGGTCAGTAACCCAACACCTTATTTCATTACCCTGACAGACCTGAAAATAGGTGGAAAAGTAGTTAAAAATCAAATGATTTCCCCCTTTGATAAATACCAATTTTCTCTGCCAAAGGGGGCCaaa                                                                                                           >  1:795436/1‑272 (MQ=255)
                                                                                                                                                           ttCAACAACGAAGCGGAAAGAGGGAAATACCCTGAAGATTGCCACACAAAGCGTCATCAAACTTTTCTGGCGACCAAAAGGTTTAGCGTATCCCTTAGGCGAGGCTCCGGCGAAACTGCGTTGCACTTCGTCAGCTGACATGGTTACGGt                                                                                                                                                                                                                                    <  2:956682/150‑1 (MQ=255)
                                                                                                                                                            tCAACAACGAAGCGGAAAGAGGGAAATACCCTGAAGATTGCCACACAAAGCGTCATCAAACTTTTCTGGCGACCAAAAGGTTTAGCGTATCCCTTAGGCGAGGCTCCGGCGAAACTGCGTTGCACTTCGTCAGCTGACATGGTTACGGTCAGTAACCCAACACCTTATTTCATTACCCTGACAGACCTGAAAATAGGTGGAAAAGTAGTTAAAAATCAAATGATTTCCCCCTTTGATAAATACCAATTTTCTCTGCCAAAGGGGGCCaaaaa                                                                                                         >  1:378413/1‑272 (MQ=255)
                                                                                                                                                                caacGAAGCGGAAAGAGGGAAATACCCTGAAGATTGCCACACAAAGCGTCATCAAACTTTTCTGGCGACCAAAAGGTTTAGCGTATCCCTTAGGCGAGGCTCCGGCGAACCTGCGTTGCACTTCGTCAGCTGACATGGTTACGGTCAGTAACCCAACACCTTATTTCATTACCCTGACAGACCTGAAAATAGGTGGAAAAGTAGTTAAAAATCAAATAATTTCCCCCTTTGATAAATACCAATTTTCTCTGCCAAAGGGGGCCAAAAATagc                                                                                                     >  1:247057/1‑272 (MQ=255)
                                                                                                                                                                           aaaGAGGGAAATACCCTGAAGATTGCCACACAAAGCGTCATCAAACTTTTCTGGCGACCACAAGGTTTAGCGTATCCCTTAGGCGAGGCTCCGGCGAAACTGCGTTGCACTTCGTCAGCTGACATGGTTACGGTCAGTAAc                                                                                                                                                                                                                             >  1:1094674/1‑141 (MQ=255)
                                                                                                                                                                            aaGAGGGAAATACCCTGAAGATTGCCACACAAAGCGTCATCAAACTTTTCTGGCGACCAAAAGGTTTAGCGTATCCCTTAGGCGAGGCTCCGGCGAAACTGCGTTGCACTTCGTCAGCTGACATGGTTACGGTCAGTAAcc                                                                                                                                                                                                                            <  2:1094673/141‑1 (MQ=255)
                                                                                                                                                                                            gAAGATTGCCACACAAAGCGTCATCAAACTTTTCTGGCGACCAAAAGGTTTAGCGTATCCCTTAGGCGAGGCTCCGGCGAAACTGCGTTGCACTTCGTCAGCTGACATGGTTACGGt                                                                                                                                                                                                                                    >  2:172416/1‑117 (MQ=255)
                                                                                                                                                                                             aaGATTGCCACACAAAGCGTCATCAAACTTTTCTGGCGACCAAAAGGTTTAGCGTATCCCTTAGGCGAGGCTCCGGCGAAACTGCGTTGCACTTCGTCAGCTGACATGGTTACGGTc                                                                                                                                                                                                                                   <  1:172416/117‑1 (MQ=255)
                                                                                                                                                                                             aaGATTGCCACACAAAGCGTCATCAAACTTTTCTGGCGACCAAAAGGTTTAGCGTATCCCTTAGGCGAGGCTCCGGCGAAACTGCGTTGCACTTCGTCAGCTGACATGGTTACGGTCAGTAACCCAACACCTTATTTCATTACCCTGACAGACCTGAAAATAGGTGGAAAAGTAGTTAAAAATCAAATGATTTCCCCCTTTGATAAATACCAATTTTCTCTGCCAAAGGGGGCCAAAAATAGCAGCGTAACGTATCGAACCATCAATGAca                                                                         >  1:1032120/1‑270 (MQ=255)
                                                                                                                                                                                              aGATTGCCACACAAAGCGTCATCAAACTTTTCTGGCGACCAAAAGGTTTAGCGTATCCCTTAGGCGAGGCTCCGGc                                                                                                                                                                                                                                                                           >  2:151856/1‑76 (MQ=255)
                                                                                                                                                                                               gATTGCCACACAAAGCGTCATCAAACTTTTCTGGCGACCAAAAGGTTTAGCGTATCCCTTAGGCGAGGCTCCGGCg                                                                                                                                                                                                                                                                          <  1:151856/76‑1 (MQ=255)
                                                                                                                                                                                                aTTGCCACACAAAGCGTCATCAAACTTTTCTGGCGACCAAAAGGTTTAGCGTATCCCTTAGGCGAGGCTCCGGCGAAACTGCGTTGCACTTCGTCAGCTGACATGGTTACGGTCAGTAACCCAACACCTTATTTCATTACCCt                                                                                                                                                                                                      >  2:681049/1‑143 (MQ=255)
                                                                                                                                                                                                aTTGCCACACAAAGCGTCATCAAACTTTTCTGGCGACCAAAAGGTTTAGCGTATCCCTTAGGCGAGGCTCCGGCGAAACTGCGTTGCACTTAGTCAGCTGACATGGTTACGGTCAGTAACCCAACACCTTATTTCATTACCCTGACAGACCTGAAAATAGGTGGAAAAGTAGTTAAAAATCAAATGATTTCCCCCTTTGATAAATACCAATTTTCTCTGCCAAAGGGGGCCaaaa                                                                                                          >  1:170859/1‑235 (MQ=255)
                                                                                                                                                                                                 ttGCCACACAAAGCGTCATCAAACTTTTCTGGCGACCAAAAGGTTTAGCGTATCCCTTAGGCGAGGCTCCGGCGAAACTGCGTTGCACTTCGTCAGCTGACATGGTTACGGTCAGTAACCCAACACCTTATTTCATTACCCTg                                                                                                                                                                                                     <  1:681050/143‑1 (MQ=255)
                                                                                                                                                                                                 ttGCCACACAAAGCGTCATCAAACTTTTCTGGCGACCAAAAGGTTTAGCGTATCCCTTAGGCGAGGCTCCGGCGAAACTGCGTTGCACTTAGTCAGCTGACATGGTTACGGTCAGTAACCCAACACCTTATTTCATTACCCTGACAGACCTGAAAATAGGTGGAAAAGTAGTTAAAAATCAAATGATTTCCCCCTTTGATAAATACCAATTTTCTCTGCCAAAGGGGGCCaaaaa                                                                                                         <  2:170859/235‑1 (MQ=255)
                                                                                                                                                                                                               gTCATCAAACTTTTCTGGCGACCAAAAGGTTTAGCGTATCCCTTAGGCGAGGCTCCGGCGAAACTGCGTTGCACTTCGTCAGCTGACATGGTTACGGTCAGTAACCCAACACCTTATTTCATTACCCTGACAGACCTGAAAATAGGTGGAAAAGTAGTTaaaaa                                                                                                                                                                  >  2:6116/1‑164 (MQ=255)
                                                                                                                                                                                                               gTCATCAAACTTTTCTGGCGACCAAAAGGTTTAGCGTATCCCTTAGGCGAGGCTCCGGCGAAACTGCGTTGCACTTCGTCAGCTGACATGGTTACGGTCAGTAACCCAACACCTTATTTCATTACCCTGACAGACCTGAAAATAGGTGGAAAAGTAGTTAAAAATCAAATGATTTCCCCCTTTGATAAATACCAATTTTCTCTGCCAAAGGGGGCCAAAAATAGCAGCGTAACg                                                                                            >  1:666079/1‑234 (MQ=255)
                                                                                                                                                                                                               gTCATCAAACTTTTCTGGCGACCAAAAGGTTTAGCGTATCCCTTAGGCGAGGCTCCGGCGAAACTGCGTTGCACTTCGTCAGCTGACATGGTTACGGTCAGTAACCCAACACCTTATTTCATTACCCTGACAGACCTGAAAATAGGTGGAAAAGTAGTTAAAAATCAAATGATTTCCCCCTTTGATAAATACCAATTTTCTCTGCCAAAGGGGGCCAAAAATAGCAGCGTAACg                                                                                            <  2:666078/234‑1 (MQ=255)
                                                                                                                                                                                                                ttttaataCTGTTCTGGCGACCAAAAGGTTTGGCTTATCCTTTAGGCGAGGCTCCGGCGAAACTGCGTTGCACTTCGTCAGCTGACATGGTTACGGTCAGTAACCCAACACCTTATTTCATTACCCTGACAGACCTGAAAATAGGTGGAAAAGTAGTTAAAAATCAAATGATTTCCCCCTTTGATAAATACCAATTTTCTCTGCCAAAGGGGGCCAAAAATAGCAGCGTAACGTATCGAACCATCAATGACTACGGGGCGGAAACGCCGcaa                                                     <  2:961407/265‑1 (MQ=255)
                                                                                                                                                                                                                tcatcaAACTTTTCTGGCGACCAAAAGGTTTAGCGTATCCCTTAGGCGAGGCTCCGGCGAAACTGCGTTGCACTTCGTCAGCTGACATGGTTACGGTCAGTAACCCAACACCTTATTTCATTACCCTGACAGACCTGAAAATAGGTGGAAAAGTAGTTAAAAAt                                                                                                                                                                 <  1:6116/164‑1 (MQ=255)
                                                                                                                                                                                                                 catcaAACTTTTCTGGCGACCAAAAGGTTTAGCGTATCCCTTAGGCGAGGCTCCGGCGAAACTGCGTTGCACTTCGTCAGCt                                                                                                                                                                                                                                                  >  2:352938/1‑82 (MQ=255)
                                                                                                                                                                                                                  atcaAACTTTTCTGGCGACCAAAAGGTTTAGCGTATCCCTTAGGCGAGGCTCCGGCGAAACTGCGTTGCACTTCGTCAGCTg                                                                                                                                                                                                                                                 <  1:352939/82‑1 (MQ=255)
                                                                                                                                                                                                                     aaaCTTTTCTGGCGACCAAAAGGTTTAGCGTATCCCTTAGGCGAGGCTCCGGCGAAACTGCGTTGCACTTCGTCAGCTGACATGGt                                                                                                                                                                                                                                          >  2:90438/1‑86 (MQ=255)
                                                                                                                                                                                                                     aaaCTTTTCTGGCGACCAAAAGGTTTAGCGTATCCCTTAGGCGAGGCTCCGGCGAAACTGCGTTGCACTTCGTCAGCTGACATGGTTACGGTCAGTAACCCAACACCTTATTTCATTACCCTGACAGACCTGAAAATAGGTGGAAAAGTAGTTAAAAATCAAATGATTTCCCCCTTTGATAAATACCAATTTTCTCTGCCAAAGGGGGCCAAAAATAGCAGCGTAACGTATCGAACCATCAATGACTACGGGGCGGAAACGCCGcaactca                                                 >  1:961408/1‑271 (MQ=255)
                                                                                                                                                                                                                      aaCTTTTCTGGCGACCAAAAGGTTTAGCGTATCCCTTAGGCGAGGCTCCGGCGAAACTGCGTTGCACTTCGTCAGCTGACATGGtt                                                                                                                                                                                                                                         <  1:90438/86‑1 (MQ=255)
                                                                                                                                                                                                                          tttCTGGCGACCAAAAGGATTAGCGTATCCCTTAGGCGAGGCTCCGGCGAAACTGCGTTGCACTTCTTCAGCTGACATGGTTACGGTCCGTAACCCAACACCTTATTTCATTACCCTGACAGACCTGAAACTAGGTGGAAAAGTAGTTAAAAATAAAATGATTTCCCCCTTTGCTAAATACCAAATTTCTCTGCCAAATGGGGCCAAAAATAGCCACGTAACGTAACGAACCACAAAAAACTACAGGGCGGAAAGGCACAAACTAAACTGt                                            >  1:954018/1‑271 (MQ=255)
                                                                                                                                                                                                                             cTGGCGACCAAAAGGTTTAGCGTATCCCTTAGGCGAGGCTCCGGCGAAACTGCGTTGCACTTCGTCAGCTGACATGGTTACGGTCAGTAACCCAACACCTTATTTCATTACCCTGACAGACCTGAAAATAGGTGGAAAAGTAGTTAAAAATCAAATGATTTCCCCCTTTGATAAATACCAATTTTCTCTGCCAAAGGGGGCCAAAAATAGCAGCGTAACGTATCGAACCATCAATGACTACGGGGCGGAAACGCCGCAACTCAACTGTAAAt                                        <  1:917803/272‑1 (MQ=255)
                                                                                                                                                                                                                               ggCGACCAAAAGGTTTAGCGTATCCCTTAGGCGAGGCTCCGGCGAAACTGCGTTGCACTTCGTCAGCTGACATGGTTACGGt                                                                                                                                                                                                                                    >  1:678497/1‑82 (MQ=255)
                                                                                                                                                                                                                                gCGACCAAAAGGTTTAGCGTATCCCTTAGGCGAGGCTCCGGCGAAACTGCGTTGCACTTCGTCAGCTGACATGGTTACGGTc                                                                                                                                                                                                                                   <  2:678496/82‑1 (MQ=255)
                                                                                                                                                                                                                                      aaaaGGTTTAGCGTATCCCTTAGGCGAGGCTCCGGCGAAACTGCGTTGCACTTCGTCAGCTGACATGGTTACGGTCAGTAACCCAACACCTTATTTCATTACCCTGACAGACCTGAAAATAGGTGGAAAAGTAGTTAAAAATCAAATGATTTCCCCCTTTGATAAATACCAATTTTCTCTGCCAAAGGGGGCCAAAAATAGCAGCGTAACGTATCGAACCATCAATGACTACGGGGCGGAAACGCCGCAACTCAACTGTAAATCGTAAGCCg                               <  1:47002/272‑1 (MQ=255)
                                                                                                                                                                                                                                        aaGGTTTAGCGTATCCCTTAGGCGAGGCTCCGGCGAAAATGCGTTGCACTTCGTCAGCTGACATg                                                                                                                                                                                                                                            >  1:689745/1‑65 (MQ=255)
                                                                                                                                                                                                                                         aGGTTTAGCGTATCCCTTAGGCGAGGCTCCGGCGAAACTGCGTTGCACTTCGTCAGCTGACATgg                                                                                                                                                                                                                                           <  2:689744/65‑1 (MQ=255)
                                                                                                                                                                                                                                                 cGTATCCCTTAGGCGAGGCTCCGGCGAAACTGCGTTGCACTTCGTCAGCTGACATGGTTACGGTCAGTAACCCAacac                                                                                                                                                                                                                      >  1:800924/1‑78 (MQ=255)
                                                                                                                                                                                                                                                  gTATCCCTTAGGCGAGGCTCCGGCGAAACTGCGTTGCACTTCGTCAGCTGACATGGTTACGGTCAGTAACCCAACAcc                                                                                                                                                                                                                     <  2:800923/78‑1 (MQ=255)
                                                                                                                                                                                                                                                       ctttcttgtagtcTCCGGGGAAATTAGGTTGCACTTCGTCTGCTGACATGGTTACTGTCAGTAACCCAACACCTAATTTCATTACCCTGTCAGACCTGAAAAGAGGTGGAAAAGTAGTTAAAAATCAAAGGATTTCCCCCTTTGATAAATACCAATTTTCTCTGCCAAAGGGGGCCAAAAATAGCAGCGTAACGTATCGAACCATCAATGACTACGGGGCGGAAACGCCGCAACTCAACTGTAAATCGTAAGCCGTCTTCAGTTAAGAGAGc              <  2:954017/260‑1 (MQ=255)
                                                                                                                                                                                                                                                       ccTTAGGCAGGTTTCCGGGGAAACTTCGTGGCACTTCTGCAGCTGACATGGTTACGGTCAGGAACCCAACACCTTATTTCATTACCCTTACAGACCTGAAAATAGGTGGAAAAGTAGTTAAAAATCAAATGATTTCCCCCTTTGATAAATACCAATTTTCTCTGCCAAAGGGGGCCAAAAATAGCAGCGTAACGTATCGAACCATCAATGACTACGGGGCGGAAACGCCGCAACTCAACTGTAAATCGTAAGCCGTCTTCAGTTAAGAGAGc              <  2:795435/272‑1 (MQ=255)
                                                                                                                                                                                                                                                        ctgaggggtgtcTCCGGCGAAACGGCGTTGCACTTCGTCAGCTGACATGGTTACGGTCAATAACCCAACACCTTATTTCATTACCCTGACAGACCTGAAAATAGGTGGAAAAGTAGTTAAAAATCAAATGATTTCCCCCTTTGATAAATACCAATTTTCTCTGCCAAAGGGGGCCAAAAATAGCAGCGTAACGTATCGAACCATCAATGACTACGGGGCGGAAACGCCGCAACTCAACTGTAAATCGTAAGCCGTCTTCAGTTAAGAGAGc              <  2:378412/260‑1 (MQ=255)
                                                                                                                                                                                                                                                        cTTAGGCGAGGCTCCGGCGAAACTGCGTTGCACTTCGTCAGCTGACATGGTTACGGTCAGTAAc                                                                                                                                                                                                                             >  1:834177/1‑64 (MQ=255)
                                                                                                                                                                                                                                                         ttAGGCGAGGCTCCGGCGAAACTGCGTTGCCCTTCGTCAGCTGACATGGTTACGGTCAGTAACCCAACACCTTATTTCATTACCCTGACAGACCTGAAAAGAGGTGGAAAAGTAGTTAAAAATCAAATGATTTCCCCCTTTGATAAATACAAATTTTCTCTGCCAAAGGGGGCCAAAAATAGCAGCGTAACGTATCGAACCATCAATGACTACGGGGCGGAAACGCCGCAACTCAACTGTAAATCGTAAGCCGTCTTCAGTTAAGAGAGc              <  1:11473/270‑1 (MQ=255)
                                                                                                                                                                                                                                                         ttAGGCGAGGCTCCGGCGAAACTGCGTTGCACTTCGTCAGCTGACATGGTTACGGTCAGTAAcc                                                                                                                                                                                                                            <  2:834176/64‑1 (MQ=255)
                                                                                                                                                                                                                                                             gCGAGGCTCCGGCGAAACTGCGTTGCACTTCGTCAGCTGACATGGTTACGGTCAGTAACCCAACACCTTATTTCATTACCCTGACAGACCTGAAAATAGGTGGAAAAGTAGTTAAAAATCAAATGATTTCCCCCTTTGATAAATACCAATTTTCTCTGCCAAAGGGGGCCAAAAATagcag                                                                                                   >  1:447117/1‑181 (MQ=255)
                                                                                                                                                                                                                                                             gCGAGGCTCCGGCGAAACTGCGTTGCACTTCGTCAGCTGACATGGTTACGGTCAGTAACCCAACACCTTATTTCATTACCCTGACAGACCTGAAAATAGGTGGAAAAGTAGTTAAAAATCAAATGATTTCCCCCTTTGATAAATACCAATTTTCTCTGCCAAAGGGGGCCAAAAATAGCAGCGTAACGTATCGAACCATCAATGACTACGGGGCGGAAACGCCGcaactcaa                                                >  2:885064/1‑232 (MQ=255)
                                                                                                                                                                                                                                                              cGAGGCTCCGGCGAAACTGCGTTGCACTTCGTCAGCTGACATGGTTACGGTCAGTAACCCAACACCTTATTTCATTACCCTGACAGACCTGAAAATAGGTGGAAAAGTAGTTAAAAATCAAATGATTTCCCCCTTTGATAAATACCAATTTTCTCTGCCAAAGGGGGCCAAAAATagcagc                                                                                                  <  2:447116/181‑1 (MQ=255)
                                                                                                                                                                                                                                                              cGAGGCTCCGGCGAAACTGCGTTGCACTTCGTCAGCTGACATGGTTACGGTCAGTAACCCAACACCTTATTTCATTACCCTGACAGACCTGAAAATAGGTGGAAAAGTAGTTAAAAATCAAATGATTTCCCCCTTTGATAAATACCAATTTTCTCTGCCAAAGGGGGCCAAAAATAGCAGCGTAACGTATCGAACCATCAATGACTACGGGGCGGAAACGCCGcaactcaac                                               <  1:885065/232‑1 (MQ=255)
                                                                                                                                                                                                                                                              cGAGGCCCCGGCGAAACTGCGTTGCACTTCGTCAGCTGACATGGTTACGGTCAGTAACCCAACACCTTATTTCATTACCCTGACAGACCTGaaaa                                                                                                                                                                                        >  2:709095/1‑95 (MQ=255)
                                                                                                                                                                                                                                                               gAGGCCCCGGCGAAACTGCGTTGCACTTCGTCAGCTGACATGGTTACGGTCAGTAACCCAACACCTTATTTCATTACCCTGACATACCTGAAAAt                                                                                                                                                                                       <  1:709096/95‑1 (MQ=255)
                                                                                                                                                                                                                                                                   cTCCGGCGAAACTGCGTTGCACTTCGTCAGCTGACATGGTTACGGTCAGTAACCCAACACCTTAt                                                                                                                                                                                                                 >  2:857927/1‑65 (MQ=255)
                                                                                                                                                                                                                                                                   cTCCGGCGAAACTGCGTTGCACTTCGTCAGCTGACATGGTTACGGTCAGTAACCCAACACCTTATTTCATTACCCTGACAGACCTGAAAATAGGTGGAAAAGTAGTTAAAAATCAAATGATTTCCCCCTTTGATAAATACCAATTTTCTCTGCCAAAGGGGGCCAAAAATAGCAGCGTAACGTATCGAACCATCAATGACTACGGGGCGGAAACGCCGCAACCCAACCGTAAAACGAAAGCCGTCATCAGGAAAAAgagagaacagtaaaaa  >  2:406760/1‑262 (MQ=255)
                                                                                                                                                                                                                                                                    tCCGGCGAAACTGCGTTGCACTTCGTCAGCTGACATGGTTACGGTCAGTAACCCAACACCTTAtt                                                                                                                                                                                                                <  1:857928/65‑1 (MQ=255)
                                                                                                                                                                                                                                                                    tCCGGCGAAACTGCGTTGCACTTCGTCAGCTGACATGGTTACGGTCAGTAACCCAACACCTTATTTCATTACCCTGACAGACCTGAAAATAGGTGGAAAAGTAGTTAAAAATCAAATg                                                                                                                                                           >  1:168420/1‑118 (MQ=255)
                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                              
CCCCAGATAATAAAAAATCAGCAGACTTTATTATTACACCACCGCTATTTGTGCTGAATCCGGCAAATGAGAATCTGTTACGCATTATGTACATTGGAGCGCCGTTGGCGAAAGACAGAGAAACCCTTTTCTTCACTAGCGTACGGGCAGTCCCTTCAACAACGAAGCGGAAAGAGGGAAATACCCTGAAGATTGCCACACAAAGCGTCATCAAACTTTTCTGGCGACCAAAAGGTTTAGCGTATCCCTTAGGCGAGGCTCCGGCGAAACTGCGTTGCACTTCGTCAGCTGACATGGTTACGGTCAGTAACCCAACACCTTATTTCATTACCCTGACAGACCTGAAAATAGGTGGAAAAGTAGTTAAAAATCAAATGATTTCCCCCTTTGATAAATACCAATTTTCTCTGCCAAAGGGGGCCAAAAATAGCAGCGTAACGTATCGAACCATCAATGACTACGGGGCGGAAACGCCGCAACTCAACTGTAAATCGTAAGCCGTCTTCAGTTAAGAGAGCGAGATGTATAGA  >  NC_000913/998695‑999224

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