Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,971,920 G→T 5.9% Q258K (CAA→AAA)  tar ← methyl‑accepting chemotaxis protein II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,971,9200GT 5.9% 363.4 / 10.6 102Q258K (CAA→AAA) tarmethyl‑accepting chemotaxis protein II
Reads supporting (aligned to +/- strand):  ref base G (50/46);  new base T (2/4);  total (52/50)
Fisher's exact test for biased strand distribution p-value = 4.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.09e-01

Only 100 of 105 total aligned reads displayed.

ATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGCCGTACCAC  >  NC_000913/1971772‑1972068
                                                                                                                                                    |                                                                                                                                                    
aTGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTg                                                                                                                                                      >  3:105428/1‑149 (MQ=255)
    tGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCata                                                                                                                                                  >  3:633837/1‑149 (MQ=255)
    tGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCata                                                                                                                                                  >  3:867800/1‑149 (MQ=255)
     ggcggcAGTTTGTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGCGTTGCCCGCCGCAATTTCACGGATACCGGCATAGATGGCATCTGAACCTTTGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCatat                                                                                                                                                 <  4:558008/149‑1 (MQ=255)
       cggcAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATAtg                                                                                                                                                <  1:654122/148‑1 (MQ=255)
            gTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACAGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGaaa                                                                                                                                           >  2:136624/1‑148 (MQ=255)
             tttCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACg                                                                                                                                         >  4:351854/1‑149 (MQ=255)
                 ttcCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTTCACATATGCAACGctct                                                                                                                                     <  4:452372/149‑1 (MQ=255)
                 ttcCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTTCACATATGCAACGctc                                                                                                                                      <  2:639244/148‑1 (MQ=255)
                   ccAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTgc                                                                                                                                   <  4:213960/149‑1 (MQ=255)
                   ccAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTg                                                                                                                                    >  2:80255/1‑148 (MQ=255)
                    cAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTTCACATATGCAACGctct                                                                                                                                     >  2:516386/1‑146 (MQ=255)
                    cAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTTCACATATGCAACGctct                                                                                                                                     <  1:516386/146‑1 (MQ=255)
                     aGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTgcgc                                                                                                                                 >  3:370409/1‑149 (MQ=255)
                      gcgcGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTgcgc                                                                                                                                 <  1:244218/148‑1 (MQ=255)
                       cgcgGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGcc                                                                                                                                <  2:328662/148‑1 (MQ=255)
                        gcgGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAg                                                                                                                              >  4:57666/1‑149 (MQ=255)
                           gATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTc                                                                                                                           >  4:452523/1‑149 (MQ=255)
                              gCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGc                                                                                                                         <  1:660176/148‑1 (MQ=255)
                               ccTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGccc                                                                                                                       <  4:867800/149‑1 (MQ=255)
                               ccTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGccc                                                                                                                       <  4:633837/149‑1 (MQ=255)
                                ctgctgTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCa                                                                                                                      >  3:766096/1‑149 (MQ=255)
                                ctgctgTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCa                                                                                                                      >  4:625468/1‑149 (MQ=255)
                                 tgctgTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCa                                                                                                                      >  1:810505/1‑148 (MQ=255)
                                  gctgTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCAt                                                                                                                     >  2:352067/1‑148 (MQ=255)
                                   ctgTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCAtt                                                                                                                    >  2:265102/1‑148 (MQ=255)
                                           tACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTgcg                                                                                                            >  2:522950/1‑148 (MQ=255)
                                                gaggaAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTc                                                                                                      <  3:93457/149‑1 (MQ=255)
                                                 aggaAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAATAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTc                                                                                                      >  1:699323/1‑148 (MQ=255)
                                                   gaAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCaa                                                                                                    <  1:741691/148‑1 (MQ=255)
                                                    aaaGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAAt                                                                                                   >  2:493091/1‑148 (MQ=255)
                                                       gATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTc                                                                                               >  4:826450/1‑149 (MQ=255)
                                                        aTCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTc                                                                                               >  2:627653/1‑148 (MQ=255)
                                                        aTCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCa                                                                                              >  4:423056/1‑149 (MQ=255)
                                                        aTCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCa                                                                                              >  4:60791/1‑149 (MQ=255)
                                                         tCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTTCACATATGCAACGctct                                                                                                                                     <  1:315710/109‑1 (MQ=255)
                                                         tCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCa                                                                                              >  2:69031/1‑148 (MQ=255)
                                                          cGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAg                                                                                             >  1:558693/1‑148 (MQ=255)
                                                          cGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAg                                                                                             <  1:522950/148‑1 (MQ=255)
                                                          cGGTGTTGCCCGCCGCAATTTCACAGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAg                                                                                             <  1:136624/148‑1 (MQ=255)
                                                              gtTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTa                                                                                         <  2:298432/148‑1 (MQ=255)
                                                                  cccGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGc                                                                                    >  4:752414/1‑149 (MQ=255)
                                                                   ccgccgCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGc                                                                                    >  1:448315/1‑148 (MQ=255)
                                                                     gccgcAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCa                                                                                  <  1:80255/148‑1 (MQ=255)
                                                                     gccgcAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAg                                                                                 <  3:423056/149‑1 (MQ=255)
                                                                      ccgcAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCACATTTCACTGCGCCCGTCAATCGTCAGGGTATTCGCCAg                                                                                 <  1:265102/148‑1 (MQ=255)
                                                                          aaTTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTgcgc                                                                                                                                 >  2:127155/1‑96 (MQ=255)
                                                                          aaTTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTgcgc                                                                                                                                 <  1:127155/96‑1 (MQ=255)
                                                                           aTTTCACGGGTACCTGCATACATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCGTTTCACTGCGCCCGTCAATGGTc                                                                                               <  1:5256/129‑1 (MQ=255)
                                                                           aTTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTc                                                                                               >  2:5256/1‑129 (MQ=255)
                                                                                 cGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCgg                                                                      <  1:69031/148‑1 (MQ=255)
                                                                                  gggTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCg                                                                    >  3:593645/1‑149 (MQ=255)
                                                                                     tACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGa                                                                   >  1:472769/1‑147 (MQ=255)
                                                                                     tACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGa                                                                   <  2:472769/147‑1 (MQ=255)
                                                                                      aCCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGAtt                                                                 >  2:292011/1‑148 (MQ=255)
                                                                                            aTAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGaa                                                          <  4:726119/149‑1 (MQ=255)
                                                                                            aTAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGa                                                           <  1:352067/148‑1 (MQ=255)
                                                                                                 tGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAAtgtg                                                      <  2:699323/148‑1 (MQ=255)
                                                                                                 tGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGa                                                     <  3:1076/149‑1 (MQ=255)
                                                                                                 tGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGa                                                     <  4:593645/149‑1 (MQ=255)
                                                                                                  ggCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGa                                                     >  2:803068/1‑148 (MQ=255)
                                                                                                   gCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGaa                                                          >  2:872292/1‑142 (MQ=255)
                                                                                                   gCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGaa                                                          <  1:872292/142‑1 (MQ=255)
                                                                                                    cATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCa                                                  <  3:351854/149‑1 (MQ=255)
                                                                                                      tCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCaat                                                >  3:846440/1‑149 (MQ=255)
                                                                                                      tCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCaa                                                 >  1:539809/1‑148 (MQ=255)
                                                                                                       cTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCaat                                                >  1:700699/1‑148 (MQ=255)
                                                                                                        tGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCaataa                                              >  4:202054/1‑149 (MQ=255)
                                                                                                          aaCCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAAtt                                            >  4:376247/1‑149 (MQ=255)
                                                                                                            ccTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAAtttt                                          >  4:574149/1‑149 (MQ=255)
                                                                                                            ccTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAAtttt                                          <  3:519275/149‑1 (MQ=255)
                                                                                                             cTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAAttttt                                         <  4:473987/149‑1 (MQ=255)
                                                                                                             cTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAAtttt                                          >  2:867631/1‑148 (MQ=255)
                                                                                                               tCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTg                                        >  2:174123/1‑148 (MQ=255)
                                                                                                               tCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTg                                        <  2:290206/148‑1 (MQ=255)
                                                                                                                      cATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCgg                                 <  1:803068/148‑1 (MQ=255)
                                                                                                                      cATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGa                                <  4:105428/149‑1 (MQ=255)
                                                                                                                       aTGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGa                                >  2:191622/1‑148 (MQ=255)
                                                                                                                       aTGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAg                               <  4:201108/149‑1 (MQ=255)
                                                                                                                        tGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAg                               >  1:277015/1‑148 (MQ=255)
                                                                                                                          aGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTa                             <  2:281160/148‑1 (MQ=255)
                                                                                                                           gTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTaa                            <  2:810505/148‑1 (MQ=255)
                                                                                                                            tGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAg                           <  1:455080/148‑1 (MQ=255)
                                                                                                                             gACGGTGTCAGTCAAAGAGCGTTGCATCTGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAAt                                                                                                   <  4:739129/75‑1 (MQ=255)
                                                                                                                             gACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAAt                                                                                                   >  3:739129/1‑75 (MQ=255)
                                                                                                                                     cagtcaAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATAc                   >  2:167877/1‑147 (MQ=255)
                                                                                                                                     cagtcaAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATAc                   <  1:167877/147‑1 (MQ=255)
                                                                                                                                     cagtcaAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACgg                 <  4:766096/149‑1 (MQ=255)
                                                                                                                                        tcaAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCg               >  1:96054/1‑148 (MQ=255)
                                                                                                                                        tcaAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGa              >  3:4324/1‑149 (MQ=255)
                                                                                                                                         caAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGa              <  1:895308/148‑1 (MQ=255)
                                                                                                                                         caAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGa              >  1:53146/1‑148 (MQ=255)
                                                                                                                                          aaaGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAAt            >  3:300646/1‑149 (MQ=255)
                                                                                                                                          aaaGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAAt            >  4:542876/1‑149 (MQ=255)
                                                                                                                                           aaGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAAt            >  1:382922/1‑148 (MQ=255)
                                                                                                                                             gagCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGc          <  2:722609/148‑1 (MQ=255)
                                                                                                                                             gagCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGc          <  2:558693/148‑1 (MQ=255)
                                                                                                                                                cGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGCCGTa      <  3:60791/149‑1 (MQ=255)
                                                                                                                                                    gCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGCCGTACCAc  <  3:376247/149‑1 (MQ=255)
                                                                                                                                                    gCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGCCGTACCAc  <  3:826450/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGCCGTACCAC  >  NC_000913/1971772‑1972068

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