Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,122,653 C→T 100% W160* (TGG→TAG)  wcaE ← putative glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,122,6530CT97.9% 930.1 / 6.0 283W160* (TGG→TAG) wcaEputative glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base C (4/2);  new base T (130/147);  total (134/149)
Fisher's exact test for biased strand distribution p-value = 4.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 283 total aligned reads displayed.

ACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCCATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCGCCAAAATCCAGCAGCGCATCGCCGGTGATCATCACGTTA  >  CP009273/2122507‑2122799
                                                                                                                                                  |                                                                                                                                                  
acaATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTAtt                                                                                                                                                 <  1:2196775/150‑1 (MQ=255)
 caATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTAttt                                                                                                                                                >  2:3263378/1‑150 (MQ=255)
  aaTTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTAtttt                                                                                                                                               >  1:3723161/1‑150 (MQ=255)
    ttCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTAtttttt                                                                                                                                             <  1:2144066/150‑1 (MQ=255)
     tCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTc                                                                                                                                            <  1:4091679/150‑1 (MQ=255)
      ccATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCa                                                                                                                                           <  2:1232750/150‑1 (MQ=255)
       cATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCaa                                                                                                                                          <  1:3097014/150‑1 (MQ=255)
        atatTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAg                                                                                                                                         <  1:4524800/150‑1 (MQ=255)
          attattGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGcc                                                                                                                                       >  2:3904230/1‑150 (MQ=255)
            tattGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCgg                                                                                                                                     <  1:363280/150‑1 (MQ=255)
            tattGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCgg                                                                                                                                     <  2:4590508/150‑1 (MQ=255)
               tggtggTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGACGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATa                                                                                                                                  >  1:4617160/1‑150 (MQ=255)
                   ggtAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAgg                                                                                                                              <  1:1700864/150‑1 (MQ=255)
                    gtAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAggg                                                                                                                             <  1:2912801/150‑1 (MQ=255)
                    gtAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAggg                                                                                                                             >  2:4340631/1‑150 (MQ=255)
                     tAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGa                                                                                                                            >  2:1666508/1‑150 (MQ=255)
                     tAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGa                                                                                                                            <  1:2992716/150‑1 (MQ=255)
                      aGATACACCACCCAATGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGAATAACCGGCTTTATACATTTTGGCTGCCAGCTCGTAGTCGGAAGAAACTATATATTCCAGGTCATAACGATATTTTTTAAATCCGGATACAGGGaa                                                                                                                           >  2:232377/1‑150 (MQ=255)
                        aTACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGaaaa                                                                                                                         >  2:4926429/1‑150 (MQ=255)
                          aCCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAat                                                                                                                       <  2:3246662/150‑1 (MQ=255)
                            cccacccaTGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAatat                                                                                                                     >  1:4207786/1‑150 (MQ=255)
                            cccacccaTGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCATTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAatat                                                                                                                     <  1:172408/150‑1 (MQ=255)
                             ccacccaTGGAAAATTCAGACACCATGCCATTGATTTTTTTAAATGCATAACCTGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATc                                                                                                                    <  2:4780580/150‑1 (MQ=255)
                             ccacccaTGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATc                                                                                                                    <  2:3811333/150‑1 (MQ=255)
                                cccaTGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGcc                                                                                                                 >  2:4975083/1‑150 (MQ=255)
                                  caTGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGATTTATACATTTTGGCTTCCTGCGCGTAGTCGTAATAAACTTTATATTCCAGGTCTTTACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCtg                                                                                                               <  2:2680568/150‑1 (MQ=255)
                                   aTGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCtga                                                                                                              >  2:4465729/1‑150 (MQ=255)
                                    tGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCtgat                                                                                                             <  1:4381631/150‑1 (MQ=255)
                                        aaaTTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACt                                                                                                         <  1:3770384/150‑1 (MQ=255)
                                          aTTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTgg                                                                                                       >  2:736263/1‑150 (MQ=255)
                                           ttCAGACACCAGGCAATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGc                                                                                                      <  2:949912/150‑1 (MQ=255)
                                            tCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCg                                                                                                     >  1:342995/1‑150 (MQ=255)
                                             cAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCgg                                                                                                    <  2:5254481/150‑1 (MQ=255)
                                              aGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCggg                                                                                                   >  2:4659034/1‑150 (MQ=255)
                                               gACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTTGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTAAATCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGc                                                                                                  >  1:4459486/1‑150 (MQ=255)
                                               gACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGc                                                                                                  <  2:1041806/150‑1 (MQ=255)
                                                 cacCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAg                                                                                                <  1:2476042/150‑1 (MQ=255)
                                                 cacCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAg                                                                                                >  1:768412/1‑150 (MQ=255)
                                                         cATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGAt                                                                                        >  1:3274111/1‑150 (MQ=255)
                                                         cATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGAt                                                                                        >  1:3876329/1‑150 (MQ=255)
                                                         cATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGAt                                                                                        <  2:1311104/150‑1 (MQ=255)
                                                          aTTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATa                                                                                       >  1:1571822/1‑150 (MQ=255)
                                                           ttGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATaa                                                                                      <  1:2111242/150‑1 (MQ=255)
                                                            tGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATaaa                                                                                     <  1:309739/150‑1 (MQ=255)
                                                             tatttttttAAATGCATAACCTTCTTTATACATTTTGGCTGCCAGCGCTTATTAGTAATAAAATTTATATTCTATTTCATAACGATATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAat                                                                                    <  2:5055846/147‑1 (MQ=255)
                                                             gAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAat                                                                                    <  1:4877384/150‑1 (MQ=255)
                                                                     tAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGc                                                                            >  2:3041958/1‑150 (MQ=255)
                                                                      aaaTGCATAACCTGATTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGcc                                                                           <  2:656719/150‑1 (MQ=255)
                                                                       aaTGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGccc                                                                          >  1:4647521/1‑150 (MQ=255)
                                                                         tGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCgg                                                                        >  2:2520166/1‑150 (MQ=255)
                                                                         tGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCgg                                                                        >  2:4163749/1‑150 (MQ=255)
                                                                            aTAACCGTCTTTATACATTTTGGCTGCCAGCGAGTATTCGTAAGAAACTTTATATTCCATGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGttt                                                                     <  2:3673596/150‑1 (MQ=255)
                                                                            aTAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGttt                                                                     <  2:4798129/150‑1 (MQ=255)
                                                                               aCCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGgc                                                                  >  2:3421946/1‑150 (MQ=255)
                                                                                ccGGCTTTATACATTTTGGCTGCGAGCGCGTAGTCGGAAGAAACTTTATATTACAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGgcg                                                                 <  2:269710/150‑1 (MQ=255)
                                                                                ccGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGgcg                                                                 >  2:1804531/1‑150 (MQ=255)
                                                                                  ggCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCt                                                               <  2:2220510/150‑1 (MQ=255)
                                                                                     tttATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACg                                                            >  2:4044969/1‑150 (MQ=255)
                                                                                      ttATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGt                                                           >  1:3147466/1‑150 (MQ=255)
                                                                                       tataCATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGtt                                                          >  2:473561/1‑150 (MQ=255)
                                                                                        ataCATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGttt                                                         <  1:1123697/150‑1 (MQ=255)
                                                                                            aTTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAAtt                                                     >  1:5137561/1‑150 (MQ=255)
                                                                                            aTTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAAtt                                                     <  1:4709171/150‑1 (MQ=255)
                                                                                             ttttGGCTGCCAGCGCGTAGTCGGAAGAATCTTTATATTCCAGTTCATAACTCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAAttt                                                    <  2:4220570/150‑1 (MQ=255)
                                                                                              tttGGCTGCCAGCGCGTAGTCGTAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAAtttt                                                   <  2:2548710/150‑1 (MQ=255)
                                                                                               ttGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTa                                                  >  2:886669/1‑150 (MQ=255)
                                                                                                tGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTAt                                                 >  2:3513379/1‑150 (MQ=255)
                                                                                                 ggCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATg                                                <  1:2582480/150‑1 (MQ=255)
                                                                                                        aGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCg                                         >  1:4234818/1‑150 (MQ=255)
                                                                                                           gcgTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCGCCa                                      >  2:4419743/1‑150 (MQ=255)
                                                                                                            cgTATTCGGAAGAAACTTTATTTTACAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCGCCaa                                     <  1:1807345/150‑1 (MQ=255)
                                                                                                             gTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCGCCaaa                                    <  1:2789542/150‑1 (MQ=255)
                                                                                                               aGTCGGAATAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCGCCAAAAt                                  <  1:5057121/150‑1 (MQ=255)
                                                                                                               aGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCGCCAAAAt                                  <  1:1308717/150‑1 (MQ=255)
                                                                                                               aGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCGCCAAAAt                                  <  2:3231756/150‑1 (MQ=255)
                                                                                                                 tCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCGCCAAAATcc                                >  1:1596580/1‑150 (MQ=255)
                                                                                                                   ggAAGAAACTTTATATTCCAGTTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCGCCAAAATCcag                              <  2:4415571/150‑1 (MQ=255)
                                                                                                                    gaagaaACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCGCCAAAATCcagc                             <  1:4910577/150‑1 (MQ=255)
                                                                                                                      agaaACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCGCCAAAATCcagcag                           >  2:2764300/1‑150 (MQ=255)
                                                                                                                        aaaCTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCGCCAAAATCCAGCAgcg                         <  1:2641407/150‑1 (MQ=255)
                                                                                                                        aaaCTTTATATTCCAGGTCATAACGCCATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCGCCAAAATCCAGCAgcg                         >  1:2031924/1‑150 (MQ=255)
                                                                                                                            tttATATTACATGTCATAACTCTATTTTTTCAAGCCGGATACATGGAAAAATATCGCCTGATGACTGGCGGGCAGGATGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCTCCAAAATCCAGCAGAGCAta                     <  2:775495/150‑2 (MQ=255)
                                                                                                                              tataTTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCGCCAAAATCCAGCAGCGCATCGc                   >  1:3948281/1‑150 (MQ=255)
                                                                                                                                  ttCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCGCCAAAATCCAGCAGCGCATCGCCGGt               >  1:1419804/1‑150 (MQ=255)
                                                                                                                                  ttCCAGGTCATAACGCCATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTAGGCGCTACGTTTAATTTTATGCCCGTCGCCAAAATCCAGCAGCGCATCGCCGGt               >  1:828256/1‑150 (MQ=255)
                                                                                                                                   tCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCGCCAAATCCAGCAGCGCATCGCCGGTGa             <  2:26093/150‑1 (MQ=255)
                                                                                                                                    ccAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCGCCAAAATCCAGCAGCGCATCGCCGGTGa             <  1:518330/150‑1 (MQ=255)
                                                                                                                                    ccAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCGCCAAAATCCAGCAGCGCATCGCCGGTGa             >  1:4567330/1‑150 (MQ=255)
                                                                                                                                    ccAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCGCCAAAATCCAGCAGCGCATCGCCGGTGa             >  2:1399921/1‑150 (MQ=255)
                                                                                                                                     cAGGTCATAACGCTATTTTTTCAAGCCGGATACAGTGAAAAATATCGCCTGATGACCGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCGCCAAAATCCAGCAGCGCATCGCCGGTGat            <  2:2974890/150‑1 (MQ=255)
                                                                                                                                     cAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCGCCAAAATCCAGCAGCGCATCGCCGGTGat            >  1:2225667/1‑150 (MQ=255)
                                                                                                                                      aGGTCATATCGCCTTTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCTGGCTGTGATTAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCGCCAAAATCCAGCAGCGCATCGCCGGTGatc           <  2:3560560/150‑1 (MQ=255)
                                                                                                                                      aGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCGCCAAAATCCAGCAGCGCATCGCCGGTGatc           <  2:5127958/150‑1 (MQ=255)
                                                                                                                                       ggTCATAACGCCATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCGCCAAAATCCAGCAGCGCATCGCCGGTGatca          >  1:3553628/1‑150 (MQ=255)
                                                                                                                                         tCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCGCCAAAATCCAGCAGCGCATCGCCGGTGatcatc        <  1:61057/150‑1 (MQ=255)
                                                                                                                                           aTAACGCTATTTTTTCCAGCCGGATACAGGGANAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTAGGCCCGTCGCCAAAATCCAGCAGCGCATCGCCGGTGAtcatcan      <  1:600375/150‑2 (MQ=255)
                                                                                                                                           aTAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCGCCAAAATCCAGAAGCGCATCGCCGGTGATCATCAc      <  1:232377/150‑1 (MQ=255)
                                                                                                                                             aaCGCTATTTTTTCAAGCCGGATACAGTGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCGCCAAAATCCAGCAGCGCATCGCCGGTGATCATCACGt    <  2:153034/150‑1 (MQ=255)
                                                                                                                                               cGCTATTTTTTCTAGCCGGATACAGGGAAATATATCGGCTGATGACTGGCGGGCAGGCTGTGATAAATTTACCAGCCCGGTTTGGCGCTACGTTTTATTTTATGCCCGTCGCCAATATCCAGCAGCGCATCGCCGGTGATCATCACGTTa  <  2:1961113/150‑1 (MQ=255)
                                                                                                                                               cGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCGCCAAAATCCAGCAGCGCATCGCCGGTGATCATCACGTTa  <  1:4968790/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
ACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCCATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCGCCAAAATCCAGCAGCGCATCGCCGGTGATCATCACGTTA  >  CP009273/2122507‑2122799

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: