Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,125,654 T→C 5.3% Q31Q (CAA→CAG)  wcaI ← putative glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,125,6540TC 5.3% 230.4 / 2.6 77Q31Q (CAA→CAG) wcaIputative glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base T (38/34);  new base C (2/2);  total (41/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

GGGTGCTCGGCTGTTTTGGCACATACAGCGGGCAGCGCCACACCGTGGCGGCCCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATGCCGGTTAACTCCGGCGAGTAGTTAATGCCGTAGACCAGTATTTTCATAATCCGGGTACTCCGGTACGCTTCTCAGCGAGAAAATAGGCGCGGCTGTTAGCAT  >  NC_000913/2125506‑2125801
                                                                                                                                                    |                                                                                                                                                   
gggTGCTCGGCTGTTTTGGCACATACAGCGGGCAGCGCCACACCGTGGCGGCCCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTAATGACCCGCACCTCATGACCt                                                                                                                                                     >  2:363686/1‑149 (MQ=255)
 ggTGCTCGGCTGTTTTGGCACATACAGCGGGCAGCGCCACACCGTGGCGGCCCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTAATGACCCGCACCTCATGACCtt                                                                                                                                                    >  2:899300/1‑149 (MQ=255)
      tCGGCTGTTTTGGCACATACAGCGGGCAGCGCCACACCGTGGCGGCCCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTAATGCCCCGCACCTCATGACCTTGTgcc                                                                                                                                               >  2:393394/1‑149 (MQ=255)
        ggCTGTTTTGGCACATACAGCGGGCAGCGCCACACCGTGGCGGCCCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCTCCCACCTGGCATTGCGGGTAGTAAGGCGGTGCGGTAATGACCCGCACCTCATGACCATGTgccgc                                                                                                                                             >  2:405720/1‑149 (MQ=255)
          cggtgttGGCACATACAGCGGGCAGCGCCACACCGGGGCGGCCCCCGCTTCTCGTTTGAAGCGCCAGGCGGAATAGTCCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTAATGACCCGCACCTCAGGACCTTGTgccgcaa                                                                                                                                           <  2:1102662/144‑3 (MQ=255)
           tGTTTTGGCACATACAGCGGGCAGCGCCACACCGTGGCGGCCCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTAATGACCCGCACCTCATGACCTTGTGCCgccag                                                                                                                                          >  1:1029493/1‑149 (MQ=255)
            gTTTTGGCACATACAGCGGGCAGCGCCACACCGTGGCGGCCCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTAATGACCCGCACCTCATGACCTTGTGCCgccagc                                                                                                                                         <  1:929115/149‑1 (MQ=255)
            gTTTTGGCACATACAGCGGGCAGCGCCACACCGTGGCGGCCCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTAATGACCCGCACCTCATGACCTTGTGCCgccagc                                                                                                                                         >  1:609934/1‑149 (MQ=255)
                tGGCAAATACAGCGGGCAGCGCCACACCGTGGCGGCCCCCTCTTCTCGTTTGTAGCGCCAGGCCGAATAGATCTCGCCCCCCGGCCAATGCTGGTAGTAAGGCCGTGCGGTAATGACCCGCACCTCATCACCTTGAGCCGCCTGCCCtt                                                                                                                                     >  1:796246/1‑149 (MQ=255)
                  gCACATACAGCGGGCAGCGCCACACCGTGGCGGCCCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTAATGACCCGCGCCTCATGACCTTGTGGCGCCCGCCATTcc                                                                                                                                   >  2:118592/1‑149 (MQ=255)
                  gCACATACAGCGGGCAGCGCCACACCGTGGCGGCCCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTcc                                                                                                                                   <  2:1165653/149‑1 (MQ=255)
                    acaTACAGCGGGCAGCGCCACACCGTGGCGGCCCCCTCTTCTCGTTTGTAGCGCCCGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTccac                                                                                                                                 >  2:58104/1‑149 (MQ=255)
                     caTACAGCGGGCAGCGCCACACCGTGGCGGCCCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTccacc                                                                                                                                <  2:593294/149‑1 (MQ=255)
                        aCAGCGGGCAGCGCCACACCGTGGCGGCCCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCCCCTGCCATTGCGGGGAGTAAGGCGGTGCGGTAATGACCCGCACCTCATGCCCCGGTGCCGCCAGCCCTTCCACCAtc                                                                                                                             >  1:955608/1‑149 (MQ=255)
                        aCAGCGGGCAGCGCCACACCGTGGCGGCCCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCAtc                                                                                                                             >  1:1136482/1‑149 (MQ=255)
                        aCAGCGGGCAGCGCCACACCGTGGCGGCCCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTAATGACCCGCACCTCATGACCCTGTGCCGCCAGCCAGTCCACCAtc                                                                                                                             >  1:960170/1‑149 (MQ=255)
                          aGCGGGCAGCGCCACACCGTGGCGGCCCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCACGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTAATGAACCGCCCCTCATGACCTTGTGCCGCCACACATTCAACCAtctc                                                                                                                           >  2:109050/1‑149 (MQ=255)
                             gggCAGCGCCACACCGTGGCGGCCCCCTCTGCTCGTTAGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCCTAGCGGGGAGTAAGGTGGTGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCCCCATCTCGcc                                                                                                                        >  2:527653/1‑149 (MQ=255)
                               gCAGCGCCACACCGTGGCGGCCCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCgg                                                                                                                      >  2:581109/1‑149 (MQ=255)
                                  gcgcCACACCGTGGCGGCCCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGgtgt                                                                                                                   <  1:1007549/149‑1 (MQ=255)
                                   cgcCACACCGTGGCGGCCCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTa                                                                                                                  >  2:582599/1‑149 (MQ=255)
                                     ccACACCGTGGCGGCCCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCCCCATCTCGCCGGTGTAtt                                                                                                                >  1:1072502/1‑149 (MQ=255)
                                         acCGTGGCGGCCCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTAAGGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTAGCCGGTGTATTTGcc                                                                                                            >  2:84571/1‑149 (MQ=255)
                                                 ggcCCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGGAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATGCCgg                                                                                                    <  1:70931/149‑1 (MQ=255)
                                                   cccccTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATGCCGGtt                                                                                                  <  1:1088341/149‑1 (MQ=255)
                                                     cccTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATGCCGGTTaa                                                                                                >  2:175494/1‑149 (MQ=255)
                                                     cccTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATGCCGGTTaa                                                                                                >  2:541046/1‑149 (MQ=255)
                                                     cccTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATGCCGGTTaa                                                                                                >  2:626861/1‑149 (MQ=255)
                                                            ctcGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATGCCGGTTAACTCCGGc                                                                                         <  1:902595/149‑1 (MQ=255)
                                                            ctcGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATGCCGGTTAACTCCGGc                                                                                         >  2:750338/1‑149 (MQ=255)
                                                            ctcGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCAGCCATTGCGGGTAGTAAGGCGGCGCGGTAATGACCCGCACCTCATGCCCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATGCCGGTTAACTCCGGc                                                                                         <  1:362362/149‑1 (MQ=255)
                                                              cGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATGCCGGTTAACTCCGGCGa                                                                                       <  1:290078/149‑1 (MQ=255)
                                                                       cgcCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTAATGACCCGCACCGCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATGCCGGTTAACTCCGGCGAGTAGTTAAt                                                                              <  2:857470/149‑1 (MQ=255)
                                                                           aGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATGCCGGTTAACTCCGGCGAGTAGTTAATGCCg                                                                          >  1:413495/1‑149 (MQ=255)
                                                                             gCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATGCCGGTTAACTCCGGCGAGTAGTTAATGCCGTa                                                                        <  2:1072502/149‑1 (MQ=255)
                                                                                gAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCAACCATCTCGCCGGTGTATTTGCCGATGCCGGTTAACTCCGGCGAGTAGTTAATGCCGTAGAc                                                                     >  1:646838/1‑149 (MQ=255)
                                                                                     gTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATGCCGGTTAACTCCGGCGAGTAGTTAATGCCGTAGACCAGTa                                                                >  1:185022/1‑149 (MQ=255)
                                                                                      ttCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATGCCGGTTAACTCCGGCGAGTAGTTAATGCCGTAGACCAGTAt                                                               >  1:22868/1‑149 (MQ=255)
                                                                                        ctcGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATGCCGGTTAACTCCGGCGAGTAGTTAATGCCGTAGACCTGTAttt                                                             <  2:337914/149‑1 (MQ=255)
                                                                                        ctcGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATGCCGGTTAACTCCGGCGAGTAGTTAATGCCGTAGACCAGTAttt                                                             >  2:465186/1‑149 (MQ=255)
                                                                                         tcGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATGCCGGTTAACTCCGGCGAGTAGTTAATGCCGTAGACCAGTAtttt                                                            <  2:675973/149‑1 (MQ=255)
                                                                                                 cTGCCATTGCGGGTAGTAAGGCGGTGCGGTAATGACCCGCACCTCATGCCCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATGCCGGTTAACTCCGGCGAGTAGTTAATGCCGTAGACCAGTATTTTCATAATcc                                                    >  1:666958/1‑149 (MQ=255)
                                                                                                   gCCATTGCGGGTAGTAAGGCGGTGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATGCCGGTTAACTCCGGCGAGTAGTTAATGCCGTAGACCAGTATTTTCATAATCCgg                                                  >  2:316028/1‑149 (MQ=255)
                                                                                                   gCCATTGCGGGTAGTAAGGCGGTGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCAGATGCCGGTTAACTCCGGCGAGTAGTTAATGCCGTAGACCAGTATTTTCATAATCCgg                                                  >  1:641194/1‑149 (MQ=255)
                                                                                                      aTTGCGGGTAGTAAGGCGGTGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATGCCGGTTAACTCCGGCGAGTAGTTAATGCCGTAGACCAGTATTTTCATAATCCGGGTa                                               >  1:838454/1‑149 (MQ=255)
                                                                                                         gCGGGTAGTAAGGCGGTGCGGTAATGCCCCGCACCTCATGACCCTGTGCCGCCAGCCATTCCACCCTCTCGCCGGTGTATTTGCCGATGCCGGTTAACTCCGGCGAGTAGTTAATGCCGTAGACCAGTATTTTCATAATCCGGGTACTc                                            <  2:955608/149‑1 (MQ=255)
                                                                                                         gCGGGTAGTAAGGCGGTGCGGTAATGACCCGCACCTCATGACCCTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATGCCGGTTAACTCCGGCGAGTAGTTAATGCCGTAGACCAGTATTTTCATAATCCGGGTACTc                                            <  2:960170/149‑1 (MQ=255)
                                                                                                           gggTAGTAAGGCGGTGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATGCCGGTTAACTCCGGCGAGTAGTTAATGCCGTAGACCAGTATTTTCATAATCCGGGTACTCCg                                          <  1:954878/149‑1 (MQ=255)
                                                                                                           gggTAGTAAGGCGGTGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAACCATTCCACCATCTCGCCGGTGTATTTGCCGATGCCGGTTAACTCCGGCGAGTAGTTAATGCCGTAGACCAGTATTTTCATAATCCGGGTACTCCg                                          <  2:49254/149‑1 (MQ=255)
                                                                                                              tagtaAGGCGGTGCGGTAATGACCCGCTCCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATGCCGGTTAACTCCGGCGAGTAGTTAATGCCGTAGACCAGTATTTTCATAATCCGGGTACTCCGGTa                                       >  1:510500/1‑149 (MQ=255)
                                                                                                              tagtaAGGCGGTGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATGCCGGTTAACTCCGGCGAGTAGTTAATGCCGTAGACCAGTATTTTCATAATCCGGGTACTCCGGTa                                       <  1:635432/149‑1 (MQ=255)
                                                                                                                 taAGGCGGTGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATGCCGGTTAACTCCGGCGAGTAGTTAATGCCGTAGACCAGTATTTTCATAATCCGGGTACTCCGGTACGc                                    <  2:185022/149‑1 (MQ=255)
                                                                                                                    ggCGGTGCGGTAATGACCCGCGCCTCAAGCCCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGAGGCCGGTTTACTCCGGCGAGTAGTTAATGCCGTAGACCAGGATTTTGATAATCCGGGTTCTCCGGTACGCTtc                                 >  1:662179/1‑149 (MQ=255)
                                                                                                                      cggtgcggtAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCCTCTCGCCGGTGTATTTGCCGATGCCGGTTAACTCCGGCGAGTAGTTAATGCCGTAGACCAGTATTTTCATAATCCGGGTACTCCGGTACGCTtctc                               <  2:163266/149‑1 (MQ=255)
                                                                                                                       tgtgaggtAATGACGCGCAGCTCATGACCTTTTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATGCCGGTTAACTCCGGCGAGTAGTTAATGCCGTAGACCAGTATTTTCATAATCCGGGTACTCCGGTACGCTTCTCa                              <  1:51063/148‑1 (MQ=255)
                                                                                                                       ggtgcggtAATGACCCGCACCGCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATGCCGGTTAACTCCGGCGAGTAGTTAATGCCGTAGACCAGTATTTTCATAATCCGGGTACTCCGGTACGCTTCTCa                              <  2:1070703/149‑1 (MQ=255)
                                                                                                                         tgcggtAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATGCCGGTTAACTCCGGCGAGTAGTTAATGCCGTAGACCAGTATTTTCATAATCCGGGTACTCCGGTACGCTTCTCAGc                            <  1:754281/149‑1 (MQ=255)
                                                                                                                         tgcggtAATGACCCGCACCGCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATGCCGGTTAACTCCGGCGAGTAGTTAATGCCGTAGACCAGTATTTTCATAATCCGGGTACTCCGGTACGCTTCTCAGc                            <  2:481671/149‑1 (MQ=255)
                                                                                                                           cggtAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGGATTTGCCGATGCCGGTTTACTCCGGCGAGTAGTTAATGCCGTAGCCCAGTATTTTCATAATCCGGGTACTCCGGTACGCTTCTCAGCga                          >  1:343862/1‑149 (MQ=255)
                                                                                                                            ggtAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATGCCGGTTAACTCCGGCGAGTAGTTAATGCCGTAGACCAGTATTTTCATAATCCGGGTACTCCGGTACGCTTCTCAGCgag                         <  2:609934/149‑1 (MQ=255)
                                                                                                                             gtAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATGGCGGTTAACTCCGGGGAGTAGTTAATGACGTAGACCTGTATTTTGATAATCAGGGTACTCCGGTACTCTTCACAGCgcga                        >  2:664367/1‑149 (MQ=255)
                                                                                                                                  gACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATGCCGGTTAACTCCGGCGAGTAGTTAATGCCGTAGACCAGTATTTTCATAATCCGGGTACTCCGGTACGCTTCTCAGCGAGAAAATa                   <  1:581109/149‑1 (MQ=255)
                                                                                                                                   aCCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTTCCGATGCCGGTTAACTCCGGCGAGTAGTTAATGCCGTAGACCAGTATTTTCATAATCCGGGTACTCCGGTACGCTTCTCAGCGAGAAAATAg                  <  1:265760/149‑1 (MQ=255)
                                                                                                                                   aCCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATGTGCCGATGCCGGTTAACTCCGGCGAGTAGTTAATGCCGTAGACCAGTATTTTCATAATCCGGGTACTCCGGTACGCTTCTCAGCGAGAAAATAg                  <  1:897076/149‑1 (MQ=255)
                                                                                                                                    cccGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATGCCGGTTAACTCCGGCGAGTAGTTAATGCCGTAGACCAGTATTTTCATAATCCGGGTACTCCGGTACGCTTCTCAGCGAGAAAATAgg                 >  2:178072/1‑149 (MQ=255)
                                                                                                                                     ccGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATGCCGGTTAACTCCGGCGAGTAGTTAATGCCGTAGACCAGTATTTTCATAATCCGGGTACTCCGGTACGCTTCTCAGCGAGAAAATAGgc                >  2:749616/1‑149 (MQ=255)
                                                                                                                                     ccGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGAAGCCGGTTAAATCAGGCGAGTAGATAATGACGTAGACCAGTATTTTCATAATCCGGGTACTCCCGTACGCTTCTCAGCGAGAAAATAGgc                >  1:160766/1‑149 (MQ=255)
                                                                                                                                         aCCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATGCCGGTTAACTCCGGCGAGTAGTTAATGCCGTAGACCAGTATTTTCATAATCCGGGTACTCCGGTACGCTTCTCAGCGAGAAAATAGGCGCgg            <  1:154507/149‑1 (MQ=255)
                                                                                                                                           cTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATGCCGGTTAACTCCGGCGAGTAGTTAATGCCGTAGACCAGTATTTTCATAATCCGGGTACTCCGGTACGCTTCTCAGCg                           >  1:1125804/1‑132 (MQ=255)
                                                                                                                                           cTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATGCCGGTTAACTCCGGCGAGTAGTTAATGCCGTAGACCAGTATTTTCATAATCCGGGTACTCCGGTACGCTTCTCAGCg                           <  2:1125804/132‑1 (MQ=255)
                                                                                                                                           cTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATGCCGGTTAACTCCGGCGAGTAGTTAATGCCGTAGACCAGTATTTTCATAATCCGGGTACTCCGGTACGCTTCTCAGCGAGAAAATAGGCGCGGCt          <  2:1136482/149‑1 (MQ=255)
                                                                                                                                           cTCATGACCTTGTGCCGCCAGCCATTCCACCATCACGCCGGTGTATTAGCCGATGCCGGTTAACTCCGGCGAGTAGTTAATGCCGTAGACCAGTATTTTCATAATCCGGGTACTCCGGTACGCTTCTCAGCGAGAAAATAGGAGCGGCt          >  1:371081/1‑149 (MQ=255)
                                                                                                                                              aTGTCCTAGTGCCGCCAGCCATTCCAGCATCTCGCCGGTGTATTAGCCTATGCCGGTTAACTCCGGCGAGTAGTTAATGCCGTAGACCAGTATTTTCTTAATCCTGGTACTCCGGTACGCTTCTCAGGGAGAAAATAGGCGCGGCTGtt       <  1:89577/149‑1 (MQ=255)
                                                                                                                                               tGACCTTGTGCCGCCAGCCATGCCACCATCTCGCCGGTGTATTTGCCGATGCCGGTTAACTCCGGCGAGTAGTTAATGCCGTAGACCAGTATTTTCATAATCCGGGCACTCCGGTACGCTTCTCAGCGAGAAAATAGGCGCGGCTGTTa      <  1:118592/149‑1 (MQ=255)
                                                                                                                                                gACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATGCCGGTTAACTCCGGCGAGTAGTTAATGCCGTAGACCAGTATTTTCATAATCCGGGTACTCCGGTACGCTTCTCAGCGAGAAAATAGGCGCGGCTGTTAg     <  1:951979/149‑1 (MQ=255)
                                                                                                                                                gACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATGCCGGTTAACTCCGGCGAGTAGTTAATGCCGTAGACCAGTATTTTCATAATCCGGGTACTCCGGTACGCTTCTCAGCGAGAAAATAGGCGCGGCTGTTAg     <  1:405720/149‑1 (MQ=255)
                                                                                                                                                gACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATGCCGGTTAACTCCGGAGAGTAGTTAATGACGTAGACCAGTATTTTCATAATCCGGGTACGCCGGTACGCTTCTCAGCGAGAAAATAGGCGCGGCTGTTAg     >  1:790248/1‑149 (MQ=255)
                                                                                                                                                   cTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATGCCGGTTAACTCCGGCGAGTAGTTAATGCCGTAGCCCAGTATTTTCATAATCCGGGTACTCCGGTACGCTTCTCAGCGAGAAAATAGGCGCGGCTGTTAGCAt  >  1:503624/1‑149 (MQ=255)
                                                                                                                                                   cTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATGCCGGTTAACTCCGGCGAGTAGTTAATGCCGTAGACCAGTATTTTCATAATCCGGGTACTCCGGTACGCTTCTCAGCGAGAAAATAGGCGCGGCTGTTAGCAt  >  2:896829/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GGGTGCTCGGCTGTTTTGGCACATACAGCGGGCAGCGCCACACCGTGGCGGCCCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATGCCGGTTAACTCCGGCGAGTAGTTAATGCCGTAGACCAGTATTTTCATAATCCGGGTACTCCGGTACGCTTCTCAGCGAGAAAATAGGCGCGGCTGTTAGCAT  >  NC_000913/2125506‑2125801

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: ‑