Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 475,972 C→G 7.3% intergenic (+21/+10) ybaY → / ← ybaZ outer membrane lipoprotein/excision repair protein, alkyltransferase‑like protein ATL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913475,9720CG 7.3% 380.0 / 15.7 109intergenic (+21/+10)ybaY/ybaZouter membrane lipoprotein/excision repair protein, alkyltransferase‑like protein ATL
Reads supporting (aligned to +/- strand):  ref base C (53/48);  new base G (4/4);  total (57/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.67e-01

Only 100 of 110 total aligned reads displayed.

ATTGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAAC  >  NC_000913/475825‑476113
                                                                                                                                                   |                                                                                                                                             
aTTGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGccc                                                                                                                                              >  3:633262/1‑149 (MQ=255)
      gCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGGGCGGcgag                                                                                                                                        <  2:90267/149‑1 (MQ=255)
       ccGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAgg                                                                                                                                       <  1:540554/149‑1 (MQ=255)
        cGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGGGCGGAGAggg                                                                                                                                      <  1:580085/149‑1 (MQ=255)
        cGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAggg                                                                                                                                      <  3:73130/149‑1 (MQ=255)
        cGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAggg                                                                                                                                      <  1:319210/149‑1 (MQ=255)
         ggTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGt                                                                                                                                     <  3:42523/149‑1 (MQ=255)
          gTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGGGCGGAGAGGGtt                                                                                                                                    <  1:337591/149‑1 (MQ=255)
          gTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGtt                                                                                                                                    >  1:53572/1‑149 (MQ=255)
           tACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTta                                                                                                                                   <  1:368327/149‑1 (MQ=255)
            aCAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTtag                                                                                                                                  >  4:621193/1‑149 (MQ=255)
              aGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTtagta                                                                                                                                <  1:311731/149‑1 (MQ=255)
              aGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTtagta                                                                                                                                <  2:913902/149‑1 (MQ=255)
               gCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTATCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGTAGAGGGTtagtac                                                                                                                               <  4:680866/149‑2 (MQ=255)
               gCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTtagtag                                                                                                                               >  3:340990/1‑149 (MQ=255)
                cAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTagtagt                                                                                                                              >  2:636954/1‑149 (MQ=255)
                    ccgccgTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCa                                                                                                                          <  1:74535/149‑1 (MQ=255)
                    ccgccgTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCa                                                                                                                          <  2:53572/149‑1 (MQ=255)
                        cgTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGa                                                                                                                      <  2:262354/149‑1 (MQ=255)
                          tGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATa                                                                                                                    >  1:485471/1‑149 (MQ=255)
                                 tCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGc                                                                                                             <  2:868189/149‑1 (MQ=255)
                                   aGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTtagta                                                                                                                                <  3:697388/128‑1 (MQ=255)
                                   aGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTtagta                                                                                                                                >  4:697388/1‑128 (MQ=255)
                                   aGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCaa                                                                                                           <  2:244619/149‑1 (MQ=255)
                                    ggCCAGCGGTGGCGCAACGACTACCGTACCGTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAg                                                                                                          <  3:853065/149‑1 (MQ=255)
                                       cAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCg                                                                                                       >  3:59519/1‑149 (MQ=255)
                                        aGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGa                                                                                                      >  1:262186/1‑149 (MQ=255)
                                         gCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGAt                                                                                                     <  1:1001609/149‑1 (MQ=255)
                                           ggTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCGGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGAttt                                                                                                   <  1:488186/149‑1 (MQ=255)
                                           ggTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCAAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGAttt                                                                                                   <  2:1039675/149‑1 (MQ=255)
                                             tGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGc                                                                                                 <  2:306349/149‑1 (MQ=255)
                                               gcgcAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGccc                                                                                               <  1:621109/149‑1 (MQ=255)
                                               gcgcAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGccc                                                                                               >  3:242587/1‑149 (MQ=255)
                                                 gcAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGc                                                                                             >  1:411187/1‑149 (MQ=255)
                                                 gcAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGc                                                                                             >  1:428818/1‑149 (MQ=255)
                                                 gcAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGc                                                                                             <  2:411187/149‑1 (MQ=255)
                                                 gcAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGc                                                                                             >  1:227651/1‑149 (MQ=255)
                                                    aCGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGTCCGCTTc                                                                                          >  3:250910/1‑149 (MQ=255)
                                                          aCCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATa                                                                                    <  3:785044/149‑1 (MQ=255)
                                                           ccGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATAc                                                                                   >  2:1047570/1‑149 (MQ=255)
                                                                 cTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATcac                                                                             >  3:83809/1‑149 (MQ=255)
                                                                 cTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATcac                                                                             <  2:1008000/149‑1 (MQ=255)
                                                                 cTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGAGATAACGCTGCAAGTCGAtt                                                                                                    >  3:119877/1‑126 (MQ=255)
                                                                 cTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGAGATAACGCTGCAAGTCGAtt                                                                                                    <  4:119877/126‑1 (MQ=255)
                                                                  ttCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATcaca                                                                            >  4:783662/1‑149 (MQ=255)
                                                                      aCTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCtt                                                                        >  3:431486/1‑149 (MQ=255)
                                                                        ttCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCt                                                                      <  2:763917/149‑1 (MQ=255)
                                                                        ttCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCt                                                                      >  3:48794/1‑149 (MQ=255)
                                                                        ttCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCt                                                                      <  2:835824/149‑1 (MQ=255)
                                                                          cACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGGGAGGGTTAGTAGTACCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCGCACCTTCTGc                                                                    >  1:764808/1‑149 (MQ=255)
                                                                          cACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGc                                                                    <  4:48794/149‑1 (MQ=255)
                                                                           aCCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCacac                                                                           >  2:777899/1‑141 (MQ=255)
                                                                           aCCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCacac                                                                           <  1:777899/141‑1 (MQ=255)
                                                                            ccAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCa                                                                  >  4:201133/1‑149 (MQ=255)
                                                                              aaCTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGt                                                                >  1:825255/1‑149 (MQ=255)
                                                                                 tCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAAt                                                             >  3:703296/1‑149 (MQ=255)
                                                                                   aGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGc                                                           >  2:450986/1‑149 (MQ=255)
                                                                                   aGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGc                                                           >  2:109541/1‑149 (MQ=255)
                                                                                     gTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGTTTTGCCCGCTTCCCGATCCCATCACACCTTCTGCCAGTAATGCCt                                                         >  2:483711/1‑149 (MQ=255)
                                                                                     gTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCt                                                         <  1:293026/149‑1 (MQ=255)
                                                                                     gTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCt                                                         >  4:259939/1‑149 (MQ=255)
                                                                                      tGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTg                                                        <  1:169081/149‑1 (MQ=255)
                                                                                             gTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGa                                                 >  3:22751/1‑149 (MQ=255)
                                                                                             gTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGa                                                 <  1:44507/149‑1 (MQ=255)
                                                                                             gTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGa                                                 <  1:113341/149‑1 (MQ=255)
                                                                                              tCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGAc                                                >  4:893034/1‑149 (MQ=255)
                                                                                              tCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGAc                                                >  1:209927/1‑149 (MQ=255)
                                                                                                 tCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCt                                             >  1:521613/1‑149 (MQ=255)
                                                                                                    gCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTa                                          >  1:1021153/1‑149 (MQ=255)
                                                                                                    gCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTa                                          >  3:810735/1‑149 (MQ=255)
                                                                                                    gCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTa                                          >  2:883795/1‑149 (MQ=255)
                                                                                                           cAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAAtccgg                                   <  3:268274/149‑1 (MQ=255)
                                                                                                                 cTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATccggtccggt                             >  3:108011/1‑149 (MQ=255)
                                                                                                                  tACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGtt                            >  3:553247/1‑149 (MQ=255)
                                                                                                                  tACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGtt                            <  1:109541/149‑1 (MQ=255)
                                                                                                                   aCGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTa                           <  3:621193/149‑1 (MQ=255)
                                                                                                                     gCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGc                         >  2:287691/1‑149 (MQ=255)
                                                                                                                     gCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGc                         >  2:515441/1‑149 (MQ=255)
                                                                                                                       aaTATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTa                           >  3:13647/1‑145 (MQ=255)
                                                                                                                       aaTATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTa                           <  4:13647/145‑1 (MQ=255)
                                                                                                                         tatTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAAGGCCTATCGCTGACGCTGTAAATCCGGTCCGGTTAGCGaaa                     >  4:685145/1‑149 (MQ=255)
                                                                                                                          atTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAAt                    >  2:878039/1‑149 (MQ=255)
                                                                                                                               gCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGc               <  4:431486/149‑1 (MQ=255)
                                                                                                                                 actactAACCCTCTCCGGGCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCg             >  4:904153/6‑149 (MQ=255)
                                                                                                                                   cGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCCGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAAtc                                      <  1:205025/122‑1 (MQ=255)
                                                                                                                                   cGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCCGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAAtc                                      >  2:205025/1‑122 (MQ=255)
                                                                                                                                    gTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTgg          >  4:902473/1‑149 (MQ=255)
                                                                                                                                    gTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTgg          >  3:439065/1‑149 (MQ=255)
                                                                                                                                     ttAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGc         <  4:343945/149‑1 (MQ=255)
                                                                                                                                      tAACCCTCTCCGGGCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCg        >  1:279139/1‑149 (MQ=255)
                                                                                                                                       aaCCCTCTCCGGGCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACGCCTTTTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGa       >  3:990820/1‑149 (MQ=255)
                                                                                                                                       aaCCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGa       <  4:22751/149‑1 (MQ=255)
                                                                                                                                        aCCCTCTCCGGGCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGAt      >  1:541438/1‑149 (MQ=255)
                                                                                                                                        aCCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGAt      >  2:244873/1‑149 (MQ=255)
                                                                                                                                         cccTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGTGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGAtt     >  3:979990/1‑149 (MQ=255)
                                                                                                                                           ctctCCGCCCGGAGAGGGATAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTaa   >  4:917463/1‑149 (MQ=255)
                                                                                                                                            tctcCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAac  <  2:485471/149‑1 (MQ=255)
                                                                                                                                            tctcCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAac  <  4:871353/149‑1 (MQ=255)
                                                                                                                                            tctcCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAac  <  2:818027/149‑1 (MQ=255)
                                                                                                                                            tctcCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAac  >  3:171032/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                             
ATTGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAAC  >  NC_000913/475825‑476113

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: