Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,020,029 G→T 9.4% E35* (GAA→TAA)  helD → DNA helicase IV

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,020,0290GT 9.4% 1013.6 / 50.5 284E35* (GAA→TAA) helDDNA helicase IV
Reads supporting (aligned to +/- strand):  ref base G (121/136);  new base T (16/11);  total (137/147)
Fisher's exact test for biased strand distribution p-value = 3.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 291 total aligned reads displayed.

GACGGTTACACTGATAAGAAATAATTTCGTGTGGATCTACAGAGTCATGGAACTGAAAGCGACAACGCTTGGAAAACGTCTGGCACAGCACCCTTACGATCGGGCGGTGATCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATACCTTATTCCTTTCAATCAATTACTGGCGATTCACTGTAAGCGCGGTCTGGTATGGGGCGAGCTGGAATTTGTACTGCCGGACGAAAAAGTGGTGCGTCTGCACGGCACCGAATGGGGCGAGACGCAGCGTTTTTACCATCATCT  >  CP009273/1019881‑1020177
                                                                                                                                                    |                                                                                                                                                    
 aCGGTTACACTGATAAGAAATAATTTCGTGTGGATCTACAGAGTCATGGAACTGAAAGCGACAACGCTTGGAAAACGTCTGGCACAGCACCCTTACGATCGGGCGGTGATCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACGa                                                                                                                                                     >  1:3217095/1‑149 (MQ=255)
   ggTTACACTTATAAGAAATAATTTCGTGTGGATCTACAGAGTCATGGAACTGAAAGCGACAACGCTTGGAAAACGTCTGGCACAGCACCCTTACGATCGGGCGGTGATCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACTAAt                                                                                                                                                   <  1:3656608/149‑1 (MQ=255)
    gTTACACTGATAAGAAATAATTTCGTGTGGATCTACAGAGTCATGGAACTGAAAGCGACAACGCTTGGAAAACGTCTGGCACAGCACCCTTACGATCGGGCGGTGATCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATa                                                                                                                                                  >  1:197327/1‑149 (MQ=255)
    gTTACACTGATAAGAAATAATTTCGTGTGGATCTACAGAGTCATGGAACTGAAAGCGACAACGCTTGGAAAACGTCTGGCACAGCACCCTTACGATCGGGCGGTGATCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATa                                                                                                                                                  <  1:714997/149‑1 (MQ=255)
    gTTACACTGATAAGAAATAATTTCGTGTGGATCTACAGAGTCATGGAACTGAAAGCGACAACGCTTGGAAAACGTCTGGCACAGCACCCTTACGATCGGGCGGTGATCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATa                                                                                                                                                  >  1:1569357/1‑149 (MQ=255)
     ttACACTTATAAGAAATAATTTCGTGTGGATCTAATGAGTCTTGGAACTGAAAGCGACAACGCTTGGAAAACGTCTGGCACAGCACCCTTACGATCGGGCGGTGATCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATAc                                                                                                                                                 <  1:1804751/149‑1 (MQ=255)
     ttACACTGATAAGAAATAATTTCGTGTGGATCTACAGAGTCATGGAACTGAAAGCGACAACGCTTGGAAAACGTCTGGCACAGCACCCTTACGATCGGGCGGTGATCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATAc                                                                                                                                                 >  2:1043288/1‑149 (MQ=255)
        cacTGATAAGAAATAATTTCGTGTGGATCTACAGAGTCATGGAACTGAAAGCGACAACGCTTGGAAAACGTCTGGCACAGCACCCTTACGATCGGGCGGTGATCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATACCtt                                                                                                                                              >  2:1507307/1‑149 (MQ=255)
        cacTGATAAGAAATAATTTCGTGTGGATCTACAGAGTCATGGAACAGAAAGCGACAACGCTTGGAAAACGTCTGGCACAGCACCCTTACGATCGGGCGGTGATCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATACCtt                                                                                                                                              >  2:3896142/1‑149 (MQ=255)
          cttaTAAGAAATAATTTCGTGTGGATCTACAGAGTCATGGAACTGAAAGCGACAACGCTTGGAAAACGTCTGGCACAGCACCCTTACGATCGGGCGGTGATCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATACCTTAt                                                                                                                                            <  1:2000708/146‑1 (MQ=255)
            gATAAGAATTAATTTCGTGTGGATCTACAGAGTCATGGAACTGAAAGCGACAACGCTTGGAAAACGTCTGGCACAGCACCCTTACGATCGGGCGGTGATCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATACCTTATTc                                                                                                                                          <  1:775362/149‑1 (MQ=255)
            gATAAGAAATAATTTCTTGTGGATCTACAGAGTCATGGAACTGAAAGCGACAACGCTTGGAAAACGTCTGGCACAGCACCCTTACGATCGGGCGGTGATCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATACCTTATTc                                                                                                                                          <  2:4763925/149‑1 (MQ=255)
            gATAAGAAATAATTTCGTGTGGATCTACAGAGTCATGGAACTGAAAGCGACAACGCTTGGAAAACGTCTGGCACAGCACCCTTACGATCGGGCGGTGATCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACGACTACCTTATTc                                                                                                                                          >  2:3456563/1‑149 (MQ=255)
             ataataATTATTTTCGTGTGTATCTACAGTGTCTTGGAACTGAAAGCGACAACGCTTGGAAAACGTCTGGCACAGCACCCTTACGATCGGGCGGTGATCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATACCTTATTcc                                                                                                                                         <  1:842180/149‑1 (MQ=255)
             aTAAGAAATAATTTCGTGTGGATCTACAGAGTCATGGAACTGAAAGCGACAACGCTTGGAAAACGTCTGGCACAGCACCCTTACGATCGGGCGGTGATCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATACCTTATTcc                                                                                                                                         <  1:1582777/149‑1 (MQ=255)
               aaGAAATAATTTCGTGTGGATCTACAGAGTCATGGAACTGAAAGCGACAACGCTTGGAAAACGTCTGGCACAGCACCCTTACGATCGGGCGGTGATCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCCCGAATACCTTATTCCtt                                                                                                                                       >  2:1073938/1‑149 (MQ=255)
               aaGAAATAATTTCGTGTGGATCTACAGAGTCATGGAACTGAAAGCGACAACGCTTGGAAAACGTCTGGCACAGCACCCTTACGATCGGGCGGTGATCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATACCTTATTCCtt                                                                                                                                       <  2:78601/149‑1 (MQ=255)
                aGAAATAATTTCGTGTGGATCTACAGAGTCATGGAACTGAAAGCGACAACGCTTGGAAAACGTCTGGCACAGCACCCTTACGATCGGGCGGTGATCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATACCTTATTCCttt                                                                                                                                      >  1:1473668/1‑149 (MQ=255)
                aGAAATAATTTCGTGTGGATCTACAGAGTCATGGAACTGAAAGCGACAACGCTTGGAAAACGTCTGGCACAGCACCCTTACGATCGGGCGGTGATCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATACCTTATTCCttt                                                                                                                                      <  2:2682041/149‑1 (MQ=255)
                 gAAATAATTTCGTGTGGATCTACAGAGTCATGGAACTGAAAGCGACAACGCTTGGAAAACGTCTGGCACAGCACCCTTACGATCGGGCGGTGATCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATACCTTATTCCTTtc                                                                                                                                     <  1:2081835/149‑1 (MQ=255)
                   aataatTTCGTGTGGATCTACAGAGTCATGGAACTGAAAGCGACAACGCTTGGAAAACGTCTGGCACAGCACCCTTACGATCGGGCGGTGATCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATACCTTATTCCTTtcaa                                                                                                                                   <  1:1200079/149‑1 (MQ=255)
                   aataatTTCGTGTGGATCTACAGAGTCATGGAACTGAAAGCGACAACGCTTGGAAAACGTCTGGCACAGCACCCTTACGATCGGGCGGTGATCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATACCTTATTCCTTtcaa                                                                                                                                   <  2:4779799/149‑1 (MQ=255)
                     tattttCGTGTGGATCTACAGAGTATTGGAACTGAAAGCGACAACGCTTGGAAAACGTCTGGCACAGCACCCTTACGATCGGGCGGTGATCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATACCTTATTCCTTtcaatc                                                                                                                                 <  1:7955/146‑1 (MQ=255)
                     taatTTCGTGTGGATCTAAATAGTCATGGAACTGAAAGCGACAACGCTTGGAAAACGTCTGGCACAGCACCCTTACGATCGGGCGGTGATCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATACCTTATTCCTTtcaatc                                                                                                                                 <  2:485986/149‑1 (MQ=255)
                             tgtgGATCTACAGAGTCATGGAACTGAAAGCGACAACGCTTGGAAAACGTCTGGCACAGCACCCTTACGATCGGGCGGTGATCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATACCTTATTCCTTTCAATCAATTACTg                                                                                                                         >  2:4282981/1‑149 (MQ=255)
                                ggATCTACAGAGTCATGGAACTGAAAGCGACAACGCTTGGAAAACGTCTGGCACAGCACCCTTACGATCGGGCGGTGATCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATACCTTATTCCTTTCAATCAATTACTGGCg                                                                                                                      >  2:451050/1‑149 (MQ=255)
                                  aTCTACAGAGTCATGGAACTGAAAGCGACAACGCTTGGAAAACGTCTGGCACAGCACCCTTACGATCGGGCGGTGATCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATACCTTATTCCTTTCAATCAATTACTGGCGAt                                                                                                                    <  1:2694199/149‑1 (MQ=255)
                                    cTACAGAGTCATGGAACTGAAAGCGACAACGCTTGGAAAACGTCTGGCACAGCACCCTTACGATCGGGCGGTGATCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATACCTTATTCCTTTCAATCAATTACTGGCGATTc                                                                                                                  <  2:2445723/149‑1 (MQ=255)
                                       cAGAGTCATGTAACTTAAAGCGACAACGCTTGGAAAACGTCTGGCACAGCACCCTTACGATCGGGCGGTGATCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACTAATACCTTATTCCTTTCAATCAATTACTGGCGATTCACt                                                                                                               <  1:1842935/149‑1 (MQ=255)
                                             cATGGAACTGAAAGCGACAACGCTTGGAAAACGTCTGGCACAGCACCCTTACGATCGGGCGGTGATCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATACCTTATTCCTTTCAATCAATTACTGGCGATTCACTGTAAgc                                                                                                         <  2:2331164/149‑1 (MQ=255)
                                               tGGAACTGAAAGCGACAACGCTTGGAAAACGTCTGGCACAGCACCCTTACGATCGGGCGGTGATCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACTAATACCTTATTCCTTTCAATCAATTACTGGCGATTCACTGTAAgcgc                                                                                                       >  2:2120005/1‑149 (MQ=255)
                                                   aCTGAAAGCGACAACGCTTGGAAAACGTCTGGCACAGCACCCTTACGATCGGGCGGTGATCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATACCTTATTCCTTTCAATCAATTACTGGCGATTCACTGTAAGCGCGGTc                                                                                                   >  2:3728009/1‑149 (MQ=255)
                                                    cTGAAAGCGACAACGCTTGGAAAACGTCTGGCACAGCACCCTTACGATCGGGCGGTGATCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATACCTTATTCCTTTCAATCAATTACTGGCGATTCACTGTAAGCGCGGTCt                                                                                                  >  2:3177354/1‑149 (MQ=255)
                                                    cTGAAAGCGACAACGCTTGGAAAACGTCTGGCACAGCACCCTTACGATCGGGCGGTGATCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATACCTTATTCCTTTCAATCAATTACTGGCGATTCACTGTAAGCGCGGTCt                                                                                                  >  1:2843625/1‑149 (MQ=255)
                                                     tGAAAGCGACAACGCTTGGAAAACGTCTGGCACAGCACCCTTACGATCGGGCGGTGATCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATACCTTATTCCTTTCAATCAATTACTGGCGATTCACTGTAAGCGCGGTCTg                                                                                                 >  2:2417933/1‑149 (MQ=255)
                                                        aaGCGACAACGCTTGGAAAACGTCTGGCACAGCACCCTTACGATCGGGCGGTGATCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATACCTTATTCCTTTCAATCAATTACTGGCGATTCACTGTAAGCGCGGTCTGGTa                                                                                              >  2:3262191/1‑149 (MQ=255)
                                                        aaGCGACAACGCTTGGAAAACGTCTGGCACAGCACCCTTACGATCGGGCGGTGATCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATACCTTATTCCTTTCAATCAATTACTGGCGATTCACTGTAAGCGCGGTCTGGTa                                                                                              <  1:1301407/149‑1 (MQ=255)
                                                            gACAACGCTTGGAAAACGTCTGGCACAGCACCCTTACGATCGGGCGGTGATCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATACCTTATTCCTTTCAATCAATTACTGGCGATTCACTGTAAGCGCGGTCTGGTATggg                                                                                          <  2:2897836/149‑1 (MQ=255)
                                                            gACAACGCTTGGAAAACGTCTGGCACAGCACCCTTACGATCGGGCGGTGATCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATACCTTATTCCTTTCAATCAATTACTGGCGATTCACTGTAAGCGCGGTCTGGTATggg                                                                                          >  1:3198601/1‑149 (MQ=255)
                                                            gACAACGCTTGGAAAACGTCTGGCACAGCACCCTTACGATCGGGCGGTGATCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATACCTTATTCCTTTCAATCAATTACTGGCGATTCACTGTAAGCGCGGTCTGGTATggg                                                                                          <  1:2355102/149‑1 (MQ=255)
                                                            gACAACGCTTGGAAAACGTCTGGCACAGCACCCTTACGATCGGGCGGTGATCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATACCTTATTCCTTTCAATCAATTACTGGCGATTCACTGTAAGCGCGGTCTGGTATggg                                                                                          >  1:2590942/1‑149 (MQ=255)
                                                               aaCGCTTGGAAAACGTCTGGCACAGCACCCTTACGATCGGGCGGTGATCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATACCTTATTCCTTTCAATCAATTACTGGCGATTCACTGTAAGCGCGGTCTGGTATGGGGCg                                                                                       <  2:466660/149‑1 (MQ=255)
                                                               aaCGCTTGGAAAACGTCTGGCACAGCACCCTTACGATCGGGCGGTGATCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATACCTTATTCCTTTCAATCAATTACTGGCGATTCACTGTAAGCGCGGTCTGGTATGGGGCg                                                                                       >  2:1034675/1‑149 (MQ=255)
                                                                 cGCTTGGAAAACGTCTGGCACAGCACCCTTACGATCGGGCGGTGATCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATACCTTATTCCTTTCAATCAATTACTGGCGATTCACTGTAAGCGCGGTCTGGTATGGGGCGAg                                                                                     <  2:4630539/149‑1 (MQ=255)
                                                                   cTTGGAAAACGTCTGGCACAGCACCCTTACGATCGGGCGGTGATCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATACCTTATTCCTTTCAATCAATTACTGGCGATTCACTGTAAGCGCGGTCTGGTATGGGGCGAGCt                                                                                   <  2:946100/149‑1 (MQ=255)
                                                                      ggAAAACGTCTGGCACAGCACCCTTACGATCGGGCGGTGATCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATACCTTATTCCTTTCAATCAATTACTGGCGATTCACTGTAAGCGCGGTCTGGTATGGGGCGAGCTGGa                                                                                >  1:3474514/1‑149 (MQ=255)
                                                                         aaaCGTCTGGCACAGCACCCTTACGATCGGGCGGTGATCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATACCTTATTCCTTTCAATCAATTACTGGCGATTCACTGTAAGCGCGGTCTGGTATGGGGCGAGCTGGAAtt                                                                             <  2:2846724/149‑1 (MQ=255)
                                                                            cGTCTGGCACAGCACCCTTACGATCGGGCGGTGATCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATACCTTATTCCTTTCAATCAATTACTGGCGATTCACTGTAAGCGCGGTCTGGTATGGGGCGAGCTGGAATTTGt                                                                          <  2:1952947/149‑1 (MQ=255)
                                                                            cGTCTGGCACAGCACCCTTACGATCGGGCGGTGATCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATACCTTATTCCTTTCAATCAATTACTGGCGATTCACTGTAAGCGCGGTCTGGTATGGGGCGAGCTGGAATTTGt                                                                          <  1:4284105/149‑1 (MQ=255)
                                                                             gTCTGGCACAGCACCCTTACGATCGGGCGGTGATCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATACCTTATTCCTTTCAATCAATTACTGGCGATTCACTGTAAGCGCGGTCTGGTATGGGGCGAGCTGGAATTTGTa                                                                         <  2:3825639/149‑1 (MQ=255)
                                                                             gTCTGGCACAGCACCCTTACGATCGGGCGGTGATCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATACCTTATTCCTTTCAATCAATTACTGGCGATTCACTGTAAGCGCGGTCTGGTATGGGGCGAGCTGGAATTTGTa                                                                         <  2:3570640/149‑1 (MQ=255)
                                                                             gTCTGGCACAGCACCCTTACGATCGGGCGGTGATCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATACCTTATTCCTTTCAATCAATTACTGGCGATTCACTGTAAGCGCGGTCTGGTATGGGGCGAGCTGGAATTTGTa                                                                         >  2:281537/1‑149 (MQ=255)
                                                                              tCTGGCACAGCACCCTTACGATCGGGCGGTGATCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATACCTTATTCCTTTCAATCAATTACTGGCGATTCACTGTAAGCGCGGTCTGGTATGGGGCGAGCTGGAATTTGTAc                                                                        <  2:781930/149‑1 (MQ=255)
                                                                              tCTGGCACAGCACCCTTACGATCGGGCGGTGATCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATACCTTATTCCTTTCAATCAATTACTGGCGATTCACTGTAAGCGCGGTCTGGTATGGGGCGAGCTGGAATTTGTAc                                                                        <  2:4375104/149‑1 (MQ=255)
                                                                               cTGGCACAGCACCCTTACGATCGGGCGGTGATCCTCAATGCCGGGATTCAAGTCTCCGGCGATCGCCACGAATACCTTATTCCGTTCAATCCATTACTGGCGAGTCACTGTAAGCGCGGTCTGGTATGGGGCGAGCTGGAAGTTGTACt                                                                       >  2:2626005/1‑149 (MQ=255)
                                                                                tGGCACAGCACCCTTACGATCGGGCGGTGATCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATACCTTATTCCTTTCAATCAATTACTGGCGATTCACTGTAAGCGCGGTCTGGTATGGGGCGAGCTGGAATTTGTACTg                                                                      >  1:1025643/1‑149 (MQ=255)
                                                                                 ggCACAGCACCCTTACGATCGGGCGGTGATCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATACCTTATTCCTTTCAATCAATTACTGGCGATTCACTGTAAGCGCGGTCTGGTATGGGGCGAGCTGGAATTTGTACTGc                                                                     >  1:2857664/1‑149 (MQ=255)
                                                                                 ggCACAGCACCCTTACGATCGGGCGGTGATCATCAATGACGGGATTAAAGTCTCCGTCGATCGCCACGAATACCTTATTCATTTCAATCAATTAAAGGCGATTCAATGTAAGCGCGGTCTGGTATGGGGCGAGCTGGAATTTGTACTGc                                                                     >  2:4166479/1‑149 (MQ=255)
                                                                                  gCACAGCACCCTTACGATCGGGCGGTGATCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATACCTTATTCCTTTCAATCAATTACTGGCGATTCACTGTAAGCGCGGTCTGGTATGGGGCGAGCTGGAATTTGTACTGcc                                                                    <  2:3064508/149‑1 (MQ=255)
                                                                                       gCACCCTTACGATCGGGCGGTGATCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATACCTTATTCCTTTCAATCAATTACTGGCGATTCACTGTAAGCGCGGTCTGGTATGGGGCGAGCTGGAATTTGTACTGCCGGACg                                                               >  1:1034190/1‑149 (MQ=255)
                                                                                        cACCCTTACGATCGGGCGGTGATCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATACCTTATTCCTTTCAATCAATTACTGGCGATTCACTGTAAGCGCGGTCTGGTATGGGGCGAGCTGGAATTTGTACTGCCGGACGa                                                              <  2:2958242/149‑1 (MQ=255)
                                                                                          cccTTACGATCGGGCGGTGATCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATACCTTATTCCTTTCAATCAATTACTGGCGATTCACTGTAAGCGCGGTCTGGTATGGGGCGAGCTGGAATTTGTACTGCCGGACGaaa                                                            >  2:3668605/1‑149 (MQ=255)
                                                                                          cccTTACGATCGGGCGGTGATCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATACCTTATTCCTTTCAATCAATTACTGGCGATTCACTGTAAGCGCGGTCTGGTATGGGGCGAGCTGGAATTTGTACTGCCGGACGaaa                                                            >  2:4783667/1‑149 (MQ=255)
                                                                                           ccTTACGATCGGGCGGTGATCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACTAATACCTTATTCCTTTCAATCAATTACTGGCGATTCACTGTAAGCGCGGTCTGGTATGGGGCGAGCTGGAATTTGTACTGCCGGACGaaaa                                                           <  1:2388476/149‑1 (MQ=255)
                                                                                            cTTACGATCGGGCGGTGATCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATACCTTATTCCTTTCAATCAATTACTGGCGATTCACTGTAAGCGCGGTCTGGTATGGGGCGAGCTGGAATTTGTACTGCCGGACGaaaaa                                                          >  1:1150309/1‑149 (MQ=255)
                                                                                              tACGATCGGGCGGTGATCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATACCTTATTCCTTTCAATCAATTACTGGCGATTCACTGTAAGCGCGGTCTGGTATGGGGCGAGCTGGAATTTGTACTGCCGGACGAAAAAgt                                                        <  2:1269487/149‑1 (MQ=255)
                                                                                                cGATCGTGCGGTGATCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATACCTTATTCCTTTCAATCAATTACTGGCGATTCACTGTAAGCGCGGTCTGGTATGGGGCGAGCTGGAATTTGTACTGCCGGACGAAAAAgtgg                                                      <  1:2527818/149‑1 (MQ=255)
                                                                                                 gATCGGGCGGTGATCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATACCTTATTCCTTTCAATCAATTACTGGCGATTCACTGTAAGCGCGGTCTGGTATGGGGCGAGCTGGAATTTGTACTGCCGGACGAAAAAgtggt                                                     >  2:1153914/1‑149 (MQ=255)
                                                                                                 gATCGGGCGGTGATCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATACCTTATTCCTTTCAATCAATTACTGGCGATTCACTGTAAGCGCGGTCTGGTATGGGGCGAGCTGGAATTTGTACTGCCGGACGAAAAAgtggt                                                     <  2:2323077/149‑1 (MQ=255)
                                                                                                  aTCGGGCGGTGATCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATACCTTATTCCTTTCAATCAATTACTGGCGATTCACTGTAAGCGCGGTCTGGTATGGGGCGAGCTGGAATTTGTACTGCCGGACGAAAAAgtggtg                                                    >  2:336786/1‑149 (MQ=255)
                                                                                                   tCGGGCGGTGATCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATACCTTATTCCTTTCAATCAATTACTGGCGATTCACTGTAAGCGCGGTCTGGTATGGGGCGAGCTGGAATTTGTACTGCCGGACGAAAAAGTGGTGc                                                   <  2:2050219/149‑1 (MQ=255)
                                                                                                      ggCGGTGATCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATACCTTATTCCTTTCAATCAATTACTGGCGATTCACTGTAAGCGCGGTCTGGTATGGGGCGAGCTGGAATTTGTACTGCCGGACGAAAAAGTGGTGCGTc                                                <  1:914526/149‑1 (MQ=255)
                                                                                                          gTGATCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATACCTTATTCCTTTCAATCAATTACTGGCGATTCACTGTAAGCGCGGTCTGGTATGGGGCGAGCTGGAATTTGTACTGCCGGACGAAAAAGTGGTGCGTCTGCa                                            >  1:3777613/1‑149 (MQ=255)
                                                                                                              tCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATACCTTATTCCTTTCAATCAATTACTGGCGATTCACTGTAAGCGCGGTCTGGTATGGGGCGAGCTGGAATTTGTACTGCCGGACGAAAAAGTGGTGCGTCTGCACGGc                                        <  1:3902280/149‑1 (MQ=255)
                                                                                                               ccTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATACCTTATTCCTTTCAATCAATTACTGGCGATTCACTGTAAGCGCGGTCTGGTATGGGGCGAGCTGGAATTTGTACTGCCGGACGAAAAAGTGGTGCGTCTGCACGGCa                                       <  2:4534323/149‑1 (MQ=255)
                                                                                                                cTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATACCTTATTCCTTTCAACCAATTACTGGCGATTCACTGTAAGCGCGGTCTGGTATGGGGCGAGCTGGAATTTGTACTGCCGGACGAAAAAGTGGTGCGTCTGCACGGCAc                                      >  2:3674539/1‑149 (MQ=255)
                                                                                                                cTCAAGGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATACCTTATTCCTTTCAATCAATTACTGGCGATTCACTGTAAGCGCGGTCTGGTATGGGGCGAGCTGGAATTTGTACTGCCGGACGAAAAAGTGGTGCGTCTGCACGGCAc                                      <  2:1708985/149‑1 (MQ=255)
                                                                                                                 tCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATACCTTATTCCTTTCAATCAATTACTGGCGATTCACTGTAAGCGCGGTCTGGTATGGGGCGAGCTGGAATTTGTACTGCCGGACGAAAAAGTGGTGCGTCTGCACGGCAcc                                     >  2:880511/1‑149 (MQ=255)
                                                                                                                   aaTGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATACCTTATTCCTTTCAATCAATTACTGGCGATTCACTGTAAGCGCGGTCTGGTATGGGGCGAGCTGGAATTTGTACTGCCGGACGAAAAGTGGTGCGTCTGCACGGCACCGaa                                  <  1:305133/149‑1 (MQ=255)
                                                                                                                    aTGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATACCTTATTCCTTTCAATCAATTACTGGCGATTCACTGTAAGCGCGGTCTGGTATGGGGCGAGCTGGAATTTGTACTGCCGGACGAAAAAGTGGTGCGTCTGCACGGCACCGaa                                  >  1:2651601/1‑149 (MQ=255)
                                                                                                                     tGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATACCTTATTCCTTTCAATCAATTACTGGCGATTCACTGTAAGCGCGGTCTGGTATGGGGCGAGCTGGAATTTGTACTGCCGGACGAAAAAGTGGTGCGTCTGCACGGCACCGAAt                                 <  2:2167483/149‑1 (MQ=255)
                                                                                                                      gCCGGGATTAAAGTCTCCGGCGATCGCCACGAATACCTTATTCCTTTCAATCAATTACTGGCGATTCACTGTAAGCGCGGTCTGGTATGGGGCGAGCTGGAATTTGTACTGCCGGACGAAAAAGTGGTGCGTCTGCACGGCACCGAATg                                <  1:3151714/149‑1 (MQ=255)
                                                                                                                       ccGGGATTAAAGTCTCCGGCGATCGCCACGAATACCTTATTCCTTTCAATCAATTACTGGCGATTCACTGTAAGCGCGGTCTGGTATGGGGCGAGCTGGAATTTGTACTGCCGGACGAAAAAGTGGTGCGTCTGCACGGCACCGAATgg                               <  1:2184071/149‑1 (MQ=255)
                                                                                                                        cGGGATTAAAGTCTCCGGCGATCGCCACGAATACCTTATTCCTTTCAATCAATTACTGGCGATTCACTGTAAGCGCGGTCTGGTATGGGGCGAGCTGGAATTTGTACTGCCGGACGAAAAAGTGGTGCGTCTGCACGGCACCGAATggg                              <  1:4100397/149‑1 (MQ=255)
                                                                                                                          ggATTAAAGTCTCCGGCGATCGCCACGAATACCTTATTCCTTTCAATCAATTACTGGCGATTCACTGTAAGCGCGGTCTGGTATGGGGCGAGCTGGAATTTGTACTGCCGGACGAAAAAGTGGTGCGTCTGCACGGCACCGAATGGGGc                            >  1:4643315/1‑149 (MQ=255)
                                                                                                                               aaaGTCTCCGGCGATCGCCACGAATACCTTATTCCTTTCAATCAATTACTGGCGATTCACTGTAAGCGCGGTCTGGTATGGGGCGAGCTGGAATTTGTACTGCCGGACGAAAAAGTGGTGCGTCTGCACGGCACCGAATGGGGCGAGAc                       <  1:4733590/149‑1 (MQ=255)
                                                                                                                                  gTCGCCGGCGATCGCCACGAATACCTTATTCCTTTCAATCAATGACTGGCGATTCACTGTAAGCGCGGTCTGGTATGGGGCGAGCTGGAATTTGTACTGCCGGACGAAAAAGTGGTGCGTCTGCACGGCACCGAATGGGGCGAGACGCa                    <  2:4203828/149‑1 (MQ=255)
                                                                                                                                      ccGGCGATCGCCACGAATACCTTATTCCTTTCAATCAATTACTGGCGATTCACTGTAAGCGCGGTCTGGTATGGGGCGAGCTGGAATTTGTACTGCCGGACGAAAAAGTGGTGCGTCTGCACGGCACCGAATGGGGCGAGACGCAGCGt                <  1:1378085/149‑1 (MQ=255)
                                                                                                                                      ccGGCGATCGCCACGAATACCTTATTCCTTTCAATCAATTACTGGCGATTCACTGTAAGCGCGGTCTGGTATGGGGCGAGCTGGAATTTGTACTGCCGGACGAAAAAGTGGTGCGTCTGCACGGCACCGAATGGGGCGAGACGCAGCGt                <  2:2126641/149‑1 (MQ=255)
                                                                                                                                       cGGCGATCGCCACGAATACCTTATTCCTTTCAATCAATTACTGGCGATTCACTGTAAGCGCGGTCTGGTATGGGGCGAGCTGGAATTTGTACTGCCGGACGAAAAAGTGGTGCGTCTGCACGGCACCGAATGGGGCGAGACGCAGCGtt               <  2:3217095/149‑1 (MQ=255)
                                                                                                                                        ggCGATCGCCACGAATACCTTATTCCTTTCAATCAATTACTGGCGATTCACTGTAAGCGCGGTCTGGTATGGGGCGAGCTGGAATTTGTACTGCCGGACGAAAAAGTGGTGCGTCTGCACGGCACCGAATGGGGCGAGACGCAGCGttt              >  2:937128/1‑149 (MQ=255)
                                                                                                                                         gCGATCGCCACGAATACCTTATTCCTTTCAATCAATTACTGGCGATTCACTGTAAGCGCGGTCTGGTATGGGGCGAGCTGGAATTTGTACTGCCGGACGAAAAAGTGGTGCGTCTGCACGGCACCGAATGGGGCGAGACGCAGCGtttt             <  2:1385405/149‑1 (MQ=255)
                                                                                                                                           gATCGCCACTAATACCTTATTCCTTTCAATCAATTACTGGCGATTCACTGTAAGCGCGGTCTGGTATGGGGCGAGCTGGAATTTGTACTGCCGGACGAAAAAGTGGTGCGTCTGCACGGCACCGAATGGGGCGAGACGCAGCGTTTTTa           >  2:3433960/1‑149 (MQ=255)
                                                                                                                                                ccACTAATACCTTATTCCTTTCAATCAATTACTGGCGATTCACTGTAAGCGCGGTCTGGTATGGGGCGAGCTGGAATTTGTACTGCCGGACGAAAAAGTGGTGCGTCTGCACGGCACCGAATGGGGCGAGACGCAGCGTTTTTACcatc      >  2:1601684/1‑149 (MQ=255)
                                                                                                                                                ccACGAATACCTTATTCCTTTCAATCAATTACTGGCGATTCACTGTAAGCGCGGTCTGGTATGGGGCGAGCTGGAATTTGTACTGCCGGACGAAAAAGTGGTGCGTCTGCACGGCACCGACTGGGGCGAGACGCAGCGTTTTTACcatc      >  1:4461073/1‑149 (MQ=255)
                                                                                                                                                ccACGAATACCTTATTCCTTTCAATCAATTACTGGCGATTCACTGTAAGCGCGGTCTGGTATGGGGCGAGCTGGAATTTGTACTGCCGGACGAAAAAGTGGTGCGTCTGCACGGCACCGAATGGGGCGAGACGCAGCGTTTTTACcatc      >  2:1054761/1‑149 (MQ=255)
                                                                                                                                                ccACGAATACCTTAGTCCTTTCAATCAATTACTGGCGATTCACTGTAAGCGCGGTCTGGTATGGGGCGAGCTGGAATTTGTACTGCCGGACGAAAAAGTGGTGCGTCTGCACGGCACCGAATGGGGCGAGACGCAGCGTTTTTACcatc      <  2:848374/149‑1 (MQ=255)
                                                                                                                                                 cACGAATACCTTATTCCTTTCAATCAATTACTGGCGATTCACTGTAAGCGCGGTCTGGTATGGGGCGAGCTGGAATTTGTACTGCCGGACGAAAAAGTGGTGCGTCTGCACGGCACCGAATGGGGCGAGACGCAGCGTTTTTACcatca     >  1:2100025/1‑149 (MQ=255)
                                                                                                                                                 cACGAATACCTTATTCCTTTCAATCAATTACTGGCGATTCACTGTAAGCGCGGTCTGGTATGGGGCGAGCTGGAATTTGTACTGCCGGACGAAAAAGTGGTGCGTCTGCACGGCACCGAATGGGGCGAGACGCAGCGTTTTTACcatca     >  2:2734879/1‑149 (MQ=255)
                                                                                                                                                    gAATACCTTATTCCTTTCAATCAATTACTGGCGATTCACTGTAAGCGCGGTCTGGTATGGGGCGAGCTGGAATTTGTACTGCCGGACGAAAAAGTGGTGCGTCTGCACGGCACCGAATGGGGCGAGACGCAGCGTTTTTACCATCATCt  >  1:4721747/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GACGGTTACACTGATAAGAAATAATTTCGTGTGGATCTACAGAGTCATGGAACTGAAAGCGACAACGCTTGGAAAACGTCTGGCACAGCACCCTTACGATCGGGCGGTGATCCTCAATGCCGGGATTAAAGTCTCCGGCGATCGCCACGAATACCTTATTCCTTTCAATCAATTACTGGCGATTCACTGTAAGCGCGGTCTGGTATGGGGCGAGCTGGAATTTGTACTGCCGGACGAAAAAGTGGTGCGTCTGCACGGCACCGAATGGGGCGAGACGCAGCGTTTTTACCATCATCT  >  CP009273/1019881‑1020177

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: