Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,306,090 G→T 11.7% Q170K (CAA→AAA)  yjcS ← metallo‑beta‑lactamase superfamily protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,306,0900GT11.7% 236.7 / 9.4 85Q170K (CAA→AAA) yjcSmetallo‑beta‑lactamase superfamily protein
Reads supporting (aligned to +/- strand):  ref base G (33/41);  new base T (9/2);  total (42/43)
Fisher's exact test for biased strand distribution p-value = 2.61e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACATTTTATCGG  >  NC_000913/4305942‑4306233
                                                                                                                                                    |                                                                                                                                               
ttttCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTg                                                                                                                                                 <  4:363095/149‑1 (MQ=255)
ttttCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGc                                                                                                                                                >  6:251212/1‑149 (MQ=255)
    cGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGa                                                                                                                                             <  5:251212/148‑1 (MQ=255)
       tGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGa                                                                                                                                          >  7:197580/1‑149 (MQ=255)
            gCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTg                                                                                                                                     <  1:132634/149‑1 (MQ=255)
            gCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTg                                                                                                                                     <  7:348969/149‑1 (MQ=255)
                cGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGa                                                                                                                                   <  1:31768/147‑1 (MQ=255)
                  tCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTa                                                                                                                               <  7:48052/149‑1 (MQ=255)
                   ccATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTa                                                                                                                               >  5:39411/1‑148 (MQ=255)
                   ccATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTa                                                                                                                               >  3:325576/1‑148 (MQ=255)
                    cATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTaaa                                                                                                                             <  1:313254/149‑1 (MQ=255)
                    cATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTaa                                                                                                                              <  1:264352/148‑1 (MQ=255)
                    cATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTa                                                                                                                               >  5:9445/1‑147 (MQ=255)
                      tAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAgg                                                                                                                           <  7:169894/149‑1 (MQ=255)
                        aaGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGAGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGt                                                                                                                          >  1:299215/1‑148 (MQ=255)
                               cAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATACTGCCTTTCACGCCACCATAGGGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCg                                                                                                                  <  7:89946/149‑1 (MQ=255)
                               cAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGc                                                                                                                   <  4:325576/148‑1 (MQ=255)
                                aGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCgg                                                                                                                 <  4:302721/149‑1 (MQ=255)
                                  gggCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCt                                                                                                               >  7:278896/1‑149 (MQ=255)
                                  gggCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTATTTGCGGACTATGCTGGAAGTAAAGGTCAAGTGCGGCt                                                                                                               >  2:329041/1‑149 (MQ=255)
                                    gCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCttt                                                                                                             >  7:58442/1‑149 (MQ=255)
                                     cAATCCCCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAAAGGCAACAATCGGTTTTTGCGGACGATTCTGGAATTAAAGGTCAAGTTCGGCtttt                                                                                                            >  1:313232/1‑149 (MQ=255)
                                     cAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCtttt                                                                                                            >  3:122821/1‑149 (MQ=255)
                                      aaTCACCTGAACTTTGTCGGATTTAATATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTg                                                                                                           <  7:144029/149‑1 (MQ=255)
                                        tCACCTGTACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCg                                                                                                         <  7:147713/149‑1 (MQ=255)
                                          aCCGGAACTTTGCCGGATTTCACATCGGCTTCAGAGATAATGCCTTTCCCGCCACCATAGTGGTCGGTGGGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCgg                                                                                                        <  1:89046/148‑1 (MQ=255)
                                           ccTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCttt                                                                                                             >  1:149459/1‑142 (MQ=255)
                                           ccTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCttt                                                                                                             <  2:149459/142‑1 (MQ=255)
                                           ccTGAACTTAGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCATACATAGTGGTCGTTGTTGCTGTGAGTGTAGATAACTGCAACAATCGGTTTTTTCGGACGTTGCTTTAAGTAATGGTCTAGTGCGGGTTTTGCGGCa                                                                                                      <  8:362654/149‑1 (MQ=255)
                                            cTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAg                                                                                                     >  4:39644/1‑149 (MQ=255)
                                            cTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAg                                                                                                     >  7:54132/1‑149 (MQ=255)
                                            cTGAACTTTGCCGGATTTAACATAGGCTTCAGATATAATGACTTTCACGCCACCATAGATGTCGGTGTGGCTGTGACTGTAGATAACGGCAACAAACCCATATATCCGACCAACCTGTAAGTAAAGGTCAAATGCGGCTTTTGCTGCat                                                                                                     >  2:182051/1‑148 (MQ=255)
                                             tGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGAAAAAAACGTATAATTAAAAAAAAAATGTAAGTAAAGGTCAAGTGAGGCTTTTTCGGCAgg                                                                                                    >  2:10500/1‑149 (MQ=255)
                                               aaCTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATCTCGGAAAACAACGTATAATTAAAAAAAAAATTGAAGTAAAGGTCAAGTTCGGCTTTTGCGGCAggcg                                                                                                  >  6:317741/1‑149 (MQ=255)
                                                 cTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAggcggc                                                                                                <  6:9445/149‑1 (MQ=255)
                                                 cTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAggcggc                                                                                                <  1:10500/149‑1 (MQ=255)
                                                  tttGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCAAAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGAAAAAAACGTTTAATTAAAAATATAATGGAAGTAAAGGTCAAGTGCGGCTTTTTAGGCAGgcggcg                                                                                               >  4:326299/1‑149 (MQ=255)
                                                  tttGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTAACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAAACGGTTATTGAGGAACATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGgcggcg                                                                                               >  4:19315/1‑149 (MQ=255)
                                                      ccGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATTCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCAc                                                                                           <  8:303111/149‑1 (MQ=255)
                                                       cGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGATTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCAc                                                                                           >  4:18426/1‑148 (MQ=255)
                                                        ggATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCa                                                                                         <  8:197580/149‑1 (MQ=255)
                                                             tAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGc                                                                                       >  5:161639/1‑146 (MQ=255)
                                                             tAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGc                                                                                       <  6:161639/146‑1 (MQ=255)
                                                              aaCATCGGCTTCAGAGATAATGCCTTTCACGCCCCCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGAAGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCATGCGGCGTCACCAGCGGGt                                                                                   >  1:259935/1‑149 (MQ=255)
                                                               aCATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTc                                                                                  >  7:341268/1‑149 (MQ=255)
                                                               aCATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTc                                                                                  <  3:18426/149‑1 (MQ=255)
                                                               aCATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTc                                                                                  >  1:6081/1‑149 (MQ=255)
                                                                 aTCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAgg                                                                                                                           <  3:354468/106‑1 (MQ=255)
                                                                 aTCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGGCACCAGCGGGTCGa                                                                                <  3:89782/149‑1 (MQ=255)
                                                                    ggCTTCACAGATAATTCCTTTCACGCCACCAAAGTGGTCGGTGTGGCTGTGACTGTAGATAACGGCAAAAATCGTATATATAGCAAGATACTGCAACTAAAGGTCAAGTGCTGCTTTTGCGGCAGTCGTCGTCACCAACGGGACGaga                                                                              >  6:161022/1‑146 (MQ=255)
                                                                      cTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGGCACCAGCGGGTCGATAACa                                                                           <  3:39644/149‑1 (MQ=255)
                                                                       ttCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACaa                                                                          >  7:393001/1‑149 (MQ=255)
                                                                         cAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACaa                                                                          <  8:54132/147‑1 (MQ=255)
                                                                         cAGAGATAATGCCTTTCACGCCACCATAGTGGGCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACaa                                                                          <  5:339469/147‑1 (MQ=255)
                                                                          agagATAATGCCTTTCACGCCACCAGAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACaataa                                                                       >  5:305011/1‑149 (MQ=255)
                                                                                  tGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGc                                                                                                                                                >  4:40727/1‑68 (MQ=255)
                                                                                  tGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGc                                                                                                                                                <  3:40727/68‑1 (MQ=255)
                                                                                     cTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCg                                                            >  3:150759/1‑149 (MQ=255)
                                                                                      tttCACGCCACAATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGc                                                           >  7:171462/1‑149 (MQ=255)
                                                                                          aCGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTTAGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCCCCAGCGGGTCGAGAACAATAATGCCTTTCTCGCGCt                                                        >  2:331042/1‑148 (MQ=255)
                                                                                          aCGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCt                                                        >  5:151509/1‑148 (MQ=255)
                                                                                            gcgaccaTAGTGGTCGGTGTTTCTGCGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGGCAAGTGCGGCTTTTGCGGCAGGGGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCt                                                        <  2:299215/143‑1 (MQ=255)
                                                                                            gCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAATGTCAAGTGCGGCTTTTGCGGCAGACGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCaa                                                     <  5:176488/149‑1 (MQ=255)
                                                                                             ccaccaTAGGGGTCGGTGTGGCTGGGAGTGTAGATAACGGCAACAATCGGGGTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGGCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCacc                                                    <  1:276206/149‑3 (MQ=255)
                                                                                                 caTAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGGCGATAACAATAATGCCTTTCTCGCCCTCaacgaac                                                <  3:326299/149‑1 (MQ=255)
                                                                                                   tAGTGGGCGGTGGGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGTAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGt                                              <  5:317741/149‑1 (MQ=255)
                                                                                                            gtgtGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTaga                                     <  6:305011/149‑1 (MQ=255)
                                                                                                                 gCTGTGAGTGTAGATAACGGAAAAAAACGTATATTTCAAAAAATTAATTAAGTAAAGGTCAAGTTCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGatat                                >  8:320336/1‑149 (MQ=255)
                                                                                                                   tgtgAGTGTAGATAACGGCAACAAACGTTTTTTTAGGACGATTCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATcc                              >  4:281946/1‑149 (MQ=255)
                                                                                                                   tgtgAGTGTAGATAACGGCAAAAAACGTATAATTAAAAAAATGATGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCACTCAACGAACGTAATGTTAGAGATATcc                              >  8:371999/1‑149 (MQ=255)
                                                                                                                    gtgAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCt                             <  1:331042/149‑1 (MQ=255)
                                                                                                                     tgAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATATCAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCATg                            <  2:313232/149‑1 (MQ=255)
                                                                                                                     tgAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTg                            <  6:151509/149‑1 (MQ=255)
                                                                                                                        gtgtAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGcc                         >  3:9697/1‑149 (MQ=255)
                                                                                                                        gtgtAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGcc                         <  1:182051/149‑1 (MQ=255)
                                                                                                                           tAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCcgcg                      >  2:335873/1‑149 (MQ=255)
                                                                                                                             gATAACGGCAACAATCGGTTTTTTCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGAGTCTCCAGCTGGTCGATAACAATAATGCATTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCgcgc                     <  5:161022/148‑1 (MQ=255)
                                                                                                                              aTAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCAc                   >  8:356377/1‑149 (MQ=255)
                                                                                                                               tAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCAcc                  >  7:397950/1‑149 (MQ=255)
                                                                                                                                   ggCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTATGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAAAAATGCCTTTCACGCCCTCAACTAACGTAATGTTAGAGAAATCCTGGCCGCGCACCTGa               >  6:18634/1‑148 (MQ=255)
                                                                                                                                      aacaaTCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGGCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACa           <  8:278896/149‑1 (MQ=255)
                                                                                                                                          aTCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACAtttt       <  3:19315/149‑1 (MQ=255)
                                                                                                                                           tCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACATTTTa      <  4:150759/149‑1 (MQ=255)
                                                                                                                                              gTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACATTTTATCg   <  8:58442/149‑1 (MQ=255)
                                                                                                                                              gTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAAAGTAATGTTAGAGATATCCTGGACGCGAACCTGATACATTTTATCg   >  3:273387/1‑149 (MQ=255)
                                                                                                                                               tttttGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACATTTTATCgg  >  8:47445/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
TTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACATTTTATCGG  >  NC_000913/4305942‑4306233

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