Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 463,179 C→A 100% I104I (ATC→ATA)  cof → thiamine pyrimidine pyrophosphate hydrolase; HMP‑PP phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273463,1790CA99.8% 1009.8 / ‑6.2 312I104I (ATC→ATA) cofthiamine pyrimidine pyrophosphate hydrolase; HMP‑PP phosphatase
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base A (144/165);  total (145/167)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
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 316 total aligned reads displayed.

GCTATCGCTGGATGCGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGA  >  CP009273/463032‑463325
                                                                                                                                                   |                                                                                                                                                  
gcTATCGCTGGATGCGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATAtt                                                                                                                                                  >  2:4171243/1‑150 (MQ=255)
gcTATCGCTGGATGCGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATAtt                                                                                                                                                  >  1:4065423/1‑150 (MQ=255)
 cTATCGCTGGATGCGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTc                                                                                                                                                 >  1:3581546/1‑150 (MQ=255)
    tCGCTGGATGCGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAAt                                                                                                                                              <  1:1000416/150‑1 (MQ=255)
    tCGCTGGATGCGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAAt                                                                                                                                              <  1:4636432/150‑1 (MQ=255)
    tCGCTGGATGCGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAAt                                                                                                                                              <  1:427102/150‑1 (MQ=255)
     cGCTGGATGCGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATg                                                                                                                                             <  2:88858/150‑1 (MQ=255)
      gCTGGATGCGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATga                                                                                                                                            >  2:2005927/1‑150 (MQ=255)
       cTGGATGCGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATgac                                                                                                                                           <  1:1575088/150‑1 (MQ=255)
       cTGGATGCGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATgac                                                                                                                                           <  2:2911589/150‑1 (MQ=255)
         ggATGCGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATgacga                                                                                                                                         <  2:3599898/150‑1 (MQ=255)
         ggATGCGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATgacga                                                                                                                                         <  2:5006505/150‑1 (MQ=255)
          gATGCGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATgacgac                                                                                                                                        >  2:4587802/1‑150 (MQ=255)
          gATGAGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATgacgac                                                                                                                                        <  1:3703013/150‑1 (MQ=255)
              cGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACggtt                                                                                                                                    <  2:2705501/150‑1 (MQ=255)
               ttATTTGTTTTCCGGCAACGGAACGCGCGTGAATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTGACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACggttg                                                                                                                                   <  2:4979553/149‑1 (MQ=255)
               gTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACggttg                                                                                                                                   >  2:3177768/1‑150 (MQ=255)
               gTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACggttg                                                                                                                                   >  2:4442557/1‑150 (MQ=255)
               gTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACggttg                                                                                                                                   >  2:2286429/1‑150 (MQ=255)
                tATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACggttgg                                                                                                                                  <  1:174781/150‑1 (MQ=255)
                  tttGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACggttggtt                                                                                                                                >  2:3272394/1‑150 (MQ=255)
                     gATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTAc                                                                                                                             >  1:1306587/1‑150 (MQ=255)
                      aTTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTAcc                                                                                                                            <  1:801301/150‑1 (MQ=255)
                      aTTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTAcc                                                                                                                            <  1:2508807/150‑1 (MQ=255)
                         aCCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCggg                                                                                                                         <  1:67761/150‑1 (MQ=255)
                         aCCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTNGGTTTACCggg                                                                                                                         <  1:2757381/150‑1 (MQ=255)
                          ccGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGa                                                                                                                        <  2:1426195/150‑1 (MQ=255)
                           cGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGaa                                                                                                                       >  2:2247184/1‑150 (MQ=255)
                                  ggAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATc                                                                                                                <  2:1325343/150‑1 (MQ=255)
                                  ggAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATc                                                                                                                >  2:838831/1‑150 (MQ=255)
                                   gAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATcc                                                                                                               <  2:40332/150‑1 (MQ=255)
                                   gAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGAGGAGTTACCTGCGGAGGTCGCGGAGCTGGTGCTGTAGCAGCAATGGGATACCCGAGCCAGCATGCGTCTCTTCATGGAGGGCGGTGGGTTTTCCGGGAAAGAGTTcc                                                                                                               >  2:4967316/1‑150 (MQ=255)
                                       gcgcgTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGc                                                                                                           <  2:158678/150‑1 (MQ=255)
                                        cgcgTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCg                                                                                                          <  2:522775/150‑1 (MQ=255)
                                          cgTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCgtt                                                                                                        <  2:4291321/150‑1 (MQ=255)
                                          cgTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCgtt                                                                                                        <  2:1647962/150‑1 (MQ=255)
                                           gTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCgttg                                                                                                       >  2:680201/1‑150 (MQ=255)
                                             gCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCgttgtt                                                                                                     >  2:4397249/1‑150 (MQ=255)
                                              cATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGttgttg                                                                                                    >  1:3356752/1‑150 (MQ=255)
                                               aTTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGc                                                                                                   <  1:1076485/150‑1 (MQ=255)
                                               aTTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGc                                                                                                   <  1:2139548/150‑1 (MQ=255)
                                               aTTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGc                                                                                                   >  1:967156/1‑150 (MQ=255)
                                                  ctctGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCCGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAgg                                                                                                >  1:497248/1‑150 (MQ=255)
                                                      ggAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCAtt                                                                                            >  1:4954202/1‑150 (MQ=255)
                                                         aGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCTATCCAGCATGCATATATTCAATGACTACTGTTGGTTTACCGGGAAAGAGATCCCTGCGTAGTTGCAGGCATTTGt                                                                                         >  2:539023/1‑150 (MQ=255)
                                                         aGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGt                                                                                         <  1:1702363/150‑1 (MQ=255)
                                                         aGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGt                                                                                         >  2:4858769/1‑150 (MQ=255)
                                                         aGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGt                                                                                         <  1:3169822/150‑1 (MQ=255)
                                                          ggTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTc                                                                                        >  1:4926528/1‑150 (MQ=255)
                                                          ggTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTc                                                                                        >  1:1790619/1‑150 (MQ=255)
                                                           gTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCt                                                                                       <  2:1118222/150‑1 (MQ=255)
                                                            tGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCta                                                                                      <  1:3075191/150‑1 (MQ=255)
                                                             gAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCtat                                                                                     >  2:1497864/1‑150 (MQ=255)
                                                               aCTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAg                                                                                   >  1:1906022/1‑150 (MQ=255)
                                                                cTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGc                                                                                  >  2:5099655/1‑150 (MQ=255)
                                                                 ttttACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCg                                                                                 >  1:3790990/1‑150 (MQ=255)
                                                                     aCATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGttt                                                                             >  1:471668/1‑150 (MQ=255)
                                                                       aTCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTc                                                                           <  1:4194723/150‑1 (MQ=255)
                                                                         cGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGt                                                                         <  1:3192437/150‑1 (MQ=255)
                                                                          gTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGtt                                                                        <  2:4152347/150‑1 (MQ=255)
                                                                           tgatgaTTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTa                                                                       <  1:3471749/150‑1 (MQ=255)
                                                                            gatgatTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTAt                                                                      <  1:4016807/150‑1 (MQ=255)
                                                                            gatgatTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTAt                                                                      >  1:2192068/1‑150 (MQ=255)
                                                                             atgatTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATc                                                                     <  1:1490679/150‑1 (MQ=255)
                                                                              tgatTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCa                                                                    >  1:1028506/1‑150 (MQ=255)
                                                                                 tttACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGAt                                                                 <  1:1205309/150‑1 (MQ=255)
                                                                                   tACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGACAGAGATCACTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATaa                                                               >  1:3328640/1‑150 (MQ=255)
                                                                                   tACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATaa                                                               >  2:3492281/1‑150 (MQ=255)
                                                                                     ccTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCCGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTCTAGCGGTTTTCGTTATCAGATAATc                                                             >  2:4699518/1‑150 (MQ=255)
                                                                                     ccTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATc                                                             <  2:2752524/150‑1 (MQ=255)
                                                                                      cTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCg                                                            >  2:2856337/1‑150 (MQ=255)
                                                                                      cTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCg                                                            <  1:4333424/150‑1 (MQ=255)
                                                                                       tGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGa                                                           <  2:4476640/150‑1 (MQ=255)
                                                                                          ggATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGt                                                        >  1:2657458/1‑150 (MQ=255)
                                                                                          ggATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGt                                                        <  2:2128045/150‑1 (MQ=255)
                                                                                          ggATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGt                                                        >  1:583605/1‑150 (MQ=255)
                                                                                            aTGTCGCGGAGCTGGTGCTGTATCAGCACTGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCa                                                      >  2:667525/1‑150 (MQ=255)
                                                                                            aTGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCa                                                      <  1:512858/150‑1 (MQ=255)
                                                                                              gTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCaaa                                                    <  2:749727/150‑1 (MQ=255)
                                                                                                                 aTCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCg                                 >  1:4501821/1‑150 (MQ=255)
                                                                                                                   cAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTc                               <  1:2005927/150‑1 (MQ=255)
                                                                                                                     gCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCAc                             <  2:1798555/150‑1 (MQ=255)
                                                                                                                       aaTGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCa                           >  1:568022/1‑150 (MQ=255)
                                                                                                                        aTGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCaa                          >  2:986736/1‑150 (MQ=255)
                                                                                                                        aTGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCaa                          >  2:2200935/1‑150 (MQ=255)
                                                                                                                        aTGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCaa                          >  2:2784630/1‑150 (MQ=255)
                                                                                                                         tGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAg                         <  2:3736752/150‑1 (MQ=255)
                                                                                                                           ggATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGAt                       >  2:2384742/1‑150 (MQ=255)
                                                                                                                             aTACCCGAGCCAGCATGCATATATTCACTGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCt                     <  2:3556038/150‑1 (MQ=255)
                                                                                                                             aTACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCt                     <  1:872324/150‑1 (MQ=255)
                                                                                                                             aTACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCt                     <  1:3967099/150‑1 (MQ=255)
                                                                                                                             aTACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTC‑AAAAAATGCCACTCGGCAGCGTCACCAAGATCTg                    <  1:4153196/150‑1 (MQ=255)
                                                                                                                               aCCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGc                   >  1:1377257/1‑150 (MQ=255)
                                                                                                                                   gAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCt               >  2:3848007/1‑150 (MQ=255)
                                                                                                                                   gAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCt               <  2:2636463/150‑1 (MQ=255)
                                                                                                                                          catgcatATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGa        >  2:3904933/1‑150 (MQ=255)
                                                                                                                                            tgcatATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATc      <  2:3458593/150‑1 (MQ=255)
                                                                                                                                              catATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCGCTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCac    <  1:838831/150‑1 (MQ=255)
                                                                                                                                               atatatTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCacg   >  1:501884/1‑150 (MQ=255)
                                                                                                                                               atatatTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCacg   >  2:2493357/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GCTATCGCTGGATGCGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGA  >  CP009273/463032‑463325

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: ‑