Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 549,681 T→G 5.8% S19A (TCC→GCC) ‡ ylbF → putative anaerobic allantoin catabolic oxamate carbamoyltransferase; DUF2877 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913549,6810TG 5.8% 356.1 / 12.7 104S19A (TCC→GCC) ‡ylbFputative anaerobic allantoin catabolic oxamate carbamoyltransferase; DUF2877 family protein
Reads supporting (aligned to +/- strand):  ref base T (49/49);  new base G (3/3);  total (52/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.17e-01

Only 100 of 106 total aligned reads displayed.

TGGCGCAGGCACCGTGCGGGCACCGCTGGCGTGCTTTGAACAGGCGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGTGCTTCGCCGTGCGCAATTCGATGCGTTATGCGGTGGATT  >  NC_000913/549535‑549829
                                                                                                                                                  |                                                                                                                                                    
tGGCGCAGGCACCGTGCGGGCACCGCTGGCGTGCTTTGAACAGGCGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTc                                                                                                                                                     >  1:658126/1‑148 (MQ=255)
    gcAGGCACCGTGCGGGCACCGCTGGCGTGCTTTGAACAGGCGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGc                                                                                                                                                >  3:456141/1‑149 (MQ=255)
     cAGGCACCGTGCGGGCACCGCTGGCGTGCTTTGAACAGGCGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGc                                                                                                                                                >  2:375446/1‑148 (MQ=255)
       ggCACCGTGCGGGCACCGCTGGCGTGCTTTGAACAGGCGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGGACTGGCgat                                                                                                                                             <  4:135602/149‑3 (MQ=255)
         cACCGTGCGGGCACCGCTGGCGTGCTTTGAACAGGCGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGGACTGGCgataa                                                                                                                                           <  4:554731/149‑5 (MQ=255)
           ccGTGCGGGCACCGCTGGCGTGCTTTGAACAGGCGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAg                                                                                                                                          >  2:178494/1‑148 (MQ=255)
             gTGCGGGCACCGCTGGCGTGCTTTGAACAGGCGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCgg                                                                                                                                       <  3:807989/149‑1 (MQ=255)
               gCGGGCACCGCTGGCGTGCTTTGAACAGGCGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCggg                                                                                                                                      >  1:376891/1‑148 (MQ=255)
                 gggCACCGCTGGCGTGCTTTGAACAGGCGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGAgt                                                                                                                                   >  4:300127/1‑149 (MQ=255)
                     aCCGCTGGCGTGCTTTGAACAGGCGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGAGTGTgg                                                                                                                               <  3:300127/149‑1 (MQ=255)
                      ccGCTGGCGTGCTTTGAACAGGCGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGc                                                                                                                                         >  3:304768/1‑138 (MQ=255)
                      ccGCTGGCGTGCTTTGAACAGGCGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGc                                                                                                                                         <  4:304768/138‑1 (MQ=255)
                       cGCTGGCGTGCTTTGAACAGGCGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGAGTGTggc                                                                                                                              <  2:350989/148‑1 (MQ=255)
                              gTGCTTTGAACAGGCGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGTGAgtgt                                                                                                                                 >  3:301296/1‑138 (MQ=255)
                              gTGCTTTGAACAGGCGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGAgtgt                                                                                                                                 <  4:301296/138‑1 (MQ=255)
                               tGCTTTGAACAGGCGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGAGTGTGGCggcgggcg                                                                                                                      <  1:514388/148‑1 (MQ=255)
                               tGCTTTGAACAGGCGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGAGTGTGGCggcgggcg                                                                                                                      <  1:666635/148‑1 (MQ=255)
                                 cTTTGAACAGGCGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGAtt                                                                                                                   <  3:519369/149‑1 (MQ=255)
                                  tttGAACAGGCGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTa                                                                                                                  >  3:329320/1‑149 (MQ=255)
                                  tttGAACAGGCGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTa                                                                                                                  >  4:641906/1‑149 (MQ=255)
                                      aaCAGGCGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAAcc                                                                                                               <  2:613023/148‑1 (MQ=255)
                                      aaCAGGCGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCt                                                                                                              <  4:70147/149‑1 (MQ=255)
                                      aaCAGGCGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCt                                                                                                              <  4:456141/149‑1 (MQ=255)
                                         aGGCGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCtga                                                                                                            <  2:575578/148‑1 (MQ=255)
                                          ggCGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCtgatg                                                                                                          >  4:667997/1‑149 (MQ=255)
                                           gcgcTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCtgatga                                                                                                         <  3:172221/149‑1 (MQ=255)
                                            cgcTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCTGATGAc                                                                                                        <  3:23392/149‑1 (MQ=255)
                                                 gAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCTGATGACGGaaa                                                                                                   >  3:766687/1‑149 (MQ=255)
                                                 gAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCTGATGACGGaa                                                                                                    >  2:908444/1‑148 (MQ=255)
                                                         ggCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCTATCCTCTGCTTGTCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGa                                                                                            <  1:652974/148‑1 (MQ=255)
                                                           cTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAAc                                                                                          <  2:658209/148‑1 (MQ=255)
                                                              aaaGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTATTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGCGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGt                                                                                       >  1:228480/1‑148 (MQ=255)
                                                              aaaGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGtt                                                                                      <  4:524599/149‑1 (MQ=255)
                                                               aaGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGtt                                                                                      >  2:882590/1‑148 (MQ=255)
                                                                  cATGGGTATGGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAAc                                                                                   >  2:615636/1‑148 (MQ=255)
                                                                   aTGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGt                                                                                 >  4:53629/1‑149 (MQ=255)
                                                                    tGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCACCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGtt                                                                                <  4:590091/149‑1 (MQ=255)
                                                                     gggTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGtt                                                                                <  1:76852/148‑1 (MQ=255)
                                                                      ggTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGc                                                                              <  3:4824/149‑1 (MQ=255)
                                                                       gTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGc                                                                              <  2:376891/148‑1 (MQ=255)
                                                                        tATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCa                                                                             >  1:510196/1‑148 (MQ=255)
                                                                        tATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCAt                                                                            >  3:321213/1‑149 (MQ=255)
                                                                          ttGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGc                                                                              <  3:463117/145‑1 (MQ=255)
                                                                          ttGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGc                                                                              >  4:463117/1‑145 (MQ=255)
                                                                             gTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAAcc                                                                                                               >  3:533536/1‑109 (MQ=255)
                                                                             gTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAAcc                                                                                                               <  4:533536/109‑1 (MQ=255)
                                                                              ttGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAg                                                                      >  4:309297/1‑149 (MQ=255)
                                                                              ttGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAg                                                                      <  3:309297/149‑1 (MQ=255)
                                                                                 aGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAggg                                                                    <  1:327936/148‑1 (MQ=255)
                                                                                  ggAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTa                                                                  >  3:894518/1‑149 (MQ=255)
                                                                                   gAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAAc                                                                                                                >  3:443764/1‑102 (MQ=255)
                                                                                   gAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAg                                                                 >  3:647000/1‑149 (MQ=255)
                                                                                   gAACGCGCAAAGACGATCATCCATCCGCTGCTTGCCAGTAGGAGCGCACCGAATTATCGGCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAACGGCGAACTGTTAACGTTGCAGCGTCAGGGTAg                                                                 <  4:435594/149‑1 (MQ=255)
                                                                                    aaCGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGt                                                                <  3:37095/149‑1 (MQ=255)
                                                                                     aCGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTTTCGCCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGAGTAACCTGATGACGGAAAGCGGCGACCTGTTAACGTTGCATCGTCAGGGTAGTg                                                               <  4:554634/149‑1 (MQ=255)
                                                                                     aCGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGt                                                                <  2:653388/148‑1 (MQ=255)
                                                                                          cAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTc                                                          <  4:321213/149‑1 (MQ=255)
                                                                                           aaTGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTc                                                          >  1:752827/1‑148 (MQ=255)
                                                                                             tGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCgg                                                        >  1:378838/1‑148 (MQ=255)
                                                                                               aCGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGcc                                                      >  2:775847/1‑148 (MQ=255)
                                                                                                  atcatcCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAAc                                                                                                                <  1:812344/87‑1 (MQ=255)
                                                                                                  atcatcCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAAc                                                                                                                >  2:812344/1‑87 (MQ=255)
                                                                                                  atcatcCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCtga                                                                                                            >  1:712691/1‑91 (MQ=255)
                                                                                                  atcatcCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCtga                                                                                                            <  2:712691/91‑1 (MQ=255)
                                                                                                  atcatcCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCg                                                   >  1:163954/1‑148 (MQ=255)
                                                                                                   tcatcCCTCCTCTGGTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGGTAGCGGGAGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCgga                                                 <  4:332659/149‑1 (MQ=255)
                                                                                                   tcatcCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGt                                                                <  3:718136/134‑1 (MQ=255)
                                                                                                   tcatcCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGt                                                                >  4:718136/1‑134 (MQ=255)
                                                                                                     atccatccTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCggag                                                <  2:412022/148‑1 (MQ=255)
                                                                                                          tccTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATgg                                           >  2:480289/1‑148 (MQ=255)
                                                                                                           ccTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATggg                                          >  1:297590/1‑148 (MQ=255)
                                                                                                            ctctGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGTg                                        >  3:850851/1‑149 (MQ=255)
                                                                                                               tGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGTGCt                                      <  1:882590/148‑1 (MQ=255)
                                                                                                                       agtagtAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGTGCTTCGCCGTgc                             <  3:878586/149‑1 (MQ=255)
                                                                                                                              gcgcACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGTGCTTCGCCGTGCGCAAtt                       <  1:615636/148‑1 (MQ=255)
                                                                                                                                  aCCGAATTATCGCCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGTGCTTCGCCGTGCGCAATTCGAt                   >  1:35824/1‑148 (MQ=255)
                                                                                                                                  aCCGAATTATCGCCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGTGCTTCGCCGTGCGCAATTCGATg                  >  4:504747/1‑149 (MQ=255)
                                                                                                                                    cGAATTATCGCCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCTCCGGAGGATGGGTGCTTCGCCGTGCGCAATTCGATGCg                >  4:69542/1‑149 (MQ=255)
                                                                                                                                    cGAATTATCGCCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGTGCTTCGCCGTGCGCAATTCGATGCg                >  4:141328/1‑149 (MQ=255)
                                                                                                                                     gAATTATCGCCAGTCCTAGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGTGCTTCGCCGTGCGCAATTCGATGCg                <  1:43598/148‑1 (MQ=255)
                                                                                                                                      aaTTATCGCCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGTGCTTCGCCGTGCGCAATTCGATGCGt               >  2:730862/1‑148 (MQ=255)
                                                                                                                                        taaccGCCAGGACTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGTGCTTCGCCGTGCGCAATTCGATGCGTTAt            >  4:441659/5‑149 (MQ=255)
                                                                                                                                         tATCGCCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGTg                                        >  2:28373/1‑120 (MQ=255)
                                                                                                                                         tATCGCCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCAGCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGTg                                        <  1:28373/120‑1 (MQ=255)
                                                                                                                                         aaccGCCAGGACTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGTGCTTCGCCGTGCGCAATTCGATGCGTTATg           >  3:161151/4‑149 (MQ=255)
                                                                                                                                          accGCCAGGACTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGAt                                             >  1:419843/3‑114 (MQ=255)
                                                                                                                                           tCGCCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCTGATGACCGAAAGCGGCGAAATGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGTGCTTCGCCGTGCGCAATTCGATGCGTTATGc          <  2:510196/148‑1 (MQ=255)
                                                                                                                                             gCCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGTGCTTCGCCGTGCGCAATTCGATGCGTTATGCGGt       >  3:605371/1‑149 (MQ=255)
                                                                                                                                              ccAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATg                                            >  2:144449/1‑111 (MQ=255)
                                                                                                                                              ccAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATg                                            <  1:144449/111‑1 (MQ=255)
                                                                                                                                              ccAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGTGCTTCGCCGTGCGCAATTCGATGCGTTATGCGGt       >  1:509535/1‑148 (MQ=255)
                                                                                                                                              ccAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGTGCTTCGCCGTGCGCAATTCGATGCGTTATGCGGt       <  1:466728/148‑1 (MQ=255)
                                                                                                                                              ccAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGTGCTTCGCCGTGCGCAATTCGATGCGTTATGCGGt       <  1:535396/148‑1 (MQ=255)
                                                                                                                                              ccAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGTGCTTCGCCGTGCGCAATTCGATGCGTTATGCGGt       >  1:578150/1‑148 (MQ=255)
                                                                                                                                              ccAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGTGCTTCGCCGTGCGCAATTCGATGCGTTATGCGGt       >  2:535396/1‑148 (MQ=255)
                                                                                                                                               cAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGTGCTTCGCCGTGCGCAATTCGATGCGTTATGCGGTgg     <  4:329320/149‑1 (MQ=255)
                                                                                                                                               cAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGTGCTTCGCCGTGCGCAATTCGATGCGTTATGCGGTg      >  2:295313/1‑148 (MQ=255)
                                                                                                                                                aGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGTGCTTCGCCGTGCGCAATTCGATGCGTTATGCGGTGGa    <  4:647000/149‑1 (MQ=255)
                                                                                                                                                  tCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGTGCTTCGCCGTGCGCAATTCGATGCGTTATGCGGTGGAtt  >  4:586235/1‑149 (MQ=255)
                                                                                                                                                  tCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGTGCTTCGCCGTGCGCAATTCGATGCGTTATGCGGTGGAtt  >  4:104583/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TGGCGCAGGCACCGTGCGGGCACCGCTGGCGTGCTTTGAACAGGCGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGTGCTTCGCCGTGCGCAATTCGATGCGTTATGCGGTGGATT  >  NC_000913/549535‑549829

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