Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,139,450 G→T 20.1% E902* (GAG→TAG)  yegE → putative diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,139,4500GT20.1% 858.6 / 141.1 293E902* (GAG→TAG) yegEputative diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base G (126/108);  new base T (26/33);  total (152/141)
Fisher's exact test for biased strand distribution p-value = 1.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 297 total aligned reads displayed.

CGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCTCCCTCCCCCTTTCCGT  >  CP009273/2139302‑2139598
                                                                                                                                                    |                                                                                                                                                    
 gAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCgag                                                                                                                                                    <  1:4189041/150‑1 (MQ=255)
     gggCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGatta                                                                                                                                                <  2:2160108/150‑1 (MQ=255)
       gcggcgATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGattatt                                                                                                                                              <  1:2741301/150‑1 (MQ=255)
        cggcgATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATtg                                                                                                                                             >  2:4022275/1‑150 (MQ=255)
         ggcgATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATtga                                                                                                                                            <  1:4087699/150‑1 (MQ=255)
         ggcgATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCACGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATtga                                                                                                                                            >  1:3110643/1‑150 (MQ=255)
          gcgATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATtgat                                                                                                                                           >  1:4580172/1‑150 (MQ=255)
           cgATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATtgatg                                                                                                                                          <  2:1136749/150‑1 (MQ=255)
               gTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGaaca                                                                                                                                      <  2:2739682/150‑1 (MQ=255)
                  gCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGaacaaac                                                                                                                                   >  1:2963615/1‑150 (MQ=255)
                   cTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCTAGATTATTGATGaacaaaca                                                                                                                                  >  2:965346/1‑150 (MQ=255)
                   cTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGaacaaaca                                                                                                                                  <  2:1079213/150‑1 (MQ=255)
                    ttGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAAAAt                                                                                                                                 <  1:699305/150‑1 (MQ=255)
                    ttGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAAAAt                                                                                                                                 <  1:3202731/150‑1 (MQ=255)
                      gatgaACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCTTGATTATTGATGAACAAACAttt                                                                                                                               <  2:1297813/150‑1 (MQ=255)
                       atgaACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTtc                                                                                                                              >  2:1774773/1‑150 (MQ=255)
                          aaCAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCATCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACAATtcgtc                                                                                                                           >  1:1779563/1‑150 (MQ=255)
                           aCAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTtcgtcg                                                                                                                          <  2:3793142/150‑1 (MQ=255)
                           aCAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTtcgtcg                                                                                                                          <  1:3816093/150‑1 (MQ=255)
                            cAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTcgtcgt                                                                                                                         <  1:1054835/150‑1 (MQ=255)
                             aGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCTAGATTATTGATGAACAAACATTTCGTCGTa                                                                                                                        >  1:1512211/1‑150 (MQ=255)
                             aGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCTAGATTATTGATGAACAAACATTTCGTCGTa                                                                                                                        >  1:1461882/1‑150 (MQ=255)
                                     atgatTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGc                                                                                                                >  1:1983103/1‑150 (MQ=255)
                                      tgatTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCg                                                                                                               <  2:142643/150‑1 (MQ=255)
                                       gatTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGa                                                                                                              >  1:558516/1‑150 (MQ=255)
                                           aaaGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCTAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCg                                                                                                          <  2:2048194/150‑1 (MQ=255)
                                           aaaGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCg                                                                                                          >  1:1935616/1‑150 (MQ=255)
                                           aaaGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCg                                                                                                          >  1:1355079/1‑150 (MQ=255)
                                           aaaGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCg                                                                                                          >  1:1934986/1‑150 (MQ=255)
                                            aaGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCgg                                                                                                         >  1:835627/1‑150 (MQ=255)
                                             agagAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGgc                                                                                                        >  2:50425/1‑150 (MQ=255)
                                             agagAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGgc                                                                                                        >  2:3755281/1‑150 (MQ=255)
                                              gagaATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCTAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGgcg                                                                                                       >  2:2109332/1‑150 (MQ=255)
                                               agaATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGgcgc                                                                                                      >  1:4440968/1‑150 (MQ=255)
                                                  aTCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTa                                                                                                   <  2:4594048/150‑1 (MQ=255)
                                                   tCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAg                                                                                                  >  1:203540/1‑150 (MQ=255)
                                                    cAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGc                                                                                                 >  1:3856391/1‑150 (MQ=255)
                                                         gatgatGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAATGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGc                                                                                            >  1:4418942/1‑150 (MQ=255)
                                                          atgatgCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGct                                                                                           >  2:3614857/1‑150 (MQ=255)
                                                          atgatgCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGct                                                                                           <  2:128670/150‑1 (MQ=255)
                                                           tgatgCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCTAGATTATTGATGAACAAACATTTCGTCGTAGCGTCAGCGATCCGGCGCTTAGCCATGctc                                                                                          <  2:4055478/150‑1 (MQ=255)
                                                             atgCTCGCCCACGGTGTCGCGTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCTAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCtt                                                                                        <  1:2951622/150‑1 (MQ=255)
                                                                 tCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCTAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGAcc                                                                                    >  2:4188520/1‑150 (MQ=255)
                                                                 tCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGAcc                                                                                    >  2:1026264/1‑150 (MQ=255)
                                                                   gCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCTAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGAccgc                                                                                  >  2:1496331/1‑150 (MQ=255)
                                                                    cccACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGAccgcc                                                                                 <  1:2290273/150‑1 (MQ=255)
                                                                      cACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCTAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCgg                                                                               >  2:4515372/1‑150 (MQ=255)
                                                                       aCGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCggg                                                                              <  2:3917660/150‑1 (MQ=255)
                                                                        cGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGt                                                                             >  2:78075/1‑150 (MQ=255)
                                                                         ggTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTa                                                                            >  1:188309/1‑150 (MQ=255)
                                                                                cTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACg                                                                     <  1:4022534/150‑1 (MQ=255)
                                                                                 ttCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGa                                                                    >  1:2635008/1‑150 (MQ=255)
                                                                                 ttCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGa                                                                    >  1:2518854/1‑150 (MQ=255)
                                                                                 ttCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGa                                                                    <  2:4277420/150‑1 (MQ=255)
                                                                                  tCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGaa                                                                   >  1:2112003/1‑150 (MQ=255)
                                                                                    gCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAAtt                                                                 >  1:619252/1‑150 (MQ=255)
                                                                                    gCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAAtt                                                                 >  2:1547988/1‑150 (MQ=255)
                                                                                    gCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAAtt                                                                 >  1:1549584/1‑150 (MQ=255)
                                                                                      cACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAAtttt                                                               <  2:3293740/150‑1 (MQ=255)
                                                                                         ggATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTcc                                                            <  1:3340636/150‑1 (MQ=255)
                                                                                           aTCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCcag                                                          <  2:2811559/150‑1 (MQ=255)
                                                                                           aTCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCcag                                                          >  2:1837001/1‑150 (MQ=255)
                                                                                           aTCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCcag                                                          <  1:3922983/150‑1 (MQ=255)
                                                                                              ccGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCcagcag                                                       >  2:650236/1‑150 (MQ=255)
                                                                                              ccGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCcagcag                                                       >  1:1216987/1‑150 (MQ=255)
                                                                                               cGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAgg                                                      <  1:2913835/150‑1 (MQ=255)
                                                                                                 gAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCTAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGcc                                                    >  2:2652773/1‑150 (MQ=255)
                                                                                                 gAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGcc                                                    >  2:3708736/1‑150 (MQ=255)
                                                                                                       cgTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCaaaa                                              <  1:3558278/150‑1 (MQ=255)
                                                                                                        gTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCTAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAg                                             >  2:2534790/1‑150 (MQ=255)
                                                                                                          aaTTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCg                                           >  1:3242275/1‑150 (MQ=255)
                                                                                                           aTTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCgg                                          <  1:2254486/150‑1 (MQ=255)
                                                                                                              tgtgGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCTAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTg                                       <  2:1098797/150‑1 (MQ=255)
                                                                                                                  gCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCTAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAg                                   >  1:172975/1‑150 (MQ=255)
                                                                                                                   cTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGt                                  >  2:565229/1‑150 (MQ=255)
                                                                                                                    tGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTa                                 <  1:1845774/150‑1 (MQ=255)
                                                                                                                     gATTTCACTTAAGCTCTGGAGTTGCGAAGGCTAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTaa                                <  1:965346/150‑1 (MQ=255)
                                                                                                                      aTTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTaaa                               >  1:2605160/1‑150 (MQ=255)
                                                                                                                        ttCACTTAAGCTCTGGAGTTGCGAAGGCTAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAgg                             <  1:3273658/150‑1 (MQ=255)
                                                                                                                        ttCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAgg                             >  1:3627872/1‑150 (MQ=255)
                                                                                                                            cTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATa                         <  2:3428549/150‑1 (MQ=255)
                                                                                                                            cTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATa                         <  2:1675953/150‑1 (MQ=255)
                                                                                                                             ttAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATaa                        >  2:4256843/1‑150 (MQ=255)
                                                                                                                              tAAGCTCTGGAGTTGCGAAGGCTAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAg                       >  2:655462/1‑150 (MQ=255)
                                                                                                                               aaGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGc                      >  1:2593316/1‑150 (MQ=255)
                                                                                                                                    ctGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCAtctc                 <  1:1150388/150‑1 (MQ=255)
                                                                                                                                    ctGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCAtctc                 <  2:2693920/150‑1 (MQ=255)
                                                                                                                                      ggAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATctccc               >  1:713629/1‑150 (MQ=255)
                                                                                                                                      ggAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATctccc               <  1:330987/150‑1 (MQ=255)
                                                                                                                                         gTTGCGAAGGCTAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCAtc                   <  2:1245553/143‑1 (MQ=255)
                                                                                                                                          ttGCGAAGGCTAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCtccctccc           <  2:3356251/150‑1 (MQ=255)
                                                                                                                                            gCGAAGGCTAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCTCCCTccccc         <  2:467871/150‑1 (MQ=255)
                                                                                                                                            gCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCATGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCTCCCTccccc         >  2:2980086/1‑150 (MQ=255)
                                                                                                                                             cGAAGGCGAGATGATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGATTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCTCCCTCCCCCt        <  2:3110643/150‑1 (MQ=255)
                                                                                                                                              gAAGGCTAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCTCCCTCCCCCtt       >  2:4426021/1‑150 (MQ=255)
                                                                                                                                              gAAGGCTAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCTCCCTCCCCCtt       >  1:1209670/1‑150 (MQ=255)
                                                                                                                                               aaGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCTCCCTCCCCCttt      >  1:3404768/1‑150 (MQ=255)
                                                                                                                                                   cGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCTCCCTCCCCCTTTCCGt  >  2:2595858/1‑150 (MQ=255)
                                                                                                                                                   cGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCTCCCTCCCCCTTTCCGt  >  1:4305092/1‑150 (MQ=255)
                                                                                                                                                   cGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCTCCCTCCCCCTTTCCGt  >  1:2688189/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCTCCCTCCCCCTTTCCGT  >  CP009273/2139302‑2139598

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: