Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009273 3,962,332 C→T R353C (CGC→TGC)  wecA → UDP‑GlcNAc:undecaprenylphosphate GlcNAc‑1‑phosphate transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,962,3320CT98.6% 233.3 / NA 70R353C (CGC→TGC) wecAUDP‑GlcNAc:undecaprenylphosphate GlcNAc‑1‑phosphate transferase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (30/39);  total (31/39)
Fisher's exact test for biased strand distribution p-value = 4.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTTCCATTGGCGTGCTGGCAGAATATTCTCATTTTGTCCCGGAGTGGGTCATGCTGGTGCTCTTTTTGCTAGCATTCTTCCTCTATGGATATTGCATTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAACGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGATATTCGTGGGTTGTTTCGTACCTTGTGGGCTG  >  CP009273/3962191‑3962475
                                                                                                                                             |                                                                                                                                               
tcTTCCATTGGCGTGTTGGCAGAATATTCTCATTTTCTCCCGGAGTGTGTCATGCTTGTGCTCTTTTTGCTAGCATTCTTCCTCTATGGATATTGCATTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAATGc                                                                                                                                               <  2:398230/143‑1 (MQ=255)
 cTTCCATTGGCGTGCTGGCAGAATATTCTCATTTTGTCCCGGAGTGGGTCATGCTGGTGCTCTTTTTGCTAGCATTCTTCCTCTATGGATATTGCATTAAGCGTGCCTGGACAGTTGCTCGCTGTATTAAGCGCGTAAAATGCa                                                                                                                                              >  2:573103/1‑144 (MQ=255)
     cATTGGCGTGCTGGCAGAATATTCTCATTTTGTCCCGGAGTGGGTCATGCTGGTGCTCTTTTTGCTAGCATTCTTCCTCTATGGATATTGCATTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAATGCAGACTGCGTa                                                                                                                                     <  1:469468/149‑1 (MQ=255)
             tGCTTGCAGAATATTCTCATTTTGTCCCGGAGTGGGGCAGGCTGGTGCTCTTTTTGCTAGCATTCTTCCTCTATGGATATTGCAGTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAATGCAGACTGCGTAGAAATCgt                                                                                                                             <  1:785054/149‑1 (MQ=255)
             tGCTGGCAGAATATTCTCATTTTGTCCCGGAGTGGGTCATGCTGGTGCTCTTTTTGCTAGCATTCTTCCTCTATGGATATTGCATTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAATGCAGACCGCGTAGaa                                                                                                                                  >  2:441202/1‑144 (MQ=255)
              gCTGGCAGAATATTCTCATTTTGTCCCGGAGTGGGTCATGCTGGTGCTCTTTTTGCTAGCATTCTTCCTCTATGGATATTGCATTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAATGCAGACTGCGTAGAAATCgtg                                                                                                                            <  1:965035/149‑1 (MQ=255)
                tGGCAGAATATTCTCATTTTGTCCCGGAGTGGGTCATGCTGGTGCTCTTTTTGCTAGCATTCTTCCTCTATGGATATTGCATTAAGCGTGCCTGGAACGTTGCTCGCTTTATTCAGCGCGTAAAACGAGGACTGCGGAGaaaac                                                                                                                               >  2:447051/1‑142 (MQ=255)
                tGGCAGAATATTCTCATTTTGTCCCGGAGTGGGTCATGCTGGTGCTCTTTTTGCTAGCATTCTTCCTCTATGGATATTGCATTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAATGCAGACTGCGTAGAAATCgtggt                                                                                                                          <  1:723904/149‑1 (MQ=255)
                tGGCAGAATATTCTCATTTTGTCCCGGAGTGGGTCATGCTGGTGCTCTTTCTGCTAGCATTCTTCCTCTATGGATATTGCATTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAATGCAGACTGCGTAGAAATCgtggt                                                                                                                          <  1:370538/149‑1 (MQ=255)
                  gCAGAATATTCTCATTTTGTCCCGGAGTGGGTCATGCTGGTGCTCTTTTTGCTAGCATTCTTCCTCTATGGATATTGCATTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAATGCAGACTGCGTAGAAATCgt                                                                                                                             <  2:368650/144‑1 (MQ=255)
                  gCAGAATATTCTCATTTTGTCCCGGAGTGGGTCATGCTGGTGCTCTTTTTGCTAGCATTCTTCCTCTATGGATATTGCATTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAATGCAGACTGCGTAGAAATCGtggtgg                                                                                                                        <  1:114925/149‑1 (MQ=255)
                    aGAATATTCTCATTTTGTCCCGGAGTGGGTCATGCTGGTGCTCTTTTTGCTAGCATTCTTCCTCTATGGATATTGCATTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGGAAAATGCGGACTGCGTAGGAATCggtg                                                                                                                           >  2:785639/1‑141 (MQ=255)
                     gatttttCTCATTTTGTCCCGGAGTGGGTCATGCTGGTGCTCTTTTTGCTAGCATTCTTCCTCTATGGATATTGCATTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAATGCAGACTGCGTAGAAATCgtggt                                                                                                                          <  2:589509/139‑1 (MQ=255)
                      aaTATTCTCATTTTGTCCCGGAGTGGGTCATGCTGGTGCTCTTTTTGCTAGCATTCTTCCTCTATGGATATTGCATTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAATGCAGACTGCGTAGAAATCGTGGTGGCAGc                                                                                                                    <  1:335288/149‑1 (MQ=255)
                            ctcATTTTGTCCCGGAGTGGGTCATGCTGGTGCTCTTTTTGCTAGCATTCTTCCTCTATGGATATTGCATTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAATGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAAt                                                                                                              <  1:236654/149‑1 (MQ=255)
                              cATTTTGTCCCGGAGTGGGTCATGCTGGTGCTCTTTTTGCTAGCATTCTTCCTCTATGGATATTGCATTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAATGCGGACTGCGTAGAAATCAGGGTGGCAGcccc                                                                                                                 >  2:383627/1‑144 (MQ=255)
                              cATTTTGTCCCGGAGTGGGTCATGCTGGTGCTCTTTTTGCTAGCATTCTTCCTCTATGGATATTGCATTAAGCGTGCCTGGAAAGTCGCTCGCTTTATTAAGCGCGTAAAATGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAAttt                                                                                                            >  1:758771/1‑149 (MQ=255)
                              cATTTTGTCCCGGAGTGGGTCATGCTGGTGCTCTTTTTGCTAGCATTCTTCCTCTATGGATATTGCATTAAGCGTGCCTGGAAAGTAGCTCGCTTTATTAAGCGCGTAAAATGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAAttt                                                                                                            >  1:758747/1‑149 (MQ=255)
                                        cGGAGTGGGTCATGCTGGTGCTCTTTTTGCTAGCATTCTTCCTCTATGGATATTGCATTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAATGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAAccg                                                                                                       >  2:441039/1‑143 (MQ=255)
                                           aGTGGGTCATGCTGGTGCTCTTTTTGCTAGCATTCTTCCTCTATGGATATTGCATTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAATGCAGACTGCGTAGAAATCGTGGTCGCAGCCCCAATTTAACCAAATAAATg                                                                                               >  1:65184/1‑149 (MQ=255)
                                              gcgtgtcactatcgCTCTTTCTGCTTAAAGTCTTCCTCTATGGATATTGCATTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAATGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAgg                                                                                            <  1:532986/136‑1 (MQ=255)
                                                  cATGCTGGTGCTCTTTTTGCTAGCATTCTTCCTCTATGGATATTGCATTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAATGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAg                                                                                             <  2:398906/144‑1 (MQ=255)
                                                   aTGCTGGTGCTCTTTTTGCTAGCATTATTCCGCTATGGATATTGCATTAAGCGCGCCTGGAAAGTTGCTCGCGTTATTAAGCGCGTAAAATGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAgg                                                                                            <  2:572976/144‑1 (MQ=255)
                                                     gCTGGTGCTCTTTTTGCTAGCATTCTTCCTCTATGGATATTGCATTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCCCGTAAAATGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGtgat                                                                                     <  1:1049109/149‑1 (MQ=255)
                                                       tGGTGCTCTTTTTGCTAGCATTCTTCCTCTATGGATATTGCATTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAATGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGtgatga                                                                                   <  1:1023462/149‑1 (MQ=255)
                                                          tGCTCTTTTTTCGAGCATTCGTCCTCTATGGAGATTGCATTAAGCGTGCCTGGAAAGTGGCTTGCTTTATTAAGCGCGTAAAATGCAGACTGCGTAGAAGTCGTGGTGGCAGCCCCAATCTAACCACATAAATGCGGATGagat                                                                                     <  2:386919/144‑1 (MQ=255)
                                                          tGCTCTTTTTGCTAGCATTCTTCCTCTATGGATATTGCATTGCGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAACATGCAGACTGCGTAGAAATCGTGGTCGCAGCCCCACTTTCACCAAATAAATGCGGATGcgat                                                                                     <  2:386958/144‑1 (MQ=255)
                                                          tGCTCTTTTTGCTAGCATTCTTCCTCTATGGATATTGCATTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAATGCAGACTGCGTAGAAATCGGGGTGGCAGCCCCAATTCAAACACATAAATGAGGATTtgat                                                                                     >  2:581759/1‑144 (MQ=255)
                                                              ctTTTTGCTAGCATTCTTCCTCTATTGATATTGCATTATGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAATGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAGTAAAGGAGGATGtgatgcgc                                                                                 >  2:175112/1‑141 (MQ=255)
                                                                 tttGCTAGCATTCTTCCTCTATGGATATTGCATTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAATGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACaa                                                                              <  2:483250/144‑1 (MQ=255)
                                                                     cTAGCATTCTTCCTCTATGGATATTGCATTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAATGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAAGGccg                                                                     <  1:1004552/149‑2 (MQ=255)
                                                                      tAGCATTCTTCCTCTATGGATATTGCATTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAATGCAGACTGCGTAGAAATCGTGGTTGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAAt                                                                         >  2:553205/1‑144 (MQ=255)
                                                                              ttcCTCTATGGATATTGCATTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAATGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCCAATACCAATTGTCGTGATGACACAACCAATGCCTGGGAAACCg                                                            <  1:1108630/149‑1 (MQ=255)
                                                                              ttcCTCTATGGATATTGCATTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAATGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGAAACCg                                                            >  1:564127/1‑149 (MQ=255)
                                                                                ccTCTATGGATATTGCATTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAATGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGaaa                                                               <  2:691084/144‑1 (MQ=255)
                                                                                 ctctATGGATATTGCATTAAGCGTGACTGGAAAGTTGCTCGCTTTATAAAGCGCGTAAAATGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAAGTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGAAACCGGcc                                                         >  1:631220/1‑149 (MQ=255)
                                                                                  tctATGGACATTGCATTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAATGCAGACTGCGTAGCAATCGTGGTGGCCGCCCCAATTCAACCACATCAATGAGGATGTGATGACACAACCAATGCCTGGGAAAcc                                                             <  2:311318/144‑1 (MQ=255)
                                                                                   ctATGGATATTGCATTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAATGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGAAACCGGCCGa                                                       >  1:125568/1‑149 (MQ=255)
                                                                                    tATGGATATTGCATTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAATGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGAAACCGGCCGaa                                                      >  1:940867/1‑149 (MQ=255)
                                                                                     aTGGATATTGCATTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAATGCAGAATGGGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGAGGGGATGACACAACCAATGCCTGGGAAACCGGCCGAAg                                                     >  1:97546/1‑149 (MQ=255)
                                                                                      tGGATATTGCATTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAATGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGAAACCGGCCGAAGa                                                    <  1:244305/149‑1 (MQ=255)
                                                                                      tGGAGATTGCATTAAGCGGGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAATGCAGACGGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCTAATAAATGAGGATGTGATGATACAACCAATGCCTGGGAAACCGGcc                                                         <  2:394183/144‑1 (MQ=255)
                                                                                       ggATATTGCATTAAGCGTGCCTGGAAAATTGCTCGCTTTATTAAGCGCGTAAAATGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATCTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGAAACCGGCCg                                                        <  2:557328/144‑1 (MQ=255)
                                                                                        gATATTGCATTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAATGCAGACTGCGTAGAAATCGAGGTGGCAGCCCAAATTTAAACAAATAAATGAGGCTGTGATGACACACCCCATCCCTGGGCACACGGCCGaccacc                                                  >  1:509490/1‑144 (MQ=255)
                                                                                         atatTGCATTAAGTGTGCCGGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAATGCAGACTGTGTAGAAATCGTGGTTGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGAAACCGGCCGaa                                                      <  2:184324/144‑1 (MQ=255)
                                                                                              gCATTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAATGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGAAACCGGCCGAAGACGc                                                 >  2:537423/1‑144 (MQ=255)
                                                                                              gCATTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAATGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGAAACCGGCCGAAGACGc                                                 >  2:626134/1‑144 (MQ=255)
                                                                                              gCATTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAATGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACAAAGGCCTGGGAAACCGGCCGAAGACGCTGaaa                                            <  1:1109966/149‑1 (MQ=255)
                                                                                               cATTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAATGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGAAACCGGCCGAAGACGCTGaaaa                                           <  1:910090/149‑1 (MQ=255)
                                                                                                   aaGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAATGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGAAATCGGCCGAAGACGCTGAAAATGaa                                       <  1:604434/149‑1 (MQ=255)
                                                                                                      cGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAATGCAGACTGCGTAGAAATCGTGTTGGCAGCTCTAATTTAACCAAATAAATGAGTATGTGATGACCCAACCAATGCCTGGGAAACCGGCCGAAGACGCTGAAAATGAACTg                                    <  1:47601/149‑1 (MQ=255)
                                                                                                      cGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAATGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGAAACCGGCCGAAGACGCTGAAAATg                                         <  2:449409/144‑1 (MQ=255)
                                                                                                           cTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAATGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGatat                               <  1:177753/149‑1 (MQ=255)
                                                                                                            tGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAATGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGAAACCGGCCGAAGACGCTGAAAATGAACTgg                                   >  2:417228/1‑144 (MQ=255)
                                                                                                                aaGTTGCTCGCTTTATTAAGCGCGTAAAATGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGATATTCGTg                          <  1:1063148/149‑1 (MQ=255)
                                                                                                                   ttGCTCGCTTTATTAAGCGCGTAAAATGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGATATTCg                            >  2:331086/1‑144 (MQ=255)
                                                                                                                    tgtttccTTTATTAAGCGCGTAAAATGCAGCGTGCGTAGAAATCGTGGTGGCAGCCCCAATTGAACCAAATAAATGTGGACGTGATGACACAACCAATGCCTGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGATATTCGt                           <  2:96071/138‑1 (MQ=255)
                                                                                                                    tGCTCGCTTTATTAAGCGCGTAAAATGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGATATTCGt                           <  2:320100/144‑1 (MQ=255)
                                                                                                                    tGCTCGCTTTATTAAGCGCGTAAAATGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGATATTCGt                           <  2:688053/144‑1 (MQ=255)
                                                                                                                    tGCTCGCTTTATTAAGCGCGTAAAATGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGATATTCGt                           <  2:688076/144‑1 (MQ=255)
                                                                                                                       tCGCTTTATTAAGCGCGTAAAATGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATGTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGATATTCGTggg                        >  2:821590/1‑144 (MQ=255)
                                                                                                                        cGCTTTATTAAGCGCGTAAAATGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGATATTCGTGGGTTGttt                  <  1:421878/149‑1 (MQ=255)
                                                                                                                          cTTTATTAAGCGCGTAAAATGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGATATTCGTGGGTTGGTTCg                >  1:163690/1‑149 (MQ=255)
                                                                                                                           tttATTAAGCGCGTAAAATGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGATATTCGTGGgttgt                    >  2:323866/1‑144 (MQ=255)
                                                                                                                                taAGCGCGTAAAATGCAGACTGCGTAGAAATCGTGGTTGCAGCACCAATTTAAAAAAAGAAAACAGGATGTGATGACACAACCAATGCCTGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGATATTCGTGGGTTGTTTCGTAccgt          >  1:713489/1‑147 (MQ=255)
                                                                                                                                   gcgcGTAAAATGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGATATTCGTGGGTTGTTTCGTACCTtgtg       >  1:980993/1‑149 (MQ=255)
                                                                                                                                     gcgTAAAATGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGATATTCGTGGGTTGTTTCGTACCtt          >  2:397578/1‑144 (MQ=255)
                                                                                                                                     gcgTAAAATGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGATATTCGTGGGTTGTTTCGTACCtt          >  2:961300/1‑144 (MQ=255)
                                                                                                                                        tAAAATGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATACCTGAGGATGTGATGACACAACCAATGCCTGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGATATTCGTGGGTTGTTTCGTACCTTGTGGGCTg  >  1:497070/1‑149 (MQ=255)
                                                                                                                                        tAAAATGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGATAGTCGTGGGTTGTTTCGTACCTTGTGGGCTg  >  1:603385/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                               
GCTTCCATTGGCGTGCTGGCAGAATATTCTCATTTTGTCCCGGAGTGGGTCATGCTGGTGCTCTTTTTGCTAGCATTCTTCCTCTATGGATATTGCATTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAACGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGATATTCGTGGGTTGTTTCGTACCTTGTGGGCTG  >  CP009273/3962191‑3962475

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: