Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 463,179 C→A 7.2% 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,1790CA 7.2% 978.1 / 31.9 258I104I (ATC→ATAcofthiamine pyrimidine pyrophosphate hydrolase; HMP‑PP phosphatase
Reads supporting (aligned to +/- strand):  ref base C (117/122);  new base A (10/9);  total (127/131)
Fisher's exact test for biased strand distribution p-value = 8.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 262 total aligned reads displayed.

GCGCTATCGCTGGATGCGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGACGA  >  CP009273/463030‑463328
                                                                                                                                                     |                                                                                                                                                     
 cgcTATCGCTGGATGCGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCt                                                                                                                                                      >  2:2217927/1‑150 (MQ=255)
 cgcTATCGCTGGATGCGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCt                                                                                                                                                      >  1:2553845/1‑150 (MQ=255)
   cTATCGCTGGATGCGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTc                                                                                                                                                    <  2:1487740/150‑1 (MQ=255)
    tATCGCTGGATGCGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCa                                                                                                                                                   >  2:4413465/1‑150 (MQ=255)
     aTCGCTGGATGCGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCTGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCaa                                                                                                                                                  <  2:2015294/150‑1 (MQ=255)
     aTCGCTGGATGCGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCaa                                                                                                                                                  >  2:4582447/1‑150 (MQ=255)
         cTGGATGCGTATTTGATTACCGTCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGAAATGGGATACCCGAGCCAGCATGCATATCTTCAATgac                                                                                                                                              <  2:2052165/150‑1 (MQ=255)
         cTGGATGCGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATgac                                                                                                                                              <  2:2861101/150‑1 (MQ=255)
         cTGGATGCGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATgac                                                                                                                                              >  2:1880744/1‑150 (MQ=255)
            gATGCGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATgacgac                                                                                                                                           <  2:414893/150‑1 (MQ=255)
            gATGCGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATgacgac                                                                                                                                           <  2:3181252/150‑1 (MQ=255)
            gATGCGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATgacgac                                                                                                                                           <  2:415632/150‑1 (MQ=255)
                cGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACggtt                                                                                                                                       <  1:1952073/150‑1 (MQ=255)
                 gTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACggttg                                                                                                                                      >  1:366738/1‑150 (MQ=255)
                       gATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTAc                                                                                                                                >  2:1278923/1‑150 (MQ=255)
                        aTTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTAcc                                                                                                                               <  1:1670859/150‑1 (MQ=255)
                          tACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCgg                                                                                                                             >  2:3676848/1‑150 (MQ=255)
                          tACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCgg                                                                                                                             <  2:4172244/150‑1 (MQ=255)
                            ccGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGCTTACCGGGa                                                                                                                           <  2:358127/150‑1 (MQ=255)
                              ggCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGaaa                                                                                                                         >  1:3546381/1‑150 (MQ=255)
                                  aCGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAgaga                                                                                                                     <  1:2698117/150‑1 (MQ=255)
                                  aCGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAgaga                                                                                                                     <  2:4063166/150‑1 (MQ=255)
                                   cGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGAt                                                                                                                    <  1:2025809/150‑1 (MQ=255)
                                   cGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGAt                                                                                                                    >  1:4360425/1‑150 (MQ=255)
                                    ggAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATc                                                                                                                   >  1:3152158/1‑150 (MQ=255)
                                    ggAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATc                                                                                                                   >  1:1233470/1‑150 (MQ=255)
                                    ggAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATc                                                                                                                   >  1:289438/1‑150 (MQ=255)
                                        cgcgcgTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTg                                                                                                               <  1:1549780/150‑1 (MQ=255)
                                        cgcgcgTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTg                                                                                                               <  1:1945831/150‑1 (MQ=255)
                                        cgcgcgTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTg                                                                                                               >  1:3775311/1‑150 (MQ=255)
                                         gcgcgTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGc                                                                                                              <  1:1018667/150‑1 (MQ=255)
                                             gTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCgttg                                                                                                          >  2:2355828/1‑150 (MQ=255)
                                                   tctcTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAg                                                                                                    <  2:326194/150‑1 (MQ=255)
                                                     tctGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGTATACCCGAGCCAGCATGCATATCTTCAATGACTACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGc                                                                                                  >  2:1153707/1‑150 (MQ=255)
                                                     tctGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGc                                                                                                  <  1:592903/150‑1 (MQ=255)
                                                     tctGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGc                                                                                                  <  1:338678/150‑1 (MQ=255)
                                                       tGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCAt                                                                                                <  1:4205905/150‑1 (MQ=255)
                                                        ggAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCAtt                                                                                               <  2:2141587/150‑1 (MQ=255)
                                                        ggAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCAtt                                                                                               <  1:1901895/150‑1 (MQ=255)
                                                         gAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCAttt                                                                                              <  2:2048735/150‑1 (MQ=255)
                                                         gAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCAttt                                                                                              <  2:3703046/150‑1 (MQ=255)
                                                           aGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGt                                                                                            <  2:1637727/150‑1 (MQ=255)
                                                           aGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGt                                                                                            >  1:1330426/1‑150 (MQ=255)
                                                            ggTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTc                                                                                           <  2:4303540/150‑1 (MQ=255)
                                                            ggTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTc                                                                                           <  1:760927/150‑1 (MQ=255)
                                                             gTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCt                                                                                          <  2:3642976/150‑1 (MQ=255)
                                                                aaCTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCtata                                                                                       >  2:743048/1‑150 (MQ=255)
                                                                    tttACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCgg                                                                                   <  1:2063858/150‑1 (MQ=255)
                                                                     ttACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGt                                                                                  >  2:2391319/1‑150 (MQ=255)
                                                                      tACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGtt                                                                                 <  2:4002122/150‑1 (MQ=255)
                                                                        cATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGtttt                                                                               >  2:2574998/1‑150 (MQ=255)
                                                                         aTCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTc                                                                              <  1:75384/150‑1 (MQ=255)
                                                                         aTCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTc                                                                              <  1:4043754/150‑1 (MQ=255)
                                                                         aTCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTc                                                                              >  2:2653246/1‑150 (MQ=255)
                                                                         aTCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTc                                                                              >  1:2969092/1‑150 (MQ=255)
                                                                             tgatgaTTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTa                                                                          <  2:2894599/150‑1 (MQ=255)
                                                                              gatgatTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTAt                                                                         <  2:2797382/150‑1 (MQ=255)
                                                                               atgatTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATc                                                                        <  1:327737/150‑1 (MQ=255)
                                                                                  atTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGa                                                                     >  1:841485/1‑150 (MQ=255)
                                                                                  atTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGa                                                                     <  1:3760463/150‑1 (MQ=255)
                                                                                        cTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCg                                                               <  1:1841277/150‑1 (MQ=255)
                                                                                         tGCGGATGTCGCTGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGa                                                              >  2:1794043/1‑150 (MQ=255)
                                                                                             gATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTc                                                          >  1:530392/1‑150 (MQ=255)
                                                                                              aTGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCa                                                         <  1:3145852/150‑1 (MQ=255)
                                                                                               tGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCCGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCAGTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCaa                                                        >  2:26820/1‑150 (MQ=255)
                                                                                               tGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATTATCGATGTCaa                                                        <  1:2396500/150‑1 (MQ=255)
                                                                                               tGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCaa                                                        >  2:1174862/1‑150 (MQ=255)
                                                                                                    cgGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCaaaaaaa                                                   >  2:2464536/1‑150 (MQ=255)
                                                                                                    cgGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCaaaaaaa                                                   >  2:1243179/1‑150 (MQ=255)
                                                                                                     ggAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAAt                                                  <  1:3713840/150‑1 (MQ=255)
                                                                                                        gCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGcc                                               <  1:2429077/150‑1 (MQ=255)
                                                                                                          tGGTGCTGTATCAGCATTGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCAc                                             <  2:1352840/150‑1 (MQ=255)
                                                                                                           ggTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACt                                            <  2:2931400/150‑1 (MQ=255)
                                                                                                            gTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTc                                           <  1:1427549/150‑1 (MQ=255)
                                                                                                                tGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCa                                       <  1:2728276/150‑1 (MQ=255)
                                                                                                                   aTCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCg                                    <  1:464730/150‑1 (MQ=255)
                                                                                                                        cAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATTATCGATGTCAAAAAAATGCCACTCGGCAGCGTCAcc                               <  2:653896/150‑1 (MQ=255)
                                                                                                                        cAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCAcc                               >  1:1111325/1‑150 (MQ=255)
                                                                                                                         aaTGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCa                              >  1:2066485/1‑150 (MQ=255)
                                                                                                                         aaTGGGATACCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCa                              >  2:4547284/1‑150 (MQ=255)
                                                                                                                           tGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAg                            >  1:4010787/1‑150 (MQ=255)
                                                                                                                             ggATAGCCGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGATAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATATTCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGAt                          <  2:4102529/150‑1 (MQ=255)
                                                                                                                             ggATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATTATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGAt                          >  2:3726968/1‑150 (MQ=255)
                                                                                                                             ggATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGAt                          >  2:3962382/1‑150 (MQ=255)
                                                                                                                                   ccGAGCCAGCATGCATATATTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTc                   <  2:1655828/150‑1 (MQ=255)
                                                                                                                                    cGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATTATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTc                   >  1:2202080/1‑150 (MQ=255)
                                                                                                                                    cGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTc                   >  1:453526/1‑150 (MQ=255)
                                                                                                                                         cAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAATAGATCCCTGCGTTGTTGCAGGAATTTGTCTATAGCGGTTTTCGTTATCAGATTATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTgg              <  2:184661/150‑1 (MQ=255)
                                                                                                                                         cAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTgg              >  1:4410258/1‑150 (MQ=255)
                                                                                                                                         cAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTgg              <  1:4507526/150‑1 (MQ=255)
                                                                                                                                          aGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGc             >  1:3321172/1‑150 (MQ=255)
                                                                                                                                          aGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGc             <  2:4578172/150‑1 (MQ=255)
                                                                                                                                          aGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGc             <  1:3952435/150‑1 (MQ=255)
                                                                                                                                           gcatgcatATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCg            >  1:1477648/1‑150 (MQ=255)
                                                                                                                                             atgcatATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAATTGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGAt          >  2:1913221/1‑150 (MQ=255)
                                                                                                                                              tgcatATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATc         <  2:716741/150‑1 (MQ=255)
                                                                                                                                               gcatATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCa        >  2:1674952/1‑150 (MQ=255)
                                                                                                                                                  tatCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCacga     >  1:938235/1‑150 (MQ=255)
                                                                                                                                                   atCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATTATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCacgac    >  2:4129148/1‑150 (MQ=255)
                                                                                                                                                     cTTCAATGACGACGGTTGGTTTACCGGGAACGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATTATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCAcgacga  >  1:704260/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCGCTATCGCTGGATGCGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACGGTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGACGA  >  CP009273/463030‑463328

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: