Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,586,883 T→C 5.3% T284T (ACA→ACG ydeQ ← putative fimbrial‑like adhesin protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,586,8830TC 5.3% 234.7 / 2.9 75T284T (ACA→ACGydeQputative fimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand):  ref base T (31/40);  new base C (2/2);  total (33/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.48e-01

CAGGTGTAGGCTAATAGCCTTCACAAAATAATTATCCGCGCAATTTTCGTGAATGGAGCGCATAGACTAAGGTCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGTTGCGCTTAAACCCAGAGGTACTTTGCTTTTATTTACTGTCCCCAGGGAGACATTTTCTCCCGTAGCCAGTATTTTTCCATTACGCATTAGAGTGACACCCACTCCCGATGCTTTTGTTGCATCAG  >  NC_000913/1586736‑1587022
                                                                                                                                                   |                                                                                                                                           
cAGGTGTAGGCTAATAGCCTTCACAAAATAATTATCCGCGCAATTTTCGTGAATGGAGCGCATAGACTAAGGTCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCt                                                                                                                                             <  7:213759/148‑1 (MQ=255)
 aGGTGTAGGCTAATAGCCTTCACAAAATAATTATCCGCGCAATTTTCGTGAATGGAGCGCATAGACTAAGGTCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGt                                                                                                                                           >  6:284743/1‑149 (MQ=255)
 aGGTGTAGGCTAATAGCCTTCACAAAATAATTATCCGCGCAATTTTCGTGAATGGAGCGCATAGACTAAGGTCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGt                                                                                                                                           <  5:81386/149‑1 (MQ=255)
   gtgtAGGCTAATAGCCTTCACAAAATAATTATCCGCGCAATTTTCGTGAATGGAGCGCATAGACTAAGGTCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGttt                                                                                                                                         <  2:34860/149‑1 (MQ=255)
      tAGGCTAATAGCCTTCACAAAATAATTATCCGCGCAATTTTCGTGAATGGAGCGCATAGACTAAGGTCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGt                                                                                                                                       <  4:230227/148‑1 (MQ=255)
      tAGGCTAATAGCCTTCACAAAATAATTATCCGCGCAATTTTCGTGAATGGAGCGCATAGACTAAGGTCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTc                                                                                                                                      >  7:210011/1‑149 (MQ=255)
        ggCTAATAGCCTTCACAAAATAATTATCCGCGCAATTTTCGTGAATGCAGCGCATAGACTAAGGGCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGACTGACTGGACGGGTCCGGCTGAAACCTTGTTTCCTGTTTGTCCa                                                                                                                                    <  1:235820/149‑1 (MQ=255)
         gCTAATAGCCTTCACAAAATAATTATCCGCGCAATTTTCGTGAATGGAGCGCATAGACTAAGGTCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCa                                                                                                                                    >  8:338165/1‑148 (MQ=255)
           tAATAGCCTTCACAAAATAATTATCCGCGCAATTTTCGTGAATGGAGCGCATAGACTAAGGTCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAg                                                                                                                                 >  1:66074/1‑149 (MQ=255)
           tAATAGCCTTCACAAAATAATTATCCGCGCAATTTTCGTGAATGGAGCGCATAGACTAAGGTCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGACTGAATGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAg                                                                                                                                 >  2:233092/1‑149 (MQ=255)
             aTAGCCTTCACAAAATAATTATCCGCTCAATTTTCGTGAATGGAGCGCATAGACTAAGGTCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAg                                                                                                                                 <  7:19564/147‑1 (MQ=255)
             aTAGCCTTCACAAAATAATTATCCGCGCAATTTTCGTGAATGGAGCGCATAGACTAAGGTCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCCGTTTGTCCATAGGt                                                                                                                               >  8:311854/1‑149 (MQ=255)
                 ccTTCACAAAATAATTATCCGCGCAATTTTCGTGAATGGAGCGCATAGACTAAGGTCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGTTgcg                                                                                                                           >  1:115609/1‑149 (MQ=255)
                  cTTCACAAAATAATTATCCGCGCAATTTTCGTGAATGGAGCGCATAGACTAAGGTCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGTTgcgc                                                                                                                          <  6:148051/149‑1 (MQ=255)
                  cTTCACAAAATAATTATCCGCGCAATTTTCGTGAATGGAGCGCATAGACTAAGGTCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGTTgcgc                                                                                                                          >  3:257690/1‑149 (MQ=255)
                   atCACAAAATAATTATCCGCGCAATTATCCTGAATGGAGCGCATAGACAAAGGTCAAAATAACAATTATTCATAGATAAAAGTGACACCAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCAAAGGTTGctct                                                                                                                         >  6:331795/2‑149 (MQ=255)
                     cacaAAATAATTATCCGCGCAATTTTCGTGAATGGAGCGCATAGACTAAGGTCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGTTGCGCtt                                                                                                                        <  2:205588/148‑1 (MQ=255)
                     cacaAAATAATTATCCGCGCAATTTTCGTGAATGGAGCGCATAGACTAAGGTCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGTTGCGCtt                                                                                                                        <  7:251719/148‑1 (MQ=255)
                      acaAAATAATTATCCGCGCAATTTTCGTGAATGGAGCGCATAGACTAAGGTCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGTTGCGCTTaa                                                                                                                      <  2:222207/149‑1 (MQ=255)
                       ctaaaTAATTATCCGCGCAATTTTCGTGAATGGAGCGCTTAGACTATGGTCAAAATCACAATTATTCATAGATTAAAGTGACCCGAATGACTGACTGGACGGTTGCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGTTGCGCTTaaa                                                                                                                     <  5:331795/147‑1 (MQ=255)
                       caAAATAATTATCCTCGCAATTTTCGTGAATGGAGCGCATAGACTAAGGTCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGTTGCGCTTaaa                                                                                                                     <  5:13422/149‑1 (MQ=255)
                        aaaaTAATTATCCGCGCAATTTTCGTGAATGGAGCGCATAGACTAAGGTCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGTTGCGCTTAAAc                                                                                                                    >  3:318173/1‑149 (MQ=255)
                        aaaaTAATTATCCGCGCAATTTTCGTGAATGGAGCGCATAGACTAAGGTCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTAGTCCATAGGTTGCGCTTAAAc                                                                                                                    >  5:268649/1‑149 (MQ=255)
                        aaaaTAATTATCCGCGCAATTTTCGTGAATGGAGCGCATAGACTAAGGTCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTCTGTCCATAGGTTGCGCTTAAAc                                                                                                                    >  5:6182/1‑149 (MQ=255)
                         aaaTAATTATCCGCGCAATTTTCGTGAATGTAGCGCATAGACTAAGGTCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGTTgcgc                                                                                                                          <  1:124363/142‑1 (MQ=255)
                         aaaTAATTATCCGCGCAATTTTCGTGAATGGAGCGCATAGACTAAGGTCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGTTgcgc                                                                                                                          >  2:124363/1‑142 (MQ=255)
                           ataatTATCCGCGCAATTTTCGTGAATGGAGCGCATAGACTAAGGTCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGTTGCGCTTAAAcc                                                                                                                   >  5:319413/1‑147 (MQ=255)
                           ataatTATCCGCGCAATTTTCGTGAATGGAGCGCATAGACTAAGGTCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGTTGCGCTTAAACCCa                                                                                                                 <  2:10922/149‑1 (MQ=255)
                            taatTATCCGCGCAATTTTCGTGAATGGAGCGCATAGACTAAGGTCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGTTGCGCTTAAACCCag                                                                                                                <  5:284743/149‑1 (MQ=255)
                              atTATCCGCGCAATTTTCGTGAATGGAGCGCATAGACTAAGGTCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGTTGCGCTTAAACCCaga                                                                                                               <  8:333193/148‑1 (MQ=255)
                               ttATCCGCGCAATTTTCGTGAATGGAGCGCATAGACTAAGGTCAAAATCACAATTATTCATAGACAAAAGTGACACCAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGTTGCGCTTAAACCCAGAgg                                                                                                             <  1:302262/149‑1 (MQ=255)
                                         aaTTTTCGTGAATGGAGCGCATAGACTAAGGTCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGTTGCGCTTAAACCCAGAGGTACTTTGCtt                                                                                                   <  7:40294/149‑1 (MQ=255)
                                         aaTTTTCGTGAATGGAGCGCATAGACTAAGGTCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGTTGCGCTTAAACCCAGAGGTACTTTGCtt                                                                                                   >  8:270511/1‑149 (MQ=255)
                                                gTGAATGGAGCGCATAGACTAAGGTCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGACTGACTGGACGGTTCCGGCTTTAACCTTGTTTCCTGTTTGTCCATAGGTTGCGCTTAAACCCAGAGGTACTTTGCTtttattta                                                                                            <  8:329545/149‑1 (MQ=255)
                                                 tGACTGGAGCGCATAGACTAAGGTCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGTTGCGCTTAAACCCAGAGGTACTTTGCTTTTATTTAc                                                                                           <  2:125687/149‑1 (MQ=255)
                                                  gAATGGAGCGCATAGACTAAGGTCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGTTGCGCTTAAACCCAGAGGTACTTTGCTTTTATTTACt                                                                                          >  4:187608/1‑149 (MQ=255)
                                                     tGGAGCGCATAGACTAAGGTCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGTTGCGCTTAAACCCAGAGGTACttt                                                                                                       <  2:21509/133‑1 (MQ=255)
                                                     tGGAGCGCATAGACTAAGGTCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGTTGCGCTTAAACCCAGAGGTACttt                                                                                                       >  1:21509/1‑133 (MQ=255)
                                                          cgcATAGACTAAGGTCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTCTGTCCATAGGTTGCGCTTAAACCCAGAGGTACTTTGCTTTTATTTACTGTCCCCAg                                                                                  <  6:6182/149‑1 (MQ=255)
                                                                gACTAAGGTCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCCGTTTGTCCATAGGTTGCGCTTAAACCCAGAGGTACTTTGCTTTTATTTACTGTCCCCAGGgaga                                                                             <  7:311854/148‑1 (MQ=255)
                                                                     aGGTCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGTTGCGCTTAAACCCAGAGGTACTTTGCTTTTATTTACTGTCCCCAGGGAGACAttt                                                                        <  6:319413/148‑1 (MQ=255)
                                                                     aGGTCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGTTGCGCTTAAACCCAGAGGTACTTTGCTTTTATTTACTGTCCCCAGGGAGACAtt                                                                         <  6:235527/147‑1 (MQ=255)
                                                                     aGGTCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGTTGCGCTTAAACCCAGAGGTACTTTGCTTTTATTTACTGTCCCCAGGGAGACAtt                                                                         <  2:300804/147‑1 (MQ=255)
                                                                        tCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGTTGCGCTTAAACCCAGAGGTACTTTGCTTTTATTTACTGTCCCCAGGGAGGCATTTtctc                                                                    >  7:153711/1‑149 (MQ=255)
                                                                        tCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGTTGCGCTTAAACCCAGAGGTACTTTGCTTTTATTTACTGTCCCCAGGGAGACATTTtctc                                                                    >  7:182122/1‑149 (MQ=255)
                                                                          aaaaTCACAATTATTCATAGATAAAAGTGACACCAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGTTGCGCTTAAACCCAGAGGTACTTTGCTTTTATTTACTGTCCCCAGGGAGACATTTTCTccc                                                                  <  2:161055/149‑1 (MQ=255)
                                                                          aaaaTCACAATTATTCATAGATAAAAGTGACACCAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGTTGCGCTTAAACCCAGAGGTACTTTGCTTTTATTTACTGTCCCCAGGGAGACATTTTCTccc                                                                  >  2:222109/1‑149 (MQ=255)
                                                                            aaTCACAATTATTCATAGATAAAAGTGACACCAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGTTGCGCTTAAACCCAGAGGTACTTTGCTTTTATTTACTGTCCCCAGGGAGAGATTTTGTCCCgg                                                                >  4:79067/1‑148 (MQ=255)
                                                                             aTCACAATTATTCATAGATAAAAGTGACACCAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGTTGCGCTTAAACCCAGAGGTACTTTGCTTTTATTTACTGTCCCCAGGGAGACATTTTCTCCCg                                                                 >  3:208882/1‑147 (MQ=255)
                                                                                acaATTATTCATAGATAAAAGTGACACCAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTcc                                                                                                                                     >  5:337167/1‑76 (MQ=255)
                                                                                acaATTATTCATAGATAAAAGTGACACCAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTcc                                                                                                                                     <  6:337167/76‑1 (MQ=255)
                                                                                          aTAGATAAAAGTGACACCAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGTTGCGTTTAAACCGAGTGGTACTTTGCTTTTATTTACTGTCCCCAGGGAGACATTTTCTCCCGTAGCCAGTAttttt                                                   <  6:268649/148‑1 (MQ=255)
                                                                                          aTAGATAAAAGTGACACCAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGTTGCGCTTAAACCCAGAGGTACTTTGCTTTTATTTACTGTCCCCAGGGAGACATTTTCTCCCGTAGCCAGTATTTTTc                                                  <  7:270511/149‑1 (MQ=255)
                                                                                             gATAAAAGTGACACCAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGTTGCGCTTAAACCCAGAGGTACTTTGCTTTTATTTACTGTCCCCAGGGAGACATTTTCTCCCGTAGCCAGTATTTTTCCa                                                <  8:210011/148‑1 (MQ=255)
                                                                                                       acacCAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGTTGCGCTTAAACCCAGAGGTACTTTGCTTTTATTTACTGTCCCCAGGGAGACATTTTCTCCCGTAGCCAGTATTTTTCCATTACGCATTa                                      <  1:233092/148‑1 (MQ=255)
                                                                                                        cacCAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGTTGCGCTTAAACCCAGAGGTACTTTGCTTTTATTTACTGTCCCCAGGGAGACATTTTCTCCCGTAGCCAGTATTTTTCCATTACGCATTaga                                    <  7:338165/149‑1 (MQ=255)
                                                                                                        cacCAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGTTGCGCTTAAACCCAGAGGTACTTTGCTTTTATTTACTGTCCCCAGGGAGACATTTTCTCCCGTAGCCAGTATTTTTCCATTACGCATTaga                                    <  3:187608/149‑1 (MQ=255)
                                                                                                        cacCAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGTTGCGCTTAAACCCAGAGGTACTTTGCTTTTATTTACTGTCCCCAGGGAGACATTTTCTCCCGTAGCCAGTATTTTTCCATTACGCATTaga                                    <  2:66074/149‑1 (MQ=255)
                                                                                                         acCAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGTTGCGCTTAAACCCAGAGGTACTTTGCTTTTATTTACTGTCCCCAGGGAGACATTTTCTCCCGTAGCCAGTATTTTTCCATTACGCATTagag                                   >  3:307804/1‑149 (MQ=255)
                                                                                                          ccAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGTTGCGCTTAAACCCAGAGGTACTTTGCTTTTATTTACTGTCCCCAGGGAGACATTTTCTCCCGTAGCCAGTATTTTTCCATTACGCATTAGAGt                                  >  1:61240/1‑149 (MQ=255)
                                                                                                           cAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCCGTTTGTCCATAGGTTGCGCTTAAACCCAGAGGTACTTTGCTTTTATTTACTGTCCCCAGGGAGACATTTTCTCCCGTAGCCAGTATTTTTc                                                  >  1:241736/1‑132 (MQ=255)
                                                                                                           cAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCCGTTTGTCCATAGGTTGCGCTTAAACCCAGAGGTACTTTGCTTTTATTTACTGTCCCCAGGGAGACATTTTCTCCCGTAGCCAGTATTTTTc                                                  <  2:241736/132‑1 (MQ=255)
                                                                                                            aaTGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGTTGCGCTTAAACCCAGAGGTACTTTGCTTTTATTTACTGTCCCCAGGGAGACATTTTCTCCCGTAGCCAGTATTTTTCCATTACGCATTAGAGTGa                                >  6:124901/1‑149 (MQ=255)
                                                                                                                actgactgGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGTTGCGCTTAAACCCAGAGGTACTTTGCTTTTATTTACTGTCCCCAGGGAGACATTTTCTCCCGTAGCCAGTATTTTTCCATTACGCATTAGAGTGACAcc                            <  4:307804/149‑1 (MQ=255)
                                                                                                                    actgGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGTTGCGCTTAAACCCAGAGGTACTTTGCTTTTATTTACTGTCCCCAGGGAGACATTTTCTCCCGTAGCCAGTATTTTTCCATTACGCATTAGAGTGACACCCAc                         <  4:208882/148‑1 (MQ=255)
                                                                                                                    actgGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGTTGCGCTTAAACCCAGAGGTACTTTGCTTTTATTTACTGTCCCCAGGGAGACATTTTCTCCCGTAGCCAGTATTTTTCCATTACGCATTAGAGTGACACCCACt                        >  8:57300/1‑149 (MQ=255)
                                                                                                                           ggTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGTTGCGCTTAAACCCAGAGGTACTTTGCTTTTATTTACTGTCCCCAGGGAGACATTTTCTCCCGTAGCCAGTATTTTTCCATTACGCATTAGAGTGACACCCACTCCCGATg                 >  6:112310/1‑149 (MQ=255)
                                                                                                                                  gCTGAAACCTTGTTTCCTGTTTGTCCATAGGATGCGCTTAAACCCAGAGGTACTTTGCTTTTATTTACTGTCCCCAGGGAGACATTTtctc                                                                    >  2:214178/1‑91 (MQ=255)
                                                                                                                                   cTGAAACCTTGTTTCCTGTTTGTCCATAGGTTGCGCTTAAACCCAGAGGTACTTTGCTTTTATTTACTGTCCCCAGGGAGACATTTTCTCCCGTAGCCAGTATTTTTCCATTACGCATTAGAGTGACACCCACTCCCGATGCTTttgtt         <  1:174296/149‑1 (MQ=255)
                                                                                                                                    tGAAACCTTGTTTCCTGTTTGTCCATAGGTTGCGCTTAAACCCAGAGGTACTTTGCTTTTATTTACTGTCCCCAGGGAGACATTTTCTCCCGTAGCCAGTATTTTTCCATTACGCATTAGAGTGACACCCACTCCCGATGCTTttgttg        >  7:111104/1‑149 (MQ=255)
                                                                                                                                        aCCTTGTTTCCTGTTTGTCCATAGGTTGCGCTTAAACCCAGAGGTACTTTGCTTTTATTTACTGTCCCCAGGGAGACATTTTCTCCCGTAGCCAGTATTTTTCCATTACGCATTAGAGTGACACCCACTCCCGATGCTTTTGTTGCATc    <  4:318173/149‑1 (MQ=255)
                                                                                                                                         ccTTGTTTCCTGTTTGTCCATAGGTTGCGCTTAAACCCAGAGGTACTTTGCtttt                                                                                                 >  8:221096/1‑55 (MQ=255)
                                                                                                                                         ccTTGTTTCCTGTTTGTCCATAGGTTGCGCTTAAACCCAGAGGTACTTTGCtttt                                                                                                 <  7:221096/55‑1 (MQ=255)
                                                                                                                                          cTTGTTTCCTGTTTGTCCATAGGTTGCGCTTAAACCCAGAGGTACTTTGCTTTTATTTACTGTCCCCAGGGAGACATTTTCTCCCGTAGCCAGTATTTTTCCATTACGCATTAGAGTGACACCCACTCCCGATGCTTTTGTTGc       <  2:127250/144‑1 (MQ=255)
                                                                                                                                          cTTGTTTCCTGTTTGTCCATAGGTTGCGCTTAAACCCAGAGGTACTTTGCTTTTATTTACTGTCCCCAGGGAGACATTTTCTCCCGTAGCCAGTATTTTTCCATTACGCATTAGAGTGACACCCACTCCCGATGCTTTTGTTGc       >  1:127250/1‑144 (MQ=255)
                                                                                                                                          cTTGTTTCCTGTTTGTCCATAGGGTGCGCTTAAACCCAGAGGTACTTTGCTTTTATTTACTGTCCCCAGGGAGACATTTTCTCCCGTAGCCAGTATTTTTCCATTACGCATTAGAGTGACACCAACTCCCGATGCTTTTGTTGCATCAg  <  8:14014/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                           
CAGGTGTAGGCTAATAGCCTTCACAAAATAATTATCCGCGCAATTTTCGTGAATGGAGCGCATAGACTAAGGTCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGTTGCGCTTAAACCCAGAGGTACTTTGCTTTTATTTACTGTCCCCAGGGAGACATTTTCTCCCGTAGCCAGTATTTTTCCATTACGCATTAGAGTGACACCCACTCCCGATGCTTTTGTTGCATCAG  >  NC_000913/1586736‑1587022

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