Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,517,528 G→A 8.5% A102V (GCA→GTA)  pdeA ← putative c‑di‑GMP phosphodiesterase PdeA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,517,5280GA 8.5% 454.4 / 22.1 152A102V (GCA→GTA) pdeAputative c‑di‑GMP phosphodiesterase PdeA
Reads supporting (aligned to +/- strand):  ref base G (80/59);  new base A (9/4);  total (89/63)
Fisher's exact test for biased strand distribution p-value = 5.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.62e-01

Only 100 of 156 total aligned reads displayed.

TTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACGGGAACGACACCGCGCCAGCCAAAGATAAACAAGATTGACAG  >  NC_000913/2517379‑2517676
                                                                                                                                                     |                                                                                                                                                    
  aGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTgcg                                                                                                                                                    <  1:2411025/150‑1 (MQ=255)
  aGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTgcg                                                                                                                                                    >  1:1767634/1‑150 (MQ=255)
     gggAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAg                                                                                                                                                 >  1:1009590/1‑150 (MQ=255)
      ggAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAgg                                                                                                                                                >  2:884484/1‑150 (MQ=255)
      ggAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAgg                                                                                                                                                >  2:1833001/1‑150 (MQ=255)
       gAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGc                                                                                                                                               <  1:724566/150‑1 (MQ=255)
        aaaTCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGca                                                                                                                                              <  2:2256750/150‑1 (MQ=255)
        aaaTCAAACAAGCACCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGca                                                                                                                                              >  2:1300707/1‑150 (MQ=255)
            cAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGcacaca                                                                                                                                          >  2:1486341/1‑150 (MQ=255)
            cAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGcacaca                                                                                                                                          >  2:2135614/1‑150 (MQ=255)
            cAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGcacaca                                                                                                                                          >  1:2027731/1‑150 (MQ=255)
            cAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGcacaca                                                                                                                                          >  1:1675533/1‑150 (MQ=255)
               agaagCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACacaaca                                                                                                                                       >  2:214634/1‑150 (MQ=255)
                gaagCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCCGCGGAATATTGCGCAGGcacaccacac                                                                                                                                      >  2:900180/1‑150 (MQ=255)
                  agCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCCAACCAAGCCAGAACAGACGATGCCAGACATATCGGCTGGTCAAGCCATAACGCCAACGCGTCACCAGCTGACAGCGTAATATTGCGCATGCACACAACAAGa                                                                                                                                    >  2:1087063/1‑150 (MQ=255)
                    cTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACa                                                                                                                                  >  2:2275460/1‑150 (MQ=255)
                        ccACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACg                                                                                                                              <  2:1583887/150‑1 (MQ=255)
                          acaAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATTGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTc                                                                                                                            <  1:2070420/150‑1 (MQ=255)
                                tACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTa                                                                                                                      >  1:1833715/1‑150 (MQ=255)
                                tACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTa                                                                                                                      >  1:1975909/1‑150 (MQ=255)
                                     gctgcATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAg                                                                                                                 <  1:376649/150‑1 (MQ=255)
                                      ctgcATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGc                                                                                                                >  1:1061895/1‑150 (MQ=255)
                                       tgcATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCa                                                                                                               <  2:976141/150‑1 (MQ=255)
                                          aTTTGATGCCGATCGGCGTCACCAAACAAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGa                                                                                                            >  2:1998811/1‑150 (MQ=255)
                                              gATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGc                                                                                                        <  1:884484/150‑1 (MQ=255)
                                               aTGCCGATCTGCTTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCTGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCa                                                                                                       <  2:523350/150‑1 (MQ=255)
                                               aTGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCa                                                                                                       >  2:2219953/1‑150 (MQ=255)
                                               aTGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCa                                                                                                       >  1:1982887/1‑150 (MQ=255)
                                                tGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAg                                                                                                      >  2:2203245/1‑150 (MQ=255)
                                                 gCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGt                                                                                                     >  2:603388/1‑150 (MQ=255)
                                                  ccGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTc                                                                                                    >  1:1678624/1‑150 (MQ=255)
                                                    gATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCTGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCt                                                                                                   <  2:2056201/149‑1 (MQ=255)
                                                          cGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTaaaaa                                                                                            >  2:415714/1‑150 (MQ=255)
                                                          cGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTaaaaa                                                                                            <  2:1982887/150‑1 (MQ=255)
                                                            tCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAag                                                                                          >  1:1319082/1‑150 (MQ=255)
                                                            tCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAag                                                                                          <  1:1833001/150‑1 (MQ=255)
                                                            tCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAag                                                                                          <  1:1724810/150‑1 (MQ=255)
                                                             cACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAaga                                                                                         >  2:651475/1‑150 (MQ=255)
                                                             cACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAaga                                                                                         >  1:2106064/1‑150 (MQ=255)
                                                             cACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCACAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAaga                                                                                         >  1:1527321/1‑150 (MQ=255)
                                                              accaaaccaaGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAagag                                                                                        <  1:35086/150‑1 (MQ=255)
                                                               ccaaaccaaGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAagaga                                                                                       >  2:804135/1‑150 (MQ=255)
                                                               ccaaaccaaGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAagaga                                                                                       >  1:2041479/1‑150 (MQ=255)
                                                                    ccaaGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGa                                                                                  >  1:709297/1‑150 (MQ=255)
                                                                         ccAGAACAGACGTTGCCAGACATATCGGCTGGACAAACCATAACGCCCACGAGTCACCAGCTGCCAGCGTAAAATTGCGCAGGAACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGtt                                                                             >  2:899277/1‑150 (MQ=255)
                                                                            gAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCca                                                                          <  1:2203245/150‑1 (MQ=255)
                                                                             aaCAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCcac                                                                         <  2:1767634/150‑1 (MQ=255)
                                                                                 gACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGAt                                                                     >  1:918809/1‑150 (MQ=255)
                                                                                   cGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCCCAGATTg                                                                   >  1:2458174/1‑150 (MQ=255)
                                                                                   cGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTg                                                                   >  2:386371/1‑150 (MQ=255)
                                                                                    gTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTgg                                                                  <  1:2219953/150‑1 (MQ=255)
                                                                                       gCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGc                                                               >  1:1029298/1‑150 (MQ=255)
                                                                                       gCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGc                                                               <  1:435655/150‑1 (MQ=255)
                                                                                        ccAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCa                                                              <  1:2275460/150‑1 (MQ=255)
                                                                                        ccAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCa                                                              <  2:713876/150‑1 (MQ=255)
                                                                                             cATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACa                                                         >  1:1706821/1‑150 (MQ=255)
                                                                                                atCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATcc                                                      >  2:1557915/1‑150 (MQ=255)
                                                                                                atCGGCTGGTCAATCCATAACGCCAACGCGTCAACAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACCTCTGCCTACCGACCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATcc                                                      >  2:460728/1‑150 (MQ=255)
                                                                                                  cGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCg                                                    <  2:1009590/150‑1 (MQ=255)
                                                                                                   ggCTGGTCAATCAATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGc                                                   >  2:338243/1‑150 (MQ=255)
                                                                                                    gCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCa                                                  <  2:1916905/150‑1 (MQ=255)
                                                                                                    gCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCa                                                  <  2:1301403/150‑1 (MQ=255)
                                                                                                       ggTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATaa                                               >  1:1271991/1‑150 (MQ=255)
                                                                                                       ggTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATaa                                               >  1:170384/1‑150 (MQ=255)
                                                                                                       ggTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATaa                                               >  1:849376/1‑150 (MQ=255)
                                                                                                           aaTCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACg                                           >  1:303423/1‑150 (MQ=255)
                                                                                                             tCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACggg                                         <  2:1678624/150‑1 (MQ=255)
                                                                                                               cATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACGGGaa                                       >  1:1267967/1‑150 (MQ=255)
                                                                                                                   aCGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACGGGAACGac                                   <  2:709297/150‑1 (MQ=255)
                                                                                                                    cGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGCGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTCGTGCAAAACATCCCGCATAAAACGGGACCGacc                                  >  2:1196601/1‑149 (MQ=255)
                                                                                                                           gcgTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACGGGAACGACACCGCgcc                           >  1:1918801/1‑150 (MQ=255)
                                                                                                                           gcgTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACGGGAACGACACCGCgcc                           <  1:1486341/150‑1 (MQ=255)
                                                                                                                              tCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACGGGAACGACACCGCgccagc                        <  1:2135614/150‑1 (MQ=255)
                                                                                                                               ccccAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACAGGAACGACACCGCgccagcc                       >  1:2382437/1‑150 (MQ=255)
                                                                                                                                gttgtgtGCCTGCGCAGTATTACGCTGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACGGGAACGACACCGCgccagcca                      >  2:2328153/7‑150 (MQ=255)
                                                                                                                                cccAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACGGGAACGACACCGCgccagcca                      >  1:1036002/1‑150 (MQ=255)
                                                                                                                                cccAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACGGGAACGACACCGCgccagcca                      >  1:1805842/1‑150 (MQ=255)
                                                                                                                                  cAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACGGGAACGACACCGCGCCAGCCaaa                    <  2:944226/150‑1 (MQ=255)
                                                                                                                                  cAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACGGGAACGACACCGCGCCAGCCaaa                    <  1:2284551/150‑1 (MQ=255)
                                                                                                                                   gtgtGCCTGCGCAATATTACGCTGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACGGGAACGACACCGCGCCAGCCAAAg                   >  1:103880/4‑150 (MQ=255)
                                                                                                                                   aGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACGGGAACGACACCGCGCCAGCCAAAg                   <  2:1911335/150‑1 (MQ=255)
                                                                                                                                    tgtGCCTGCGCAATATTACGCTGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACGGGAACGACACCGCGCCAGCCAAAGa                  >  1:422671/3‑150 (MQ=255)
                                                                                                                                    tgtGCCTGCGCAATATTACGCTGGCACACAACACGACAAACGTCTAGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACGGGAACGACACcgcg                             >  1:1978726/3‑139 (MQ=255)
                                                                                                                                    tgtGCCTGCGCAATATTACGCTGGCACACAACACGACAAACGTCTAGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACGGGAACGACACcgcg                             <  2:1978726/137‑1 (MQ=255)
                                                                                                                                    gCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACGGGAACGACACCGCGCCAGCCAAAGa                  >  2:307851/1‑150 (MQ=255)
                                                                                                                                    gCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACGGGAACGACACCGCGCCAGCCAAAGa                  <  2:933022/150‑1 (MQ=255)
                                                                                                                                     cTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACGGGAACGACACCGCGCCAGCCAAAGAt                 <  1:1998811/150‑1 (MQ=255)
                                                                                                                                      tGCCTGCGCAATATTACGCTGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACGGGAACGACACCGCGCCAGCCAAAGATa                >  1:1607142/1‑150 (MQ=255)
                                                                                                                                       gCCTGCGCAATATTACGCTGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACGGGAACGACACCGCGCCAGCCAAAGATaa               >  2:637402/1‑150 (MQ=255)
                                                                                                                                       gCCTGCGCAATATTACGCTGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACGGGAACGACACCGCGCCAGCCAAAGATaa               >  2:768909/1‑150 (MQ=255)
                                                                                                                                         cAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACGGGAACGACACCGCGCCAGCCAAAGATAAAc             <  2:170384/150‑1 (MQ=255)
                                                                                                                                         cAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTATAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACGGGAACGACACCGCGCCAGCCAAAGATAAAc             <  1:34938/150‑1 (MQ=255)
                                                                                                                                          aGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACGGGAACGACACCGCGCCAGCCAAAGATAAACa            >  2:1446957/1‑150 (MQ=255)
                                                                                                                                            cGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGCCCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACGGGAACGACACCGCGCCAGCCAAAGATAAACAAg          >  1:1607330/1‑150 (MQ=255)
                                                                                                                                            cGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACGGGAACGACACCGCGCCAGCCAAAGATAAACAAg          <  2:1198590/150‑1 (MQ=255)
                                                                                                                                             gTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACGGGAACGACACCGCGCCAGCCAAAGATAAACAAGa         <  2:1721653/150‑1 (MQ=255)
                                                                                                                                                atatTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACGGGAACGACACCGCGCCAGCCAAAGATAAACAAGATTg      >  1:1153273/1‑150 (MQ=255)
                                                                                                                                                atatTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACGGGAACGACACCGCGCCAGCCAAAGATAAACAAGATTg      <  2:2106064/150‑1 (MQ=255)
                                                                                                                                                  atTGCGCAGGCACACAACACGGCAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACGGGAACGACACCGCGCCAGCCAAAGATAAACAAGATTGAc    >  2:2411111/1‑150 (MQ=255)
                                                                                                                                                    tGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACGGGAACGACACCGCGCCAGCCAAAGATAAACAAGATTGACAg  >  2:2136934/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACGGGAACGACACCGCGCCAGCCAAAGATAAACAAGATTGACAG  >  NC_000913/2517379‑2517676

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