Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,969,436 G→T 5.5% A136S (GCG→TCG)  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,4360GT 5.5% 313.9 / 4.0 119A136S (GCG→TCG) wzzEEntobacterial Common Antigen (ECA) polysaccharide chain length modulation protein
Reads supporting (aligned to +/- strand):  ref base G (64/47);  new base T (4/3);  total (69/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01

Only 100 of 121 total aligned reads displayed.

AACATGGCTTCTGCCGACCAACCATCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATGTTGCTTTTGCCAGCCAGC  >  NC_000913/3969289‑3969583
                                                                                                                                                    |                                                                                                                                                   
aaCATGGCTTCTGCCGACCAACCATCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCGGCAAACCGACTATTACCAACAGCGGATGGTGGGCAACAGCAAGGCCGATgc                                                                                                                                                    >  1:1066270/1‑149 (MQ=255)
      gctgctCCCGACCAACCATCGGTCATGGACGAAGCCTATAAAGAGTTGTTTATGCAACTGGCCTCGTGGGAGACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGCTGGTGGGCAACAGCAAAGCCGATGCGGCGtt                                                                                                                                              <  1:1410565/149‑1 (MQ=255)
        ttCTGCCGACCAACCATCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAACGCCGATGCGGCGTTGc                                                                                                                                            >  1:638046/1‑149 (MQ=255)
        ttCTGCCGACCAACCATCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGc                                                                                                                                            >  2:47434/1‑149 (MQ=255)
           tGCCGACCAACCATCGGTCATGCACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGGGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGGGGGCAACAGCAACGCCGCATCGGCTTTGCTgt                                                                                                                                         <  2:760730/149‑2 (MQ=255)
             ccGACCAACCATCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGAt                                                                                                                                       >  2:55028/1‑149 (MQ=255)
              cGACCAACCATCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCTCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATg                                                                                                                                      >  2:1268787/1‑149 (MQ=255)
               gACCAACCATCGGTCATGGACGAAGCCTATAAAGACTTTCTTATGCAACTGGCCTAGTGGGATACCAGCAGAGAGTTATGGATGCAAACCGACTATTAAAAACAGCGGATGGTGGACAACCGCAAAGCCGATGCGGCGTCGCTGGATGa                                                                                                                                     >  2:396883/1‑149 (MQ=255)
                accaaccaTCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAATCCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGaa                                                                                                                                    >  2:591922/1‑149 (MQ=255)
                accaaccaTCGCGAATGGAGGAAGCCTATAAAGAGTTTGTTTTGCATCTGGCCGCGTGGGATACCCGCAGCGAGTTCTGGCTGCAAACAGACTATTACAAACCGCGGATGGTGGGCAACAGCAACGCCGCATCGGCTTTGCTgttgccc                                                                                                                                    <  2:252966/149‑7 (MQ=255)
                 ccaaccaTCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGaaa                                                                                                                                   >  2:921791/1‑149 (MQ=255)
                 ccaaccaTCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGGTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGaaa                                                                                                                                   <  2:901604/149‑1 (MQ=255)
                   acccaTCGGTCATGGACGCAGCCGATAAAGAGTGTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGCGTTCTGGCTGCAAACCGACTATTACAACCAGCGGATGGAGGGCAACAGCCAAGCCTATGCGGCGTAGCTGGATGAAATg                                                                                                                                 >  1:852873/3‑149 (MQ=255)
                   aaccaTCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACAGACTATTACAAACAGCGGATGGTGCGCAACAGCAAAGCCGATGAGGCGTTGCTCGATCAAATg                                                                                                                                 >  1:465315/1‑149 (MQ=255)
                     ccaTCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGACTCGTGGGAAACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGAAGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGAt                                                                                                                               <  1:232682/149‑1 (MQ=255)
                         cGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTaac                                                                                                                           <  2:787626/149‑1 (MQ=255)
                          ggTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTaaca                                                                                                                          <  1:921791/149‑1 (MQ=255)
                           gTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTaacaa                                                                                                                         >  1:1129743/1‑149 (MQ=255)
                              aTGGACGAAGCCTATAAAGAGGTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACAt                                                                                                                      <  1:1182895/149‑1 (MQ=255)
                               tGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGCTACCCGCAGAGAGTTCTGTCTGCCAGCCGGCTATTTCAAACAGCGGATGGTGGGCAACAGCAAAGAGGAAGCGGCGTTGCTGGATGATATGATTAACAACATc                                                                                                                     >  2:315743/1‑149 (MQ=255)
                                  aCGAAGCGTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTAGAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAg                                                                                                                  <  1:586109/149‑1 (MQ=255)
                                  aCGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAg                                                                                                                  >  1:5996/1‑149 (MQ=255)
                                   cGAAGCCTATAATGAGTTTGTTATGCAACGGGCCTCGTGGGATACCCGCAGCGAGGGCTGGCTGCAAACCGACTATTACAAACAGCGGATGGGGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAgg                                                                                                                 <  2:1391915/149‑2 (MQ=255)
                                   cGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGAGACACGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGGGGGCAACAGCAAAGCCGATCCGGCGGTGCTGGTTGAAATGATTACCAACCTCCAGt                                                                                                                 >  1:1079832/1‑149 (MQ=255)
                                    gAAGCCTAAAAAGAGGGGGCTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAACCAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGtt                                                                                                                <  1:1001243/149‑1 (MQ=255)
                                     taGCCTATAAAGAGTTTGTTATGCAACAGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGttt                                                                                                               <  2:229765/148‑1 (MQ=255)
                                        ccTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTCAAACAGCGGATGGTGGGCAAAAGCAACGCCGATTCGGCGTTGCTCGATGAAATGATTAACAACAACCAGTGTATcc                                                                                                           >  1:1493788/1‑149 (MQ=255)
                                         cTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATcc                                                                                                           >  2:1484608/1‑149 (MQ=255)
                                           ttaAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATcccc                                                                                                         <  2:1050138/148‑1 (MQ=255)
                                           ataAAGAGTTTGTTATGCAACTGGCCTCGTGGGTGACCCGGAGAGCGGTCAGGCTGCAAACCGCCGATTACAAACAGCGGAGGGTGGGCAACAGCAAAGCCGCTGCGGCGTTGCTGGATGAAATGAGTAACAACATCCTGTTTATcccc                                                                                                         >  2:1026870/1‑149 (MQ=255)
                                           ataAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAACGCCGCATCGGCTTTGCt                                                                                                                                           <  1:880769/115‑1 (MQ=255)
                                           ataAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATcccc                                                                                                         >  1:225381/1‑149 (MQ=255)
                                             aaaGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACACCATCCAGTTTATCCCCgg                                                                                                       >  1:1142173/1‑149 (MQ=255)
                                             aaaGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCgg                                                                                                       >  2:1089437/1‑149 (MQ=255)
                                             aaaGAGTTTGTTATGCAACTGGCATCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACAGACTATTACAAAAAGCTGAGGGTGGGCAACAGCAAAGCCTATGCGGCGTTGCTGGATGAAAGGATTAACAACCTCCAGTTTATACCCgg                                                                                                       >  2:113429/1‑149 (MQ=255)
                                              aaGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAATCCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGga                                                                                                      >  1:201198/1‑149 (MQ=255)
                                              aaGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGga                                                                                                      >  2:61641/1‑149 (MQ=255)
                                                 agTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTCGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGAc                                                                                                   >  2:278244/1‑149 (MQ=255)
                                                    ttGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACttt                                                                                                <  1:200559/149‑1 (MQ=255)
                                                     tGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTTTTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACAGCCCGTTTATCCCCGGAGACTTTa                                                                                               >  2:501122/1‑149 (MQ=255)
                                                     tGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTa                                                                                               >  1:1423618/1‑149 (MQ=255)
                                                     cgTTATGCCCCTGGCCTCGTGGGACACCCGCAGCGAGTCCTGTCTGCAAACCGACTATTACACACAGCGGATGGTTGGCCACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTa                                                                                               <  1:1073351/148‑1 (MQ=255)
                                                             aaCTGGCCTCGTGGGATTCCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCACCAGCAAAGCCGAGGCGGCGTTGCTGGATGAAATGATTCACAACGTCCAGTGTATCCCCGGAGACTTTCCCcgcgcg                                                                                       >  2:392957/1‑149 (MQ=255)
                                                                      cGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGAc                                                                              >  1:552000/1‑149 (MQ=255)
                                                                       gTGGGATACCCGCAGAGAGTTCTGGCTGCATACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTATCAAAATCCAGTTTATCCCAGGAGACTTTACCCGCGCGGTCAATGACa                                                                             <  2:1451849/149‑1 (MQ=255)
                                                                          ggATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCg                                                                          >  2:1244935/1‑149 (MQ=255)
                                                                           gATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGt                                                                         >  1:1510446/1‑149 (MQ=255)
                                                                            aTACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTg                                                                        >  2:915005/1‑149 (MQ=255)
                                                                            aTACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTg                                                                        >  2:1206893/1‑149 (MQ=255)
                                                                             tACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGa                                                                       <  1:1390136/149‑1 (MQ=255)
                                                                                ccGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGc                                                                    >  2:514478/1‑149 (MQ=255)
                                                                                    agagagTTCTGGCTGCAAACCGACTTTTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTAt                                                                <  1:224490/149‑1 (MQ=255)
                                                                                    agagagTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTAt                                                                <  2:1510446/149‑1 (MQ=255)
                                                                                        agTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATGCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGcc                                                            <  1:1376757/149‑1 (MQ=255)
                                                                                          ttCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGa                                                          >  1:1238448/1‑149 (MQ=255)
                                                                                            cTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGaaa                                                        <  1:728252/149‑1 (MQ=255)
                                                                                            cTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGaaa                                                        <  2:1182851/149‑1 (MQ=255)
                                                                                             tGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCTGCGTTGCTTGTTGAAATGATTAAATACATCAAGTTTATCCACGGATACTTTACCCGCGCGTTCAATGACAGATTGAAGCTTATTGCAGAAAc                                                       >  1:1017927/1‑149 (MQ=255)
                                                                                               gCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACcg                                                     >  2:1347745/1‑149 (MQ=255)
                                                                                                     aaCCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATACCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTg                                               >  2:1543273/1‑149 (MQ=255)
                                                                                                       ccGACTATTACAAACAGCGGATGGTTGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGGCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGAc                                             <  1:1501094/149‑1 (MQ=255)
                                                                                                       ccGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGAc                                             >  2:487889/1‑149 (MQ=255)
                                                                                                        cGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACg                                            >  2:703563/1‑149 (MQ=255)
                                                                                                         gACTATTACAAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACg                                            >  2:576859/1‑149 (MQ=255)
                                                                                                            tATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCtaa                                        >  2:262902/1‑149 (MQ=255)
                                                                                                             aTTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCtaat                                       <  1:1347745/149‑1 (MQ=255)
                                                                                                              ttACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTTTCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTAATTGCCGAAACCGCGCCTGACGCtaata                                      >  1:855547/1‑149 (MQ=255)
                                                                                                                 cAAACAGCGGATGGGGGGCAACCTCAAAGCCGATGCCGCGTTGCTGGATGTAATGATTAACAACCTCCATGTTATCCCCGGTGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAAcc                                   <  1:1484608/149‑1 (MQ=255)
                                                                                                                   aaaCAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGCCTTTCCCCGCGCGGGCATGGCCAGCGTGAAGCTTATTGCCGAACCCGCGCATGCGGGTAATAcccg                                  >  1:1448412/1‑148 (MQ=255)
                                                                                                                   aaaCAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGaaa                                                        <  2:1141464/127‑1 (MQ=255)
                                                                                                                   aaaCAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGaaa                                                        >  1:1141464/1‑127 (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)
                                                                                                                       aGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTa                              <  1:98447/149‑1 (MQ=255)
                                                                                                                        gCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCATTGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTAc                             <  1:557588/149‑1 (MQ=255)
                                                                                                                          ggATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGGCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGt                           <  2:771385/149‑1 (MQ=255)
                                                                                                                           gATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGCAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTc                          >  1:1339345/1‑149 (MQ=255)
                                                                                                                              ggtggGCAACATCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCAAATAACCTGTTACGTCAGt                       <  1:940735/149‑1 (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)
                                                                                                                                  ggCAACAGCAATGCCGATGCGGCGTGTCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTTCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATGt                   <  2:397541/149‑1 (MQ=255)
                                                                                                                                  ggCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATGt                   <  1:61641/149‑1 (MQ=255)
                                                                                                                                   gCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCa                         <  1:796517/142‑1 (MQ=255)
                                                                                                                                   gCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCa                         >  2:796517/1‑142 (MQ=255)
                                                                                                                                   gCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATGtt                  <  2:201198/149‑1 (MQ=255)
                                                                                                                                   gCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATGtt                  <  1:915005/149‑1 (MQ=255)
                                                                                                                                   gCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATGtt                  >  2:1324212/1‑149 (MQ=255)
                                                                                                                                   gCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATGtt                  <  1:1206893/149‑1 (MQ=255)
                                                                                                                                     tccAGCAACGCCGCATCGGCTTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCGGACGCTAATACCCGGGTACGTCAGTATGTTga                >  1:1529311/3‑148 (MQ=255)
                                                                                                                                       cAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATGTTGCtt              <  2:62851/149‑1 (MQ=255)
                                                                                                                                        aGCAACGCCGCATCGGCTTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATGTTGCttt             >  2:951201/1‑149 (MQ=255)
                                                                                                                                         gCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGAGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATGTTGCtttt            >  1:1417452/1‑149 (MQ=255)
                                                                                                                                          cAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGCCGCTAATAACCTGTTACGTCAGTATGTTGCTTTTg           >  1:1319465/1‑149 (MQ=255)
                                                                                                                                           aaaGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATGTTGCTTTTgc          <  2:5996/149‑1 (MQ=255)
                                                                                                                                            aaGCAGATGCGGCGTGGCTGGATGAAATGATTAACAACATCCAGTTTATCCCAGGAGACTTTACCCGCGCGGTCAATGACAGCGTTAAGCGTTTTGCAGAAACCGCGCCATAGGCTTATTACCTGGTAGGTGAGTAGTTTGAGTTTgcc         >  2:247219/1‑149 (MQ=255)
                                                                                                                                                   tGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATGTTGCTTTTGCcagccagc  <  1:1244935/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AACATGGCTTCTGCCGACCAACCATCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATGTTGCTTTTGCCAGCCAGC  >  NC_000913/3969289‑3969583

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: