Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,802,110 T→A 6.5% pseudogene (436/1178 nt) ygaY → pseudogene, major facilitator transporter superfamily;putative transport; Not classified; putative transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,802,1100TA 6.5% 896.4 / 40.2 247pseudogene (436/1178 nt)ygaYpseudogene, major facilitator transporter superfamily;putative transport; Not classified; putative transport protein
Reads supporting (aligned to +/- strand):  ref base T (137/94);  new base A (9/7);  total (146/101)
Fisher's exact test for biased strand distribution p-value = 7.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

Only 100 of 252 total aligned reads displayed.

TGGCGATGATGATCCTCGGTACGGCATTAACCGGTTTATTCTCAGTCGTGGCACAAATTCTGGTTCCGCTGGCAGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTCTTTTGGGTTGCTTCGGTGTTAATGGCACTGATGGCGCTGGCATTATGGCGTGGTCTGCCACAAATGAAATCAGAAACCCACCTCAACTACCCA  >  CP009273/2801963‑2802258
                                                                                                                                                   |                                                                                                                                                    
tGGCGATGATGATCCTCGGTACGGCATTAACCGGTTTATTCTCAGTCGTGGCACAAATTCTGGTTCCGCTGGCAGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTg                                                                                                                                                     >  1:2394888/1‑149 (MQ=255)
tGGCGATGATGATCCTCGGTACGGCATTAACCGGTTTATTCTCAGTCGTGGCACAAATTCTGGTTCCGCTGGCAGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTg                                                                                                                                                     >  2:1934247/1‑149 (MQ=255)
 ggCGATGATGATCCTCGGTACGGCATTAACCGGTTTATTCTCAGTCGTGGCACAAATTCTGGTTCCGCTGGCAGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGc                                                                                                                                                    >  2:3253270/1‑149 (MQ=255)
 ggCGATGATGATCCTCGGTACGGCATTAACCGGTTTATTCTCAGTCGTGGCACAAAGTCTGGTTCCGCTGGCAGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGc                                                                                                                                                    >  2:2814378/1‑149 (MQ=255)
   cGATGATGATCCTCGGTACGGCATTAACCGGTTTATTCTCAGTCGTGGCACAAATTCTGGTTCCGCTGGCAGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTg                                                                                                                                                  <  1:4005993/149‑1 (MQ=255)
     atgatgatCCTCGGTACGGCATTAACCGGTTTATTCTCAGTCGTGGCACAAATTCTGGTTCCGCTGGCAGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTAGCTGGc                                                                                                                                                >  1:1609023/1‑149 (MQ=255)
       gatgatCCTCGGTACGGCATTAACCGGTTTATTCTCAGTCGTGGCACAAATTCTGGTTCCGCTGGCAGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTGGCAc                                                                                                                                              >  1:3966025/1‑149 (MQ=255)
       gatgatCCTCGGTACGGCATTAACCGGTTTATTCTCAGTCGTGGCACAAATTCTGGTTCCGCTGGCAGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTGGCAc                                                                                                                                              <  1:1486576/149‑1 (MQ=255)
        atgatCCTCGGTACGGCATTAACCGGTTTATTCTCAGTCGTGGCACAAATTCTGGTTCCGCTGGCAGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACg                                                                                                                                             >  2:2807568/1‑149 (MQ=255)
        atgatCCTCGGTACGGCATTAACCGGTTTATTCTCAGTCGTGGCACAAATTCTGGTTCCGCTGGCAGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACg                                                                                                                                             >  1:4618218/1‑149 (MQ=255)
         tgatCCTCGGTACGGCATTAACCGGTTTATTCTCAGTCGTGGCACAAATTCTGGTTCCGCTGGCAGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACgg                                                                                                                                            >  1:3561793/1‑149 (MQ=255)
          gatCCTCGGTACGGCATTAACCGGTTTATTCTCAGTCGTGGCACAAATTCTGGTTCCGCTGGCAGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGa                                                                                                                                           >  1:2615219/1‑149 (MQ=255)
           atCCTCGGTACGGCATTAACCGGTTTATTCTCAGTCGTGGCACAAATTCTGGTTCCGCTGGCAGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGAc                                                                                                                                          >  1:2905457/1‑149 (MQ=255)
               tCGGTACGGCATTAACCGGTTTATTCTCAGTCGTGGCACAAATTCTGGTTCCGCTGGCAGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGtt                                                                                                                                      <  2:967595/149‑1 (MQ=255)
                cGGTACGGCATTAACCGGTTTATTCTCAGTCGTGGCACAAATTCTGGTTCCGCTGGCAGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTg                                                                                                                                     >  2:785565/1‑149 (MQ=255)
                  gTACGGCATTAACCGGTTTATTCTCAGTCGTGGCACAAATTCTGGTTCCGCTGGCAGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGcc                                                                                                                                   >  1:4597542/1‑149 (MQ=255)
                   tACGGCATTAACCGGTTTATTCTCAGTCGTGGCACAAATTCTGGTTCCGCTGGCAGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCg                                                                                                                                  >  2:3583625/1‑149 (MQ=255)
                    aCGGCATTAACCGGTTTATTCTCAGTCGTGGCACAAATTCTGGTTCCGCTGGCAGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCgg                                                                                                                                 >  1:4478198/1‑149 (MQ=255)
                       gCATTAACCGGTTTATTCTCAGTCGTGGCACAAATTCTGGTTCCGCTGGCAGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCGGAtt                                                                                                                              >  1:2026707/1‑149 (MQ=255)
                        cATTAACCGGTTTATTCTCAGTCGTGGCACAAATTCTGGTTCCGCTGGCAGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTg                                                                                                                             >  1:3301814/1‑149 (MQ=255)
                         aTTAACCGGTTTATTCTCAGTCGTGGCACAAATTCTGGTTCCGCTGGCAGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGc                                                                                                                            >  1:3177430/1‑149 (MQ=255)
                          ttAACCGGTTTATTCTCAGTCGTGGCACAAATTCTGGTTCCGCTGGCAGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCt                                                                                                                           >  2:2514764/1‑149 (MQ=255)
                             aCCGGTTTATTCTCAGTCGTGGCACAAATTCTGGTTCCGCTGGCAGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGc                                                                                                                        >  1:896668/1‑149 (MQ=255)
                             aCCGGTTTATTCTCAGTCGTGGCACAAATTCTGGTTCCGCTGGCAGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGc                                                                                                                        >  2:501177/1‑149 (MQ=255)
                                 gTTTATTCTCAGTCGTGGCACAAATTCTGGTTCCGCTGGCAGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAAt                                                                                                                    >  1:1133597/1‑149 (MQ=255)
                                 gTTTATTCTCAGTCGTGGCACAAATTCTGGTTCCGCTGGCAGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAAt                                                                                                                    <  1:2460237/149‑1 (MQ=255)
                                    tagtctcAGTCGTGGCACAAATTCTGGTTCCGCTGGCAGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAAtctc                                                                                                                 <  1:4154030/146‑1 (MQ=255)
                                    tagtctcAGTCGTGGCACAAATTCTGGTTCCGCTGGCAGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAAtctc                                                                                                                 <  2:2378758/146‑1 (MQ=255)
                                    tATTCTCAGTCGTGGCACAAATTCTGGTTCCGCTGGCAGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAAtctc                                                                                                                 >  2:1374991/1‑149 (MQ=255)
                                     aTTCTCAGTCGTGGCACAAATTCTGGTTCCGCTGGCAGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTcg                                                                                                                >  2:1700040/1‑149 (MQ=255)
                                      ttCTCAGTCGTGGCACAAATTCTGGTTCCGCTGGCAGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTcgg                                                                                                               <  1:4506407/149‑1 (MQ=255)
                                      ttCTCAGTCGTGGCACAAATTCTGGTTCCGCTGGCAGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTcgg                                                                                                               >  1:4189812/1‑149 (MQ=255)
                                       tctcAGTCGTGGCACAAATTCTGGTTCCGCTGGCAGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTcggc                                                                                                              >  2:1834083/1‑149 (MQ=255)
                                          cAGTCGTGGCACAAATTCTGGTTCCGCTGGCAGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCggcggc                                                                                                           >  2:3014352/1‑149 (MQ=255)
                                          cAGTCGTGGCACAAATTCTGGTTCCGCTGGCAGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTAGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCggcggc                                                                                                           <  1:2992494/149‑1 (MQ=255)
                                            gTCGTGGCACAAATTCTGGTTCCGCTGGCAGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTg                                                                                                         >  1:4263807/1‑149 (MQ=255)
                                               gTGGCACAAATTCTGGTTCCGCTGGCAGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTAGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGgcg                                                                                                      >  2:4370813/1‑149 (MQ=255)
                                                      aaaTTCTGGTTCCGCTGGCAGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTc                                                                                               <  1:103429/149‑1 (MQ=255)
                                                      aaaTTCTGGTTCCGCTGGCAGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTc                                                                                               >  1:60243/1‑149 (MQ=255)
                                                       aaTTCTGGTTCCGCTGGCAGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTCt                                                                                              <  2:1515763/149‑1 (MQ=255)
                                                       aaTTCTGGTTCCGCTGGCAGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTCt                                                                                              >  1:931852/1‑149 (MQ=255)
                                                        agtCTGGTTCCGCTGGCAGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTAGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTCtt                                                                                             <  1:516831/147‑1 (MQ=255)
                                                        aTTCTGGTTCCGCTGGCAGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTCtt                                                                                             <  1:1807366/149‑1 (MQ=255)
                                                        aTTCTGGTTCCGCTGGCAGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTCtt                                                                                             >  2:3840432/1‑149 (MQ=255)
                                                        aTTCTGGTTCCGCTGGCAGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTCtt                                                                                             >  1:1851523/1‑149 (MQ=255)
                                                         ttCTGGTTCCGCTGGCAGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTCttt                                                                                            <  1:3606282/149‑1 (MQ=255)
                                                              gTTCCGCTGGCAGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTCTTTTGGGt                                                                                       <  1:1786777/149‑1 (MQ=255)
                                                              gTTCCGCTGGCAGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTCTTTTGGGt                                                                                       >  2:2242317/1‑149 (MQ=255)
                                                               ttCCGCTGGCAGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTCTTTTGGGtt                                                                                      >  1:2942945/1‑149 (MQ=255)
                                                               ttCCGCTGGCAGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTCTTTTGGGtt                                                                                      <  2:2013873/149‑1 (MQ=255)
                                                                 ccGCTGGCAGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTCTTTTGGGTTGc                                                                                    <  2:3168059/149‑1 (MQ=255)
                                                                  cGCTGGCAGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTCTTTTGGGTTGCt                                                                                   <  2:2640033/149‑1 (MQ=255)
                                                                     tGGCAGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTCTTTTGGGTTGCTTCg                                                                                >  1:4125965/1‑149 (MQ=255)
                                                                         aGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTCTTTTGGGTTGCTTCGgtgt                                                                            <  1:4364050/149‑1 (MQ=255)
                                                                         aGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTCTTTTGGGTTGCTTCGgtgt                                                                            >  1:2893777/1‑149 (MQ=255)
                                                                         aGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTCTTTTGGGTTGCTTCGgtgt                                                                            <  2:2207318/149‑1 (MQ=255)
                                                                         aGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTCTTTTGGGTTGCTTCGgtgt                                                                            <  2:3438356/149‑1 (MQ=255)
                                                                         aGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTCTTTTGGGTTGCTTCGgtgt                                                                            >  2:4283513/1‑149 (MQ=255)
                                                                             aCGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTCTTTTGGGTTGCTTCGGTGTTAAt                                                                        >  2:4079976/1‑149 (MQ=255)
                                                                             aCGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTCTTTTGGGTTGCTTCGGTGTTAAt                                                                        >  2:3389927/1‑149 (MQ=255)
                                                                               gCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTAGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTCTTTTGGGTTGCTTCGGTGTTAATgg                                                                      >  1:1370968/1‑149 (MQ=255)
                                                                                cTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTCTTTTGGGTTGCTTCGGTGTTAATGGc                                                                     >  1:183561/1‑149 (MQ=255)
                                                                                  ggCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTCTTTTGGGTTGCTTCGGTGTTAATGGCAc                                                                   >  1:878563/1‑149 (MQ=255)
                                                                                   gCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTCTTTTGGGTTGCTTCGGTGTTAATGGCACt                                                                  >  2:3540756/1‑149 (MQ=255)
                                                                                     ttCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTCTTTTGGGTTGCTTCGGTGTTAATGGCACTGa                                                                <  1:3909294/149‑1 (MQ=255)
                                                                                       cACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTCTTTTGGGTTGCTTCGGTGTTAATGGCACTGATg                                                              <  2:4238557/149‑1 (MQ=255)
                                                                                           ggATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTCTTTTGGGTTGCTTCGGTGTTAATGGCACTGATGgcgc                                                          >  1:3450312/1‑149 (MQ=255)
                                                                                           ggATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTCTTTTGGGTTGCTTCGGTGTTAATGGCACTGATGgcgc                                                          >  1:2731053/1‑149 (MQ=255)
                                                                                             aTAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTCTTTTGGGTTGCTTCGGTGTTAATGGCACTGATGGCGCTg                                                        <  2:3872672/149‑1 (MQ=255)
                                                                                             aTAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTCTTTTGGGTTGCTTCGGTGTTAATGGCACTGATGGCGCTg                                                        <  2:4049573/149‑1 (MQ=255)
                                                                                             aTAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTCTTTTGGGTTGCTTCGGTGTTAATGGCACTGATGGCGCTg                                                        <  2:748452/149‑1 (MQ=255)
                                                                                              tAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTCTTTTGGGTTGCTTCGGTGTTAATGGCACTGATGGCGCTgg                                                       >  1:2096941/1‑149 (MQ=255)
                                                                                                 aCGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTCTTTTGGGTTGCTTCGGTGTTAATGGCACTGATGGCGCTGGCAt                                                    <  1:215042/149‑1 (MQ=255)
                                                                                                   gcgGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTCTTTTGGGTTGCTTCGGTGTTAATGGCACTGATGGCGCTGGCATTa                                                  <  1:4375667/149‑1 (MQ=255)
                                                                                                     ggCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTCTTTTGGGTTGCTTCGGTGTTAATGGCACTGATGGCGCTGGCATTATg                                                <  1:2392188/149‑1 (MQ=255)
                                                                                                        aaaGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTCTTTTGGGTTGCTTCGGTGTTAATGGCACTGATGGCGCTGGCATTATGGCg                                             >  2:293551/1‑149 (MQ=255)
                                                                                                         aaGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTCTTTTGGGTTGCTTCGGTGTTAATGGCACTGATGGCGCTGGCATTATGGCGt                                            <  2:334760/149‑1 (MQ=255)
                                                                                                          aGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTCTTTTGGGTTGCTTCGGTGTTAATGGCACTGATGGCGCTGGCATTATGGCGTg                                           <  2:2706563/149‑1 (MQ=255)
                                                                                                                 ggCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTCTTTTGGGTTGCTTCGGTGTTAATGGCACTGATGGCGCTGGCATTATGGCGTGGTCTGcc                                    >  2:4221637/1‑149 (MQ=255)
                                                                                                                       attatGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTCTTTTGGGTTGCTTCGGTGTTAATGGCACTGATGGCGCTGGCATTATGGCGTGGTCTGCCACAAAt                              <  1:774686/149‑1 (MQ=255)
                                                                                                                       attatGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTCTTTTGGGTTGCTTCGGTGTTAATGGCACTGATGGCGCTGGCATTATGGCGTGGTCTGCCACAAAt                              <  2:1612952/149‑1 (MQ=255)
                                                                                                                        ttatGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTCTTTTGGGTTGCTTCGGTGTTAATGGCACTGATGGCGCTGGCATTATGGCGTGGTCTGCCACAAATg                             >  1:3705726/1‑149 (MQ=255)
                                                                                                                        ttatGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTCTTTTGGGTTGCTTCGGTGTTAATGGCACTGATGGCGCTGGCATTATGGCGTGGTCTGCCACAAATg                             <  1:679807/149‑1 (MQ=255)
                                                                                                                        ttatGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTCTTTTGGGTTGCTTCGGTGTTAATGGCACTGATGGCGCTGGCATTATGGCGTGGTCTGCCACAAATg                             >  2:4070949/1‑149 (MQ=255)
                                                                                                                        ttatGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTCTTTTGGGTTGCTTCGGTGTTAATGGCACTGATGGCGCTGGCATTATGGCGTGGTCTGCCACAAATg                             >  2:1739537/1‑149 (MQ=255)
                                                                                                                        ttatGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTCTTTTGGGTTGCTTCGGTGTTAATGGCACTGATGGCGCTGGCATTATGGCGTGGTCTGCCACAAATg                             <  2:3509425/149‑1 (MQ=255)
                                                                                                                              gCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTCTTTTGGGTTGCTTCGGTGTTAATGGCACTGATGGCGCTGGCATTATGGCGTGGTCTGCCACAAATGAAATCa                       >  1:1181823/1‑149 (MQ=255)
                                                                                                                                ggTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTATTTTGGGTTGCTTCGGTGTTAATGGCACTGATGGCGCTGGCATTATGGCGAGGTCTGCCACAAATGAAATCAGa                     >  2:1248219/1‑149 (MQ=255)
                                                                                                                                  tCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTCTTTTGGGTTGCTTCGGTGTTAATGGCACTGATGGCGCTGGCATTATGGCGTGGTCTGCCACAAATGAAATCAGaaa                   >  2:1895542/1‑149 (MQ=255)
                                                                                                                                  tCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTCTTTTGGGTTGCTTCGGTGTTAATGGCACTGATGGCGCTGGCATTATGGCGTGGTCTGCCACAAATGAAATCAGaaa                   <  1:1374991/149‑1 (MQ=255)
                                                                                                                                   ctgctgTTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTCTTTTGGGTTGCTTCGGTGTTAATGGCACTGATGGCGCTGGCATTATGGCGTGGTCTGCCACAAATGAAATCAGAAAc                  <  2:617756/149‑1 (MQ=255)
                                                                                                                                        gTTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTCTTTTGGGTTGCTTCGGTGTTAATGGCACTGATGGCGCTGGCATTATGGCGTGGTCTGCCACAAATGAAATCAGAAACCCAcc             >  1:4079554/1‑149 (MQ=255)
                                                                                                                                          tGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTCTTTTGGGTTGCTTCGGTGTTAATGGCACTGATGGCGCTGGCATTATGGCGTGGTCTGCCACAAATGAAATCAGAAACCCACCTc           >  1:1273794/1‑149 (MQ=255)
                                                                                                                                           ggggATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTCTTTTGGGTTGCTTCGGTGTTAATGGCACTGATGGCGCTGGCATTATGGCGTGGTCTGCCACAAATGAAATCAGAAACCCACCTCa          >  1:4397736/1‑149 (MQ=255)
                                                                                                                                            gggATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTCTTTTGGGTTGCTTCGGTGTTAATGGCACTGATGGCGCTGGCATTATGGCGTGGTCTGCCACAAATGAAATCAGAAACCCACCTCaa         <  2:4448660/149‑1 (MQ=255)
                                                                                                                                               aTCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTCTTTTGGGTTGCTTCGGTGTTAATGGCACTGATGGCGCTGGCATTATGGCGTGGTCTGCCACAAATGAAATCAGAAACCCACCTCAACTa      >  2:4000423/1‑149 (MQ=255)
                                                                                                                                                tCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTCTTTTGGGTTGCTTCGGTGTTAATGGCACTGATGGCGCTGGCATTATGGCGTGGTCTGCCACAAATGAAATCAGAAACCCACCTCAACTAc     <  1:2454298/149‑1 (MQ=255)
                                                                                                                                                  ttGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTCTTTTGGGTTGCTTCGGTGTTAATGGCACTGATGGCGCTGGCATTATGGCGTGGTCTGCCACAAATGAAATCAGAAACCCACCTCAACTAccc   >  1:1974844/1‑149 (MQ=255)
                                                                                                                                                  ttGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTCTTTTGGGTTGCTTCGGTGTTAATGGCACTGATGGCGCTGGCATTATGGCGTGGTCTGCCACAAATGAAATCAGAAACCCACCTCAACTAccc   <  2:4597542/149‑1 (MQ=255)
                                                                                                                                                  ttGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTCTTTTGGGTTGCTTCGGTGTTAATGGCACTGATGGCGCTGGCATTATGGCGTGGTCTGCCACAAATGAAATCAGAAACCCACCTCAACTAccc   >  2:645991/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TGGCGATGATGATCCTCGGTACGGCATTAACCGGTTTATTCTCAGTCGTGGCACAAATTCTGGTTCCGCTGGCAGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTCTTTTGGGTTGCTTCGGTGTTAATGGCACTGATGGCGCTGGCATTATGGCGTGGTCTGCCACAAATGAAATCAGAAACCCACCTCAACTACCCA  >  CP009273/2801963‑2802258

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: