Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,692,573 C→A 38.5% D289Y (GAT→TAT)  malI ← transcriptional repressor of Mal regulon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,692,5730CA38.5% 358.5 / 290.3 286D289Y (GAT→TAT) malItranscriptional repressor of Mal regulon
Reads supporting (aligned to +/- strand):  ref base C (81/93);  new base A (60/50);  total (141/145)
Fisher's exact test for biased strand distribution p-value = 2.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 291 total aligned reads displayed.

GAGCCGGGCGGGAATAATAAGATTGCGTGAATGCGTCTCTTCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGCAAACCAAACCATGCCCCCATCGCAATAGTTTCGTTATAGCAAACCACGGCACTGATGGTCGGGTTATGA  >  CP009273/1692424‑1692721
                                                                                                                                                     |                                                                                                                                                    
gAGCCGGGCGGGAATAATAAGATTGCGTGAATGCGTCTCTTCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATc                                                                                                                                                      >  1:1279634/1‑150 (MQ=255)
gAGCCGGGCGGGAATAATAAGATTGCGTGAATGCGTCTCTTCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATc                                                                                                                                                      >  1:731109/1‑150 (MQ=255)
gAGCCGGGCGGGAATAATAAGATTGCGTGAATGCGTCTCTTCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATc                                                                                                                                                      <  1:2101912/150‑1 (MQ=255)
 aGCCGGGCGGGAATAATAAGATTGCGTGAATGCGTCTCTTCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCg                                                                                                                                                     <  2:3405251/150‑1 (MQ=255)
  gCCGGGCGGGAATAATAAGATTGCGTGAATGCGTCTCTTCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCgg                                                                                                                                                    <  1:646134/150‑1 (MQ=255)
    cgggcgggAATAATAAGATTGCGTGAATGCGTCTCTTCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTa                                                                                                                                                  <  2:1573901/150‑1 (MQ=255)
    cgggcgggAATAATAAGAGTGCGTGAATGCGTCTCTTCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTGCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCACGTGTGGTTGGTGTCGCATCGGTa                                                                                                                                                  <  2:4066471/150‑1 (MQ=255)
     gggcgggAATAATAAGATTGCGTGAATGCGTCTCTTCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTaa                                                                                                                                                 >  1:1824053/1‑150 (MQ=255)
       gcgggAATAATAAGATTGCGTGAATGCGTCTCTTCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATAGGTAAAt                                                                                                                                               <  1:3283183/150‑1 (MQ=255)
         gggAATAATAAGATTGCGTGAATGCGTCTCTTCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATAGGTAAATgc                                                                                                                                             <  2:3364088/150‑1 (MQ=255)
            aataatAAGATTGCGTGAATGCGTCTCTTCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATAGGTAAATgccgc                                                                                                                                          >  2:491310/1‑150 (MQ=255)
             ataataAGATTGCGTGAATGCGTCTCTTCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATgccgcc                                                                                                                                         <  1:5126274/150‑1 (MQ=255)
             ataataAGATTGCGTGAATGCGTCTCTTCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATgccgcc                                                                                                                                         <  1:3190398/150‑1 (MQ=255)
             ataataAGATTGCGTGAATGCGTCTCTTCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATgccgcc                                                                                                                                         <  2:2055082/150‑1 (MQ=255)
             ataataAGATTGCGTGAATGCGTCTCTTCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATgccgcc                                                                                                                                         <  2:2571246/150‑1 (MQ=255)
              taataaGATTGCGTGAATGCGTCTCTTCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCa                                                                                                                                        >  1:498272/1‑150 (MQ=255)
              taataaGATTGCGTGAATGCGTCTCTTCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATAGGTAAATGCCGCCa                                                                                                                                        <  2:1789731/150‑1 (MQ=255)
                ataaGATTGCGTGAATGCGTCTCTTCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGc                                                                                                                                      >  1:3565363/1‑150 (MQ=255)
                ataaGATTGCGTGAATGCGTCTCTTCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATAGGTAAATGCCGCCAGc                                                                                                                                      >  2:3060635/1‑150 (MQ=255)
                 taaGATTGCGTGAATGCGTCTCTTCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCg                                                                                                                                     >  1:388001/1‑150 (MQ=255)
                  aaGATTGCGTGAATGCGTCTCTTCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATAGGTAAATGCCGCCAGCGa                                                                                                                                    >  1:632801/1‑150 (MQ=255)
                   aGATTGCGTGAATGCGTCTCTTCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATAGGTAAATGCCGCCAGCGaa                                                                                                                                   >  1:4618139/1‑150 (MQ=255)
                         cGTGAATGCGTCTCTTCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGt                                                                                                                             >  1:3245254/1‑150 (MQ=255)
                            gAATGCGTCTCTTCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACAAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGTTTGGTGTCTCATCGGTAAATGACGCCAGCGAAACCTGTTGc                                                                                                                          >  1:4893581/1‑150 (MQ=255)
                               tGCGTCTCTTCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATAGGTAAATGCCGCCAGCGAAACCTGTTGCTCa                                                                                                                       <  1:1938210/150‑1 (MQ=255)
                                 cGTCTCTTCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATAGGTAAATGCCGCCAGCGAAACCTGTTGCTCaaa                                                                                                                     >  2:3802358/1‑150 (MQ=255)
                                    ctctTCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATAGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTa                                                                                                                  <  2:1164186/150‑1 (MQ=255)
                                      ctTCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAAc                                                                                                                >  2:1936056/1‑150 (MQ=255)
                                       ttCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACg                                                                                                               >  1:4369346/1‑150 (MQ=255)
                                        tCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGa                                                                                                              <  2:315956/150‑1 (MQ=255)
                                        tCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATAGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGa                                                                                                              >  2:4043998/1‑150 (MQ=255)
                                         cATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGAt                                                                                                             >  2:4535666/1‑150 (MQ=255)
                                         cATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGAt                                                                                                             >  1:992544/1‑150 (MQ=255)
                                                gATTTTTTGCAACATGCGATCCGCAAGTGTGATACCAATTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTTTGTTTGGTTTCGCATCGTTAAATGCCGCCAGCGAAACCTGTTGCTCATAGTAACGATCGACTcc                                                                                                      <  2:231861/150‑1 (MQ=255)
                                                 aTTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATAGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCg                                                                                                     >  1:2786305/1‑150 (MQ=255)
                                                 aTTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATAGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCg                                                                                                     <  2:4702003/150‑1 (MQ=255)
                                                    ttttGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCtt                                                                                                  <  2:1043584/150‑1 (MQ=255)
                                                    ttttGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCtt                                                                                                  >  1:2881955/1‑150 (MQ=255)
                                                       tGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTcc                                                                                               >  2:4360346/1‑150 (MQ=255)
                                                        gCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATAGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTccc                                                                                              <  2:52182/150‑1 (MQ=255)
                                                            catGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCt                                                                                          >  2:76920/1‑150 (MQ=255)
                                                             atGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCtt                                                                                         >  1:4703849/1‑150 (MQ=255)
                                                                cGATCCGCAAGTGTGAGACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATAGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTgc                                                                                      <  2:4907456/150‑1 (MQ=255)
                                                                 gATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTgcc                                                                                     >  2:4743360/1‑150 (MQ=255)
                                                                 gATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATAGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTgcc                                                                                     >  1:1409846/1‑150 (MQ=255)
                                                                   tCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTgccgc                                                                                   >  2:1994594/1‑150 (MQ=255)
                                                                    ccGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTgccgcc                                                                                  >  2:4585863/1‑150 (MQ=255)
                                                                           tgtgATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCttt                                                                           <  2:2237961/150‑1 (MQ=255)
                                                                           tgtgATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATAGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCttt                                                                           >  1:2681386/1‑150 (MQ=255)
                                                                           tgtgATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATAGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCttt                                                                           >  1:896507/1‑150 (MQ=255)
                                                                             tgATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGa                                                                         <  1:751492/150‑1 (MQ=255)
                                                                              gATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAg                                                                        <  2:270711/150‑1 (MQ=255)
                                                                                tACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGca                                                                      <  2:3594376/150‑1 (MQ=255)
                                                                                   cAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATAGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGcaaac                                                                   >  1:1000330/1‑150 (MQ=255)
                                                                                    aaGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGcaaacc                                                                  <  2:4851574/150‑1 (MQ=255)
                                                                                     aGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGcaaacca                                                                 >  1:418010/1‑150 (MQ=255)
                                                                                     aGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATAGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGcaaacca                                                                 <  2:383360/150‑1 (MQ=255)
                                                                                      gTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGGGGTTGGTGTCGAATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACCATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGcaaaccaa                                                                <  1:3047124/150‑1 (MQ=255)
                                                                                       ttCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGcaaaccaaa                                                               >  1:4003468/1‑150 (MQ=255)
                                                                                        tCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGcaaaccaaac                                                              >  2:1661163/1‑150 (MQ=255)
                                                                                           cgcgCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGCAAACCAAACCAt                                                           <  2:3245254/150‑1 (MQ=255)
                                                                                           cgcgCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGCAAACCAAACCAt                                                           <  2:2499481/150‑1 (MQ=255)
                                                                                               cTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGCAAACCAAACCATGccc                                                       <  2:3718972/150‑1 (MQ=255)
                                                                                               cTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGCAAACCAAACCATGccc                                                       >  1:2381424/1‑150 (MQ=255)
                                                                                                tGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTTAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGCAAACCAAACCATGcccc                                                      <  1:4284664/150‑1 (MQ=255)
                                                                                                 ggCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGCAAACCAAACCATGccccc                                                     <  2:3511912/150‑1 (MQ=255)
                                                                                                   cGTGCTGGCCCAGGTAACAGGTATTTCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCGGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGTCTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGCAAACCAAACCATGCCCCCAt                                                   <  2:288051/150‑1 (MQ=255)
                                                                                                   cGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGCAAACCAAACCATGCCCCCAt                                                   >  1:1798382/1‑150 (MQ=255)
                                                                                                         ggCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGCAAACCAAACCATGCCCCCATCGCAAt                                             <  1:1582618/150‑1 (MQ=255)
                                                                                                             cAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGCAAACCAAACCATGCCCCCATCGCAATAGtt                                         <  1:1704686/150‑1 (MQ=255)
                                                                                                               ggTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGCAAACCAAACCATGCCCCCATCGCAATAGTTTc                                       <  2:4413944/150‑1 (MQ=255)
                                                                                                               ggTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATAGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGCAAACCAAACCATGCCCCCATCGCAATAGTTTc                                       >  2:1411210/1‑150 (MQ=255)
                                                                                                                gTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGCAAACCAAACCATGCCCCCATCGCAATAGTTTCg                                      >  1:1211298/1‑150 (MQ=255)
                                                                                                                 tAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGCAAACCAAACCATGCCCCCATCGCAATAGTTTCGt                                     >  2:2134320/1‑150 (MQ=255)
                                                                                                                 tAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATAGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGCAAACCAAACCATGCCCCCATCGCAATAGTTTCGt                                     <  2:1351578/150‑1 (MQ=255)
                                                                                                                 tAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATAGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGCAAACCAAACCATGCCCCCATCGCAATAGTTTCGt                                     <  2:1351545/150‑1 (MQ=255)
                                                                                                                   aCAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGCAAACCAAACCATGCCCCCATCGCAATAGTTTCGTta                                   >  1:2143280/1‑150 (MQ=255)
                                                                                                                   aCAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGCAAACCAAACCATGCCCCCATCGCAATAGTTTCGTta                                   <  2:2306321/150‑1 (MQ=255)
                                                                                                                    cAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGCAAACCAAACCATGCCCCCATCGCAATAGTTTCGTtat                                  >  2:4250490/1‑150 (MQ=255)
                                                                                                                     aGGTATATCATCAAGTGTGGTTGGTGTCGCATAGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGCAAACCAAACCATGCCCCCATCGCAATAGTTTCGTtata                                 >  2:2409168/1‑150 (MQ=255)
                                                                                                                         atatCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGCAAACCAAACCATGCCCCCATCGCAATAGTTTCGTTATAGCaa                             <  2:3661592/150‑1 (MQ=255)
                                                                                                                         atatCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGCAAACCAAACCATGCCCCCATCGCAATAGTTTCGTTATAGCaa                             <  1:1884840/150‑1 (MQ=255)
                                                                                                                         atatCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGCAAACCAAACCATGCCCCCATCGCAATAGTTTCGTTATAGCaa                             <  1:3333497/150‑1 (MQ=255)
                                                                                                                         atatCATCAAGTGTGGTTGGTGTCGCATAGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGCAAACCAAACCATGCCCCCATCGCAATAGTTTCGTTATAGCaa                             >  2:2673309/1‑150 (MQ=255)
                                                                                                                          tatCATCAAGTGTGGTTGGTGTCGCATAGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGCAAACCAAACCATGCCCCCATCGCAATAGTTTCGTTATAGCaaa                            <  1:1080222/150‑1 (MQ=255)
                                                                                                                           atcatcAAGTGTGGTTGGTGTCGCATCGTTAAATGCGGCCAGCGAAACCTGTTGCTCAAATTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGCAAACCAAACCATGCCCCCATCGCAATAGTTTCGTTATAGCAAAc                           <  1:22565/150‑1 (MQ=255)
                                                                                                                            tcatcaAGTGTGGTTGGTGTCGCATAGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGCAAACCAAACCATGCCCCCATCGCAATAGTTTCGTTATAGCAAAcc                          >  2:5068667/1‑150 (MQ=255)
                                                                                                                              atcaAGTGTGGTTGGTGTCGCATAGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGCAAACCAAACCATGCCCCCATCGCAATAGTTTCGTTATAGCAAACCAc                        >  1:4897161/1‑150 (MQ=255)
                                                                                                                              atcaAGTGTGGTTGGTGTCGCATAGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGCAAACCAAACCATGCCCCCATCGCAATAGTTTCGTTATAGCAAACCAc                        >  1:3992735/1‑150 (MQ=255)
                                                                                                                                    tgtgGTTGGTGTCGCATAGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGCAAACCAAACCATGCCCCCATCGCAATAGTTTCGTTATAGCAAACCACGGCACt                  >  1:534325/1‑150 (MQ=255)
                                                                                                                                         ttggtGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGCAAACCAAACCATGCCCCCATCGCAATAGTTTCGTTATAGCAAACCACGGCACTGATgg             >  2:4442714/1‑150 (MQ=255)
                                                                                                                                         ttggtGTCGCATAGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGCAAACCAAACCATGCCCCCATCGCAATAGTTTCGTTATAGCAAACCACGGCACTGATgg             <  2:889485/150‑1 (MQ=255)
                                                                                                                                            gtgtCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAACGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGCAAACCAAACCATGCCCCCATCGCAATAGTTTCGTTATAGCAAACCACGGCACTGATGGTCg          <  2:2456536/150‑1 (MQ=255)
                                                                                                                                             tgtCGCATAGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGCAAACCAAACCATGCCCCCATCGCAATAGTTTCGTTATAGCAAACCACGGCACTGATGGTCgg         <  1:5084101/150‑1 (MQ=255)
                                                                                                                                              gtCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGCAAACCAAACCATGCCCCCATCGCAATAGTTTCGTTATAGCAAACCACGGCACTGATGGTCggg        >  1:3963703/1‑150 (MQ=255)
                                                                                                                                                cGCATAGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGCAAACCAAACCATGCCCCCATCGCAATAGTTTCGTTATAGCAAACCACGGCACTGATGGTCGGGtt      >  2:2789182/1‑150 (MQ=255)
                                                                                                                                                 gCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGCAAACCAAACCATGCCCCCATCGCAATAGTTTCGTTATAGCAAACCACGGCACTGATGGTCGGGTTa     >  2:4369696/1‑150 (MQ=255)
                                                                                                                                                  cATAGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGCAAACCAAACCATGCCCCCATCGCAATAGTTTCGTTATAGCAAACCACGGCACTGATGGTCGGGTTAt    >  2:226330/1‑150 (MQ=255)
                                                                                                                                                  cATAGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGCAAACCAAACCATGCCCCCATCGCAATAGTTTCGTTATAGCAAACCACGGCACTGATGGTCGGGTTAt    <  2:557307/150‑1 (MQ=255)
                                                                                                                                                   aTCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGCAAACCAAACCATGCCCCCATCGCAATAGTTTCGTTATAGCAAACCACGGCACTGATGGTCGGGTTATg   >  1:3751552/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GAGCCGGGCGGGAATAATAAGATTGCGTGAATGCGTCTCTTCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGCAAACCAAACCATGCCCCCATCGCAATAGTTTCGTTATAGCAAACCACGGCACTGATGGTCGGGTTATGA  >  CP009273/1692424‑1692721

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: