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,1790CA96.2% 644.1 / 14.8 212I104I (ATC→ATA) cofthiamine pyrimidine pyrophosphate hydrolase; HMP‑PP phosphatase
Reads supporting (aligned to +/- strand):  ref base C (3/5);  new base A (86/117);  total (89/123)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 219 total aligned reads displayed.

GCTATCGCTGGATGCGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCAC  >  CP009273/463032‑463323
                                                                                                                                                   |                                                                                                                                                
  tATCGCTGGATGCGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATAtt                                                                                                                                                >  1:3116098/1‑148 (MQ=255)
  tATCGCTGGATGCGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATAtt                                                                                                                                                >  2:3385320/1‑148 (MQ=255)
    tCGCTGGATGCGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCa                                                                                                                                              >  1:346600/1‑148 (MQ=255)
    tCGCTGGATGCGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCa                                                                                                                                              <  1:2307629/148‑1 (MQ=255)
     cGCTGGATGCGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCaa                                                                                                                                             <  2:1940542/148‑1 (MQ=255)
          gATGCGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATgacg                                                                                                                                        >  2:3064407/1‑148 (MQ=255)
             gCGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGacgacg                                                                                                                                     <  2:1465033/148‑1 (MQ=255)
              cGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACgg                                                                                                                                    <  1:808760/148‑1 (MQ=255)
              cGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACgg                                                                                                                                    <  1:2191634/148‑1 (MQ=255)
                 aTTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCCTGCATATATTCAATGACGACggttg                                                                                                                                 >  1:2016054/1‑148 (MQ=255)
                 aTTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACggttg                                                                                                                                 <  1:581084/148‑1 (MQ=255)
                  tttGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACggttgg                                                                                                                                <  2:1540330/148‑1 (MQ=255)
                  tttGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACggttgg                                                                                                                                <  2:1541327/148‑1 (MQ=255)
                   ttGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACggttggt                                                                                                                               <  2:2553147/148‑1 (MQ=255)
                      aTTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTa                                                                                                                            <  2:263964/148‑1 (MQ=255)
                      aTTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTa                                                                                                                            <  1:2476932/148‑1 (MQ=255)
                        tACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTAcc                                                                                                                          >  1:1871257/1‑148 (MQ=255)
                        tACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTAcc                                                                                                                          <  2:2983841/148‑1 (MQ=255)
                        tACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTAcc                                                                                                                          >  2:1770885/1‑148 (MQ=255)
                         aCCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCg                                                                                                                         <  2:2220658/148‑1 (MQ=255)
                         aCCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCg                                                                                                                         <  2:3093453/148‑1 (MQ=255)
                         aCCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCg                                                                                                                         <  1:2539006/148‑1 (MQ=255)
                            ggCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGa                                                                                                                      <  2:1572611/148‑1 (MQ=255)
                             gCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGaa                                                                                                                     >  2:1819432/1‑148 (MQ=255)
                               aaCGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAag                                                                                                                   <  1:1241373/148‑1 (MQ=255)
                               aaCGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAag                                                                                                                   <  2:2105155/148‑1 (MQ=255)
                                 cGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAagag                                                                                                                 <  1:60372/148‑1 (MQ=255)
                                    aaCGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATc                                                                                                              >  1:1698376/1‑148 (MQ=255)
                                      cgcgcgTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATccc                                                                                                            <  1:2465063/148‑1 (MQ=255)
                                      cgcgcgTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATccc                                                                                                            >  2:2730133/1‑148 (MQ=255)
                                           gTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCgt                                                                                                       <  2:1645613/148‑1 (MQ=255)
                                           gTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCgt                                                                                                       <  1:1153359/148‑1 (MQ=255)
                                              cATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCgttgt                                                                                                    <  2:466016/148‑1 (MQ=255)
                                               aTTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCgttgtt                                                                                                   >  1:701378/1‑148 (MQ=255)
                                                ttCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGttgttg                                                                                                  <  1:3007816/148‑1 (MQ=255)
                                                ttCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGttgttg                                                                                                  <  2:2326417/148‑1 (MQ=255)
                                                ttCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGttgttg                                                                                                  >  1:479982/1‑148 (MQ=255)
                                                 tctcTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGc                                                                                                 >  1:146686/1‑148 (MQ=255)
                                                  ctctGGAAGGTGAACTTTTACATCGTGATGATTTACGTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGNGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCa                                                                                                >  1:1807640/1‑148 (MQ=255)
                                                    ctGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAgg                                                                                              <  1:2606965/148‑1 (MQ=255)
                                                    ctGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAgg                                                                                              <  2:2988448/148‑1 (MQ=255)
                                                     tGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGc                                                                                             <  1:2451744/148‑1 (MQ=255)
                                                          ggTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTg                                                                                        >  1:1972572/1‑148 (MQ=255)
                                                          ggTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTg                                                                                        >  2:701423/1‑148 (MQ=255)
                                                            tGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTc                                                                                      <  1:1384547/148‑1 (MQ=255)
                                                             gAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCt                                                                                     >  2:1034401/1‑148 (MQ=255)
                                                                 ttttACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAg                                                                                 <  1:3327195/148‑1 (MQ=255)
                                                                   ttACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCg                                                                               <  2:42804/148‑1 (MQ=255)
                                                                    tACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCgg                                                                              <  2:1730075/148‑1 (MQ=255)
                                                                    tACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCgg                                                                              <  1:1695514/148‑1 (MQ=255)
                                                                       aTCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGttt                                                                           <  2:3443883/148‑1 (MQ=255)
                                                                        tCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGtttt                                                                          <  2:1871257/148‑1 (MQ=255)
                                                                          gTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCg                                                                        <  1:3214126/148‑1 (MQ=255)
                                                                          gTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCg                                                                        <  2:2468748/148‑1 (MQ=255)
                                                                          gTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCg                                                                        <  1:2225926/148‑1 (MQ=255)
                                                                            gatgatTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCAAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGtt                                                                      <  1:3520794/148‑1 (MQ=255)
                                                                             atgatTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCCGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTa                                                                     >  1:1632916/1‑148 (MQ=255)
                                                                             atgatTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTa                                                                     <  2:3160387/148‑1 (MQ=255)
                                                                               gatTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGACAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATc                                                                   >  1:2692117/1‑148 (MQ=255)
                                                                                 tttACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAg                                                                 <  2:1310654/148‑1 (MQ=255)
                                                                                    aCCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATa                                                              <  1:49535/148‑1 (MQ=255)
                                                                                    aCCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATa                                                              <  2:3372578/148‑1 (MQ=255)
                                                                                    aCCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATa                                                              <  1:2114529/148‑1 (MQ=255)
                                                                                     ccTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATaa                                                             <  2:554240/148‑1 (MQ=255)
                                                                                      cTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAAt                                                            <  2:762621/148‑1 (MQ=255)
                                                                                        gCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCg                                                          <  2:968330/148‑1 (MQ=255)
                                                                                         cGGTTGTCGCGGAGCTGGTGCTGTATGAGCAATGGGATACCCGAGCCAGCATGCATATTTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGa                                                         <  1:3561682/148‑1 (MQ=255)
                                                                                         cGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATaa                                                             <  1:1212437/144‑1 (MQ=255)
                                                                                         cGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATaa                                                             >  2:1212437/1‑144 (MQ=255)
                                                                                         cGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGa                                                         <  2:2162548/148‑1 (MQ=255)
                                                                                            aTGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGt                                                      >  1:1176107/1‑148 (MQ=255)
                                                                                            aTGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGt                                                      >  1:860844/1‑148 (MQ=255)
                                                                                             tGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTc                                                     >  1:1724843/1‑148 (MQ=255)
                                                                                                    gAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCaaaaaaa                                              <  2:2313382/148‑1 (MQ=255)
                                                                                                     aGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAAt                                             >  1:3127463/1‑148 (MQ=255)
                                                                                                      gCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATg                                            <  2:2413152/148‑1 (MQ=255)
                                                                                                        tGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGcc                                          <  2:610555/148‑1 (MQ=255)
                                                                                                        tGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGcc                                          >  2:1595883/1‑148 (MQ=255)
                                                                                                         ggTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCa                                         <  2:1233462/148‑1 (MQ=255)
                                                                                                          gTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCAc                                        >  2:2017885/1‑148 (MQ=255)
                                                                                                           tGCTTTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACt                                       <  1:1841897/148‑1 (MQ=255)
                                                                                                           tGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACt                                       <  1:39255/148‑1 (MQ=255)
                                                                                                               gTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGc                                   >  1:2863289/1‑148 (MQ=255)
                                                                                                                   cAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCg                               >  2:1396909/1‑148 (MQ=255)
                                                                                                                    aGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGt                              <  2:1784097/148‑1 (MQ=255)
                                                                                                                     gCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTc                             <  2:28783/148‑1 (MQ=255)
                                                                                                                      cAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCa                            <  1:122579/148‑1 (MQ=255)
                                                                                                                       aaTGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCAc                           <  2:192032/148‑1 (MQ=255)
                                                                                                                           ggATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAg                       >  2:3466772/1‑148 (MQ=255)
                                                                                                                             aTACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGAt                     >  1:2952167/1‑148 (MQ=255)
                                                                                                                              tACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATc                    <  1:1093556/148‑1 (MQ=255)
                                                                                                                                 ccGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGc                 <  2:2642721/148‑1 (MQ=255)
                                                                                                                                     gCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCt             <  1:2460117/148‑1 (MQ=255)
                                                                                                                                         gcatgcatATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTgg         >  2:7011/1‑148 (MQ=255)
                                                                                                                                            tgcatATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGa      >  2:71025/1‑148 (MQ=255)
                                                                                                                                             gcatATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGAt     >  1:2287578/1‑148 (MQ=255)
                                                                                                                                             gcatATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGAt     <  1:1662072/148‑1 (MQ=255)
                                                                                                                                              catATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATCGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATc    >  2:1807474/1‑148 (MQ=255)
                                                                                                                                               atatatTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCa   >  2:2355603/1‑148 (MQ=255)
                                                                                                                                                tataTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCac  >  2:321307/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
GCTATCGCTGGATGCGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCAC  >  CP009273/463032‑463323

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