Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,000,101 G→T 5.7% A296S (GCC→TCC)  elfC → putative outer membrane fimbrial subunit export usher protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,000,1010GT 5.7% 246.0 / 5.2 70A296S (GCC→TCC) elfCputative outer membrane fimbrial subunit export usher protein
Reads supporting (aligned to +/- strand):  ref base G (37/29);  new base T (2/2);  total (39/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.00e-01

GAACTGACCGTAGGTGATGATTATACTGCAGGCGATTTTTTCGATAGTGTCAGCTTTCGTGGTGTGCAGCTGGCGTCAGATGACAACATGCTGCCAGACAGCCTGAAAGGGTTTGCGCCTGTGGTGCGTGGTATCGCCAAAAGCAATGCCCAGATAACGATTAAGCAAAATGGTTACACCATTTACCAAACTTATGTATCGCCTGGTGCTTTTGAAATTAGTGATCTCTATTCCACGTCGTCGAGCGGTGATTTGTTAGTTGAAATCAAAGAAGCGGACGGCAGCGTCAATAGCTA  >  NC_000913/999954‑1000249
                                                                                                                                                   |                                                                                                                                                    
gAACTGACCGTAGGTGATGATTATACTGCAGGCGATTTTTTCGATAGTGTCAGCTTTCGTGGTGTGCAGCTGGCGTCAGATGACAACATGCTGCCAGACAGCCTGAAAGGGTTTGCGCCTGTGGTGCGTGGTATCGCCAAAAGCAATGc                                                                                                                                                     >  1:304266/1‑149 (MQ=255)
 aaCTGACCGTAGGTGATGATTATACTGCAGGCGATTTTTTCGATAGTGTCAGCTTTCGTGGTGTGCAGCTGGCGTCAGATGACAACATGCTGCCAGACAGCCTGAAAGGGTTTGCGCCTGTGGTGCGTGGTATCGCCAAAAGCAATGcc                                                                                                                                                    >  2:159798/1‑149 (MQ=255)
    tGACCGTAGGTGATGATTATACTGCAGGCGATTTTTTCGGTAGTGTCAGCTTTCGTAGTGTGAAGCTGGCGTCACATGACAACATGCTGCCAGACAGCCTGAAAGGGTTTGCGCCTGCGGTGCGTGGTATCGCCTAAAGCTATGCCCAg                                                                                                                                                 >  3:350428/1‑149 (MQ=255)
              tgatgaTTATACTGCAGGCGATTTTTTCGATAGTGTCAGCTTTCGTGGTGTGCAGCTGGCGTCAGATGACAACATGCTGCCAGACAGCCTGAAAGGGTTTGCGCCTGTGGTGCGTGGTATCGCCAAAAGCAATGCCCAGATAACGATTa                                                                                                                                       >  4:249665/1‑149 (MQ=255)
               gatgatTATACTGCAGGCGATTTTTTCGATAGTGTCAGCTTTCGTGGTGTGCAGCTGGCGTCAGATGACAACATGCTGCCAGACAGCCTGAAAGGGTTTGCGCCTGTGGTGCGTGGTATCGCCAAAAGCAATGCCCAGATAACGATTaa                                                                                                                                      <  4:433301/149‑1 (MQ=255)
               gatgatTATACTGCAGGCGATTTTTTCGATAGTGTCAGCTTTCGTGGTGTGCAGCTGGCGTCAGATGACAACATGCTGCCAGACAGCCTGAAAGGGTTTGCGCCTGTGGTGCGTGGTATCGCCAAAAGCAATGCCCAGATAACGATTaa                                                                                                                                      <  4:565478/149‑1 (MQ=255)
               gatgatTATACTGCAGGCGATTTTTTCGATAGTGTCAGCTTTCGTGGTGTGCAGCTGGCGTCAGATGACAACATGCTGCCAGACAGCCTGAAAGGGTTTGCGCCTGTGGTGCGTGGTATCGCCAAAAGCAATGCCCAGATAACGATTaa                                                                                                                                      <  3:448131/149‑1 (MQ=255)
                  gatTATACTGCAGGCGATTTTTTCGATAGTGTCAGCTTTCGTGGTGTGCAGCTGGCGTCAGATGACAACATGCTGCCAGACAGCCTGAAAGGGTTTGCGCCTGTGGTGCGTGGTATCGCCAAAAGCAATGCCCAGATAACGATTAAGCa                                                                                                                                   >  4:264039/1‑149 (MQ=255)
                  gatTATACTGCAGGCGATTTTTTCGATAGTGTCAGCTTTCGTGGTGTGCAGCTGGCGTCAGATGACAACATGCTGCCAGACAGCCTGAAAGGGTTTGCGCCTGTGGTGCGTGGTATCGCCAAAAGCAATGCCCAGATAACGATTAAGCa                                                                                                                                   >  1:460272/1‑149 (MQ=255)
                    ttATACTGCAGGCGATTTTTTCGATAGTGTCAGCTTTCGTGGTGTGCAGCTGGCGTCAGATGACAACATGCTGCCAGACAGCCTGAAAGGGTTTGCGCCTGTGGTGCGTGGTATCGCCAAAAGCAATGCCCAGATAACGATTAAGCaaa                                                                                                                                 >  1:337269/1‑149 (MQ=255)
                    ttATACTGCAGGCGATTTTTTCGATAGTGTCAGCTTTCGTGGTGTGCAGCTGGCGTCAGATGACAACATGCTGCCAGACAGCCTGAAAGGGTTTGCGCCTGTGGTGCGTGGTATCGCCAAAAGCAATGCCCAGATAACGATTAAGCaaa                                                                                                                                 >  3:491009/1‑149 (MQ=255)
                      ataCTGCAGGCGATTTTTTCGATAGTGTCAGCTTTCGTGGTGTGCAGCTGGCGTCAGATGACAACATGCTGCCAGACAGCCTGAAAGGGTTTGCGCCTGTGGTGCGTGGTATCGCCAAAAGCAATGCCCAGATAACGATTAAGCAAAAt                                                                                                                               <  1:391925/149‑1 (MQ=255)
                                 gATTTTTTCGATAGTGTCAGCTTTCGTGGTGTGCAGCTGGCGTCAGATGACAACATGCTGCCAGACAGCCTGAAAGGGTTTGCGCCTGTGGTGCGTGGTATCGCCAAAAGCAATGCCCAGATAACGATTAAGCAAAATGGTTACACCAt                                                                                                                    >  3:263475/1‑149 (MQ=255)
                                 gATTTTTTCGATAGTGTCAGCTTTCGTGGTGTGCAGCTGGCGTCAGATGACAACATGCTGCCAGACAGCCTGAAAGGGTTTGCGCCTGTGGTGCGTGGTATCGCCAAAAGCAATGCCCAGATAACGATTAAGCAAAATGGTTACACCAt                                                                                                                    >  1:223370/1‑149 (MQ=255)
                                  aTTTTTTCGATAGTGTCAGCTTTCGTGGTGTGCAGCTGGCGTCAGATGACAACATGCTGCCAGACAGCCTGAAAGGGTTTGCGCCTGTGGTGCGTGGTATCGCCAAAAGCAATGCCCAGATAACGATTAAGCAAAATGGTTACACCAtt                                                                                                                   <  2:390407/149‑1 (MQ=255)
                                   ttatttCGATAGTGTCAGCTTTTGTGGTGTGCAGGTGGCGTCAGATGCCAAAACGCTGCCAGACAGCCTGAGAGGCTTTGCGCCTGAGGTGCGAGGTAACGTCAAAAGCAATGCCCAGATAACGATTAAGCAAAATGGTTACACCAttt                                                                                                                  >  4:476518/4‑149 (MQ=255)
                                              gtgtCAGCTTTCGTGGTGTGCAGCTGGCGTCAGATGACAACATGCTGCCAGACAGCCTGAAAGGGTTTGCGCCTGTGGTGCGTGGTATCGCCAAAAGCAATGCCCAGATAACGATTAAGCAAAATGGTTACACCATTTACCAAACTTAt                                                                                                       <  4:550490/149‑1 (MQ=255)
                                                 tCAGCTTTCGTGGTGTGCAGCTGGCGTCAGATGACAACATGCTGCCAGACAGCCTGAAAGGGTTTGCGCCTGTGGTGCGTGGTATCGCCAAAAGCAATGCCCAGATAACGATTAAGCAAAATGGTTACACCATTTACCAAACTTAt                                                                                                       <  1:489824/146‑1 (MQ=255)
                                                 tCAGCTTTCGTGGTGTGCAGCTGGCGTCAGATGACAACATGCTGCCAGACAGCCTGAAAGGGTTTGCGCCTGTGGTGCGTGGTATCGCCAAAAGCAATGCCCAGATAACGATTAAGCAAAATGGTTACACCATTTACCAAACTTAt                                                                                                       >  2:489824/1‑146 (MQ=255)
                                                     cTTTCGTGGTGCGCAGCTGGCGTCAGATGATAACATGCTGCCAGACAGCCTGAAAGGGTTTGCGCCTGTGGTGCTTGGTCTCGCCAAAAGCAATGCCCAGATAACGATTAAGCAAAATAGTTACACCATTTTCTAAACTTATGTATCGc                                                                                                <  4:350428/149‑1 (MQ=255)
                                                      tttCGTGGTGTGCAGCTGGCGTCAGATGACAACATGCTGCCAGACAGCCTGAAAGGGTTTGCGCCTGTGGTGCGTGGTATCGCCAAAAGCAATGCCCAGATAACGATTAAGCAAAATGGTTACACCATTTACCAAACTTATGTATCGcc                                                                                               >  1:90558/1‑149 (MQ=255)
                                                          gtggtgTGCAGCTGGCGTCAGATGACAACATGCTGCCAGACAGCCTGAAAGGGTTTGCGCCTGTGGTGCGTGGTATCGCCAAAAGCAATGCCCAGATAACGATTAAGCAAAATGGTTACACCATTTACCAAACTTATGTATCGCCTGGt                                                                                           <  2:266923/149‑1 (MQ=255)
                                                          gtggtgTGCAGCTGGCGTCAGATGACAACATGCTGCCAGACAGCCTGAAAGGGTTTGCGCCTGTGGTGCGTGGTATCGCCAAAAGCAATGCCCAGATAACGATTAAGCAAAATGGTTACACCATTTACCAAACTTATGTATCGCCTGGt                                                                                           >  4:147140/1‑149 (MQ=255)
                                                             gtgtGCAGCTGGCGTCAGATGACAACATGCTGCCAGACAGCCTGAAAGGGTTTGCGCCTGTGGTGCGTGGTATCGCCAAAAGCAATGCCCAGATAACGATTAAGCAAAATGGTTACACCATTTACCAAACTTATGTATCGCCTGGTGCt                                                                                        <  1:54932/149‑1 (MQ=255)
                                                              tgtgCAGCTGGCGTCAGATGACAACATGCCGCCAGACAGCCTGAAAGGGTTTGCGCCTGTGGTGCGTGGTATCGCCAAAAGCAATGCCCAGATAACGATTAAGCAAAATGGTTACACCATTTACCAAACTTATGTATCGCCTGGTGCtt                                                                                       >  3:85634/1‑149 (MQ=255)
                                                               gtgCAGCTGGCGTCAGATGACAACATGCTGCCAGACAGCCTGAAAGGGTTTGCGCCTGTGGTGCGTGGTATCGCCAAAAGCAATGCCCAGATAACGATTAAGCAAAATGGTTACACCATTTACCAAACTTATGTATCGCCTGGTGCttt                                                                                      <  1:159798/149‑1 (MQ=255)
                                                                tgCAGCTGGCGTCAGATGACAACATGCTGCCAGACAGCCTGAAAGGGTTTGCGCCTGTGGTGCGTGGTATCGCCAAAAGCAATGCCCAGATAACGATTAAGCAAAATGGTTACACCATTTACCAAACTTATGTATCGCCTGGTGCtttt                                                                                     <  1:99243/149‑1 (MQ=255)
                                                                  cAGCTGGCGTCAGATGACAACATGCTGCCAGACAGCCTGAAAGGGTTTGCGCCTGTGGTGCGTGGTATCGCCAAAAGCAATGCCCAGATAACGATTAAGCAAAATGGTTACACCATTTACCAAACTTATGTATCGCCTGGTGCTTTTGa                                                                                   >  2:312288/1‑149 (MQ=255)
                                                                  cAGCTGGCGTCAGATGACAACATGCTGCCAGACAGCCTGAAAGGGTTTGCGCCTGTGGTGCGTGGTATCGCCAAAAGCAATGCCCAGATAACGATTAAGCAAAATGGTTACACCATTTACCAAACTTATGTATCGCCTGGTGCTTTTGa                                                                                   >  4:7984/1‑149 (MQ=255)
                                                                       ggCGTCAGATGACAACATGCTGCCAGACAGCCTGAAAGGGTTTGCGCCTGTGGTGCGTGGTATCGCCAAAAGCAATGCCCAGATAACGATTAAGCAAAATGGTTACACCATTTACCAAACTTATGTATCGCCTGGTGCTTTTGAAATTa                                                                              <  2:143013/149‑1 (MQ=255)
                                                                         cGTCAGATGACAACATGCTGCCAGACAGCCTGAAAGGGTTTGCGCCTGTGGTGCGTGGTATCGCCAAAAGCAATGCCCAGATAACGATTAAGCAAAATGGTTACACCATTTACCAAACTTATGTATCGCCTGGTGCTTTTGAAATTAGt                                                                            <  1:490744/149‑1 (MQ=255)
                                                                          gTCAGATGACAACATGCTGCCAGACAGCCTGAAAGGGTTTGCGCCTGTGGTGCGTGGTATCGCCAAAAGCAATGCCCAGATAACGATTAAGCAAAATGGTTACACCATTTACCAAACTTATGTATCGCCTGGTGCTTTTGAAATTAg                                                                             <  3:61230/147‑1 (MQ=255)
                                                                          gTCAGATGACAACATGCTGCCAGACAGCCTGAAAGGGTTTGCGCCTGTGGTGCGTGGTATCGCCAAAAGCAATGCCCAGATAACGATTAAGCAAAATGGTTACACCATTTACCAAACTTATGTATCGCCTGGTGCTTTTGAAATTAg                                                                             >  4:61230/1‑147 (MQ=255)
                                                                             aGATGACAACATGCTGCCAGACAGCCTGAAAGGGTTTGCGCCTGTGGTGCGTGGTATCGCCAAAAGCAATGCCCAGATAACGATTAAGCAAAATGGTTACACCATTTACCAAACTTATGTATCGCCTGGTGCTTTTGAAATTAGTGAtc                                                                        >  4:297421/1‑149 (MQ=255)
                                                                              gATGACAACATGCTGCCAGACAGCCTGAAAGGGTTTGCGCCTGTGGTGCGTGGTATCGCCAAAAGCAATTCCCAGATAACGATTAAGCAAAATGGTTACACCATTTACCAAACTTATGt                                                                                                     >  3:56082/1‑119 (MQ=255)
                                                                              gATGACAACATGCTGCCAGACAGCCTGAAAGGGTTTGCGCCTGTGGTGCGTGGTATCGCCAAAAGCAATTCCCAGATAACGATTAAGCAAAATGGTTACACCATTTACCAAACTTATGt                                                                                                     >  1:426207/1‑119 (MQ=255)
                                                                              gATGACAACATGCTGCCAGACAGCCTGAAAGGGTTTGCGCCTGTGGTGCGTGGTATCGCCAAAAGCAATTCCCAGATAACGATTAAGCAAAATGGTTACACCATTTACCAAACTTATGt                                                                                                     <  4:56082/119‑1 (MQ=255)
                                                                              gATGACAACATGCTGCCAGACAGCCTGAAAGGGTTTGCGCCTGTGGTGCGTGGTATCGCCAAAAGCAATTCCCAGATAACGATTAAGCAAAATGGTTACACCATTTACCAAACTTATGt                                                                                                     <  2:426207/119‑1 (MQ=255)
                                                                              gATGACAACATGCTGCCAGACAGCCTGAAAGGGTTTGCGCCTGTGGTGCGTGGTATCGCCAAAAGCAATGCCCAGATAACGATTAAGCAAAATGGTTACACCATTTACCAAACTTATGTATCGCCTGGTGCTTTTGAAATTAGTGAtct                                                                       >  2:86851/1‑149 (MQ=255)
                                                                                 gACAACATGCTGCCAGACAGCCTGAAAGGGTTTGCGCCTGTGGTGCGTGGTATCGCCAAAAGCAATGCCCAGATAACGATTAAGCAAAATGGTTACACCATTTACCAAACTTATGTATCGCCTGGTGCTTTTGAAATTAGTGATCTCTa                                                                    >  1:556972/1‑149 (MQ=255)
                                                                                 gACAACATGCTGCCAGACAGCCTGAAAGGGTTTGCGCCTGTGGTGCGTGGTATCGCCAAAAGCAATGCCCAGATAACGATTAAGCAAAATGGTTACACCATTTACCAAACTTATGTATCGCCTGGTGCTTTTGAAATTAGTGATCTCTa                                                                    >  3:389395/1‑149 (MQ=255)
                                                                                    aacaTGCTGCCAGACAGCCTGAAAGGGTTTGCGCCTGTGGTGCGTGGTATCGCCAAAAGCAATGCCCAGATAACGATTAAGCAAAATGGTTACACCATTTACCAAACTTATGTATCGCCTGGTGCTTTTGAAATTAGTGATCTCTATTc                                                                 >  1:389745/1‑149 (MQ=255)
                                                                                    aacaTGCTGCCAGACAGCCTGAAAGGGTTTGCGCCTGTGGTGCGTGGTATCGCCAAAAGCAATGCCCAGATAACGATTAAGCAAAATGGTTACACCATTTACCAAACTTATGTATCGCCTGGTGCTTTTGAAATTAGTGATCTCTATTc                                                                 >  2:246080/1‑149 (MQ=255)
                                                                                    aacaTGCTGACAGACAGCCTGAAAGGGTTTGCGCCTGTGGTGCGTGGTATCGCCAAAAGCAATGCCCAGATAACGATTAAGCAAAATGGTTACACCATTTACCAAACTTATGTATCGCCTGGTGCTTTTGAAATTAGTGATCTCTATTc                                                                 >  1:502880/1‑149 (MQ=255)
                                                                                             ccAGACAGCCTGAAAGGGTTTGCGCCTGTGGTGCGTGGTATCGCCAAAAGCAATGCCCAGATAACGATTAAGCAAAATGGTTACACCATTTACCAAACTTATGTATCGCCTGGTGCTTTTGAAATTAGTGATCTCTATTCCACgtcgtc                                                        >  2:445114/1‑149 (MQ=255)
                                                                                             ccAGACAGCCTGAAAGGGTTTGCGCCTGTGGTGCGTGGTATCGCCAAAAGCAATGCCCAGATAACGATTAAGCAAAATGGTTACACCATTTACCAAACTTATGTATCGCCTGGTGCTTTTGAAATTAGTGATCTCTATTCCACgtcgtc                                                        <  3:297421/149‑1 (MQ=255)
                                                                                                    gCCTGAAAGGGTTTGCGCCTGTGGTGCGTGGTATCGCCAAAAGCAATGCCCAGATAACGATTAAGCAAAATGGTTACACCATTTACCAAACTTATGTATCGCCTGGTGCTTTTGAAATTAGTGATCTCTATTCCACGTCGTCGAGCGGt                                                 <  4:403216/149‑1 (MQ=255)
                                                                                                       tGAAAGGGTTTGCGCCTGTGGTGCGTGGTATCGCCAAAAGCAATGCCCAGATAACGATTAAGCAAAATGGTTACACCATTTACCAAACTTATGTATCGCCTGGTGCTTTTGAAATTAGTGATCTCTATTCCACGTCGTCGAGCGGTGAt                                              >  1:149152/1‑149 (MQ=255)
                                                                                                        gAAAGGGTTTGCGCCTGTGGTGCGTGGTATCGCCAAAAGCAATGCCCAGATAACGATTAAGCAAAATGGTTACACCATTTACCAAACTTATGTATCGCCTGGTGCTTTTGAAATTAGTGATCTCTATTcc                                                                <  1:202142/130‑1 (MQ=255)
                                                                                                        gAAAGGGTTTGCGCCTGTGGTGCGTGGTATCGCCAAAAGCAATGCCCAGATAACGATTAAGCAAAATGGTTACACCATTTACCAAACTTATGTATCGCCTGGTGCTTTTGAAATTAGTGATCTCTATTcc                                                                >  2:202142/1‑130 (MQ=255)
                                                                                                           aGGGTTTGCGCCTGTGGTGCGTGGTATCGCCAAAAGCAATGCCCAGATAACGATTAAGCAAAATGGTTACACCATTTACCAAACTTATGTATCGCCTGGTGCTTTTGAAATTAGTGATCTCTATTCCACGTCGTCGAGCGGTGATTTGt                                          <  1:450808/149‑1 (MQ=255)
                                                                                                           aGGGTTTGCGCCTGTGGTGCGTGGTATCGCCAAAAGCAATGCCCAGATAACGATTAAGCAAAATGGTTACACCATTTACCAAACTTATGTATCGCCTGGTGCTTTTGAAATTAGTGATCTCTATTCCACGTCGTCGAGCGGTGATTTGt                                          <  2:518236/149‑1 (MQ=255)
                                                                                                            gggTTTGCGCCTGTGGTGCGTGGTATCGCCAAAAGCAATGCCCAGATAACGATTAAGCAAAATGGTTACACCATTTGCCAAACTTATGTATCGCCTGGTGCTTTTGAAATTAGTGATCTCTATTCCACGTCGTCGAGCGGTGATTTGtt                                         <  4:462989/149‑1 (MQ=255)
                                                                                                             ggTTTGCGCCTGTGGTGCGTGGTATCGCCAAAAGCAATGCCCAGATAACGATTAAGCAAAATGGTTACACCATTTACCAAACTTATGTATCGCCTGGTGCTTTTGAAATTAGTGATCTCTATTCCACGTCGTCGAGCGGTGATTTGTTa                                        <  2:502880/149‑1 (MQ=255)
                                                                                                               tttGCGCCTGTGGTGCGTGGTATCGCCAAAAGCAATGCCCAGATAACGATTAAGCAAAATGGTTACACCATTTACCAAACTTATGTATCGCCTGGTGCTTTTGAAATTAGTGATCTCTATTCCACGTCGTCGAGCGGTGATTTGTTAGt                                      >  1:351138/1‑149 (MQ=255)
                                                                                                               tttGCGCCTGTGGTGCGTGGTATCGCCAAAAGCAATGCCCAGATAACGATTAAGCAAAATGGTTACACCATTTACCAAACTTATGTATCGCCTGGTGCTTTTGAAATTAGTGATCTCTATTCCACGTCGTCGAGCGGTGATTTGTTAGt                                      >  1:379316/1‑149 (MQ=255)
                                                                                                                       tgtgGTGCGTGGTATCGCCAAAAGCAATGCCCAGATAACGATTAAGCAAAATGGTTACACCATTTACCAAACTTATGTATCGCCTGGTGCTTTTGAAATTAGTGATCTCTATTCCACGTCGTCGAGCGGTGATTTGTTAGTTGAAATCa                              >  1:579032/1‑149 (MQ=255)
                                                                                                                       tgtgGTGCGTGGTATCGCCAAAAGCAATGCCCAGATAACGATTAAGCAAAATGGTTACACCATTTACCAAACTTATGTATCGCCTGGTGCTTTTGAAATTAGTGATCTCTATTCCACGTCGTCGAGCGGTGATTTGTTAGTTGAAATCa                              >  1:322228/1‑149 (MQ=255)
                                                                                                                         tggtgCGTGGTATCGCCAAAAGCAATGCCCAGATAATGATTAAGCAAAATGGTTACACCATTTACCAAACTTATGTATCGCCTGGTGCTTTTGAAATTAGTGATCTCTATTCCACGTCGTCGAGCGGTGATTTGTTAGTTGAAATCaaa                            <  3:208618/149‑1 (MQ=255)
                                                                                                                         tggtgCGTGGTATCGCCAAAAGCAATGCCCAGATAACGATTAAGCAAAATGGTTACACCATTTACCAAACTTATGTATCGCCTGGTGCTTTTGAAATTAGTGATCTCTATTCCACGTCGTCGAGCGGTGATTTGTTAGTTGAAATCaaa                            <  1:275924/149‑1 (MQ=255)
                                                                                                                              cGTGGTATCGCCAAAAGCAATGCCCAGATAACGATTAAGCAAAATGGTTACACCATTTACCAAACTTATGTATCGCCTGGTGCTTTTGAAATTAGTGATCTCTATTCCACGTCGTCGAGCGGTGATTTGTTAGTTGAAATCAAAGAAGc                       <  4:491009/149‑1 (MQ=255)
                                                                                                                               gTGGTATCGCCAAAAGCAATGCCCAGATAACGATTAAGCAAAATGGTTACACCATTTACCAAACTTATGTATCGCCTGGTGCTTTTGAAATTAGTGATCTCTATTCCACGTCGTCGAGCGGTGATTTGTTAGTTGAAATCAAAGAAGCg                      <  1:475046/149‑1 (MQ=255)
                                                                                                                                  gtttCGCCAAAAGCAATGCCCAGATAACGATTAAGCAAAATGGTTACACCATTTACCAAACTTATGTATCGCCTGGTGCTTTTGAAATTAGTGATCTCTATTCCACGTCGTCGAGCGGTGATTTGTTAGTTGAAATCAAAGAGGCGAAc                   >  3:377308/4‑149 (MQ=255)
                                                                                                                                   tATCGCCAAAAGCAATGCCCAGATAACGATTAAGCAAAATGGTTACACCATTTACCAAACTTATGTATCGCCTGGTGCTTTTGAAATTAGTGATCTCTATTCCACGTCGTCGAGCGGTGATTTGTTAGTTGAAATCAAAGAAGCGGACg                  >  3:85821/1‑149 (MQ=255)
                                                                                                                                   tATCGCCAAAAGCAATGCCCAGATAACGATTAAGCAAAATGGTTACACCATTTACCAAACTTATGTATCGCCTGGTGCTTTTGAAATTAGTGATCTCTATTCCACGTCGTCGAGCGGTGATTTGTTAGTTGAAATCAAAGAAGCGGACg                  >  1:463749/1‑149 (MQ=255)
                                                                                                                                    aTCGCCAAAAGCAATGCCCAAATAACGATTAAGCAAAATGGTTACACCATTTACCAAACTTATGTATCGCCTGGTGCTTTTGAAATTAGTGATCTCTATTCCACGTCGTCGAGCGGTGATTTGTTAGTTGAAATCAAAGAAGCGGACgg                 >  1:225200/1‑149 (MQ=255)
                                                                                                                                     tCGCCAAAAGCAATGCCCAGATAACGATTAAGCAAAATGGTTACACCATTTACCAAACTTATGTATCGCCTGGTGCTTTTGAAATTAGTGATCTCTATTCCACGTCGTCGAGCGGTGATTTGTTAGTTGAAATCAAAGAAGCGGACGGc                >  3:17673/1‑149 (MQ=255)
                                                                                                                                          aaaaGCAATGCCCAGATAACGATTAAGCAAAATGGTTACACCATTTACCAAACTTATGTATCGCCTGGTGCTTTTGAAATTAGTGATCTCTATTCCACGTCGTCGAGCGGTGATTTGTTAGTTGAAATCAAAGAAGCGGACGGCAGCGt           <  1:445114/149‑1 (MQ=255)
                                                                                                                                               cAATGCCCAGATAACGATTAAGCAAAATGGTTACACCATTTACCAAACTTATGTATCGCCTGGTGCTTTTGAAATTAGTGATCTCTATTCCACGTCGTCGAGCGGTGATTTGTTAGTTGAAATCAAAGAAGCGGACGGCAGCGTCAATa      <  1:312288/149‑1 (MQ=255)
                                                                                                                                               cAATGCCCAGATAACGATTAAGCAAAATGGTTACACCATTTACCAAACTTATGTATCGCCTGGTGCTTTTGAAATTAGTGATCTCTATTCCACGTCGTCGAGCGGTGATTTGTTAGTTGAAATCAAAGAAGCGGACGGCAGCGTCAATa      <  3:7984/149‑1 (MQ=255)
                                                                                                                                                   gCCCAGATAACGATTAAGCAAAATGGTTACACCATTTACCAAACTTATGTATCGCCTGGTGCTTTTGAAATTAGTGATCTCTATTCCACGTCGTCGAGCGGTGATTTGTTAGTTGAAATCAAAGAAGCGGACGGCAGCGTCAATAGCTa  >  2:436844/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GAACTGACCGTAGGTGATGATTATACTGCAGGCGATTTTTTCGATAGTGTCAGCTTTCGTGGTGTGCAGCTGGCGTCAGATGACAACATGCTGCCAGACAGCCTGAAAGGGTTTGCGCCTGTGGTGCGTGGTATCGCCAAAAGCAATGCCCAGATAACGATTAAGCAAAATGGTTACACCATTTACCAAACTTATGTATCGCCTGGTGCTTTTGAAATTAGTGATCTCTATTCCACGTCGTCGAGCGGTGATTTGTTAGTTGAAATCAAAGAAGCGGACGGCAGCGTCAATAGCTA  >  NC_000913/999954‑1000249

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