Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,738,088 T→A 8.2% F134Y (TTC→TAC)  yiaM → 2,3‑diketo‑L‑gulonate TRAP transporter small permease protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,738,0880TA 8.2% 630.6 / 39.0 206F134Y (TTC→TAC) yiaM2,3‑diketo‑L‑gulonate TRAP transporter small permease protein
Reads supporting (aligned to +/- strand):  ref base T (85/104);  new base A (3/14);  total (88/118)
Fisher's exact test for biased strand distribution p-value = 3.90e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01

Only 100 of 216 total aligned reads displayed.

ATTCCTTAATCTTATTTATTTGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTG  >  CP009273/3737941‑3738235
                                                                                                                                                   |                                                                                                                                                   
aTTCCTTAATCTTATTTATTTGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCAtt                                                                                                                                                     >  2:2204541/1‑148 (MQ=255)
aTTCCTTAATCTTATTTATTTGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCAtt                                                                                                                                                     <  2:1391508/148‑1 (MQ=255)
aTTCCTTAATCTTATTTATTTGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCAtt                                                                                                                                                     >  1:742214/1‑148 (MQ=255)
aTTCCTTAATCTTATTTATTTGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCAtt                                                                                                                                                     >  1:1361282/1‑148 (MQ=255)
  tCCTTAATCTTATTTATTTGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCt                                                                                                                                                   <  1:3152002/148‑1 (MQ=255)
   ccTTAATCTTATTTATTTGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCtt                                                                                                                                                  <  2:860457/148‑1 (MQ=255)
   ccTTAATCTTATTTATTTGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCtt                                                                                                                                                  <  2:419432/148‑1 (MQ=255)
        aTCTTATTTATTTGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAAc                                                                                                                                             <  2:2305353/148‑1 (MQ=255)
          cTTATTTATTTGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTg                                                                                                                                           <  2:1989364/148‑1 (MQ=255)
          cTTATTTATTTGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTg                                                                                                                                           >  1:1168186/1‑148 (MQ=255)
            tatttattTGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCg                                                                                                                                         <  1:30979/148‑1 (MQ=255)
              tttatttGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTc                                                                                                                                       >  2:1673084/1‑148 (MQ=255)
               ttatttGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCa                                                                                                                                      >  1:3012795/1‑148 (MQ=255)
                tatttGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATACTTTGAACTGCGTCAt                                                                                                                                     >  2:2414105/1‑148 (MQ=255)
                 atttGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCAtt                                                                                                                                    >  2:65413/1‑148 (MQ=255)
                   ttGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTta                                                                                                                                  >  2:1867103/1‑148 (MQ=255)
                     gtgGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATACTTTGAACTGCGTCATTtata                                                                                                                                <  2:2297106/148‑1 (MQ=255)
                      tgGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTatat                                                                                                                               >  1:620231/1‑148 (MQ=255)
                         cgcACGGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCaa                                                                                                                            <  1:1236793/148‑1 (MQ=255)
                             cTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATACTTTGAACTGCGTCATTTATATCAACtca                                                                                                                        <  1:1538602/148‑1 (MQ=255)
                               ggCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACtcatc                                                                                                                      >  2:2951391/1‑148 (MQ=255)
                                     gggCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGt                                                                                                                >  2:2483728/1‑148 (MQ=255)
                                       gcgcAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAg                                                                                                              <  1:4291/148‑1 (MQ=255)
                                          cAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCaa                                                                                                           <  2:328473/148‑1 (MQ=255)
                                          cAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCaa                                                                                                           >  1:1558336/1‑148 (MQ=255)
                                          cAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCaa                                                                                                           <  1:2474499/148‑1 (MQ=255)
                                          cAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCaa                                                                                                           <  1:1733814/148‑1 (MQ=255)
                                          cAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCaa                                                                                                           <  1:1666904/148‑1 (MQ=255)
                                              gCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATACTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATtct                                                                                                       <  2:1219271/148‑1 (MQ=255)
                                                 aaaaaCAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTa                                                                                                    >  1:56608/1‑148 (MQ=255)
                                                  aaaaCAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTaa                                                                                                   <  2:284847/148‑1 (MQ=255)
                                                    aaCAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACg                                                                                                 <  1:2327307/148‑1 (MQ=255)
                                                    aaCAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACg                                                                                                 <  1:811864/148‑1 (MQ=255)
                                                        aTTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGc                                                                                             >  2:1014894/1‑148 (MQ=255)
                                                         ttCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGcc                                                                                            >  1:2851400/1‑148 (MQ=255)
                                                             ggACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCAGTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCg                                                                                        >  2:1009462/1‑148 (MQ=255)
                                                              gACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGc                                                                                       >  2:497476/1‑148 (MQ=255)
                                                               aCTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCa                                                                                      <  2:2904646/148‑1 (MQ=255)
                                                                  ggAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAgg                                                                                   >  1:2546106/1‑148 (MQ=255)
                                                                   gAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGa                                                                                  >  1:2064122/1‑148 (MQ=255)
                                                                   gAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATACTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGa                                                                                  >  2:3129294/1‑148 (MQ=255)
                                                                    aGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAg                                                                                 >  1:514196/1‑148 (MQ=255)
                                                                     gCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGc                                                                                <  2:2064064/148‑1 (MQ=255)
                                                                     gCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGc                                                                                <  2:1493508/148‑1 (MQ=255)
                                                                       gATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCACGGAGCtt                                                                              >  2:1706762/1‑148 (MQ=255)
                                                                        attattCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTa                                                                             <  1:2894651/148‑1 (MQ=255)
                                                                         ttattCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTaa                                                                            <  1:683892/148‑1 (MQ=255)
                                                                         ttattCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATACTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTaa                                                                            <  1:2819018/148‑1 (MQ=255)
                                                                               aCCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATgg                                                                      <  1:3161491/148‑1 (MQ=255)
                                                                                    tCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCtgtg                                                                 >  1:1653847/1‑148 (MQ=255)
                                                                                    tCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCtgtg                                                                 >  1:2639112/1‑148 (MQ=255)
                                                                                    tCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCtgtg                                                                 <  1:2295685/148‑1 (MQ=255)
                                                                                    tCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATACTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCtgtg                                                                 <  2:2748180/148‑1 (MQ=255)
                                                                                    tCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATACTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCtgtg                                                                 <  2:33114/148‑1 (MQ=255)
                                                                                     ccTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGc                                                                >  2:366690/1‑148 (MQ=255)
                                                                                     ccTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGc                                                                >  1:2271239/1‑148 (MQ=255)
                                                                                     ccTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGc                                                                >  1:609981/1‑148 (MQ=255)
                                                                                       tCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTg                                                              <  1:2257649/148‑1 (MQ=255)
                                                                                          gTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGAtt                                                           <  2:1363752/148‑1 (MQ=255)
                                                                                               cccATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCt                                                      >  2:68976/1‑148 (MQ=255)
                                                                                                ccATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTg                                                     >  2:1579043/1‑148 (MQ=255)
                                                                                                ccATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTg                                                     <  2:785886/148‑1 (MQ=255)
                                                                                                ccATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTg                                                     <  2:2085642/148‑1 (MQ=255)
                                                                                                   tCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGc                                                  >  1:771377/1‑148 (MQ=255)
                                                                                                   tCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGc                                                  >  2:858789/1‑148 (MQ=255)
                                                                                                          gATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGctgtctgt                                           >  2:468656/1‑148 (MQ=255)
                                                                                                          gATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGctgtctgt                                           <  2:2222700/148‑1 (MQ=255)
                                                                                                          gATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGctgtctgt                                           <  2:922616/148‑1 (MQ=255)
                                                                                                          gATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATACTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGctgtctgt                                           <  1:2375331/148‑1 (MQ=255)
                                                                                                           aTGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGtt                                          >  1:1090667/1‑148 (MQ=255)
                                                                                                            tGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTg                                         >  2:1373448/1‑148 (MQ=255)
                                                                                                             gTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTgg                                        <  1:912168/148‑1 (MQ=255)
                                                                                                               aTGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGggt                                      >  1:2990315/1‑148 (MQ=255)
                                                                                                                 gCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGggtgg                                    <  2:3012795/148‑1 (MQ=255)
                                                                                                                  ccGCGTGTCTGCCCCCCAGCCTCGTTATCGCATTCTTTTAACTGCGTCATTTATATCAACTCATCACGCGGAGCAATTCTCTAACTTCGCCACCGCCAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGggtggt                                   <  1:3277134/148‑1 (MQ=255)
                                                                                                                   cGCTTGTCTGCCCACCAGCCTCGTTATCGCATACTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTa                                  <  1:2414105/148‑1 (MQ=255)
                                                                                                                    gCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTAt                                 >  1:3312254/1‑148 (MQ=255)
                                                                                                                       tGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGc                              <  2:2546106/148‑1 (MQ=255)
                                                                                                                          cTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTAt                           <  2:3055250/148‑1 (MQ=255)
                                                                                                                              ccaccaGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGt                       <  1:1014894/148‑1 (MQ=255)
                                                                                                                               caccaGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGtt                      <  2:1620185/148‑1 (MQ=255)
                                                                                                                               caccaGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGtt                      >  2:1414115/1‑148 (MQ=255)
                                                                                                                                   aGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTCACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGcc                  <  2:1848586/148‑1 (MQ=255)
                                                                                                                                    gCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCt                 <  2:3104800/148‑1 (MQ=255)
                                                                                                                                     ccTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCTGTTTGCCTa                >  2:3061398/1‑148 (MQ=255)
                                                                                                                                     ccTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTa                <  2:2442808/148‑1 (MQ=255)
                                                                                                                                      cTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTAt               <  1:2726407/148‑1 (MQ=255)
                                                                                                                                        cGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTg             >  1:276605/1‑148 (MQ=255)
                                                                                                                                         gTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGc            <  1:1709445/148‑1 (MQ=255)
                                                                                                                                          ttATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGcc           >  2:323505/1‑148 (MQ=255)
                                                                                                                                           tATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCt          >  1:387640/1‑148 (MQ=255)
                                                                                                                                            aTCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTg         >  1:127300/1‑148 (MQ=255)
                                                                                                                                             tCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCTAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTgg        >  2:2780241/1‑148 (MQ=255)
                                                                                                                                              cGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTggg       <  1:34787/148‑1 (MQ=255)
                                                                                                                                               gCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGc      >  2:2178312/1‑148 (MQ=255)
                                                                                                                                                cATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGcc     >  2:3009611/1‑147 (MQ=255)
                                                                                                                                                cATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCa     <  1:1437804/148‑1 (MQ=255)
                                                                                                                                                cATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCa     >  1:250882/1‑148 (MQ=255)
                                                                                                                                                  ttCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACt   >  2:609855/1‑148 (MQ=255)
                                                                                                                                                   tCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACtg  >  2:1730603/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ATTCCTTAATCTTATTTATTTGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTG  >  CP009273/3737941‑3738235

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: