Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,709,841 T→A 31.2% intergenic (‑33/‑272) grcA ← / → ung autonomous glycyl radical cofactor/uracil‑DNA‑glycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,709,8410TA31.2% 374.8 / 249.6 250intergenic (‑33/‑272)grcA/ungautonomous glycyl radical cofactor/uracil‑DNA‑glycosylase
Reads supporting (aligned to +/- strand):  ref base T (76/96);  new base A (35/43);  total (111/139)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 254 total aligned reads displayed.

GCATAACCTGCTTTTGCAACGATGCAACGCGCTTCGCCTTTTTCGCTGTCCAGCAGCCAGAAAGAGTTCAGCAGATCGTCGTTAGCGGCTTTAGTAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTGTATATATACTCCTCAAACACCCTTGAATCTTTGATTTAAATCAATAAAAACCACACATCAAGTATGGTCGCAAATGGATTTTATTGTTTTACATCAACTTATGC  >  CP009273/2709693‑2709989
                                                                                                                                                    |                                                                                                                                                    
gCATAACCTGCTTTTGCAACGATGCAACGCGCTTCGCCTTTTTCGCTGTCCAGCAGCCAGAAAGATTTCAGCAGATCGTCGTTAGCGGCTTTAGTAATCTGGATACCTGTAATAATGTGATGCCTCCCCGGCAAAATTATTTGATTTGtt                                                                                                                                                     >  2:1092062/1‑150 (MQ=255)
gCATAACCTGCTTTTGCAACGATGCAACGCGCTTCGCCTTTTTCGCTGTCCAGCAGCCAGAAAGAGTTCAGCAGATCGTCGTTAGCGGCTTTAGTAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGtt                                                                                                                                                     >  1:930781/1‑150 (MQ=255)
 cATACCCTGCTTTTGCAACGATGCAACGCGCTTCGCCTTTTTCGCTGTCCAGCAGCCAGAAAGAGTTCAGCAGATCGTCGTTAGCGGCTTTAGTAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGTTc                                                                                                                                                    <  2:3075520/150‑1 (MQ=255)
 cATAACCTGCTTTTGCAACGATGCAACGCGCTTCGCCTTTTTCGCTGTCCAGCAGCCAGAAAGAGTTCAGCAGATCGTCGTTAGCGGCTTTAGTAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGTTc                                                                                                                                                    >  1:1890974/1‑150 (MQ=255)
 cATAACCTGCTTTTGCAACGATGCAACGCGCTTCGCCTTTTTCGCTGTCCAGCAGCCAGAAAGAGTTCAGCAGATCGTCGTTAGCGGCTTTAGTAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGTTc                                                                                                                                                    <  2:2500586/150‑1 (MQ=255)
  aTAACCTGCTTTTGCAACGATGCAACGCGCTTCGCCTTTTTCGCTGTCCAGCAGCCAGAAAGAGTTCAGCAGATCGTCGTTAGCGGCTTTAGTAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGTTCa                                                                                                                                                   <  1:4370802/150‑1 (MQ=255)
       cTGCTTTTGCAACGATGCAACGCGCTTCGCCTTTTTCGCTGTCCAGCAGCCAGAAAGAGTTCAGCAGATCGTCGTTAGCGGCTTTAGTAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGATCAGCCTg                                                                                                                                              <  2:1850384/150‑1 (MQ=255)
                cAACGATGCAACGCGCTTCGCCTTTTTCGCTGTCCAGCAGCCAGAAAGAGTTCAGCAGATCGTCGTTAGCGGCTTTAGTAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCa                                                                                                                                     >  2:1478898/1‑150 (MQ=255)
                   cGATGCAACGCGCTTCGCCTTTTTCGCTGTCCAGCAGCCAGAAAGAGTTCAGCAGATCGTCGTTAGCGGCTTTAGTAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAAtt                                                                                                                                  <  1:1038535/150‑1 (MQ=255)
                    gATGCAACGCGCTTCGCCTTTTTCGCTGTCCAGCAGCCAGAAAGAGTTCAGCAGATCGTCGTTAGCGGCTTTAGTAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGATCAGCCTGTCGCGGCCAATTg                                                                                                                                 <  2:468164/150‑1 (MQ=255)
                    gATGCAACGCGCTTCGCCTTTTTCGCTGTCCAGCAGCCAGAAAGAGTTCAGCAGATCGTCGTTAGCGGCTTTAGTAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGATCAGCCTGTCGCGGCCAATTg                                                                                                                                 <  2:572776/150‑1 (MQ=255)
                       gCAACGCGCTTCGCCTTTTTCGCTGTCCAGCAGCCAGAAAGAGTTCAGCAGATCGTCGTTAGCGGCTTTAGTAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTa                                                                                                                              >  2:2533425/1‑150 (MQ=255)
                         aaCGCGCTTCGCCTTTTTCGCTGTCCAGCAGCCAGAAAGAGTTCAGCAGATCGTCGTTAGCGGCTTTAGTAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTaaa                                                                                                                            <  1:1642608/150‑1 (MQ=255)
                         aaCGCGCGTCGCCTTTTTCGCTGTCCAGCAGCCAGAAAGAGTTCAGCAGATCGTCGTTAGCGGCTTTAGTAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGATCAGCCTGTCGCGGCCAATTGGTaaa                                                                                                                            <  2:4615434/150‑1 (MQ=255)
                          aCGCGCTTCGCCTTTTTCGCTGTCCAGCAGCCAGAAAGAGTTCAGCAGATCGTCGTTAGCGGCTTTAGTAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTaaaa                                                                                                                           <  1:2977693/150‑1 (MQ=255)
                           cgcgCTTCGCCTTTTTCGCTGTCCAGCAGCCAGAAAGAGTTCAGCAGATCGTCGTTAGCGGCTTTAGTAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAAc                                                                                                                          <  2:2018817/150‑1 (MQ=255)
                            gcgcTTCGCCTTTTTCGCTGTCCAGCAGCCAGAAAGAGTTCAGCAGATCGTCGTTAGCGGCTTTAGTAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGATCAGCCTGTCGCGGCCAATTGGTAAAAcc                                                                                                                         <  2:109597/150‑1 (MQ=255)
                             cgcTTCGCCTTTTTCGCTGTCCAGCAGCCAGAAAGAGTTCAGCAGATCGTCGTTAGCGGCTTTAGTAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCa                                                                                                                        >  2:825779/1‑150 (MQ=255)
                               cTTCGCCTTTTTCGCTGTCCAGCAGCCAGAAAGAGTTCAGCAGATCGTCGTTAGCGGCTTTAGTAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCAtt                                                                                                                      <  1:1039369/150‑1 (MQ=255)
                               cTTCGCCTTTTTCGCTGTCCAGCAGCCAGAAAGAGTTCAGCAGATCGTCGTTAGCGGCTTTAGTAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCAtt                                                                                                                      >  1:4134974/1‑150 (MQ=255)
                                ttCGCCTTTTTCGCTGTCCAGCAGCCAGAAAGAGTTCAGCAGATCGTCGTTAGCGGCTTTAGTAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCAttg                                                                                                                     >  2:3524889/1‑150 (MQ=255)
                                ttCGCCTTTTTCGCTGTCCAGCAGCCAGAAAGAGTTCAGCAGATCGTCGTTAGCGGCTTTAGTAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCAttg                                                                                                                     <  2:2350419/150‑1 (MQ=255)
                                ttCGCCTTTTTCGCTGTCCAGCAGCCAGAAAGAGTTCAGCAGATCGTCGTTAGCGGCTTTAGTAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCAttg                                                                                                                     <  2:1455186/150‑1 (MQ=255)
                                 tCGCCTTTTTCGCTGTCCAGCAGCCAGAAAGAGTTCAGCAGATCGTCGTTAGCGGCTTTAGTAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCAttgt                                                                                                                    >  1:3936084/1‑150 (MQ=255)
                                     cTTTTTCGCTGTCCAGCAGCCAGAAAGAGTTCAGCAGATCGTCGTTAGCGGCTTTAGTAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCt                                                                                                                <  2:1824041/150‑1 (MQ=255)
                                      tttttCGCTGTCCAGCAGCCAGAAAGAGTTCAGCAGATCGTCGTTAGCGGCTTTAGTAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCtt                                                                                                               <  2:4514636/150‑1 (MQ=255)
                                      tttttCGCTGTCCAGCAGCCAGAAAGAGTTCAGCAGATCGTCGTTAGCGGCTTTAGTAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCtt                                                                                                               >  2:3731421/1‑150 (MQ=255)
                                        tttCGCTGTCCAGCAGCCAGAAAGAGTTCAGCAGATCGTCGTTAGCGGCTTTAGTAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGATCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGa                                                                                                             >  2:866067/1‑150 (MQ=255)
                                         ttCGCTGTCCAGCAGCCAGAAAGAGTTCAGCAGATCGTCGTTAGCGGCTTTAGTAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAg                                                                                                            >  1:2907697/1‑150 (MQ=255)
                                         ttCGCTGTCCAGCAGCCAGAAAGAGTTCAGCAGATCGTCGTTAGCGGCTTTAGTAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGATCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAg                                                                                                            <  2:2847622/150‑1 (MQ=255)
                                             cTGTCCAGCAGCCAGAAAGAGTTCAGCAGATCGTCGTTAGCGGCTTTAGTAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTGta                                                                                                        >  2:2739286/1‑150 (MQ=255)
                                              tGTCCAGCAGCCAGAAAGAGTTCAGCAGATCGTCGTTAGCGGCTTTAGTAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTGtat                                                                                                       <  1:2659770/150‑1 (MQ=255)
                                              tGTCCAGCAGCCAGAAAGAGTTCAGCAGATCGTCGTTAGCGGCTTTAGTAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGATCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTGtat                                                                                                       >  2:545845/1‑150 (MQ=255)
                                               gTCCAGCAGCCAGAAAGAGTTCAGCAGATCGTCGTTAGCGGCTTTAGTAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTGtata                                                                                                      >  1:2823900/1‑150 (MQ=255)
                                               gTCCAGCAGCCAGAAAGAGTTCAGCAGATCGTCGTTAGCGGCTTTAGTAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTGtata                                                                                                      <  1:1776519/150‑1 (MQ=255)
                                                  cagcagCCAGAAAGAGTTCAGCAGATCGTCGTTAGCGGCTTTAGTAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTGtatatat                                                                                                   >  2:2186329/1‑150 (MQ=255)
                                                  cagcagCCAGAAAGAGTTCAGCAGATCGTCGTTAGCGGCTTTAGTAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTGtatatat                                                                                                   >  2:11973/1‑150 (MQ=255)
                                                     cagcCAGAAAGAGTTCAGCAGATCGTCGTTAGCGGCTTTAGTAATCTGGATACCTGTAATCATGTGTTGCCTCCCCGGCAAAATTATTTGATTTGATCAGCCTGTCGCGGCCAATTGGTACAACCATTGTTGCTTGAGTGTATTTTTAct                                                                                                <  1:82819/150‑1 (MQ=255)
                                                         cAGAAAGAGTTCAGCAGATCGTCGTTAGCGGCTTTAGTAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGATCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTGTATATATActcctc                                                                                            >  1:3896469/1‑150 (MQ=255)
                                                          aGAAAGAGTTCAGCAGATCGTCGTTAGCGGCTTTAGTAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTGTATATATACTCCTCa                                                                                           >  2:2968591/1‑150 (MQ=255)
                                                          aGAAAGAGTTCAGCAGATCGTCGTTAGCGGCTTTAGTAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGATCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTGTATATATACTCCTCa                                                                                           <  1:3545039/150‑1 (MQ=255)
                                                           gAAAGAGTTCAGCAGATCGTCGTTAGCGGCTTTAGTAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTGTATATATACTCCTCaa                                                                                          >  1:3353535/1‑150 (MQ=255)
                                                           gAAAGAGTTCAGCAGATCGTCGTTAGCGGCTTTAGTAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTGTATATATACTCCTCaa                                                                                          <  1:4567318/150‑1 (MQ=255)
                                                           gAAAGAGTTCAGCAGATCGTCGTTAGCGGCTTTAGTAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGATCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTGTATATATACTCCTCaa                                                                                          >  2:2256175/1‑150 (MQ=255)
                                                              agagTTCAGCAGATCGTCGTTAGCGGCTTTAGTAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTGTATATATACTCCTCAAaca                                                                                       >  2:3948298/1‑150 (MQ=255)
                                                                agTTCAGCAGATCGTCGTTAGCGGCTTTAGTAATCTGGATACCTGTAATCATGTGACGCCTCCCCGGCAAAATTATTTGATTTGATCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTGTATATATACTCCTCAAACAcc                                                                                     >  1:3336881/1‑150 (MQ=255)
                                                                  ttCAGCAGATCGTCGTTAGCGGCTTTAGTAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTGTATATATACTCCTCAAACACCCt                                                                                   <  1:3760573/150‑1 (MQ=255)
                                                                  ttCAGCAGATCGTCGTTAGCGGCTTTAGTAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGATCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTGTATATATACTCCTCAAACACCCt                                                                                   <  2:1801286/150‑1 (MQ=255)
                                                                   tCAGCAGATCGTCGTTAGCGGCTTTAGTAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTGTATATATACTCCTCAAACACCCtt                                                                                  <  2:1318046/150‑1 (MQ=255)
                                                                     agcagATCGTCGTTAGCGGCTTTAGTAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTGTATATATACTCCTCAAACACCCTTGa                                                                                >  2:1711298/1‑150 (MQ=255)
                                                                     agcagATCGTCGTTAGCGGCTTTAGTAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTGTATATATACTCCTCAAACACCCTTGa                                                                                <  2:3737949/150‑1 (MQ=255)
                                                                        agATCGTCGTTAGCGGCTTTAGTAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTGTATATATACTCCTCAAACACCCTTGAATc                                                                             >  1:3297438/1‑150 (MQ=255)
                                                                        agATCGTCGTTAGCGGCTTTAGTAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTGTATATATACTCCTCAAACACCCTTGAATc                                                                             >  2:2619010/1‑150 (MQ=255)
                                                                         gATCGTCGTTAGCGGCTTTAGTAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGATCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTGTATATATACTCCTCAAACACCCTTGAATCt                                                                            <  2:71111/150‑1 (MQ=255)
                                                                            cgtcgtTAGCGGCTTTAGTAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTGTATATATACTCCTCAAACACCCTTGAATCTTTg                                                                         >  2:187352/1‑150 (MQ=255)
                                                                            cgtcgtTAGCGGCTTTAGTAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTGTATATATACTCCTCAAACACCCTTGAATCTTTg                                                                         >  1:1969510/1‑150 (MQ=255)
                                                                             gtcgtTAGCGGCTTTAGTAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTGTATATATACTCCTCAAACACCCTTGAATCTTTGa                                                                        <  1:1391104/150‑1 (MQ=255)
                                                                             gtcgtTAGCGGCTTTAGTAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGATCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTGTATATATACTCCTCAAACACCCTTGAATCTTTGa                                                                        >  2:3147225/1‑150 (MQ=255)
                                                                               cgtTAGCGGCTTTAGTAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTGTATATATACTCCTCAAACACCCTTGAATCTTTGAtt                                                                      <  2:3833252/150‑1 (MQ=255)
                                                                               cgtTAGCGGCTTTAGTAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGATCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTGTATATATACTCCTCAAACACCCTTGAATCTTTGAtt                                                                      <  1:4610409/150‑1 (MQ=255)
                                                                                gtTAGCGGCTTTAGTAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGATCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTGTATATATACTCCTCAAACACCCTTGAATCTTTGAttt                                                                     <  1:2983047/150‑1 (MQ=255)
                                                                                gtTAGCGGCTTTAGTAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGATCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTGTATATATACTCCTCAAACACCCTTGAATCTTTGAttt                                                                     <  1:2741000/150‑1 (MQ=255)
                                                                                 ttAGCGGCTTTAGTAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTGTATATATACTCCTCAAACACCCTTGAATCTTTGATTTa                                                                    >  2:3400796/1‑150 (MQ=255)
                                                                                  tAGCGGCTTTAGTAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGATCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTGTATATATACTCCTCAAACACCCTTGAATCTTTGATTTaa                                                                   >  2:2397903/1‑150 (MQ=255)
                                                                                   aGCGGCTTTAGTAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTGTATATATACTCCTCAAACACCCTTGAATCTTTGATTTaaa                                                                  >  1:2100166/1‑150 (MQ=255)
                                                                                          ttAGTAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGATCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTGTATATATACTCCTCAAACACCCTTGAATCTTTGATTTAAATCAATaa                                                           >  1:408906/1‑150 (MQ=255)
                                                                                              tAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGATCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTGTATATATACTCCTCAAACACCCTTGAATCTTTGATTTAAATCAATAAAAAc                                                       >  2:3280988/1‑150 (MQ=255)
                                                                                               aaTCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGATCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTGTATATATACTCCTCAAACACCCTTGAATCTTTGATTTAAATCAATAAAAAcc                                                      <  2:406773/150‑1 (MQ=255)
                                                                                                  cTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGATCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTGTATATATACTCCTCAAACACCCTTGAATCTTTGATTTAAATCAATAAAAACcaca                                                   <  1:4196917/150‑1 (MQ=255)
                                                                                                   tGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTGTATATATACTCCTCAAACACCCTTGAATCTTTGATTTAAATCAATAAAAACcacac                                                  <  1:696725/150‑1 (MQ=255)
                                                                                                   tGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTGTATATATACTCCTCAAACACCCTTGAATCTTTGATTTAAATCAATAAAAACcacac                                                  <  1:1394750/150‑1 (MQ=255)
                                                                                                     gATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTGTATATATACTCCTCAAACACCCTTGAATCTTTGATTTAAATCAATAAAAACCACACAt                                                <  1:515326/150‑1 (MQ=255)
                                                                                                     gATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTGTATATATACTCCTCAAACACCCTTGAATCTTTGATTTAAATCAATAAAAACCACACAt                                                <  1:3342057/150‑1 (MQ=255)
                                                                                                         ccTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTGTATATATACTCCTCAAACACCCTTGAATCTTTGATTTAAATCAATAAAAACCACACATCAAg                                            <  1:3814343/150‑1 (MQ=255)
                                                                                                           tGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTGTATATATACTCCTCAAACACCCTTGAATCTTTGATTTAAATCAATAAAAACCACACATCAAGTa                                          >  1:3705829/1‑150 (MQ=255)
                                                                                                             tAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGACGCGGCCAATAGGTAAAACCATTGTTGCTAGAGTGTATATATACTCCTCAAACACCCTTGAATCTTTGATTTAAATCAATAAAAACCACACATCAAGAATg                                        >  2:3420671/1‑150 (MQ=255)
                                                                                                               aTCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGATCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTGTATATATACTCCTCAAACACCCTTGAATCTTTGATTTAAATCAATAAAAACCACACATCAAGTATGGt                                      <  1:2439573/150‑1 (MQ=255)
                                                                                                                    gtgATGCCTCCCCGGCAAAATTATTTGATTTGATCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTGTATATATACTCCTCAAACACCCTTGAATCTTTGATTTAAATCAATAAAAACCACACATCAAGTATGGTCGCaa                                 >  2:241146/1‑150 (MQ=255)
                                                                                                                      gATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTGTATATATACTCCTCAAACACCCTTGAATCTTTGATTTAAATCAATAAAAACCACACATCAAGTATGGTCGCAAAt                               >  1:1872690/1‑150 (MQ=255)
                                                                                                                      gATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTGTATATATACTCCTCAAACACCCTTGAATCTTTGATTTAAATCAATAAAAACCACACATCAAGTATGGTCGCAAAt                               <  1:3663858/150‑1 (MQ=255)
                                                                                                                         gCCTCCCCGGCAAAATTATTTGATTTGATCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTGTATATATACTCCTCAAACACCCTTGAATCTTTGATTTAAATCAATAAAAACCACACATCAAGTATGGTCGCAAATGGa                            >  1:572040/1‑150 (MQ=255)
                                                                                                                          ccTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTGTATATATACTCCTCAAACACCCTTGAATCTTTGATTTAAATCAATAAAAACCACACATCAAGTATGGTCGCAAATGGAt                           >  1:3024810/1‑150 (MQ=255)
                                                                                                                          ccTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTGTATATATACTCCTCAAACACCCTTGAATCTTTGATTTAAATCAATAAAAACCACACATCAAGTATGGTCGCAAATGGAt                           >  1:4310805/1‑150 (MQ=255)
                                                                                                                           cTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTGTATATATACTCCTCAAACACCCTTGAATCTTTGATTTAAATCAATAAAAACCACACATCAAGTATGGTCGCAAATGGAtt                          >  1:3618156/1‑150 (MQ=255)
                                                                                                                             ccccGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTGTATATATACTCCTCAAACACCCTTGAATCTTTGATTTAAATCAATAAAAACCACACATCAAGTATGGTCGCAAATGGAtttt                        >  2:3710953/1‑150 (MQ=255)
                                                                                                                             ccccGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTGTATATATACTCCTCAAACACCCTTGAATCTTTGATTTAAATCAATAAAAACCACACATCAAGTATGGTCGCAAATGGAtttt                        <  1:3638519/150‑1 (MQ=255)
                                                                                                                              cccGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATGGGTAAAACCATTGTTGCTTGAGTGTATATATACTCCTCAAACACCCTTGAATCTTTGATTTAAATCAATAAAAACCACACATCAAGTATGGTCGCAAATGGATTTTa                       <  2:4015116/150‑1 (MQ=255)
                                                                                                                                 ggCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTGTATATATACTCCTCAAACACCCTTGAATCTTTGATTTAAATCAATAAAAACCACACATCAAGTATGGTCGCAAATGGATTTTATTg                    >  1:1412381/1‑150 (MQ=255)
                                                                                                                                    aaaaTTATTTGATTTGATCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTGTATATATACTCCTCAAACACCCTTGAATCTTTGATTTAAATCAATAAAAACCACACATCAAGTATGGTCGCAAATGGATTTTATTGttt                 >  1:2270090/1‑150 (MQ=255)
                                                                                                                                        ttattTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTGTATATATACTCCTCAAACACCCTTGAATCTTTGATTTAAATCAATAAAAACCACACATCAAGTATGGTCGCAAATGGATTTTATTGTTTTACa             <  2:4579424/150‑1 (MQ=255)
                                                                                                                                        ttattTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTGTATATATACTCCTCAAACACCCTTGAATCTTTGATTTAAATCAATAAAAACCACACATCAAGTATGGTCGCAAATGGATTTTATTGTTTTACa             <  2:1814939/150‑1 (MQ=255)
                                                                                                                                        ttattTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTGTATATATACTCCTCAAACACCCTTGAATCTTTGATTTAAATCAATAAAAACCACACATCAAGTATGGTCGCAAATGGATTTTATTGTTTTACa             >  2:3915583/1‑150 (MQ=255)
                                                                                                                                         tattTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTGTATATATACTCCTCAAACACCCTTGAATCTTTGATTTAAATCAATAAAAACCACACATCAAGTATGGTCGCAAATGGATTTTATTGTTTTACAt            <  2:4177117/150‑1 (MQ=255)
                                                                                                                                         tattTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTGTATATATACTCCTCAAACACCCTTGAATCTTTGATTTAAATCAATAAAAACCACACATCAAGTATGGTCGCAAATGGATTTTATTGTTTTACAt            <  2:541916/150‑1 (MQ=255)
                                                                                                                                         tattTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTGTATATATACTCCTCAAACACCCTTGAATCTTTGATTTAAATCAATAAAAACCACACATCAAGTATGGTCGCAAATGGATTTTATTGTTTTACAt            <  1:3927422/150‑1 (MQ=255)
                                                                                                                                          atttgatttgTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTGTATATATACTCCTCAAACACCCTTGAATCTTTGATTTAAATCAATAAAAACCACACATCAAGTATGGTCGCAAATGGATTTTATTGTTTTACATc           <  1:228193/150‑1 (MQ=255)
                                                                                                                                            ttgatttgatCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTGTATATATACTCCTCAAACACCCTTGAATCTTTGATTTAAATCAATAAAAACCACACATCAAGTATGGTCGCAAATGGATTTTATTGTTTTACATCaa         <  1:1753208/150‑1 (MQ=255)
                                                                                                                                               atttgATCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTGTATATATACTCCTCAAACACCCTTGAATCTTTGATTTAAATCAATAAAAACCACACATCAAGTATGGTCGCAAATGGATTTTATTGTTTTACATCAACtt      >  1:2603014/1‑150 (MQ=255)
                                                                                                                                                tttgTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTGTATATATACTCCTCAAACACCCTTGAATCTTTGATTTAAATCAATAAAAACCACACATCAAGTATGGTCGCAAATGGATTTTATTGTTTTACATCAACTTa     <  2:668820/150‑1 (MQ=255)
                                                                                                                                                   gTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTGTATATATACTCCTCAAACACCCTTGAATCTTTGATTTAAATCAATAAAAACCACACATCAAGTATGGTCGCAAATGGATTTTATTGTTTTACATCAACTTATGc  >  1:1813081/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCATAACCTGCTTTTGCAACGATGCAACGCGCTTCGCCTTTTTCGCTGTCCAGCAGCCAGAAAGAGTTCAGCAGATCGTCGTTAGCGGCTTTAGTAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTGTATATATACTCCTCAAACACCCTTGAATCTTTGATTTAAATCAATAAAAACCACACATCAAGTATGGTCGCAAATGGATTTTATTGTTTTACATCAACTTATGC  >  CP009273/2709693‑2709989

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: