Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,242,509 C→T 93.0% 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,5090CT93.0% 1030.8 / 50.6 316L26F (CTC→TTC) ubiCchorismate‑‑pyruvate lyase
Reads supporting (aligned to +/- strand):  ref base C (15/7);  new base T (140/153);  total (156/160)
Fisher's exact test for biased strand distribution p-value = 7.80e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

Only 100 of 316 total aligned reads displayed.

ATTTGCTTTACATCTCCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGCTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAAGAGTCTCGTT  >  CP009273/4242362‑4242653
                                                                                                                                                   |                                                                                                                                                
 tttGCTTTACATCTCCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCg                                                                                                                                               <  1:1141989/150‑1 (MQ=255)
 tttGCTTTACATCTCCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCg                                                                                                                                               >  2:1840489/1‑150 (MQ=255)
  ttGCTTTACATCTCCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGa                                                                                                                                              >  1:63300/1‑150 (MQ=255)
   tGCTTTACATCTCCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAAGGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGAc                                                                                                                                             <  2:1878085/150‑1 (MQ=255)
     cTTTACATCTCCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTg                                                                                                                                           <  1:1676427/150‑1 (MQ=255)
      tttACATCTCCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTgg                                                                                                                                          <  1:44204/150‑1 (MQ=255)
        tACATCTCCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCt                                                                                                                                        >  2:4555278/1‑150 (MQ=255)
        tACATCTCCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGCTCGACTGGCt                                                                                                                                        >  2:2745560/1‑150 (MQ=255)
               cccGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTg                                                                                                                                 >  2:3602868/1‑150 (MQ=255)
                  gTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTggag                                                                                                                              <  1:2982472/150‑1 (MQ=255)
                   tAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTggagg                                                                                                                             <  1:4035695/150‑1 (MQ=255)
                       cacTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGCTCGACTGGCTGTTGCTGGAGGATTc                                                                                                                         >  2:1299742/1‑150 (MQ=255)
                         cTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCa                                                                                                                       <  2:3430849/150‑1 (MQ=255)
                         cTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCa                                                                                                                       >  1:393842/1‑150 (MQ=255)
                         cTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCa                                                                                                                       >  1:3845262/1‑150 (MQ=255)
                          ttttCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCAt                                                                                                                      <  1:23160/150‑1 (MQ=255)
                           tttCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATg                                                                                                                     <  1:2415106/150‑1 (MQ=255)
                           tttCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATg                                                                                                                     <  2:3205272/150‑1 (MQ=255)
                            ttCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGa                                                                                                                    <  1:2719943/150‑1 (MQ=255)
                             tCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGAc                                                                                                                   <  1:428317/150‑1 (MQ=255)
                                  gATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAAc                                                                                                              <  1:2082942/150‑1 (MQ=255)
                                   aTACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACg                                                                                                             >  1:1277215/1‑150 (MQ=255)
                                       aaTGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGtttt                                                                                                         >  1:3011748/1‑150 (MQ=255)
                                        aTGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTg                                                                                                        <  2:1017230/150‑1 (MQ=255)
                                         tGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGa                                                                                                       >  1:2493460/1‑150 (MQ=255)
                                           ccTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAAc                                                                                                     <  2:4393038/150‑1 (MQ=255)
                                            cTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAAca                                                                                                    <  1:317572/150‑1 (MQ=255)
                                             tttACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAAcag                                                                                                   <  1:3725762/150‑1 (MQ=255)
                                               tACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGACcagca                                                                                                 >  2:3090063/1‑150 (MQ=255)
                                               tACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAAcagca                                                                                                 <  2:4500591/150‑1 (MQ=255)
                                               tACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAAcagca                                                                                                 >  2:334446/1‑150 (MQ=255)
                                                aCGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAAcagcag                                                                                                >  1:4163557/1‑150 (MQ=255)
                                                    tATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGaa                                                                                            <  2:574397/150‑1 (MQ=255)
                                                    tATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGaa                                                                                            >  1:3862916/1‑150 (MQ=255)
                                                    tATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGaa                                                                                            >  2:2766156/1‑150 (MQ=255)
                                                    tATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGCTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGaa                                                                                            >  2:786337/1‑150 (MQ=255)
                                                       gTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGCTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGaaaaa                                                                                         <  2:1137923/150‑1 (MQ=255)
                                                        tAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAAc                                                                                        <  2:1908517/150‑1 (MQ=255)
                                                           cGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGt                                                                                     >  2:1952695/1‑150 (MQ=255)
                                                           cGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGCTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGt                                                                                     >  2:573542/1‑150 (MQ=255)
                                                            ggAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTa                                                                                    >  1:3540347/1‑150 (MQ=255)
                                                            ggAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTa                                                                                    <  1:2010122/150‑1 (MQ=255)
                                                            ggAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTa                                                                                    >  1:839771/1‑150 (MQ=255)
                                                             gagaGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTaa                                                                                   >  2:2414292/1‑150 (MQ=255)
                                                              agagTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAg                                                                                  >  2:1419440/1‑150 (MQ=255)
                                                               gagTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGc                                                                                 <  1:2766156/150‑1 (MQ=255)
                                                                agTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCg                                                                                >  1:472138/1‑150 (MQ=255)
                                                                agTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCg                                                                                >  2:2482000/1‑150 (MQ=255)
                                                                    cGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGAc                                                                            >  1:2260855/1‑150 (MQ=255)
                                                                    cGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGAc                                                                            >  2:1798172/1‑150 (MQ=255)
                                                                     ggCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACg                                                                           >  2:148574/1‑150 (MQ=255)
                                                                      gCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCACCTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACga                                                                          <  2:945562/150‑1 (MQ=255)
                                                                        aTGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACgatg                                                                        >  1:1900486/1‑150 (MQ=255)
                                                                            cacaCCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATcc                                                                    <  2:4592779/150‑1 (MQ=255)
                                                                              cacCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCcgc                                                                  <  2:2942269/150‑1 (MQ=255)
                                                                                ccccGCGTTAACGCACCTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCAGCGa                                                                >  2:2069753/1‑150 (MQ=255)
                                                                                   cgcgTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAgg                                                             >  2:4271795/1‑150 (MQ=255)
                                                                                    gcgTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAggg                                                            <  1:77288/150‑1 (MQ=255)
                                                                                     cgTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGt                                                           <  1:1870848/150‑1 (MQ=255)
                                                                                     cgTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGt                                                           >  2:2140339/1‑150 (MQ=255)
                                                                                      gTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGtt                                                          >  1:4503338/1‑150 (MQ=255)
                                                                                       ttAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGttt                                                         >  1:2493632/1‑150 (MQ=255)
                                                                                         aaCGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGt                                                       >  2:219173/1‑150 (MQ=255)
                                                                                            gCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGa                                                    <  1:344258/150‑1 (MQ=255)
                                                                                             cAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAg                                                   >  2:1208015/1‑150 (MQ=255)
                                                                                             cAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAg                                                   <  2:2754869/150‑1 (MQ=255)
                                                                                                  gcgtgcgCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGaa                                              <  1:1776565/150‑1 (MQ=255)
                                                                                                   cgtgcgCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAAt                                             <  2:4174239/150‑1 (MQ=255)
                                                                                                    gtgcgCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATg                                            <  1:1438128/150‑1 (MQ=255)
                                                                                                       cgctgcgctATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGaaa                                         >  2:2696520/1‑150 (MQ=255)
                                                                                                          tgcgctATTGTTAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTTTTGCTTGAGGATTCCATGACAAAACGTTTTGAACAGCATGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATcc                                      <  1:1785577/150‑1 (MQ=255)
                                                                                                            cgctATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATcccc                                    <  2:2284013/150‑1 (MQ=255)
                                                                                                             gctATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGCTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCg                                   >  1:2829839/1‑150 (MQ=255)
                                                                                                              ctATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCga                                  <  2:3154751/150‑1 (MQ=255)
                                                                                                                aTTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCgaag                                >  1:541055/1‑150 (MQ=255)
                                                                                                                aTTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCgaag                                <  2:1740437/150‑1 (MQ=255)
                                                                                                                  tGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCgaagaa                              >  1:4204611/1‑150 (MQ=255)
                                                                                                                  tGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCgaagaa                              <  1:2107540/150‑1 (MQ=255)
                                                                                                                   gTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAAc                             >  1:691606/1‑150 (MQ=255)
                                                                                                                   gTAAAGAGATCCCTGCCCTGGATCCGCAACTGCTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAAc                             >  1:7506/1‑150 (MQ=255)
                                                                                                                    tAAAGAGATCCCTTCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGATGATTCCATGACAAAACTTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACt                            <  1:2460251/150‑1 (MQ=255)
                                                                                                                    tAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACt                            <  1:3347740/150‑1 (MQ=255)
                                                                                                                       agagATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGcc                         <  2:3733175/150‑1 (MQ=255)
                                                                                                                          gaTCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCgct                      >  1:167926/1‑150 (MQ=255)
                                                                                                                          gaTCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCgct                      >  2:760732/1‑150 (MQ=255)
                                                                                                                          gaTCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCgct                      <  1:3661826/150‑1 (MQ=255)
                                                                                                                          gaTCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTACTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCgct                      <  2:2234513/150‑1 (MQ=255)
                                                                                                                            tCCCTGCCCTGGATCCGCAACTGAACGACTGGATGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAATCGGTAAGCGTGACGATGATCCGCGAAGGGGTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCgctgc                    >  2:4205534/1‑150 (MQ=255)
                                                                                                                                tgccctgGATCCGCAACTGCTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGcc                >  1:1007074/1‑150 (MQ=255)
                                                                                                                                 gccctgGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCg               >  1:896537/1‑150 (MQ=255)
                                                                                                                                   cctgGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGaa             >  2:4009376/1‑150 (MQ=255)
                                                                                                                                    ctgGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGaaa            <  2:38264/150‑1 (MQ=255)
                                                                                                                                     tgGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAag           >  2:1672803/1‑150 (MQ=255)
                                                                                                                                       gATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAagag         <  1:719155/150‑1 (MQ=255)
                                                                                                                                         tCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAAGAGtc       >  2:1444667/1‑150 (MQ=255)
                                                                                                                                           cGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAAGAGtctc     <  2:2856876/150‑1 (MQ=255)
                                                                                                                                           cGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAAGAGtctc     >  1:805496/1‑150 (MQ=255)
                                                                                                                                           cGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAAGAGtctc     <  2:4289971/150‑1 (MQ=255)
                                                                                                                                           cGCAACTGTTCGACGGGCTGTTGCTGGTGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGGAGAATTAAATCCCCGTAGAACTGCCGCTGCTGCCGAAAGAGtctc     <  2:3867280/150‑1 (MQ=255)
                                                                                                                                              aaCTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAAGAGTCTCGtt  <  1:4191057/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
ATTTGCTTTACATCTCCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGCTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAAGAGTCTCGTT  >  CP009273/4242362‑4242653

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: