Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,969,429 A→C 5.8% K133N (AAA→AAC)  wzzE → Entobacterial Common Antigen (ECA) polysaccharide chain length modulation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,969,4290AC 5.8% 316.2 / 3.6 122K133N (AAA→AAC) wzzEEntobacterial Common Antigen (ECA) polysaccharide chain length modulation protein
Reads supporting (aligned to +/- strand):  ref base A (68/44);  new base C (5/3);  total (74/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 125 total aligned reads displayed.

TCCGTTCAAACATGGCTTCTGCCGACCAACCATCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTAC‑AAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATGTTGCTTTTGCC  >  NC_000913/3969281‑3969576
                                                                                                                                                     |                                                                                                                                                   
    ttCAAACATGGCTTCTGCCGACCAACCATCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTAC‑AAACAGCGGATGGTGGGCAACAGCAAAGCCg                                                                                                                                                 >  1:771385/1‑149 (MQ=255)
        aaCATGGCTTCTGCCGACCAACCATCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCGGCAAACCGACTATTAC‑CAACAGCGGATGGTGGGCAACAGCAAGGCCGATgc                                                                                                                                             >  1:1066270/1‑149 (MQ=255)
           aTGGCTTCTGCCGACCAACCATCGGTCTTGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTAC‑AAACAGCGGATGGTGGTCAACAGCAAAGCCGATgcggc                                                                                                                                          >  1:985285/1‑149 (MQ=255)
              gctgctCCCGACCAACCATCGGTCATGGACGAAGCCTATAAAGAGTTGTTTATGCAACTGGCCTCGTGGGAGACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTAC‑AAACAGCGGCTGGTGGGCAACAGCAAAGCCGATGCGGCGtt                                                                                                                                       <  1:1410565/149‑1 (MQ=255)
                ttCTGCCGACCAACCATCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTAC‑AAACAGCGGATGGTGGGCAACAGCAACGCCGATGCGGCGTTGc                                                                                                                                     >  1:638046/1‑149 (MQ=255)
                ttCTGCCGACCAACCATCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTAC‑AAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGc                                                                                                                                     >  2:47434/1‑149 (MQ=255)
                   tGCCGACCAACCATCGGTCATGCACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGGGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTAC‑AAACAGCGGATGGGGGGCAACAGCAACGCCGCATCGGCTTTGCTgt                                                                                                                                  <  2:760730/149‑2 (MQ=255)
                     ccGACCAACCATCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTAC‑AAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGAt                                                                                                                                >  2:55028/1‑149 (MQ=255)
                      cGACCAACCATCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCTCAGAGAGTTCTGGCTGCAAACCGACTATTAC‑AAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATg                                                                                                                               >  2:1268787/1‑149 (MQ=255)
                       gACCAACCATCGGTCATGGACGAAGCCTATAAAGACTTTCTTATGCAACTGGCCTAGTGGGATACCAGCAGAGAGTTATGGATGCAAACCGACTATTAA‑AAACAGCGGATGGTGGACAACCGCAAAGCCGATGCGGCGTCGCTGGATGa                                                                                                                              >  2:396883/1‑149 (MQ=255)
                        accaaccaTCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAATCCGACTATTAC‑AAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGaa                                                                                                                             >  2:591922/1‑149 (MQ=255)
                        accaaccaTCGCGAATGGAGGAAGCCTATAAAGAGTTTGTTTTGCATCTGGCCGCGTGGGATACCCGCAGCGAGTTCTGGCTGCAAACAGACTATTAC‑AAACCGCGGATGGTGGGCAACAGCAACGCCGCATCGGCTTTGCTgttgccc                                                                                                                             <  2:252966/149‑7 (MQ=255)
                         ccaaccaTCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTAC‑AAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGaaa                                                                                                                            <  1:1168270/149‑1 (MQ=255)
                         ccaaccaTCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTAC‑AAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGaaa                                                                                                                            >  2:921791/1‑149 (MQ=255)
                         ccaaccaTCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGGTCTGGCTGCAAACCGACTATTAC‑AAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGaaa                                                                                                                            <  2:901604/149‑1 (MQ=255)
                           acccaTCGGTCATGGACGCAGCCGATAAAGAGTGTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGCGTTCTGGCTGCAAACCGACTATTAC‑AACCAGCGGATGGAGGGCAACAGCCAAGCCTATGCGGCGTAGCTGGATGAAATg                                                                                                                          >  1:852873/3‑149 (MQ=255)
                           aaccaTCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACAGACTATTAC‑AAACAGCGGATGGTGCGCAACAGCAAAGCCGATGAGGCGTTGCTCGATCAAATg                                                                                                                          >  1:465315/1‑149 (MQ=255)
                             ccaTCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGACTCGTGGGAAACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTAC‑AAACAGCGGAAGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGAt                                                                                                                        <  1:232682/149‑1 (MQ=255)
                                 cGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTAC‑AAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTaac                                                                                                                    <  2:787626/149‑1 (MQ=255)
                                  ggTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTAC‑AAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTaaca                                                                                                                   <  1:921791/149‑1 (MQ=255)
                                   gTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTAC‑AAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTaacaa                                                                                                                  >  1:1129743/1‑149 (MQ=255)
                                    tCATGGACGAAGCCTATAAATAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTAA‑AAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTaacaac                                                                                                                 >  1:648507/1‑149 (MQ=255)
                                       tGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTAC‑AAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATc                                                                                                              >  2:146160/1‑149 (MQ=255)
                                          aCGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTAC‑AAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAg                                                                                                           >  1:5996/1‑149 (MQ=255)
                                           cGAAGCCTATAATGAGTTTGTTATGCAACGGGCCTCGTGGGATACCCGCAGCGAGGGCTGGCTGCAAACCGACTATTAC‑AAACAGCGGATGGGGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAgg                                                                                                          <  2:1391915/149‑2 (MQ=255)
                                           cGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGAGACACGCAGAGAGTTCTGGCTGCAAACCGACTATTAC‑AAACAGCGGATGGGGGGCAACAGCAAAGCCGATCCGGCGGTGCTGGTTGAAATGATTACCAACCTCCAGt                                                                                                          >  1:1079832/1‑149 (MQ=255)
                                            gAAGCCTAAAAAGAGGGGGCTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTAC‑AACCAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGtt                                                                                                         <  1:1001243/149‑1 (MQ=255)
                                             taGCCTATAAAGAGTTTGTTATGCAACAGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTAC‑AAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGttt                                                                                                        <  2:229765/148‑1 (MQ=255)
                                                ccTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATT‑C‑AAACAGCGGATGGTGGGCAAAAGCAACGCCGATTCGGCGTTGCTCGATGAAATGATTAACAACAACCAGTGTATcc                                                                                                    >  1:1493788/1‑149 (MQ=255)
                                                   ttaAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTAC‑AAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATcccc                                                                                                  <  2:1050138/148‑1 (MQ=255)
                                                   ataAAGAGTTTGTTATGCAACTGGCCTCGTGGGTGACCCGGAGAGCGGTCAGGCTGCAAACCGCCGATTAC‑AAACAGCGGAGGGTGGGCAACAGCAAAGCCGCTGCGGCGTTGCTGGATGAAATGAGTAACAACATCCTGTTTATcccc                                                                                                  >  2:1026870/1‑149 (MQ=255)
                                                   ataAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTAC‑AAACAGCGGATGGTGGGCAACAGCAACGCCGCATCGGCTTTGCt                                                                                                                                    <  1:880769/115‑1 (MQ=255)
                                                    taAAGAGTTTGTTATGCAACTGGCCACGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTAC‑AAACAGCGGATGGTGGGCAACAGCAAAGCCGCTGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCg                                                                                                 >  1:245307/1‑149 (MQ=255)
                                                     aaaGAGTTTGTTATGCCACTGGCCTCGTGGGATACCCGCAGAGCGTTCTGTCTGCAACCCGACCATTAC‑ACCCAGCGGAAGGTGGGCACCGGCAAAGCCGATGCGGCGTTTCTGGATGAAATGATTAACACCATCACGGTGATCCCCgg                                                                                                >  1:107812/1‑149 (MQ=255)
                                                     aaaGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTAC‑AAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACACCATCCAGTTTATCCCCgg                                                                                                >  1:1142173/1‑149 (MQ=255)
                                                     aaaGAGTTTGTTATGCAACTGGCATCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACAGACTATTAC‑AAAAAGCTGAGGGTGGGCAACAGCAAAGCCTATGCGGCGTTGCTGGATGAAAGGATTAACAACCTCCAGTTTATACCCgg                                                                                                >  2:113429/1‑149 (MQ=255)
                                                      aaGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAATCCGACTATTAC‑AAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGga                                                                                               >  1:201198/1‑149 (MQ=255)
                                                      aaGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTAC‑AAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGga                                                                                               >  2:61641/1‑149 (MQ=255)
                                                        gagTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTAC‑AAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAAGAACATCCAGTTTATCCCCGgaga                                                                                             >  1:1263748/1‑149 (MQ=255)
                                                         agTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTAC‑AAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTCGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGAc                                                                                            >  2:278244/1‑149 (MQ=255)
                                                           tttGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTAC‑AAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACtt                                                                                          <  2:1281844/149‑1 (MQ=255)
                                                            ttGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTAC‑AAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACttt                                                                                         <  1:200559/149‑1 (MQ=255)
                                                            ttGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTAC‑AAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACttt                                                                                         <  2:316899/149‑1 (MQ=255)
                                                             tGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTTTTAC‑AAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACAGCCCGTTTATCCCCGGAGACTTTa                                                                                        >  2:501122/1‑149 (MQ=255)
                                                             tGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTAC‑AAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTa                                                                                        >  1:1423618/1‑149 (MQ=255)
                                                             cgTTATGCCCCTGGCCTCGTGGGACACCCGCAGCGAGTCCTGTCTGCAAACCGACTATTAC‑ACACAGCGGATGGTTGGCCACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTa                                                                                        <  1:1073351/148‑1 (MQ=255)
                                                                    cAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTAC‑AAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCcgcgc                                                                                 <  2:665828/149‑1 (MQ=255)
                                                                     aaCTGGCCTCGTGGGATTCCCGCAGAGAGTTCTGGCTGCAAACCGACTATTAC‑AAACAGCGGATGGTGGGCACCAGCAAAGCCGAGGCGGCGTTGCTGGATGAAATGATTCACAACGTCCAGTGTATCCCCGGAGACTTTCCCcgcgcg                                                                                >  2:392957/1‑149 (MQ=255)
                                                                              cGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTAC‑AAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGAc                                                                       >  1:552000/1‑149 (MQ=255)
                                                                               gTGGGATACCCGCAGAGAGTTCTGGCTGCATACCGACTATTAC‑AAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTATCAAAATCCAGTTTATCCCAGGAGACTTTACCCGCGCGGTCAATGACa                                                                      <  2:1451849/149‑1 (MQ=255)
                                                                                  ggATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTAC‑AAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCg                                                                   >  2:1244935/1‑149 (MQ=255)
                                                                                   gATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTAC‑AAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGt                                                                  >  1:1510446/1‑149 (MQ=255)
                                                                                    aTACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTAC‑AAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTg                                                                 >  2:915005/1‑149 (MQ=255)
                                                                                     tACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTAC‑AAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGa                                                                <  1:1390136/149‑1 (MQ=255)
                                                                                       cccGCAGAGAGTTCTGGCTGCAAACCGACTATTAC‑AAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAg                                                              <  1:871466/149‑1 (MQ=255)
                                                                                            agagagTTCTGGCTGCAAACCGACTTTTAC‑AAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTAt                                                         <  1:224490/149‑1 (MQ=255)
                                                                                            agagagTTCTGGCTGCAAACCGACTATTAC‑AAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTAt                                                         <  2:1510446/149‑1 (MQ=255)
                                                                                                agTTCTGGCTGCAAACCGACTATTAC‑AAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATGCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGcc                                                     <  1:1376757/149‑1 (MQ=255)
                                                                                                  ttCTGGCTGCAAACCGACTATTAC‑AAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGa                                                   >  1:1238448/1‑149 (MQ=255)
                                                                                                    cTGGCTGCAAACCGACTATTAC‑AAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGaaa                                                 <  2:1182851/149‑1 (MQ=255)
                                                                                                     tGGCTGCAAACCGACTATTAC‑AAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCTGCGTTGCTTGTTGAAATGATTAAATACATCAAGTTTATCCACGGATACTTTACCCGCGCGTTCAATGACAGATTGAAGCTTATTGCAGAAAc                                                >  1:1017927/1‑149 (MQ=255)
                                                                                                          gCAAACCGACTATTAC‑AAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATAGCCGAAACCgcgc                                           >  1:1560193/1‑149 (MQ=255)
                                                                                                             aaCCGACTATTAC‑AAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATACCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTg                                        >  2:1543273/1‑149 (MQ=255)
                                                                                                               ccGACTATTAC‑AAACAGCGGATGGTTGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGGCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGAc                                      <  1:1501094/149‑1 (MQ=255)
                                                                                                               ccGACTATTAC‑AAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGAc                                      >  2:487889/1‑149 (MQ=255)
                                                                                                                cGACTATTAC‑AAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACg                                     >  2:703563/1‑149 (MQ=255)
                                                                                                                 gACTATTACAAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACg                                     >  2:576859/1‑149 (MQ=255)
                                                                                                                    tATTAC‑AAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCtaa                                 >  2:262902/1‑149 (MQ=255)
                                                                                                                     aTTAC‑AAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCtaat                                <  1:1347745/149‑1 (MQ=255)
                                                                                                                         c‑AAACAGCGGATGGGGGGCAACCTCAAAGCCGATGCCGCGTTGCTGGATGTAATGATTAACAACCTCCATGTTATCCCCGGTGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAAcc                            <  1:1484608/149‑1 (MQ=255)
                                                                                                                           aaaCAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGCCTTTCCCCGCGCGGGCATGGCCAGCGTGAAGCTTATTGCCGAACCCGCGCATGCGGGTAATAcccg                           >  1:1448412/1‑148 (MQ=255)
                                                                                                                           aaaCAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGaaa                                                 >  1:1141464/1‑127 (MQ=255)
                                                                                                                           aaaCAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGaaa                                                 <  2:1141464/127‑1 (MQ=255)
                                                                                                                           aaaCAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCTCCTGACGCTAATAACCt                           >  1:764578/1‑149 (MQ=255)
                                                                                                                           aaaCAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCt                           >  1:1213498/1‑149 (MQ=255)
                                                                                                                           aaaCAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGACAATGACAGCGTGAAGCTTATTGCCGACACCGCGCCTGACGCTAATAACCt                           >  1:1076866/1‑149 (MQ=255)
                                                                                                                           aaaCAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTATATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCTCCTGACGCTAATAACCt                           >  1:295380/1‑149 (MQ=255)
                                                                                                                              cAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGtt                        <  2:552000/149‑1 (MQ=255)
                                                                                                                               aGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTa                       <  1:98447/149‑1 (MQ=255)
                                                                                                                                gCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCATTGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTAc                      <  1:557588/149‑1 (MQ=255)
                                                                                                                                gCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTAc                      <  2:427862/149‑1 (MQ=255)
                                                                                                                                  ggATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGt                    >  2:1151397/1‑149 (MQ=255)
                                                                                                                                  ggATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGGCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGt                    <  2:771385/149‑1 (MQ=255)
                                                                                                                                   gATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGCAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTc                   >  1:1339345/1‑149 (MQ=255)
                                                                                                                                   gATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGCAACCGCGCCTGACGCTAATAACCTGTTACGTc                   >  1:580660/1‑149 (MQ=255)
                                                                                                                                      ggtggGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGt                <  1:278244/149‑1 (MQ=255)
                                                                                                                                      ggtggGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGt                <  2:1423618/149‑1 (MQ=255)
                                                                                                                                        tggGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTAt              >  2:356524/1‑149 (MQ=255)
                                                                                                                                         gggCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATg             >  1:113810/1‑149 (MQ=255)
                                                                                                                                          ggCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATGt            <  1:61641/149‑1 (MQ=255)
                                                                                                                                           gCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCa                  >  2:796517/1‑142 (MQ=255)
                                                                                                                                           gCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCa                  <  1:796517/142‑1 (MQ=255)
                                                                                                                                           gCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATGtt           <  1:915005/149‑1 (MQ=255)
                                                                                                                                           gCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATGtt           <  2:201198/149‑1 (MQ=255)
                                                                                                                                           gCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATGtt           >  2:1324212/1‑149 (MQ=255)
                                                                                                                                           gCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATGtt           <  1:1206893/149‑1 (MQ=255)
                                                                                                                                               cAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATGTTGCtt       <  2:62851/149‑1 (MQ=255)
                                                                                                                                                aGCAACGCCGCATCGGCTTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATGTTGCttt      >  2:951201/1‑149 (MQ=255)
                                                                                                                                                 gCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGAGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATGTTGCtttt     >  1:1417452/1‑149 (MQ=255)
                                                                                                                                                    aaGCAGATGCGGCGTGGCTGGATGAAATGATTAACAACATCCAGTTTATCCCAGGAGACTTTACCCGCGCGGTCAATGACAGCGTTAAGCGTTTTGCAGAAACCGCGCCATAGGCTTATTACCTGGTAGGTGAGTAGTTTGAGTTTgcc  >  2:247219/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TCCGTTCAAACATGGCTTCTGCCGACCAACCATCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTAC‑AAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATGTTGCTTTTGCC  >  NC_000913/3969281‑3969576

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