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,1790CA97.5% 741.5 / 9.0 237I104I (ATC→ATA) cofthiamine pyrimidine pyrophosphate hydrolase; HMP‑PP phosphatase
Reads supporting (aligned to +/- strand):  ref base C (5/1);  new base A (117/114);  total (122/115)
Fisher's exact test for biased strand distribution p-value = 2.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 241 total aligned reads displayed.

TATCGCTGGATGCGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCA  >  CP009273/463034‑463322
                                                                                                                                                 |                                                                                                                                               
tATCGCTGGATGCGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATAtt                                                                                                                                               <  1:3764680/148‑1 (MQ=255)
 aTCGCTGGATGCGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTc                                                                                                                                              >  1:334897/1‑148 (MQ=255)
 aTCGCTGGATGCGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTc                                                                                                                                              >  2:601789/1‑148 (MQ=255)
 aTCGCTGGATGCGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTc                                                                                                                                              >  2:602015/1‑148 (MQ=255)
  tCGCTGGATGCGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCa                                                                                                                                             <  1:913945/148‑1 (MQ=255)
   cGCTGGATGCGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCaa                                                                                                                                            <  1:1012807/148‑1 (MQ=255)
     cTGGATGCTTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTGTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATg                                                                                                                                          <  1:2496441/148‑1 (MQ=255)
     cTGGATGCGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATg                                                                                                                                          >  2:1489980/1‑148 (MQ=255)
         aTGCGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATgacga                                                                                                                                      <  1:299694/148‑1 (MQ=255)
           gCGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGacgacg                                                                                                                                    >  2:3751857/1‑148 (MQ=255)
            cGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACgg                                                                                                                                   <  1:686318/148‑1 (MQ=255)
            cGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACgg                                                                                                                                   <  1:938965/148‑1 (MQ=255)
               aTTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACggttg                                                                                                                                >  1:68609/1‑148 (MQ=255)
               aTTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACggttg                                                                                                                                >  2:2832032/1‑148 (MQ=255)
                tttGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACggttgg                                                                                                                               >  1:2681456/1‑148 (MQ=255)
                tttGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACggttgg                                                                                                                               >  2:2283084/1‑148 (MQ=255)
                 ttGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACggttggt                                                                                                                              >  1:421209/1‑148 (MQ=255)
                     ttACCGGCAACGGAACGCTCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTTGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTAc                                                                                                                          <  2:163922/148‑1 (MQ=255)
                      tACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGAGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTAcc                                                                                                                         >  2:744744/1‑148 (MQ=255)
                        ccGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCgg                                                                                                                       >  2:3673994/1‑148 (MQ=255)
                         cGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCggg                                                                                                                      >  2:1213378/1‑148 (MQ=255)
                            cAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCCGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGaaa                                                                                                                   >  1:156223/1‑148 (MQ=255)
                            cAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGaaa                                                                                                                   >  1:2584467/1‑148 (MQ=255)
                              aCGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAaga                                                                                                                 >  2:648006/1‑148 (MQ=255)
                                 gAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGAt                                                                                                              >  1:2888014/1‑148 (MQ=255)
                                 gAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGAt                                                                                                              >  1:1450010/1‑148 (MQ=255)
                                        cgTGCATTCTCTGGAAGGTTAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCg                                                                                                       <  2:3016414/148‑1 (MQ=255)
                                         gTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCgt                                                                                                      <  1:368496/148‑1 (MQ=255)
                                         gTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCgt                                                                                                      <  2:2872111/148‑1 (MQ=255)
                                         gTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCgt                                                                                                      <  2:2871804/148‑1 (MQ=255)
                                          tGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCgtt                                                                                                     >  2:1856680/1‑148 (MQ=255)
                                          tGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCgtt                                                                                                     >  1:3671548/1‑148 (MQ=255)
                                            cATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCgttgt                                                                                                   >  1:1797504/1‑148 (MQ=255)
                                              ttCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGttgttg                                                                                                 <  1:1746513/148‑1 (MQ=255)
                                               tctcTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGc                                                                                                >  2:2686949/1‑148 (MQ=255)
                                                   tGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGc                                                                                            <  2:914035/148‑1 (MQ=255)
                                                     gAAGGTGAACTTTTACATCGTGATGATTTACCTGCTGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCAt                                                                                          <  1:2161798/148‑1 (MQ=255)
                                                     gAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCAt                                                                                          <  1:3590805/148‑1 (MQ=255)
                                                     gAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGGTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCAt                                                                                          >  2:3272846/1‑148 (MQ=255)
                                                         gTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGt                                                                                      <  1:2372806/148‑1 (MQ=255)
                                                              cTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAATGAGATCCCTGCGTTGTTGCAGGCATTTGTCtata                                                                                 >  2:2379803/1‑148 (MQ=255)
                                                              cTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCtata                                                                                 >  2:3279065/1‑148 (MQ=255)
                                                                tttACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGc                                                                               >  1:3436852/1‑148 (MQ=255)
                                                                 ttACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCg                                                                              <  2:2399733/148‑1 (MQ=255)
                                                                  tACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCgg                                                                             <  2:3550683/148‑1 (MQ=255)
                                                                    cATCGTGATGATTTACCTGAGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGtt                                                                           <  1:810678/148‑1 (MQ=255)
                                                                       cGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTc                                                                        >  1:48650/1‑148 (MQ=255)
                                                                        gTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCg                                                                       <  2:518454/148‑1 (MQ=255)
                                                                          gatgatTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGtt                                                                     <  1:1630950/148‑1 (MQ=255)
                                                                          gatgatTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGtt                                                                     <  1:1381915/148‑1 (MQ=255)
                                                                           atgatTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTa                                                                    >  1:3793191/1‑148 (MQ=255)
                                                                           atgatTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTa                                                                    <  2:2004676/148‑1 (MQ=255)
                                                                            tgatTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTAt                                                                   <  2:2778569/148‑1 (MQ=255)
                                                                             gatTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTTCCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATc                                                                  >  1:2519119/1‑148 (MQ=255)
                                                                               tttACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAg                                                                <  2:461090/148‑1 (MQ=255)
                                                                                ttACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGa                                                               >  1:3366573/1‑148 (MQ=255)
                                                                                  aCCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATa                                                             <  1:1426667/148‑1 (MQ=255)
                                                                                  aCCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATa                                                             <  1:1434158/148‑1 (MQ=255)
                                                                                   ccTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATaa                                                            >  1:1273472/1‑148 (MQ=255)
                                                                                   ccTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATaa                                                            <  1:2695514/148‑1 (MQ=255)
                                                                                    cTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTCTAGCGGTTTTCGTTATCAGATAAt                                                           <  1:3780440/148‑1 (MQ=255)
                                                                                    cTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAAt                                                           <  2:952760/148‑1 (MQ=255)
                                                                                     tGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATc                                                          >  2:1693361/1‑148 (MQ=255)
                                                                                     tGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATc                                                          >  1:282260/1‑148 (MQ=255)
                                                                                      gCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCg                                                         >  1:197932/1‑148 (MQ=255)
                                                                                      gCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCg                                                         >  1:2126602/1‑148 (MQ=255)
                                                                                      gCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATAACCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCg                                                         <  2:2516819/148‑1 (MQ=255)
                                                                                       cGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGa                                                        >  1:2255339/1‑148 (MQ=255)
                                                                                       cGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGa                                                        >  1:2344355/1‑148 (MQ=255)
                                                                                       cGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGa                                                        <  2:3514147/148‑1 (MQ=255)
                                                                                           tGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTc                                                    >  2:328301/1‑148 (MQ=255)
                                                                                             tCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCaa                                                  >  1:3295835/1‑148 (MQ=255)
                                                                                                 ggAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCATGCGTTCATGCAGGCATGTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCaaaaaa                                              >  1:1534086/1‑148 (MQ=255)
                                                                                                  gAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCaaaataa                                             <  1:1969053/148‑1 (MQ=255)
                                                                                                     cTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGc                                          <  1:2832032/148‑1 (MQ=255)
                                                                                                      tGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGcc                                         <  2:2208855/148‑1 (MQ=255)
                                                                                                       ggTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCa                                        <  1:1145554/148‑1 (MQ=255)
                                                                                                        gTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCAc                                       >  2:1858426/1‑148 (MQ=255)
                                                                                                         tGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACt                                      <  2:2097214/148‑1 (MQ=255)
                                                                                                         tGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACt                                      <  2:1344038/148‑1 (MQ=255)
                                                                                                           cTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCg                                    >  1:2550364/1‑148 (MQ=255)
                                                                                                             gTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGc                                  >  2:139374/1‑148 (MQ=255)
                                                                                                                tCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGc                               >  1:2736885/1‑148 (MQ=255)
                                                                                                                 cAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCg                              >  2:1529719/1‑148 (MQ=255)
                                                                                                                 cAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCg                              >  1:2632184/1‑148 (MQ=255)
                                                                                                                   gCAATGGGATACCCGAGCCCGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTc                            >  2:1904739/1‑148 (MQ=255)
                                                                                                                       tGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCa                        <  1:3149425/148‑1 (MQ=255)
                                                                                                                         ggATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAg                      <  1:1322359/148‑1 (MQ=255)
                                                                                                                         ggATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAg                      >  2:3182060/1‑148 (MQ=255)
                                                                                                                          gATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGa                     >  1:34098/1‑148 (MQ=255)
                                                                                                                           aTACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGAt                    <  2:72827/148‑1 (MQ=255)
                                                                                                                           aTAACCGAGCCATCATTCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCATGCATTTGTCTATAGCGGTTTTCTTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGAt                    <  2:1718392/148‑1 (MQ=255)
                                                                                                                               ccGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGc                <  2:3790975/148‑1 (MQ=255)
                                                                                                                                    ccAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCtg           <  1:3639921/148‑1 (MQ=255)
                                                                                                                                     cAGCATGCATAGATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCtgt          <  2:3426435/148‑1 (MQ=255)
                                                                                                                                        catgcatATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGc       >  2:2967559/1‑148 (MQ=255)
                                                                                                                                          tgcatATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGa     >  2:1258598/1‑148 (MQ=255)
                                                                                                                                          tgaatatatTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGa     <  2:2701380/145‑1 (MQ=255)
                                                                                                                                           gcatATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGAt    >  2:1777414/1‑148 (MQ=255)
                                                                                                                                             atatatTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCa  <  2:2888014/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
TATCGCTGGATGCGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCA  >  CP009273/463034‑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: ‑