Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,194,733 G→T 10.6% V56V (GTG→GTT ygiM → SH3 domain protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,194,7330GT10.6% 1151.8 / 76.9 313V56V (GTG→GTTygiMSH3 domain protein
Reads supporting (aligned to +/- strand):  ref base G (164/115);  new base T (18/16);  total (182/131)
Fisher's exact test for biased strand distribution p-value = 5.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 317 total aligned reads displayed.

TCGGATTAACTTTACTCGCACTTAGCGCGACTGCCGTCTCACACGCTGAAGAAACGCGCTATGTTTCCGACGAACTGAATACCTGGGTCCGCAGCGGTCCGGGAGATCATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGAGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCC  >  CP009273/3194585‑3194878
                                                                                                                                                    |                                                                                                                                                 
tCGGATTAACTTTACTCGCACTTAGCGCGACTGCCGTCTCACACGCTGAAGAAACGCGCTATGTTTCCGACGAACTGAATACCTGGGTCCGCAGCGGTCCGGGAGATCATTATCGCCTCGTGGGCACGGGTCACGCCGGCGAGGACGTg                                                                                                                                                   >  2:3650184/1‑149 (MQ=255)
  ggATTAACTTTACTCGCACTTAGCGCGACTGCCGTCTCACACGCTGAAGAAACGCGCTATGTTTCCGACGAACTGAATACCTGGGTCCGCAGCGGTCCGGGAGATCATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTGAc                                                                                                                                                 <  1:1721129/149‑1 (MQ=255)
   gATTAACTTTACTCGCACTTAGCGCGACTGCCGTCTCACACGCTGAAGAAACGCGCTATGTTTCCGACGAACTGAATACCTGGGTCCGCAGCGGTCCGGGAGATCATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTTAcc                                                                                                                                                <  1:1466368/149‑1 (MQ=255)
   gATTAACTTTACTCGCACTTAGCGCGACTGCCGTCTCACACGCTGAAGAAACGCGCTATGTTTCCGACGAACTGAATACCTGGGTCCGCAGCGGTCCGGGAGATCATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTGAcc                                                                                                                                                >  1:1072935/1‑149 (MQ=255)
   gATTAACTTTACTCGCACTTAGCGCGACTGCCGTCTCACACGCTGAAGAAACGCGCTATGTTTCCGACGAACTGAATACCTGGGTCCGCAGCGGTCCGGGAGATCATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTGAcc                                                                                                                                                >  1:4916302/1‑149 (MQ=255)
    aTTAACTTTACTCGCACTTAGCGCGACTGCCGTCTCACACGCTGAAGAAACGCGCTATGTTTCCGACGAACTGAATACCTGGGTCCGCAGCGGTCCGGGAGATCATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTGACCt                                                                                                                                               >  2:1093460/1‑149 (MQ=255)
         cTTTACTCGCACTTAGCGCGACTGCCGTCTCACACGCTGAAGAAACGCGCTATGTTTCCGACGAACTGAATACCTGGGTCCGCAGCGGTCCGGGAGATCATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTGACCttatta                                                                                                                                          >  1:1358954/1‑149 (MQ=255)
            tACTCGCACTTAGCGCGACTGCCGTCTCACACGCTGAAGAAACGCGCTATGTTTCCGACGAACTGAATACCTGGGTCCGCAGCGGTCCGGGAGATCATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTGACCTTATTACaa                                                                                                                                       >  1:3982396/1‑149 (MQ=255)
             aCTCGCACTTAGCGCGACTGCCGTCTCACACGCTGAAGAAACGCGCTATGTTTCCGACGAACTGAATACCTGGGTCCGCAGCGGTCCGGGAGATCATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTGACCTTATTACaaa                                                                                                                                      >  1:2458307/1‑149 (MQ=255)
              cTCGCACTTAGCGCGACTGCCGTCTCACACGCTGAAGAAACGCGCTATGTTTCCGACGAACTGAATACCTGGGTCCGCAGCGGTCCGGGAGATCATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTGACCTTATTACAAAc                                                                                                                                     <  1:1406370/149‑1 (MQ=255)
              cTCGCACTTAGCGCGACTGCCGTCTCACACGCTGAAGAAACGCGCTATGTTTCCGACGAACTGAATACCTGGGTCCGCAGCGGTCCGGGAGATCATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTGACCTTATTACAAAc                                                                                                                                     <  1:2004721/149‑1 (MQ=255)
               tCGCACTTAGCGCGACTGCCGTCTCACACGCTGAAGAAACGCGCTATGTTTCCGACGAACTGAATACCTGGGTCCGCAGCGGTCCGGGAGATCATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTGACCTTATTACAAACt                                                                                                                                    >  1:869955/1‑149 (MQ=255)
                 gCACTTAGCGCGACTGCCGTCTCACACGCTGAAGAAACGCGCTATGTTTCCGACGAACTGAATACCTGGGTCCGCAGCGGTCCGGGAGATCATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTTACCTTATTACAAACTGa                                                                                                                                  <  2:1004058/149‑1 (MQ=255)
                 gCACTTAGCGCGACTGCCGTCTCACACGCTGAAGAAACGCGCTATGTTTCCGACGAACTGAATACCTGGGTCCGCAGCGGTCCGGGAGATCATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTGACCTTATTACAAACTGa                                                                                                                                  <  1:2192015/149‑1 (MQ=255)
                 gCACTTAGCGCGACTGCCGTCTCACACGCTGAAGAAACGCGCGATGTTTCCGACGAACTGAATACCTGGGTCCGCAGCGGTCCGGGAGATCATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTTACCTTATTACAAACTGa                                                                                                                                  <  1:273880/149‑1 (MQ=255)
                        gcgcGACTGCCGTCTCACACGCTGAAGAAACGCGCTATGTTTCCGACGAACTGAATACCTGGGTCCGCAGCGGTCCGGGAGATCATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTGACCTTATTACAAACTGACGCCAac                                                                                                                           >  2:4249350/1‑149 (MQ=255)
                           cgACTGCCGTCTCACACGCTGAAGAAACGCGCTATGTTTCCGACGAACTGAATACCTGGGTCCGCAGCGGTCCGGGAGATCATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTGACCTTATTACAAACTGACGCCAACAcc                                                                                                                        <  2:2097668/149‑1 (MQ=255)
                             aCTGCCGTCTCACACGCTGAAGAAACGCGCTATGTTTCCGACGAACTGAATACCTGGGTCCGCAGCGGTCCGGGAGATCATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTTACCTTATTACAAACTGACGCCAACACCaa                                                                                                                      <  2:1281635/149‑1 (MQ=255)
                               tGCCGTCTCACACGCTGAAGAAACGCGCTATGTTTCCGACGAACTGAATACCTGGGTCCGCAGCGGTCCGGGAGATCATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAAtt                                                                                                                    <  1:5020021/149‑1 (MQ=255)
                               tGCCGTCTCACACGCTGAAGAAACGCGCTATGTTTCCGACGAACTGAATACCTGGGTCCGCAGCGGTCCGGGAGATCATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAAtt                                                                                                                    <  1:4724196/149‑1 (MQ=255)
                                  cGTCTCACACGCTGAAGAAACGCGCTATGTTTCCGACGAACTGAATACCTGGGTCCGCAGCGGTCCGGGAGATCATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAATTATg                                                                                                                 <  1:3110665/149‑1 (MQ=255)
                                  cGTCTCACACGCTGAAGAAACGCGCTATGTTTCCGACGAACTGAATACCTGGGTCCGCAGCGGTCCGGGAGATCATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAATTATg                                                                                                                 >  1:3904469/1‑149 (MQ=255)
                                   gTCTCACACGCTGAAGAAACGCGCTATGTTTCCGACGAACTGAATACCTGGGTCCGCAGCGGTCCGGGAGATCATTATCGCCTCGTGGGCCCGGTTAACGCCGGCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAATTATGc                                                                                                                >  2:4704670/1‑149 (MQ=255)
                                     ctcACACGCTGAAGAAACGCGCTATGTTTCCGACGAACTGAATACCTGGGTCCGCAGCGGTCCGGGAGATCATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAATTATGccc                                                                                                              <  1:771610/149‑1 (MQ=255)
                                     ctcACACGCTGAAGAAACGCGCTATGTTTCCGACGAACTGAATACCTGGGTCCGCAGCGGTCCGGGAGATCATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAATTATGccc                                                                                                              <  1:5243466/149‑1 (MQ=255)
                                      tcACACGCTGAAGAAACGCGCTATGTTTCCGACGAACTGAATACCTGGGTCCGCAGCGGTCCGGGAGATCATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAATTATGCCCa                                                                                                             >  1:4132539/1‑149 (MQ=255)
                                      tcACACGCTGAAGAAACGCGCTATGTTTCCGACGAACTGAATACCTGGGTCCGCAGCGGTCCGGGAGATCATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAATTATGCCCa                                                                                                             <  1:351172/149‑1 (MQ=255)
                                       cacaCGCTGAAGAAACGCGCTATGTTTCCGACGAACTGAATACCTGGGTCCGCAGCGGTCCGGGAGATCATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAATTATGCCCAg                                                                                                            >  2:970710/1‑149 (MQ=255)
                                       cacaCGCTGAAGAAACGCGCTATGTTTCCGACGAACTGAATACCTGGGTCCGCAGCGGTCCGGGAGATCATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAATTATGCCCAg                                                                                                            >  2:4651000/1‑149 (MQ=255)
                                         cacGCTGAAGAAACGCGCTATGTTTCCGACGAACTGAATACCTGGGTCCGCAGCGGTCCGGGAGATCATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGt                                                                                                          <  1:113871/149‑1 (MQ=255)
                                            gCTGAAGAAACGCGCTATGTTTCCGACGAACTGAATACCTGGGTCCGCAGCGGTCCGGGAGATCATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGaa                                                                                                       >  2:1468867/1‑149 (MQ=255)
                                                aagaaACGCGCTATGTTTCCGACGAACTGAATACCTGGGTCCGCAGCGGTCCGGGAGATCATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGAc                                                                                                   <  2:1744632/149‑1 (MQ=255)
                                                aagaaACGCGCTATGTTTCCGACGAACTGAATACCTGGGTCCGCAGCGGTCCGGGAGATCATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGAc                                                                                                   <  1:1265380/149‑1 (MQ=255)
                                                 agaaACGCGCTATGTTTCCGACGAACTGAATACCTGGGTCCGCAGCGGTCCGGGAGATCATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACa                                                                                                  <  1:3594509/149‑1 (MQ=255)
                                                  gaaACGCGCTATGTTTCCGACGAACTGAATACCTGGGTCCGCAGCGGTCCGGGAGATCATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAg                                                                                                 >  2:3343706/1‑149 (MQ=255)
                                                     acgctatATGTTTCCGACGAACTGAATACCTGGTTCCGCAGCGGTCCGGGAGATCATTATCGACTCGTGGGCACGGTTAACGCCGGCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGctc                                                                                              >  1:3059507/7‑149 (MQ=255)
                                                     aCGCGCTATGTTTCCGACGAACTGAATACCTGGGTCCGCAGCGGTCCGGGAGATCATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGctc                                                                                              <  2:4093909/149‑1 (MQ=255)
                                                     aCGCGCTATGTTTCCGACGAACTGAATACCTGGGTCCGCAGCGGTCCGGGAGATCATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGctc                                                                                              >  1:435538/1‑149 (MQ=255)
                                                       gcgcTATGTTTCCGACGAACTGAATACCTGGGTCCGCAGCGGTCCGGGAGATCATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTg                                                                                            <  1:1670531/149‑1 (MQ=255)
                                                       gcgcTATGTTTCCGACGAACTGAATACCTGGGTCCGCAGCGGTCCGGGAGATCATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTg                                                                                            <  2:455827/149‑1 (MQ=255)
                                                           tATGTTTCCGACGAACTGAATACCTGGGTCCGCAGCGGTCCGGGAGATCATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTTACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCg                                                                                        >  2:4440248/1‑149 (MQ=255)
                                                           tATGTTTCCGACGAACTGAATACCTGGGTCCGCAGCGGTCCGGGAGATCATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCg                                                                                        >  2:2704281/1‑149 (MQ=255)
                                                             tGTTTCCGACGAACTGAATACCTGGGTCCGCAGCGGTCCGGGAGATCATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTa                                                                                      >  2:1237598/1‑149 (MQ=255)
                                                                 tCCGACGAACTGAATACCTGGGTCCGCAGCGGTCCGGGAGATCATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGc                                                                                  <  2:4479957/149‑1 (MQ=255)
                                                                    gacgaACTGAATACCTGGGTCCGCAGCGGTCCGGGAGATCATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTg                                                                               >  2:4788375/1‑149 (MQ=255)
                                                                     acgaACTGAATACCTGGGTCCGCAGCGGTCCGGGAGATCATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTgg                                                                              >  1:4625924/1‑149 (MQ=255)
                                                                     acgaACTGAATACCTGGGTCCGCAGCGGTCCGGGAGATCATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTgg                                                                              <  2:3186560/149‑1 (MQ=255)
                                                                      cgaACTGAATACCTGGGTCCGCAGCGGTCCGGGAGATCATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGa                                                                             >  2:5072925/1‑149 (MQ=255)
                                                                       gaACTGAATACCTGGGTCCGCAGCGGTCCGGGAGATCATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGAt                                                                            <  1:4390706/149‑1 (MQ=255)
                                                                       gaACTGAATACCTGGGTCCGCAGCGGTCCGGGAGATCATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGAt                                                                            >  1:1282558/1‑149 (MQ=255)
                                                                        aaCTGAATACCTGGGTCCGCAGCGGTCCGGGAGATCATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATc                                                                           >  2:4456485/1‑149 (MQ=255)
                                                                         aCTGAATACCTGGGTCCGCAGCGGTCCGGGAGATCATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTTACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATcc                                                                          <  1:4304261/149‑1 (MQ=255)
                                                                         aCTGAATACCTGGGTCCGCAGCGGTCCGGGAGATCATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATcc                                                                          >  2:2789843/1‑149 (MQ=255)
                                                                         aCTGAATACCTGGGTCCGCAGCGGTCCGGGAGATCATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATcc                                                                          >  2:4402269/1‑149 (MQ=255)
                                                                            gAATACCTGGGTCCGCAGCGGTCCGGGAGATCATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGt                                                                       >  1:3232691/1‑149 (MQ=255)
                                                                            gAATACCTGGGTCCGCAGCGGTCCGGGAGATCATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGt                                                                       >  1:4074913/1‑149 (MQ=255)
                                                                              ctACCTGGGTCCGCAGCGGTCCGGGAGATCATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTg                                                                     >  2:50838/2‑149 (MQ=255)
                                                                              aTACCTGGGTCCGCAGCGGTCCGGGAGATCATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTg                                                                     >  1:4172258/1‑149 (MQ=255)
                                                                              aTACCTGGGTCCGCAGCGGTCCGGGAGATCATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTg                                                                     <  1:3317502/149‑1 (MQ=255)
                                                                               tACCTGGGTCCGCAGCGGTCCGGGAGATCATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGa                                                                    >  1:1033834/1‑149 (MQ=255)
                                                                                 ccTGGGTCCGCAGCGGTCCGGGAGATCATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGaaa                                                                  >  2:1872155/1‑149 (MQ=255)
                                                                                        ccGCAGCGGTCCGGGAGATCATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTa                                                           >  1:4683905/1‑149 (MQ=255)
                                                                                          gCAGCGGTCCGGGAGATCATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGc                                                         <  1:2317011/149‑1 (MQ=255)
                                                                                          gCAGCGGTCCGGGAGATCATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGc                                                         >  1:120226/1‑149 (MQ=255)
                                                                                          gCAGCGGTCCGGGAGATCATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGc                                                         <  1:2880384/149‑1 (MQ=255)
                                                                                                 tCCGGGAGATCATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGAGc                                                  >  1:3772641/1‑149 (MQ=255)
                                                                                                     ggAGATCATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGAGCCAAg                                              >  2:361499/1‑149 (MQ=255)
                                                                                                      gagaTCATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGAGCCAAGc                                             <  2:4132539/149‑1 (MQ=255)
                                                                                                       agaTCATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTTACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGAGCCAAGcc                                            <  1:2117217/149‑1 (MQ=255)
                                                                                                       agaTCATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGAGCCAAGcc                                            >  2:2225367/1‑149 (MQ=255)
                                                                                                       agaTCATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGAGCCAAGcc                                            >  1:2540902/1‑149 (MQ=255)
                                                                                                       agaTCATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGAGCCAAGcc                                            >  2:1528340/1‑149 (MQ=255)
                                                                                                       agaTCATTATCGCCTCGTGGGCACGGTTAACGCCAGCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGAGCCAAGcc                                            >  2:733532/1‑149 (MQ=255)
                                                                                                         aTCATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGAGCCAAGCCTg                                          >  1:4047070/1‑149 (MQ=255)
                                                                                                         aTCATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGAGCCAAGCCTg                                          <  1:5100492/149‑1 (MQ=255)
                                                                                                           cATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTTACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGAGCCAAGCCTgcg                                        >  2:4924617/1‑149 (MQ=255)
                                                                                                                tCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGAGCCAAGCCTGCGCTccc                                   <  2:2526748/149‑1 (MQ=255)
                                                                                                                tCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGAGCCAAGCCTGCGCTccc                                   <  1:3396796/149‑1 (MQ=255)
                                                                                                                tCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGAGCCAAGCCTGCGCTccc                                   >  2:4188525/1‑149 (MQ=255)
                                                                                                                 cGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGAGCCAAGCCTGCGCTCCCg                                  >  2:3681926/1‑149 (MQ=255)
                                                                                                                   ccTCGTGGGCACGGTTAACGCCGGCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGAGCCAAGCCTGCGCTCCCgtg                                <  2:3541087/149‑1 (MQ=255)
                                                                                                                    cTCGTGGGCACGGTTAACGCCGGCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGAGCCAAGCCTGCGCTCCCgtgt                               >  1:2873172/1‑149 (MQ=255)
                                                                                                                       gTGGGCACGGTTAACGCCGGCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGAGCCAAGCCTGCGCTCCCGTGTGcc                            <  2:1144075/149‑1 (MQ=255)
                                                                                                                        tGGGCACGGTTAACGCCGGCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGAGCCAAGCCTGCGCTCCCGTGTGCCa                           >  2:3471714/1‑149 (MQ=255)
                                                                                                                          ggCACGGTTAACGCCGGCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGAGCCAAGCCTGCGCTCCCGTGTGCCAGa                         >  2:4153478/1‑149 (MQ=255)
                                                                                                                           gCACGGTTAACGCCGGCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGAGCCAAGCCTGCGCTCCCGTGTGCCAGAt                        <  2:3904469/149‑1 (MQ=255)
                                                                                                                           gCACGGTTAACGCCGGCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGAGCCAAGCCTGCGCTCCCGTGTGCCAGAt                        >  1:5158977/1‑149 (MQ=255)
                                                                                                                               ggTTAACGCCGGCGAGGAAGTTACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGAGCCAAGCCTGCGCTCCCGTGTGCCAGATCTgg                    <  2:2383502/149‑1 (MQ=255)
                                                                                                                                gTTAACGCCGGCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGAGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGa                   <  1:2934765/149‑1 (MQ=255)
                                                                                                                                  tAACGCCGGCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGAGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGaaa                 <  1:810894/149‑1 (MQ=255)
                                                                                                                                  tAACGCCGGCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGAGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGaaa                 >  2:2997450/1‑149 (MQ=255)
                                                                                                                                   aaCGCCGGCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGAGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGaaaa                >  2:5043367/1‑149 (MQ=255)
                                                                                                                                       ccGGCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGAGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAg            <  2:4822980/149‑1 (MQ=255)
                                                                                                                                        cGGCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGAGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAgg           <  1:1310142/149‑1 (MQ=255)
                                                                                                                                         ggCGAGGAAGTTACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGAGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGt          >  2:1661232/1‑149 (MQ=255)
                                                                                                                                         ggCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGAGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGt          <  1:2109653/149‑1 (MQ=255)
                                                                                                                                               gAAGTTACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGAGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAAc    >  1:3928219/1‑149 (MQ=255)
                                                                                                                                               gAAGTGACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGAGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAAc    <  2:575454/149‑1 (MQ=255)
                                                                                                                                                 aGTGACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGAGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAAccc  >  2:1441177/1‑149 (MQ=255)
                                                                                                                                                 aGTGACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGAGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAAccc  <  2:5161282/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
TCGGATTAACTTTACTCGCACTTAGCGCGACTGCCGTCTCACACGCTGAAGAAACGCGCTATGTTTCCGACGAACTGAATACCTGGGTCCGCAGCGGTCCGGGAGATCATTATCGCCTCGTGGGCACGGTTAACGCCGGCGAGGAAGTGACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGAGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCC  >  CP009273/3194585‑3194878

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: