Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,242,509 C→T 100% L26F (CTC→TTC)  ubiC → chorismate‑‑pyruvate lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,242,5090CT98.2% 936.9 / 3.6 282L26F (CTC→TTC) ubiCchorismate‑‑pyruvate lyase
Reads supporting (aligned to +/- strand):  ref base C (4/1);  new base T (135/141);  total (140/142)
Fisher's exact test for biased strand distribution p-value = 2.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 282 total aligned reads displayed.

ATTTGCTTTACATCTCCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGCTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAAGAGTCTCGTTA  >  CP009273/4242362‑4242654
                                                                                                                                                   |                                                                                                                                                 
 tttGCTTTACATCTCCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCg                                                                                                                                                >  2:4512189/1‑150 (MQ=255)
 tttGCTTTACATCTCCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCg                                                                                                                                                >  1:4597871/1‑150 (MQ=255)
     cTTTACATCTCCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTg                                                                                                                                            >  2:3169578/1‑150 (MQ=255)
         aCATCTCCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTg                                                                                                                                        >  1:838603/1‑150 (MQ=255)
         aCATCTCCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTg                                                                                                                                        >  2:1290114/1‑150 (MQ=255)
         aCATCTCCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTg                                                                                                                                        >  1:2606559/1‑150 (MQ=255)
          cATCTCCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGt                                                                                                                                       >  2:4883714/1‑150 (MQ=255)
           aTCTCCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGtt                                                                                                                                      <  1:2704908/150‑1 (MQ=255)
             ctcCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGc                                                                                                                                    >  2:1951263/1‑150 (MQ=255)
               cccGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTg                                                                                                                                  <  2:92375/150‑1 (MQ=255)
               cccGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTg                                                                                                                                  <  2:3038728/150‑1 (MQ=255)
               cccGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTg                                                                                                                                  <  2:3131488/150‑1 (MQ=255)
                ccGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTgg                                                                                                                                 <  2:2238027/150‑1 (MQ=255)
                 cGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTgga                                                                                                                                >  1:2933099/1‑150 (MQ=255)
                     aaCACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGAt                                                                                                                            <  1:3789727/150‑1 (MQ=255)
                        acTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTcc                                                                                                                         <  1:1677579/150‑1 (MQ=255)
                         cTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCa                                                                                                                        <  2:740472/150‑1 (MQ=255)
                          ttttCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCAt                                                                                                                       <  2:2075929/150‑1 (MQ=255)
                          ttttCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCAt                                                                                                                       >  1:4736625/1‑150 (MQ=255)
                           tttCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGTCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATg                                                                                                                      >  2:2607270/1‑150 (MQ=255)
                           tttCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATg                                                                                                                      <  2:4217430/150‑1 (MQ=255)
                             tCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGAc                                                                                                                    <  1:228415/150‑1 (MQ=255)
                                  gATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAAc                                                                                                               >  2:1277883/1‑150 (MQ=255)
                                       aaTGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGtttt                                                                                                          <  1:2191006/150‑1 (MQ=255)
                                       aaTGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGtttt                                                                                                          <  1:4393126/150‑1 (MQ=255)
                                        aTGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTg                                                                                                         >  2:1015468/1‑150 (MQ=255)
                                         tGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGa                                                                                                        <  1:4859787/150‑1 (MQ=255)
                                          gCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGaa                                                                                                       >  2:1364152/1‑150 (MQ=255)
                                          gCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGaa                                                                                                       >  1:1088413/1‑150 (MQ=255)
                                             tttACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAAcag                                                                                                    <  2:4283578/150‑1 (MQ=255)
                                                  gTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAggg                                                                                               <  2:1053808/150‑1 (MQ=255)
                                                   ttATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGa                                                                                              <  1:1277883/150‑1 (MQ=255)
                                                   ttATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGa                                                                                              >  2:574450/1‑150 (MQ=255)
                                                    tATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGaa                                                                                             <  2:2452435/150‑1 (MQ=255)
                                                      tGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGaaaa                                                                                           >  1:302276/1‑150 (MQ=255)
                                                        tAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAAc                                                                                         >  1:1482089/1‑150 (MQ=255)
                                                        tAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAAc                                                                                         <  2:3536512/150‑1 (MQ=255)
                                                         aaCGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTTTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACg                                                                                        <  1:1390918/150‑1 (MQ=255)
                                                          aCGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACgg                                                                                       <  1:72333/150‑1 (MQ=255)
                                                             gagaGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTaa                                                                                    <  2:602539/150‑1 (MQ=255)
                                                                agTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCg                                                                                 >  1:5112602/1‑150 (MQ=255)
                                                                agTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCg                                                                                 <  2:2018641/150‑1 (MQ=255)
                                                                 gTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGt                                                                                >  1:4934439/1‑150 (MQ=255)
                                                                  ttCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTg                                                                               >  1:5150653/1‑150 (MQ=255)
                                                                    cGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTCAGCGTGAc                                                                             >  2:66851/1‑150 (MQ=255)
                                                                    cGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGAc                                                                             <  1:4494012/150‑1 (MQ=255)
                                                                     ggCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACg                                                                            >  1:2328936/1‑150 (MQ=255)
                                                                     ggCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACg                                                                            <  1:1709177/150‑1 (MQ=255)
                                                                       cATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACgat                                                                          >  1:179597/1‑150 (MQ=255)
                                                                          gTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGAGCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACgatgat                                                                       <  1:1914679/150‑1 (MQ=255)
                                                                           tCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATc                                                                      <  2:1681228/150‑1 (MQ=255)
                                                                           tCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATc                                                                      >  2:1137504/1‑150 (MQ=255)
                                                                            cacaCCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATcc                                                                     <  1:2008467/150‑1 (MQ=255)
                                                                               acCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCcgcg                                                                  >  2:3811118/1‑150 (MQ=255)
                                                                                ccccGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGa                                                                 >  2:4116339/1‑150 (MQ=255)
                                                                                  ccGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAg                                                               <  2:1435827/150‑1 (MQ=255)
                                                                                   cgcgTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAgg                                                              <  2:2752032/150‑1 (MQ=255)
                                                                                     cgTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGt                                                            >  1:1047999/1‑150 (MQ=255)
                                                                                            gCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGTTGATCCGCGAAGGGTTTGTCGa                                                     <  2:1517586/150‑1 (MQ=255)
                                                                                            gCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGa                                                     <  2:505456/150‑1 (MQ=255)
                                                                                              taCTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGc                                                   >  2:3345052/2‑150 (MQ=255)
                                                                                              aaCTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGc                                                   >  1:1151596/1‑150 (MQ=255)
                                                                                                 tgcgtgcgCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGa                                                >  1:429404/1‑150 (MQ=255)
                                                                                                 tgcgtgcgCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGa                                                >  2:2548047/1‑150 (MQ=255)
                                                                                                  gcgtgcgCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGaa                                               >  1:1536786/1‑150 (MQ=255)
                                                                                                    gtgcgCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTTTTCGACTGGCTGTTGCTGGAGGATTCCTTGACAAAACGTTTTGAACAGCAGGGAAAAACTGTAAGCGTGACGATGATCCGCGAAGGTTTTGTCGAGCAGATTg                                             <  1:4125826/150‑1 (MQ=255)
                                                                                                     tgcgctgcgcTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGa                                            <  2:3980680/150‑1 (MQ=255)
                                                                                                      gcgctgcgctATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGaa                                           <  1:4354135/150‑1 (MQ=255)
                                                                                                       cgctgcgctATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGaaa                                          >  1:3946929/1‑150 (MQ=255)
                                                                                                       cgctgcgctATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGaaa                                          >  2:4141636/1‑150 (MQ=255)
                                                                                                       cgctgcgctATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGaaa                                          >  2:4806530/1‑150 (MQ=255)
                                                                                                       cgctgcgctATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGaaa                                          >  2:392456/1‑150 (MQ=255)
                                                                                                        gctgcgctATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAAt                                         >  2:2666497/1‑150 (MQ=255)
                                                                                                        gctgcgctATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAAt                                         >  2:2913282/1‑150 (MQ=255)
                                                                                                                aTTGTAAAGAGATCCCTGCCCTGGATCCGCAACAGTTCGACTGGCTCTAGCTGTAGGATTCCATGACAAAACGTTTTGAACAGCAGTGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGATTGTCGAGAAGAATGAAATCCCCgaag                                 >  1:2753191/1‑150 (MQ=255)
                                                                                                                     aaaGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTg                            <  1:4358217/150‑1 (MQ=255)
                                                                                                                         agaTCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCgc                        >  2:2469495/1‑150 (MQ=255)
                                                                                                                           aTCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCgctg                      <  2:368671/150‑1 (MQ=255)
                                                                                                                               ctgccctgGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCtgctgc                  <  1:1598655/150‑1 (MQ=255)
                                                                                                                                tgccctgGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGcc                 >  2:4353405/1‑150 (MQ=255)
                                                                                                                                 gccctgGATCCGCTACTGTTCGACTGGCTGTTGCTGGAGTATTCCATGAAAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCTCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCg                <  1:3661228/150‑1 (MQ=255)
                                                                                                                                   cctgGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGac              >  2:4861873/1‑149 (MQ=255)
                                                                                                                                    ctgGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAATCGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGaaa             >  2:93204/1‑150 (MQ=255)
                                                                                                                                     tgGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAag            >  1:862144/1‑150 (MQ=255)
                                                                                                                                      ggATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAaga           <  2:4451998/150‑1 (MQ=255)
                                                                                                                                      ggATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAaga           <  2:1881709/150‑1 (MQ=255)
                                                                                                                                       gATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAagag          <  2:179597/150‑1 (MQ=255)
                                                                                                                                        aTCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAAGAGt         >  1:3200680/1‑150 (MQ=255)
                                                                                                                                         tCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAAGAGtc        <  1:4187581/150‑1 (MQ=255)
                                                                                                                                         tCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAAGAGtc        <  1:512212/150‑1 (MQ=255)
                                                                                                                                           cGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAAGAGtctc      <  2:4249989/150‑1 (MQ=255)
                                                                                                                                            gCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAAGAGTCTCg     <  1:5028307/150‑1 (MQ=255)
                                                                                                                                            gCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAAGAGTCTCg     <  1:3304993/150‑1 (MQ=255)
                                                                                                                                              aaCTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAAGAGTCTCGtt   <  2:3455717/150‑1 (MQ=255)
                                                                                                                                              aaCTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAAGAGTCTCGtt   <  1:3345052/150‑1 (MQ=255)
                                                                                                                                              aaCTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAAGAGTCTCGtt   <  2:1159102/150‑1 (MQ=255)
                                                                                                                                              aaCTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAAGAGTCTCGtt   >  1:2855619/1‑150 (MQ=255)
                                                                                                                                              aaCTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAAGAGTCTCGtt   >  2:2255836/1‑150 (MQ=255)
                                                                                                                                               aCTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAAGAGTCTCGTTa  >  2:4597199/1‑150 (MQ=255)
                                                                                                                                               aCTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAAGAGTCTCGTTa  >  2:853847/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
ATTTGCTTTACATCTCCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGCTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAAGAGTCTCGTTA  >  CP009273/4242362‑4242654

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