Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,572,483 G→T 5.4% intergenic (‑33/+6) eutM ← / ← eutD ethanolamine utilization protein, putative carboxysome structural protein/phosphate acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,572,4830GT 5.4% 383.9 / 9.5 111intergenic (‑33/+6)eutM/eutDethanolamine utilization protein, putative carboxysome structural protein/phosphate acetyltransferase
Reads supporting (aligned to +/- strand):  ref base G (53/52);  new base T (3/3);  total (56/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.80e-01

Only 100 of 114 total aligned reads displayed.

CGCCGCCAATCTGCTTAACACCAACCAGCTTGACGCGCGCAGCTTTCACCATCGCGTCAGAAGCCTCAATCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGGCACAGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCACGAGAGAGATCGTGCATCGGCGCGGCAAGTCCTTGTATC  >  NC_000913/2572336‑2572630
                                                                                                                                                   |                                                                                                                                                   
cgccgccAATCTGCTTAACACCAACCAGCTTGACGCGCGCAGCTTTCACCATCGCGTCAGAAGCCTCAATCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACggg                                                                                                                                                     <  4:490092/148‑1 (MQ=255)
cgccgccAATCTGCTTAACACCAACCAGCTTGACGCGCGCAGCTTTCACCATCGCGTCAGAAGCCTCAATCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGa                                                                                                                                                    <  2:810331/149‑1 (MQ=255)
  ccgccAATCTGCTTAACACCAACCAGCTTGACGCGCGCAGCTTTCACCATCGCGTCAGAAGCCTCAATCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGAcc                                                                                                                                                  <  2:332128/149‑1 (MQ=255)
  ccgccAATCTGCTTAACACCAACCAGCTTGACGCGCGCAGCTTTCACCATCGCGTCAGAAGCCTCAATCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGAc                                                                                                                                                   >  3:295727/1‑148 (MQ=255)
    gccAATCTGCTTAACACCAACCAGCTTGACGCGCGCAGCTTTCACCATCGCGTCAGAAGCCTCAATCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCg                                                                                                                                                 <  3:945801/148‑1 (MQ=255)
        aTCTGCTTAACACCAACCAGCTTGACGCGCGCAGCTTTCACCATCGCGTCAGAAGCCTCAATCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGttcat                                                                                                                                            <  1:97144/149‑1 (MQ=255)
         tCTGCTTAACCCCAACCAGCTTGACGCGCGCAGCTTTCACCATCGCGTCAGAAGCCTCAATCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGttcat                                                                                                                                            >  3:846172/1‑148 (MQ=255)
         tCTGCTTAACACCAACCAGCTTGACGCGCGCAGCTTTCACCATCGCGTCAGAAGCCTCAATCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGTCCCGttct                                                                                                                                             >  4:294204/1‑146 (MQ=255)
         tCTGCTTAACACCAACCAGCTTGACGCGCGCAGCTTTCACCATCGCGTCAGAAGCCTCAATCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGTCCCGttct                                                                                                                                             <  3:294204/147‑2 (MQ=255)
           tGCTTAACACCAACCAGCTTGACGCGCGCAGCTTTCACCATCGCGTCAGAAGCCTCAATCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGttcattca                                                                                                                                         >  2:349828/1‑149 (MQ=255)
            gCTTAACACCAACCAGCTTGACGCGCGCAGCTTTCACCATCGCGTCAGAAGCCTCAATCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGttcattca                                                                                                                                         <  3:544967/148‑1 (MQ=255)
             cTTAACACCAACCAGCTTGACGCGCGCAGCTTTCACCATCGCGTCAGAAGCCTCAATCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGTCCCGttctggacc                                                                                                                                        <  4:129280/148‑7 (MQ=255)
             cTTAACACCAACCAGCTTGACGCGCGCAGCTTTCACCATCGCGTCAGAAGCCTCAATCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCaa                                                                                                                                        <  4:354535/148‑1 (MQ=255)
               tAACACCAACCAGCTTGACGCGCGCAGCTTTCACCATCGCGTCAGAAGCCTCAATCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAAcc                                                                                                                                      >  3:280384/1‑148 (MQ=255)
               tAACACCAACCAGCTTGACGCGCGCAGCTTTCACCATCGCGTCAGAAGCCTCAATCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAAcc                                                                                                                                      >  3:166069/1‑148 (MQ=255)
                aaCACCAACCAGCTTGACGCGCGCAGCTTTCACCATCGCGTCAGAAGCCTCAATCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAg                                                                                                                                    >  2:988422/1‑149 (MQ=255)
                 acacCAACCAGCTTGACGCGCGCAGCTTTCACCATCGCGTCAGAAGCCTCAATCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAgt                                                                                                                                   <  1:159533/149‑1 (MQ=255)
                 acacCAACCAGCTTGACGCGCGCAGCTTTCACCATCGCGTCAGAAGCCTCAATCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAg                                                                                                                                    <  3:257248/148‑1 (MQ=255)
                  cacCAACCAGCTTGACGCGCGCAGCTTTCACCATCGCGTCAGAAGCCTCAATCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAgtg                                                                                                                                  >  2:746895/1‑149 (MQ=255)
                  cacCAACCAGCTTGACGCGCGCAGCTTTCACCATCGCGTCAGAAGCCTCAATCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAgt                                                                                                                                   >  4:926194/1‑148 (MQ=255)
                   accaaccaGCTTGACGCGCGCAGCTTTCACCATCGCGTCAGAAGCCTCAATCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAgtgt                                                                                                                                 >  2:761612/1‑149 (MQ=255)
                   accaaccaGCTTGACGCGCGCAGCTTTCACCATCGCGTCAGAAGCCTCAATCAGTGCAACCAGGCCCCGGGTTTCGATCACTCCTAATGCTTCCATTGTGCTATCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAgtg                                                                                                                                  >  3:911641/1‑148 (MQ=255)
                    ccaaccaGCTTGACGCGCGCAGCTTTCACCATCGCGTCAGAAGCCTCAATCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAgtgt                                                                                                                                 >  4:382388/1‑148 (MQ=255)
                     caaccaGCTTGACGCGCGCAGCTTTCACCATCGCGTCAGAAGCCTCAATCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGtt                                                                                                                                <  3:1022618/148‑1 (MQ=255)
                     caaccaGCTTGACGCGCGCAGCTTTCACCATCGCGTCAGAAGCCTCAATCAGTGCAACCAGGCCCCAGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGttt                                                                                                                               <  2:596929/149‑1 (MQ=255)
                       accaGCTTGACGCGCGCAGCTTTCACCATCGCGTCAGAAGCCTCAATCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTg                                                                                                                              >  4:625891/1‑148 (MQ=255)
                        ccaGCTTGACGCGCGCAGCTTTCACCATCGCGTCAGAAGCCTCAATCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTa                                                                                                                            >  2:710518/1‑149 (MQ=255)
                          aGCTTGACGCGCGCAGCTTTCACCATCGCGTCAGAAGCCTCAATCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTaa                                                                                                                           >  3:104904/1‑148 (MQ=255)
                              tGACGCGCGCAGCTTTCACCATCGCGTCAGAAGCCTCAATCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGTCCCGttctgga                                                                                                                                          <  1:19322/129‑5 (MQ=255)
                              tGACGCGCGCAGCTTTCACCATCGCGTCAGAAGCCTCAATCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGTCCCGttctgga                                                                                                                                          >  2:19322/1‑125 (MQ=255)
                              tGACGCGCGCAGCTTTCACCATCGCGTCAGAAGCCTCAATCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTg                                                                                                                       >  4:65283/1‑148 (MQ=255)
                                aCGCGCGCAGCTTTCACCATCGCGTCAGAAGCCTCAATCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCt                                                                                                                     <  3:65283/148‑1 (MQ=255)
                                aCGCGCGCAGCTTTCACCATCGCGTCAGAAGCCTCAATCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCt                                                                                                                     <  3:926194/148‑1 (MQ=255)
                                         gCTTTCACCATCGCGTCAGAAGCCTCAATCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGtt                                                                                                            <  4:93567/148‑1 (MQ=255)
                                          cTTTCACCATCGCGTCAGAAGCCTCAATCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCa                                                                                                          <  1:588924/149‑1 (MQ=255)
                                            ttCACCATCGCGTCAGAAGCCTCAATCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCAc                                                                                                         >  3:742574/1‑148 (MQ=255)
                                            ttCACCATCGCGTCAGAAGCCTCAATCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCAc                                                                                                         >  4:946121/1‑148 (MQ=255)
                                                ccATCGCGTCAGAAGCCTCAATCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTtct                                                                                                     >  4:37298/1‑148 (MQ=255)
                                                ccATCGCGCCAGAAGCCTCAATCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCtt                                                                                                                    >  1:163872/1‑133 (MQ=255)
                                                ccATCGCGCCAGAAGCCTCAATCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCtt                                                                                                                    <  2:163872/133‑1 (MQ=255)
                                                 cATCGCGTCAGAAGCCTCAATCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTtctgt                                                                                                   <  1:54170/149‑1 (MQ=255)
                                                  aTCGCGTCAGAAGCCTCAATCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTtctgt                                                                                                   >  3:853953/1‑148 (MQ=255)
                                                  aTCGCGTCAGAAGCCTCAATCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTtctgt                                                                                                   <  3:524777/148‑1 (MQ=255)
                                                   tCGCGTCAGAAGCCTCAATCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTtctgtct                                                                                                 >  2:435714/1‑149 (MQ=255)
                                                   tCGCGTCAGAAGCCTCAATCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTtctgtc                                                                                                  >  4:387826/1‑148 (MQ=255)
                                                      cgTCAGAAGCCTCAATCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGa                                                                                               >  4:583792/1‑148 (MQ=255)
                                                        tCAGAAGCCTCAATCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGAcgc                                                                                            <  2:244473/149‑1 (MQ=255)
                                                            aaGCCTCAATCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGGc                                                                                         >  4:734468/1‑148 (MQ=255)
                                                             aGCCTCAATCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGGcac                                                                                       <  1:445327/149‑1 (MQ=255)
                                                                 tCAATCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGGCACAGCt                                                                                   >  1:1002295/1‑149 (MQ=255)
                                                                  cAATCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACtt                                                                                                       <  4:118038/128‑1 (MQ=255)
                                                                  cAATCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACtt                                                                                                       >  3:118038/1‑128 (MQ=255)
                                                                   aaTCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCATTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGGCACAGCTGc                                                                                 <  1:761612/149‑1 (MQ=255)
                                                                     tCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGGCACAGCTGcc                                                                                >  3:815839/1‑148 (MQ=255)
                                                                        gTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGGCACAGCTGccacca                                                                            >  2:48168/1‑149 (MQ=255)
                                                                             aCCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGGCACAGCTGCCACCAgcgc                                                                        >  3:605900/1‑148 (MQ=255)
                                                                                 ggCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGGCACAGCTGCCACCAGCGCCAGc                                                                    >  4:433104/1‑148 (MQ=255)
                                                                                     ccGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGGca                                                                                        >  3:858121/1‑124 (MQ=255)
                                                                                     ccGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGGca                                                                                        <  4:858121/124‑1 (MQ=255)
                                                                                          tttCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTtctg                                                                                                    >  2:456372/1‑107 (MQ=255)
                                                                                           ttCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGGCACAGCTGCCACCAGCGCCAGCTCGATAATTTc                                                         <  1:349828/149‑1 (MQ=255)
                                                                                           ttCGATCACTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGGCACAGCTGCCACCAGCGCCAGCTCGATAAttt                                                          <  4:911641/148‑1 (MQ=255)
                                                                                            tCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGGCACAGCTGCCACCAGCGCCAGCTCGATAATTTc                                                         <  3:964768/148‑1 (MQ=255)
                                                                                             cGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGGCACAGCTGCCACCAGCGCCAGCTCGATAATTTcc                                                        <  3:547934/148‑1 (MQ=255)
                                                                                             cGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGGCACAGCTGCCACCAGCGCCAGCTCGATAATTTCCt                                                       <  2:204907/149‑1 (MQ=255)
                                                                                              gATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGGCACAGCTGCCACCAGCGCCAGCTCGATAATTTCCTg                                                      >  2:477246/1‑149 (MQ=255)
                                                                                                   ttCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGGCACAGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACg                                                  <  4:367459/148‑1 (MQ=255)
                                                                                                        aaTGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGGCACAGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTACa                                             >  4:1005846/1‑148 (MQ=255)
                                                                                                         aTGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGGCACAGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTACaa                                            >  3:308553/1‑148 (MQ=255)
                                                                                                            cTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGGCACAGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCa                                         >  3:377112/1‑148 (MQ=255)
                                                                                                            cTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTAGCGGTTCACTTCTGTCTGACGCGGCACAGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCAc                                        >  1:813924/1‑149 (MQ=255)
                                                                                                              tCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGGCACAGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCACga                                      <  1:48168/149‑1 (MQ=255)
                                                                                                              tCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGGCACAGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCACg                                       >  4:614359/1‑148 (MQ=255)
                                                                                                               ccATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGGCACAGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCACgag                                     >  1:53087/1‑149 (MQ=255)
                                                                                                                cATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGGCACAGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCACgaga                                    <  2:360368/149‑1 (MQ=255)
                                                                                                                cATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTCTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGGCACAGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCACgaga                                    <  1:187410/149‑1 (MQ=255)
                                                                                                                 aTTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGGCACAGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCACgagag                                   <  1:499563/149‑1 (MQ=255)
                                                                                                                    gtgCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGGCACAGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCACgagagag                                 <  3:583792/148‑1 (MQ=255)
                                                                                                                     tgCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGGCACAGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCACGAGAGAGAt                               <  1:477246/149‑1 (MQ=255)
                                                                                                                       cTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGGCACAGCTGCCACCAGCGCCAGCTCGATAATTTCCTTCACGCTACAACCACGAGAGAGATc                              >  3:525905/1‑148 (MQ=255)
                                                                                                                        tttCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGGCACAGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCACGAGAGAGATCg                             <  3:820868/148‑1 (MQ=255)
                                                                                                                           ccTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGGCACAGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCACGAGAGAGATCGTGc                          <  4:605900/148‑1 (MQ=255)
                                                                                                                            ctctTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGGCACAGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCACGAGAGAGATCGTGCAt                        >  1:812366/1‑149 (MQ=255)
                                                                                                                             tctTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGGCACAGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGc                                                 >  1:762682/1‑123 (MQ=255)
                                                                                                                             tctTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGGCACAGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCACGAGAGAGATCGt                            <  1:850757/144‑1 (MQ=255)
                                                                                                                             tctTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGGCACAGCTGCCACCAGCGCCAGCTCGATAATTTCCCGCACGc                                                 <  2:762682/123‑1 (MQ=255)
                                                                                                                             tctTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGAGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGGCACAGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCACGAGAGAGATCGt                            >  2:850757/1‑144 (MQ=255)
                                                                                                                                ttATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGGCACAGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCACGAGAGAGATCGTGCATCGgc                    >  2:709296/1‑149 (MQ=255)
                                                                                                                                  aTCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGGCACAGCTGc                                                                                 <  1:831096/86‑1 (MQ=255)
                                                                                                                                  aTCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGGCACAGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCACGAGAGAGATCGTGCATCGgcgc                  >  2:371918/1‑149 (MQ=255)
                                                                                                                                  aTCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGGCACAGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCACGAGAGAGATCGTGCATCGgcgc                  <  1:842161/149‑1 (MQ=255)
                                                                                                                                  aTCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGGCACAGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCACGAGAGAGATCGTGCATCGgcgc                  <  1:188633/149‑1 (MQ=255)
                                                                                                                                   tCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGGCACAGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCACGAGAGAGATCGTGCATCGgcgcc                 <  2:911415/149‑2 (MQ=255)
                                                                                                                                        gTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGGCACAGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCACGAGAGAGATCGTGCATCGGCGCGGCaa             <  3:384669/148‑1 (MQ=255)
                                                                                                                                        gTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGGCACAGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCACGAGAGAGATCGTGCATCGGCGCGGCaa             >  4:927071/1‑148 (MQ=255)
                                                                                                                                             gAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGGCACAGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCACGAGAGAGATCGTGCATCGGCGCGGCAAGTCCtt       <  1:543989/149‑1 (MQ=255)
                                                                                                                                             gAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGGCACAGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCACGAGAGAGATCGTGCATCGGCGCGGCAAGTCCt        <  4:801814/148‑1 (MQ=255)
                                                                                                                                               aCGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGGCACAGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCACGAGAGAGATCGTGCATCGGCGCGGCAAGTCCTTGt     <  2:812366/149‑1 (MQ=255)
                                                                                                                                                 gggACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGGCACAGCTGCCACCAGCGCCAGCTCGATAAGTTCCTGCACGCTACAACCACGAGAGAGATCGTGCATCGGCGCGGCAAGTCCTTGTa    <  3:980812/148‑1 (MQ=255)
                                                                                                                                                   gACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGGCACAGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCACGAGAGAGATCGTGCATCGGCGCGGCAAGTCCTTGTATc  <  3:47080/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CGCCGCCAATCTGCTTAACACCAACCAGCTTGACGCGCGCAGCTTTCACCATCGCGTCAGAAGCCTCAATCAGTGCAACCAGGCCCCGGGTTTCGATCATTCCTAATGCTTCCATTGTGCTTTCCTCTTTATCAGGGTCCAGAACGGGACCGTTCATTCAACCAGTGTTTGTAAACTGCTTTCGCGGTTCACTTCTGTCTGACGCGGCACAGCTGCCACCAGCGCCAGCTCGATAATTTCCTGCACGCTACAACCACGAGAGAGATCGTGCATCGGCGCGGCAAGTCCTTGTATC  >  NC_000913/2572336‑2572630

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