Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 819,862 G→T 12.2% P31Q (CCG→CAG)  ybhP ← endo/exonuclease/phosphatase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273819,8620GT12.2% 1262.9 / 104.4 370P31Q (CCG→CAG) ybhPendo/exonuclease/phosphatase family protein
Reads supporting (aligned to +/- strand):  ref base G (144/181);  new base T (13/32);  total (157/213)
Fisher's exact test for biased strand distribution p-value = 5.44e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

Only 100 of 374 total aligned reads displayed.

AATCGCTCCACATAGTGTCGGCGAGAAACTCGTAGTGCGAGGTATCGGGCCAGTTTTCCACATGCAGCGGATGAACTTCGTGCGCGCCCATCACTTCCTGCAGGCAAACAATATCGGCGCTGACGGTACGCACGGCGTCGCGAAGTTCCGGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAA  >  CP009273/819713‑820010
                                                                                                                                                     |                                                                                                                                                    
aaTCGCTCCACATAGTGTCGGCGAGAAACTCGTAGTGCGAGGTATCGGGCCAGTTTTCCACATGCAGCGGATGAACTTCGTGCGCGCCCATCACTTCCTGCAGGCAAACAATATCGGCGCTGACGGTACGCACGGCGTCGCGAAGTTCCg                                                                                                                                                      >  1:4164119/1‑150 (MQ=255)
   cGCTCCACATAGTGTCGGCGAGAAACTCGTAGTGCGAGGTATCGGGCCAGTTTTCCACATGCAGCGGATGAACTTCGTGCGCGCCCATCACTTCCTGCAGGCAAACAATATCGGCGCTGACGGTACGCACGGCGTCGCGAAGTTCCGGCa                                                                                                                                                   <  1:71139/150‑1 (MQ=255)
   cGCTCCACATAGTGTCGGCGAGAAACTCGTAGTGCGAGGTATCGGGCCAGTTTTCCACATGCAGCGGATGAACTTCGTGCGCGCCCATCACTTCCTGCAGGCAAACAATATCGGCGCTGACGGTACGCACGGCGTCGCGAAGTTCCGGCa                                                                                                                                                   >  1:402920/1‑150 (MQ=255)
       ccagttAGTGTCGGCGAGAAACTAGTAGTGCGAGGTATCGGGCCAGTTTTCAACATGCTGCGGATGAACTTCGTGAGCGCCCATCACTTCCTGCAGGCAAACAATTTCGGCGATGACGGTACGCACGGCGTCGCGAAGTTCCGGCAAAAt                                                                                                                                               <  1:1723964/145‑1 (MQ=255)
       ccACATAGTGTCGGCGAGAAACTCGTAGTGCGAGGTATCGGGCCAGTTTTCCACATGCAGCGGATGAACTTCGTGCGCGCCCATCACTTCCTGCAGGCAAACAATATCGGCGCTGACGGTACGCACGGCGTCGCGAAGTTCCGGCAAAAt                                                                                                                                               <  1:874613/150‑1 (MQ=255)
        cacaTAGTGTCGGCGAGAAACTCGTAGTGCGAGGTATCGGGCCAGTTTTCCACATGCAGCGGATGAACTTCGTGCGCGCCCATCACTTCCTGCAGGCAAACAATATCGGCGCTGACGGTACGCACGGCGTCGCGAAGTTCCGGCAAAATg                                                                                                                                              <  2:1893585/150‑1 (MQ=255)
             aGTGTCGGCGAGAAACTCGTAGTGCGAGGTATCGGGCCAGTTTTCCACATGCAGCGGATGAACTTCGTGCGCGCCCATCACTTCCTGCAGGCAAACAATATCGGCGCTGACGGTACGCACGGCGTCGCGAAGTTCCGGCAAAATGAAGCg                                                                                                                                         >  1:3590520/1‑150 (MQ=255)
              gtgtCGGCGAGAAACTCGTAGTGCGAGGTATCGGGCCAGTTTTCCACATGCAGCGGATGAACTTCGTGCGCGCCCATCACTTCCTGCAGGCAAACAATATCGGCGCTGACGGTACGCACGGCGTCGCGAAGTTCCGGCAAAATGAAGCGt                                                                                                                                        <  1:3857630/150‑1 (MQ=255)
               tgtCGGCGAGAAACTCGTAGTGCGAGGTATCGGGCCAGTTTTCCACATGCAGCGGATGAACTTCGTGCGCGCCCATCACTTCCTGCAGGCAAACAATATCGGCGCTGACGGTACGCACGGCGTCGCGAAGTTCCGGCAAAATGAAGCGTc                                                                                                                                       >  2:3561648/1‑150 (MQ=255)
               tgtCGGCGAGAAACTCGTAGTGCGAGGTATCGGGCCAGTTTTCCACATGCAGCGGATGAACTTCGTGCGCGCCCATCACTTCCTGCAGGCAAACAATATCGGCGCTGACGGTACGCACGGCGTCGCGAAGTTCCGGCAAAATGAAGCGTc                                                                                                                                       >  2:2817758/1‑150 (MQ=255)
                  cGGCGAGAAACTCGTAGTGCGAGGTATCGGGCCAGTTTTCCACATGCAGCGGATGAACTTCGTGCGCGCCCATCACTTCCTGCAGGCAAACAATATCGGCGCTGACGGTACGCACGGCGTCGCGAAGTTCCGGCAAAATGAAGCGTCGGt                                                                                                                                    >  2:667921/1‑150 (MQ=255)
                   ggCGAGAAACTCGTAGTGCGAGGTATCGGGCCAGTTTTCCACATGCAGCGGATGAACTTCGTGCGCGCCCATCACTTCCTGCAGGCAAACAATATCGGCGCTGACGGTACGCACGGCGTCGCGAAGTTCCGGCAAAATGAAGCGTCGGtt                                                                                                                                   <  2:4491715/150‑1 (MQ=255)
                      gagaAACTCGTAGTGCGAGGTATCGGGCCAGTTTTCCACATGCAGCGGATGAACTTCGTGCGCGCCCATCACTTCCTGCAGGCAAACAATATCGGCGCTGACGGTACGCACGGCGTCGCGAAGTTCCGGCAAAATGAAGCGTCGGTTaaa                                                                                                                                >  1:1020558/1‑150 (MQ=255)
                       agaAACTCGTAGTGCGAGGTATCGGGCCAGTTTTCCCCATGCAGCGGATGCACTTCGTGCGCGCCCATCACTTCCTGCAGGCAAACAATATCGGCGCTGACGGTACGCACGGCGTCGCGAAGTTCCGGCAAAATGAAGCGTCGGTTAAAc                                                                                                                               <  1:4932855/150‑1 (MQ=255)
                        gaAACTCGTAGTGCGAGGTATCGGGCCAGTTTTCCACATGCAGCGGATGAACTTCGTGCGCGCCCATCACTTCCTGCAGGCAAACAATATCGGCGCTGACGGTACGCACGGCGTCGCGAAGTTCCGGCAAAATGAAGCGTCGGTTAAAcg                                                                                                                              <  1:2487584/150‑1 (MQ=255)
                        gaAACTCGTAGTGCGAGGTATCGGGCCAGTTTTCCACATGCAGCGGATGAACTTCGTGCGCGCCCATCACTTCCTGCAGGCAAACAATATCGGCGCTGACGGTACGCACGGCGTCGCGAAGTTCCGGCAAAATGAAGCGTCGGTTAAAcg                                                                                                                              >  2:5075035/1‑150 (MQ=255)
                          aaCTCGTAGTGCGAGGTATCGGGCCAGTTTTCCACATGCAGCGGATGAACTTCGTGCGCGCCCATCACTTCCTGCAGGCAAACAATATCGGCGCTGACGGTACGCACGGCGTCGCGAAGTTCCTGCAAAATGAAGCGTCGGTTAAAcgcg                                                                                                                            >  1:2557230/1‑150 (MQ=255)
                            cTCGTAGTGCGAGGTATCGGGCCAGTTTTCCACATGCAGCGGATGAACTTCGTGCGCGCCCATCACTTCCTGCAGGCAAACAATATCGGCGCTGACGGTACGCACGGCGTCGCGAAGTTCCGGCAAAATGAAGCGTCGGTTAAACGCGGt                                                                                                                          <  1:2817758/150‑1 (MQ=255)
                            cTCGTAGTGCGAGGTATCGGGCCAGTTTTCCACATGCAGCGGATGAACTTCGTGCGCGCCCATCACTTCCTGCAGGCAAACAATATCGGCGCTGACGGTACGCACGGCGTCGCGAAGTTCCGGCAAAATGAAGCGTCGGTTAAACGCGGt                                                                                                                          <  2:1453106/150‑1 (MQ=255)
                                   tGCGAGGTATCGGGCCAGTTTTCCACATGCAGCGGATGAACTTCGTGCGCGCCCATCACTTCCTGCAGGCAAACAATATCGGCGCTGACGGTACGCACGGCGTCGCGAAGTTCCGGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCt                                                                                                                   >  1:2073028/1‑150 (MQ=255)
                                       aGGTATCGGGCCAGTTTTCCACATGCAGCGGATGAACTTCGTGCGCGCCCATCACTTCCTGCAGGCAAACAATATCGGCGCTGACGGTACGCACGGCGTCGCGAAGTTCCGGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTtgt                                                                                                               <  2:3592728/150‑1 (MQ=255)
                                        ggTATCGGGCCAGTTTTCCACATGCAGCGGATGAACTTCGTGCGCGCCCATCACTTCCTGCAGGCAAACAATATCGGCGCTGACGGTACGCACGGCGTCGCGAAGTTCCGGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTtgtg                                                                                                              >  2:3738490/1‑150 (MQ=255)
                                           aTCGGGCCAGTTTTCCACATGCAGCGGATGAACTTCGTGCGCGCCCATCACTTCCTGCAGGCAAACAATATCGGCGCTGACGGTACGCACGGCGTCGCGAAGTTCCTGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAat                                                                                                           >  1:1979734/1‑150 (MQ=255)
                                            tCGGGCCAGTTTTCCACATGCAGCGGATGAACTTCGTGCGCGCCCATCACTTCCTGCAGGCAAACAATATCGGCGCTGACGGTACGCACGGCGTCGCGAAGTTCCGGCAAAATGAAGCGTCGGTTAAACGCGGGAAAGCCTTTGTGAata                                                                                                          >  2:1409799/1‑150 (MQ=255)
                                                  cAGTTTTCCACATGCAGCGGATGAACTTCGTGCGCGCCCATCACTTCCTGCAGGCAAACAATAACGGCGCTGACGGTACGCACGGCGTCGCGAAGTTCCGGCAAAATGAATCGTCGGTTAAACGCGGTAAAGCCTTGTGGAATATTGATg                                                                                                    >  2:3101413/1‑150 (MQ=255)
                                                   aGTTTTCCACATGCAGCGGATGAACTTCGTGCGCGCCCATCACTTCCTGCAGGCAAACAATATCGGCGCTGACGGTACGCACGGCGTCGCGAAGTTCCGGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAATATTGATgg                                                                                                   >  1:3790312/1‑150 (MQ=255)
                                                    gTTTTCCACATGCAGCGGATGAACTTCGTGCGCGCCCATCACTTCCTGCAGGCAAACAATATCGGCGCTGACGGTACGCACGGCGTCGCGAAGTTCCGGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGt                                                                                                  >  2:4210706/1‑150 (MQ=255)
                                                    gTTTTCCACATGCAGCGGATGAACTTCGTGCGCGCCCATCACTTCCTGCAGGCAAACAATATCGGCGCTGACGGTACGCACGGCGTCGCGAAGTTCCGGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGt                                                                                                  <  2:1948625/150‑1 (MQ=255)
                                                          cacaTGCAGCGGATGAACTTCGTGCGCGCCCATCACTTCCTGCAGGCAAACAATATCGGCGCTGACGGTACGCACGGCGTCGCGAAGTTCCTGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCAc                                                                                            <  1:3713419/150‑1 (MQ=255)
                                                          cacaTGCAGCGGATGAACTTCGTGCGCGCCCATCACTTCCTGCAGGCAAACAATATCGGCGCTGACGGTACGCACGGCGTCGCGAAGTTCCGGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCAc                                                                                            <  1:3161380/150‑1 (MQ=255)
                                                          cacaTGCAGCGGATGAACGTCGTGCGCGCCCATCACTTCCTGCAGGCAACCACTAGCGGCGCTGCCGGTACGCACGGCGTCGCGAAGTTCCGGCAACATGCAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGACTATTGATGGTGAGCAc                                                                                            <  2:2717715/150‑1 (MQ=255)
                                                           acaTGCAGCGGATGAACTTCGTGCGCGCCCATCACTTCCTGCAGGCAAACAATATCGGCGCTGACGGTACGCACGGCGTCGCGAAGTTCCGGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCAcc                                                                                           <  2:396766/150‑1 (MQ=255)
                                                            caTGCAGCGGATGAACTTCGTGCGCGCCCATCACTTCCTGCAGGCAAACAATATCGGCGCTGACGGTACGCACGGCGTCGCGAAGTTCCGGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCt                                                                                          <  2:1229632/150‑1 (MQ=255)
                                                             aTGCAGCGGATGAACTTCGTGCGCGCCCATCACTTCCTGCAGGCAAACAATATCGGCGCTGACGGTACGCACGGCGTCGCGAAGTTCCGGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCtt                                                                                         <  1:4638666/150‑1 (MQ=255)
                                                             aTGCAGCGGATGAACTTCGTGCGCGCCCATCACTTCCTGCAGGCAAACAATATCGGCGCTGACGGTACGCACGGCGTCGCGAAGTTCCGGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCtt                                                                                         <  1:1387173/150‑1 (MQ=255)
                                                               gCAGCGGATGAACTTCGTGCGCGCCCATCACTTCCTGCAGGCAAACAATATCGGCGCTGACGGTACGCACGGCGTCGCGAAGTTCCGGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGa                                                                                       >  2:104454/1‑150 (MQ=255)
                                                                 aGCGGATGAACTTCGTGCGCGCCCATCACTTCCTGCAGGCAAACAATATCGGCGCTGACGGTACGCACGGCGTCGCGAAGTTCCGGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAAc                                                                                     >  2:153482/1‑150 (MQ=255)
                                                                    ggATGAACTTCGTGCGCGCCCATCACTTCCTGCAGGCAAACAATATCGGCGCTGACGGTACGCACGGCGTCGCGAAGTTCCGGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGaa                                                                                  >  2:257290/1‑150 (MQ=255)
                                                                      aTGAACTTCGTGCGCGCCCATCACTTCCTGCAGGCAAACAATATCGGCGCTGACGGTACGCACGGCGTCGCGAAGTTCCGGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGaaaa                                                                                <  1:4226178/150‑1 (MQ=255)
                                                                        gAACTTCGTGCGCGCCCATCACTTCCTGCAGGCAAACAATATCGGCGCTGACGGTACGCACGGCGTCGCGAAGTTCCGGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAAtt                                                                              >  1:3562500/1‑150 (MQ=255)
                                                                         aaCTTCGTGCGCGCCCATCACTTCCTGCAGGCAAACAATATCGGCGCTGACGGTACGCACGGCGTCGCGAAGTTCCGGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAAttg                                                                             <  2:2175205/150‑1 (MQ=255)
                                                                         aaCTTCGTGCGCGCCCATCACTTCCTGCAGGCAAACAATATCGGCGCTGACGGTACGCACGGCGTCGCGAAGGTCCGGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAAttg                                                                             <  2:4675421/150‑1 (MQ=255)
                                                                            ttCGTGCGCGCCCATCACTTCCTGCAGGCAAACAATATCGGCGCTGACGGTACGCACGGCGTCGCGAAGTTCCGGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAAttgttg                                                                          >  2:4274341/1‑150 (MQ=255)
                                                                            ttCGTGCGCGCCCATCACTTCCTGCAGGCAAACAATATCGGCGCTGACGGTACGCACGGCGTCGCGAAGTTCCGGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAAttgtgg                                                                          >  2:4602156/1‑148 (MQ=255)
                                                                              cGTGCGCGCCCATCACTTCCTGCAGGCAAACAATATCGGCGCTGACGGTACGCACGGCGTCGCGAAGTTCCTGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTtgtg                                                                        <  1:895378/150‑1 (MQ=255)
                                                                              cGTGCGCGCCCATCACTTCCTGCAGGCAAACAATATCGGCGCTGACGGTACGCACGGCGTCGCGAAGTTCCGGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTtgtg                                                                        >  2:4954338/1‑150 (MQ=255)
                                                                                  cgcgcCCATCACTTCCTGCAGGCAAACAATATCGGCGCTGACGGTACGCACGGCGTCGCGAAGTTCCGGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTg                                                                    >  2:10921/1‑150 (MQ=255)
                                                                                   gcgcCCATCACTTCCTGCAGGCAAACAATATCGGCGCTGACGGTACGCACGGCGTCGCGAAGTTCCGGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGa                                                                   >  2:860112/1‑150 (MQ=255)
                                                                                    cgcCCATCACTTCCTGCAGGCAAACAATATCGGCGCTGACGGTACGCACGGCGTCGCGAAGTTCCGGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGAt                                                                  <  2:755413/150‑1 (MQ=255)
                                                                                    cgcCCATCACTTCCTGCAGGCAAACAATATCGGCGCTGACGGTACGCACGGCGTCGCGAAGTTCCGGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGAt                                                                  >  2:2947722/1‑150 (MQ=255)
                                                                                     gcCCATCACTTCCTGCAGGCAAACAATATCGGCGCTGACGGTACGCACGGCGTCGCGAAGTTCCGGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGACTATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATc                                                                 >  1:4885834/1‑150 (MQ=255)
                                                                                      cccATCACTTCCTGCAGGCAAACAATATCGGCGCTGACGGTACGCACGGCGTCGCGAAGTTCCGGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCg                                                                >  1:3501471/1‑150 (MQ=255)
                                                                                        cATCACTTCCTGCAGGCAAACAATATCGGCGCTGACGGTACGCACGGCGTCGCGAAGTTCCGGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCggg                                                              <  2:2068949/150‑1 (MQ=255)
                                                                                         aTCACTTCCTGCAGGCAAACAATATCGGCGCTGACGGTACGCACGGCGTCGCGAAGTTCCGGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGc                                                             <  2:4468873/150‑1 (MQ=255)
                                                                                           cACTTCCTGCAGGCAAACAATATCGGCGCTGACGGTACGCACGGCGTCGCGAAGTTCCGGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCAt                                                           >  1:325765/1‑150 (MQ=255)
                                                                                            aCTTCCTGCAGGCAAACAATATCGGCGCTGACGGTACGCACGGCGTCGCGAAGTTCCGGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATa                                                          >  1:458000/1‑150 (MQ=255)
                                                                                               tCCTGCAGGCAAACAATATCGGCGCTGACGGTACGCACGGCGTCGCGAAGTTCCGGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAAtt                                                       <  1:781255/150‑1 (MQ=255)
                                                                                                    cAGGCAAACAATATCGGCGCTGACGGTACGCACGGCGTCGCGAAGTTCCGGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCt                                                  >  1:4542142/1‑150 (MQ=255)
                                                                                                    cAGGCAAACAAGCTCGGCGCTGACGGTACGCACGGCGTCGCGAAGGTCCGGCAAAATGAAGCGTCGGTGAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCt                                                  <  2:3344019/150‑1 (MQ=255)
                                                                                                     aGGCAAACAATATCGGCGCTGACGGTACGCACGGCGTCGCGAAGTTCCGGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTg                                                 <  1:1263892/150‑1 (MQ=255)
                                                                                                     aGGCAAACAATATCGGCGCTGACGGTACGCACGGCGTCGCGAAGTTCCGGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTg                                                 <  2:4038064/150‑1 (MQ=255)
                                                                                                     aGGCAAACAATATCGGCGCTGACGGTACGCACGGCGTCGCGAAGTTCCGGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTg                                                 <  1:320877/150‑1 (MQ=255)
                                                                                                       gCAAACAATATCGGCGCTGACGGTACGCACGGCGTCGCGAAGTTCCGGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTc                                               <  1:510375/150‑1 (MQ=255)
                                                                                                         aaaCAATATCGGCGCTGACGGTACGCACGGCGTCGCGAAGTTCCGGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCtt                                             <  2:2575278/150‑1 (MQ=255)
                                                                                                          aaCAATATCGGCGCTGACGGTACGCACGGCGTCGCGAAGTTCCGGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCttt                                            <  1:2371939/150‑1 (MQ=255)
                                                                                                             aaTATCGGCGCTGACGGTACGCACGGCGTCGCGAAGTTCCGGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGcc                                         >  1:7388/1‑150 (MQ=255)
                                                                                                               tatCGGCGCTGACGGTACGCACGGCGTCGCGAAGTTCCGGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCctc                                       >  1:1895388/1‑150 (MQ=255)
                                                                                                                atCGGCGCTGACGGTACGCGCGGCGTCGCGAAGTTCCGGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCctct                                      >  2:4405767/1‑150 (MQ=255)
                                                                                                                atCGGCGCTGACGGTACGCACGGCGTCGCGAAGTTCCGGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCctct                                      <  1:1670829/150‑1 (MQ=255)
                                                                                                                 tCGGCGCTGACGGTACGCACGGCGTCGCGAAGTTCCGGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCtt                                     >  1:4385588/1‑150 (MQ=255)
                                                                                                                 tCGGCGCTGACGGTACGCACGGCGTCGCGAAGTTCCGGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCtt                                     <  1:585621/150‑1 (MQ=255)
                                                                                                                   ggCGCTGACGGTACGCACGGCGTCGCGAAGTTCCTGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTAt                                   <  2:2114947/150‑1 (MQ=255)
                                                                                                                      gcTGACGGTACGCACGGCGTCGCGAAGTTCCGGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTAtctc                                <  2:3456292/150‑1 (MQ=255)
                                                                                                                       cTGACGGTACGCAAGGCGTCGCGAAGTTCCTGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCa                               <  2:2234713/150‑1 (MQ=255)
                                                                                                                        tGACGGTACGCACGGCGTCGCGAAGTTCCGGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCAt                              <  2:3215415/150‑1 (MQ=255)
                                                                                                                         gACGGTACGCACGGCGTCGCGAAGTTCCGGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCAtt                             <  2:962809/150‑1 (MQ=255)
                                                                                                                           cGGTACGCACGGCGTCGCGAAGTTCCTGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGa                           >  2:553574/1‑150 (MQ=255)
                                                                                                                           cGGTACGCACGGCGTCGCGAAGTTCCGGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGa                           >  2:3038214/1‑150 (MQ=255)
                                                                                                                           cGGTACGCACGGCGTCGCGAAGTTCCGGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGa                           <  1:1252406/150‑1 (MQ=255)
                                                                                                                           cGGTACGCACGGCGTCGCGAAGTTCCGGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGa                           <  2:3317172/150‑1 (MQ=255)
                                                                                                                            ggTACGCACGGCGTCGCGAAGTTCCGGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGGTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCAGTGaa                          >  1:3896678/1‑150 (MQ=255)
                                                                                                                                 gCACGGCGTCGCGAAGTTCCGGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAgt                     >  2:4800626/1‑150 (MQ=255)
                                                                                                                                  cACGGCGTCGCGAAGTTCCGGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAgtg                    <  1:1185840/150‑1 (MQ=255)
                                                                                                                                  cACGGCGTCGCGAAGTTCCGGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAgtg                    <  1:996419/150‑1 (MQ=255)
                                                                                                                                   aCGGCGTCGCGAAGTTCCGGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAgtgt                   <  2:1261886/150‑1 (MQ=255)
                                                                                                                                   aCGGCGTCGCGAAGTTCCGGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAgtgt                   >  1:4471111/1‑150 (MQ=255)
                                                                                                                                    cGGCGTCGCGAAGTTCCTGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTa                  <  1:4289908/150‑1 (MQ=255)
                                                                                                                                    cGGCGTCGCGAAGTTCCGGCAAAATGAAGCGTCGGTTAAACGGGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTa                  >  1:1567665/1‑150 (MQ=255)
                                                                                                                                     ggCGTCGCGAAGTTCCGGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAg                 <  1:3013362/150‑1 (MQ=255)
                                                                                                                                       cGTCGCGAAGTTCCGGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTc               >  2:2149384/1‑150 (MQ=255)
                                                                                                                                          cgcgAAGTTCCGGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGc            >  1:2121003/1‑150 (MQ=255)
                                                                                                                                            cgAAGTTCCGGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGt          <  1:2061965/150‑1 (MQ=255)
                                                                                                                                              aaGTTCCGGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTca        >  2:4532625/1‑150 (MQ=255)
                                                                                                                                              aaGTTCCGGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTca        >  2:1686582/1‑150 (MQ=255)
                                                                                                                                               aGTTCCGGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTcac       >  2:1750205/1‑150 (MQ=255)
                                                                                                                                                gTTCCGGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTGGTCGGCGTcaca      <  2:394289/150‑1 (MQ=255)
                                                                                                                                                gTTCCGGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTcaca      <  2:4086188/150‑1 (MQ=255)
                                                                                                                                                 ttCCGGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACaa     >  1:4690760/1‑150 (MQ=255)
                                                                                                                                                   ccGGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACaaaa   >  1:2689316/1‑150 (MQ=255)
                                                                                                                                                    cGGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACaaaaa  >  2:3271700/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AATCGCTCCACATAGTGTCGGCGAGAAACTCGTAGTGCGAGGTATCGGGCCAGTTTTCCACATGCAGCGGATGAACTTCGTGCGCGCCCATCACTTCCTGCAGGCAAACAATATCGGCGCTGACGGTACGCACGGCGTCGCGAAGTTCCGGCAAAATGAAGCGTCGGTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAA  >  CP009273/819713‑820010

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: