Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 463,179 C→A 87.3% 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,1790CA87.3% 527.0 / 71.5 219I104I (ATC→ATAcofthiamine pyrimidine pyrophosphate hydrolase; HMP‑PP phosphatase
Reads supporting (aligned to +/- strand):  ref base C (10/18);  new base A (84/107);  total (94/125)
Fisher's exact test for biased strand distribution p-value = 5.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 220 total aligned reads displayed.

CTATCGCTGGATGCGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCA  >  CP009273/463033‑463322
                                                                                                                                                  |                                                                                                                                               
   tCGCTGGATGCGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCCGCATGCATATATTCa                                                                                                                                             >  2:2070715/1‑148 (MQ=255)
   tCGCTGGATGCGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCa                                                                                                                                             <  2:2557924/148‑1 (MQ=255)
    cGCTGGATGCGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCaa                                                                                                                                            <  1:2694325/148‑1 (MQ=255)
    cGCTGGATGCGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCaa                                                                                                                                            <  1:919692/148‑1 (MQ=255)
      cTGGATGCGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATg                                                                                                                                          >  2:909056/1‑148 (MQ=255)
      cTGGATGCGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATg                                                                                                                                          >  2:361385/1‑148 (MQ=255)
      cTGGATGCGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATg                                                                                                                                          <  2:2733877/148‑1 (MQ=255)
       tGGATGCGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATga                                                                                                                                         <  2:1103602/148‑1 (MQ=255)
        ggATGCGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATgac                                                                                                                                        >  1:2961755/1‑148 (MQ=255)
         gATGCGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATgacg                                                                                                                                       <  2:830041/148‑1 (MQ=255)
         gATGCGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATgacg                                                                                                                                       >  1:1195751/1‑148 (MQ=255)
            gCGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGacgacg                                                                                                                                    >  1:2031995/1‑148 (MQ=255)
             cGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACgg                                                                                                                                   <  2:394682/148‑1 (MQ=255)
              gTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACggt                                                                                                                                  <  1:1446375/148‑1 (MQ=255)
               tATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACggtt                                                                                                                                 <  1:717587/148‑1 (MQ=255)
                aTTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACggttg                                                                                                                                >  1:1723517/1‑148 (MQ=255)
                aTTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACggttg                                                                                                                                <  1:3198986/148‑1 (MQ=255)
                aTTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACggttg                                                                                                                                <  2:86906/148‑1 (MQ=255)
                aTTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACggttg                                                                                                                                >  2:2825891/1‑148 (MQ=255)
                       tACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTAcc                                                                                                                         >  2:1401306/1‑148 (MQ=255)
                              aaCGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAag                                                                                                                  <  2:1444753/148‑1 (MQ=255)
                               aCGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAaga                                                                                                                 >  2:2810964/1‑148 (MQ=255)
                                  gAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGAt                                                                                                              >  1:533288/1‑148 (MQ=255)
                                  gAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGAt                                                                                                              >  2:2563609/1‑148 (MQ=255)
                                   aaCGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATc                                                                                                             >  1:890514/1‑148 (MQ=255)
                                    tcgcgcgTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATcc                                                                                                            <  1:347410/147‑1 (MQ=255)
                                    aCGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATcc                                                                                                            <  2:1228134/148‑1 (MQ=255)
                                         cgTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCg                                                                                                       <  1:404500/148‑1 (MQ=255)
                                         cgTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCg                                                                                                       <  2:2031995/148‑1 (MQ=255)
                                             cATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCgttgt                                                                                                   <  1:2652610/148‑1 (MQ=255)
                                             cATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATACAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCgttgt                                                                                                   >  1:2041702/1‑148 (MQ=255)
                                              aTTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCgttgtt                                                                                                  >  2:3003809/1‑148 (MQ=255)
                                                tctcTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCCGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGc                                                                                                >  2:1992112/1‑148 (MQ=255)
                                                tctcTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGc                                                                                                <  1:3147696/148‑1 (MQ=255)
                                                 ctctGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCa                                                                                               <  1:1396619/148‑1 (MQ=255)
                                                 ctctGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCa                                                                                               <  1:766762/148‑1 (MQ=255)
                                                  tctGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTAAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAg                                                                                              >  2:172849/1‑148 (MQ=255)
                                                     ggAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCa                                                                                           >  1:2771568/1‑148 (MQ=255)
                                                     ggAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCa                                                                                           >  2:2104095/1‑148 (MQ=255)
                                                     ggAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCa                                                                                           >  1:511231/1‑148 (MQ=255)
                                                       aaGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCAtt                                                                                         <  1:1897457/148‑1 (MQ=255)
                                                         ggTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTg                                                                                       >  2:1500035/1‑148 (MQ=255)
                                                          gTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCCGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGt                                                                                      >  2:2382707/1‑148 (MQ=255)
                                                          gTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGt                                                                                      <  1:2430977/148‑1 (MQ=255)
                                                            gAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCt                                                                                    >  1:2183675/1‑148 (MQ=255)
                                                              aCTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCtat                                                                                  <  1:122374/148‑1 (MQ=255)
                                                              aCTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCtat                                                                                  <  1:122564/148‑1 (MQ=255)
                                                                ttttACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAg                                                                                >  2:1782280/1‑148 (MQ=255)
                                                                ttttACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAg                                                                                >  2:1780991/1‑148 (MQ=255)
                                                                    aCATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTTTATCAGCTATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGt                                                                            <  2:555825/148‑1 (MQ=255)
                                                                    aCATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGt                                                                            <  2:2306851/148‑1 (MQ=255)
                                                                       tCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGtttt                                                                         <  1:2628686/148‑1 (MQ=255)
                                                                          tgatgaTTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGt                                                                      <  1:949837/148‑1 (MQ=255)
                                                                          tgatgaTTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGt                                                                      <  1:2825891/148‑1 (MQ=255)
                                                                          tgatgaTTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGt                                                                      <  2:3020653/148‑1 (MQ=255)
                                                                            atgatTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTa                                                                    >  1:1280462/1‑148 (MQ=255)
                                                                                tttACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAg                                                                <  2:2547268/148‑1 (MQ=255)
                                                                                 ttACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCCGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGa                                                               >  2:1575329/1‑148 (MQ=255)
                                                                                 ttACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGa                                                               >  2:1320224/1‑148 (MQ=255)
                                                                                   aCCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATa                                                             <  1:1666845/148‑1 (MQ=255)
                                                                                    ccTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATaa                                                            >  1:1817426/1‑148 (MQ=255)
                                                                                    ccTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATaa                                                            <  1:2196627/148‑1 (MQ=255)
                                                                                       gCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCg                                                         <  2:1234214/148‑1 (MQ=255)
                                                                                       gCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCg                                                         <  2:1467563/148‑1 (MQ=255)
                                                                                          gATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATg                                                      <  1:3254072/148‑1 (MQ=255)
                                                                                           aTGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGt                                                     >  1:2102739/1‑148 (MQ=255)
                                                                                                   gAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTaaaaaaaa                                             >  2:1427217/1‑148 (MQ=255)
                                                                                                   gAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCaaaaaaa                                             <  1:1059272/148‑1 (MQ=255)
                                                                                                    aGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAAt                                            <  1:2563609/148‑1 (MQ=255)
                                                                                                        ggTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCa                                        <  2:1729517/148‑1 (MQ=255)
                                                                                                          tGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACt                                      <  2:3304027/148‑1 (MQ=255)
                                                                                                             tGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCgg                                   <  2:2902703/148‑1 (MQ=255)
                                                                                                             tGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCgg                                   <  1:961921/148‑1 (MQ=255)
                                                                                                             tGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCgg                                   <  2:1035933/148‑1 (MQ=255)
                                                                                                              gTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGc                                  >  2:1257753/1‑148 (MQ=255)
                                                                                                                          ggATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAg                      >  1:2519560/1‑148 (MQ=255)
                                                                                                                          ggATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAg                      >  1:2335894/1‑148 (MQ=255)
                                                                                                                          ggATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAg                      <  2:3183050/148‑1 (MQ=255)
                                                                                                                            aTACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGAt                    >  1:2550299/1‑148 (MQ=255)
                                                                                                                              aCCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCt                  >  2:1313042/1‑148 (MQ=255)
                                                                                                                                ccGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGTGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGc                <  1:1611405/148‑1 (MQ=255)
                                                                                                                                ccGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGc                >  2:2685284/1‑148 (MQ=255)
                                                                                                                                 cGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCt               >  2:2335658/1‑148 (MQ=255)
                                                                                                                                   aGCCAGCATGCATATCGTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTc             <  2:1890612/148‑1 (MQ=255)
                                                                                                                                   aGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTc             <  2:1723091/148‑1 (MQ=255)
                                                                                                                                    gCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCt            <  1:153024/148‑1 (MQ=255)
                                                                                                                                     ccAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCtg           <  1:196459/148‑1 (MQ=255)
                                                                                                                                       aGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCtgtg         <  2:1197186/148‑1 (MQ=255)
                                                                                                                                       aGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCtgtg         >  2:3253994/1‑148 (MQ=255)
                                                                                                                                         catgcatATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGc       >  2:1953016/1‑148 (MQ=255)
                                                                                                                                          atgcatATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCg      <  2:246111/148‑1 (MQ=255)
                                                                                                                                           tgcatATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGa     >  1:2355290/1‑148 (MQ=255)
                                                                                                                                           tgcatATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGa     <  1:2899893/148‑1 (MQ=255)
                                                                                                                                            gcatATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATTTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGAt    <  2:164100/148‑1 (MQ=255)
                                                                                                                                             catATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATc   <  2:234673/148‑1 (MQ=255)
                                                                                                                                             catATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATc   <  2:1668023/148‑1 (MQ=255)
                                                                                                                                             catATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATc   >  2:718417/1‑148 (MQ=255)
                                                                                                                                             catATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATc   >  1:230075/1‑148 (MQ=255)
                                                                                                                                             catATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATc   >  1:1165348/1‑148 (MQ=255)
                                                                                                                                              atatatTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCa  >  1:2892835/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
CTATCGCTGGATGCGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCA  >  CP009273/463033‑463322

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: