Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,605,542 T→A 8.4% G447G (GGA→GGT yjjI ← DUF3029 family protein, putative glycine radical enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,605,5420TA 8.4% 575.6 / 35.6 191G447G (GGA→GGTyjjIDUF3029 family protein, putative glycine radical enzyme
Reads supporting (aligned to +/- strand):  ref base T (84/91);  new base A (9/7);  total (93/98)
Fisher's exact test for biased strand distribution p-value = 6.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

Only 100 of 197 total aligned reads displayed.

GTGTTGCGTGCGGCTTCTTCGCCCAGCCAGGTGGTGTTGGTGCGTGAACCTTCGGCGCGATATTTTTCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGC  >  CP009273/4605396‑4605689
                                                                                                                                                  |                                                                                                                                                   
   ttGCGTGCGGCTTCTTCGCCCAGCCAGGTGGTGTTGGTGCGTGAACCTTCGGCGCGATATTTTTCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCgg                                                                                                                                                 <  2:1601393/148‑1 (MQ=255)
   ttGCGTGCGGCTTCTTCGCCCAGCCAGGTGGTGTTGGTGCGTGAACCTTCGGCGCGATATTTTTCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCgg                                                                                                                                                 <  2:632663/148‑1 (MQ=255)
      cgtgcgGCTTCTTCGCCCAGCCAGGTGGTGTTGGTGCGTGAACCTTCGGCGCGATATTTTTCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCtt                                                                                                                                              <  1:313172/148‑1 (MQ=255)
          cgGCTTCTTCGCCCAGCCAGGTGGTGTTGGTGCGTGAACCTTCGGCGCGATATTTTTCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTaaa                                                                                                                                          <  2:83776/148‑1 (MQ=255)
          cgGCTTCTTCGCCCAGCCAGGTGGTGTTGGTGCGTGAACCTTCGGCGCGATATTTTTCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTaaa                                                                                                                                          <  1:2931284/148‑1 (MQ=255)
          cgGCTTCTTCGCCCAGCCAGGTGGTGTTGGTGCGTGAACCTTCGGCGCGATATTTTTCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTaaa                                                                                                                                          <  2:2788987/148‑1 (MQ=255)
          cgGCTTCTTCGCCCAGCCAGGTGGTGTTGGTGCGTGAACCTTCGGCGCGATATTTTTCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTaaa                                                                                                                                          >  1:3298270/1‑148 (MQ=255)
               tcttcGCCCAGCCAGGTGGTGTTGGTGCGTGAACCTTCGGCGCGATATTTTTCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAATGTCAc                                                                                                                                     >  2:3119372/1‑148 (MQ=255)
               tcttcGCCCAGCCAGGTGGTGTTGGTGCGTGAACCTTCGGCGCGATATTTTTCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCAc                                                                                                                                     <  2:1232318/148‑1 (MQ=255)
                       cagccagGTGGTGTTGGTGCGTGAACCTTCGGCGCGATATTTTTCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCaa                                                                                                                             <  2:1729587/148‑1 (MQ=255)
                       cagccagGTGGTGTTGGTGCGTGAACCTTCGGCGCGATATTTTTCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCaa                                                                                                                             >  1:2398862/1‑148 (MQ=255)
                       cagccagGTGGTGTTGGTGCGTGAACCTTCGGCGCGATATTTTTCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATACCGGCTTTAAAGGCACCGAGGCaa                                                                                                                             <  1:2062141/148‑1 (MQ=255)
                         gccagGTGGTGTTGGTGCGTGAACCTTCGGCGCGATATTTTTCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAg                                                                                                                           >  1:885920/1‑148 (MQ=255)
                           cagGTGGTGTTGGTGCGTGAACCTTCGGCGCGATATTTTTCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCt                                                                                                                         <  2:1363482/148‑1 (MQ=255)
                            agGTGGTGTTGGTGCGTGAACCTTCGGCGCGATATTTTTCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTg                                                                                                                        <  1:2468288/148‑1 (MQ=255)
                            agGTGGTGTTGGTGCGTGAACCTTCGGCGCGATATTTTTCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATACCGGCTTTATAGGCACCGAGGCAAAGCTg                                                                                                                        <  2:1008791/148‑1 (MQ=255)
                              gtggtgTTGGTGCGTGAACCTTCGGCGCGATATTTTTCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTa                                                                                                                      >  1:2832690/1‑148 (MQ=255)
                              gtggtgTTGGTGCGTGAACCTTCGGCGCGATATTTTTCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCGTTAAAGGCACCGAGGCAAAGCTGTa                                                                                                                      >  1:1889640/1‑148 (MQ=255)
                                  tgtTGGTGCGTGAACCTTCGGCGCGATATTTTTCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAg                                                                                                                  >  1:489208/1‑148 (MQ=255)
                                          cGTGAACCTTCGGCGCGATATTTTTCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCg                                                                                                          >  2:2099155/1‑148 (MQ=255)
                                            tGAACCTTCGGCGCGATATTTTTCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCggg                                                                                                        >  2:1821981/1‑148 (MQ=255)
                                             gAACCTTCGGCGCGATATTTTTCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGt                                                                                                       >  1:105416/1‑148 (MQ=255)
                                               aCCTTCGGCGCGATATTTTTCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTa                                                                                                     <  2:1869977/148‑1 (MQ=255)
                                                ccTTCGGCGCGATATTTTTCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTAc                                                                                                    <  1:2768232/148‑1 (MQ=255)
                                                ccTTCGGCGCGATATTTTTCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCGGCGGGTTAc                                                                                                    >  1:1810576/1‑148 (MQ=255)
                                                  ttCGGCGCGATATTTTTCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGt                                                                                                  <  2:409237/148‑1 (MQ=255)
                                                      gcgcGATATTTTTCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGa                                                                                              <  1:580432/148‑1 (MQ=255)
                                                        gcgATATTTTTCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATg                                                                                            <  1:1998151/148‑1 (MQ=255)
                                                           atatTTTTCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGtt                                                                                         >  2:3237071/1‑148 (MQ=255)
                                                            tatTTTTCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGttt                                                                                        >  1:1542544/1‑148 (MQ=255)
                                                              tttttCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTtcg                                                                                      <  2:1755109/148‑1 (MQ=255)
                                                              ttatacTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTtcg                                                                                      >  2:3186306/6‑148 (MQ=255)
                                                               ttttCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTtcgt                                                                                     >  1:2005981/1‑148 (MQ=255)
                                                               ttttCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTtcgt                                                                                     >  1:1816619/1‑148 (MQ=255)
                                                                 ttCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTtcgtcg                                                                                   >  2:2123553/1‑148 (MQ=255)
                                                                     aaaTCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCg                                                                               >  1:819455/1‑148 (MQ=255)
                                                                       aTCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTc                                                                             >  2:2329305/1‑148 (MQ=255)
                                                                        tCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCa                                                                            >  2:2438756/1‑148 (MQ=255)
                                                                          cGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGa                                                                          <  1:781046/148‑1 (MQ=255)
                                                                          cGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGa                                                                          <  2:193967/148‑1 (MQ=255)
                                                                          cGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGGTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGa                                                                          <  1:3082299/148‑1 (MQ=255)
                                                                             cAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATg                                                                       >  2:1398013/1‑148 (MQ=255)
                                                                             cAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATg                                                                       <  1:1642394/148‑1 (MQ=255)
                                                                              aaaCGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGt                                                                      <  1:856419/148‑1 (MQ=255)
                                                                                aCGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCg                                                                    <  1:341552/148‑1 (MQ=255)
                                                                                   cACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTg                                                                 <  2:2036522/148‑1 (MQ=255)
                                                                                      cATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATg                                                              >  1:1019457/1‑148 (MQ=255)
                                                                                       atatAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGc                                                             >  1:3086873/1‑148 (MQ=255)
                                                                                        tataACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGcc                                                            >  1:602133/1‑148 (MQ=255)
                                                                                         ataACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCg                                                           >  1:541705/1‑148 (MQ=255)
                                                                                         ataACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCg                                                           >  2:502866/1‑148 (MQ=255)
                                                                                          taACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCgg                                                          >  2:180155/1‑148 (MQ=255)
                                                                                            aCCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGTAAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTAACAGTGCCGGCGGGTTGCGTTTGATGTTATCGTCGAGCGTCAGAATGACGCAGATGCCGGaa                                                        >  2:1643930/1‑148 (MQ=255)
                                                                                            aCCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGaa                                                        >  1:1740910/1‑148 (MQ=255)
                                                                                             ccGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGaat                                                       <  2:514632/148‑1 (MQ=255)
                                                                                             ccGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGaat                                                       >  1:221700/1‑148 (MQ=255)
                                                                                              cGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGaata                                                      <  1:1250559/148‑1 (MQ=255)
                                                                                               ggTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGACACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGaataa                                                     <  1:1128988/148‑1 (MQ=255)
                                                                                                 tAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGaataata                                                   <  2:607043/148‑1 (MQ=255)
                                                                                                 tAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGaataata                                                   <  2:2383299/148‑1 (MQ=255)
                                                                                                    cgcgAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAataataata                                                <  1:3237062/148‑1 (MQ=255)
                                                                                                    cgcgAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAataataata                                                <  1:1030815/148‑1 (MQ=255)
                                                                                                       gAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGc                                             <  1:289147/148‑1 (MQ=255)
                                                                                                         aCCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCat                                           >  2:2572627/1‑148 (MQ=255)
                                                                                                         aCCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCat                                           >  2:783503/1‑148 (MQ=255)
                                                                                                         aCCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCat                                           <  1:46600/148‑1 (MQ=255)
                                                                                                          ccAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCatg                                          <  1:1861748/148‑1 (MQ=255)
                                                                                                           cAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCACAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGCAATAATAATAAGCatga                                         <  2:1727449/148‑1 (MQ=255)
                                                                                                           cAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCatga                                         <  1:195390/148‑1 (MQ=255)
                                                                                                             gATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCatgatg                                       <  2:2702715/148‑1 (MQ=255)
                                                                                                               tCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATTCg                                     >  2:1720892/1‑148 (MQ=255)
                                                                                                               tCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCg                                     <  1:1473336/148‑1 (MQ=255)
                                                                                                               tCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCg                                     <  1:438282/148‑1 (MQ=255)
                                                                                                               tCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCg                                     >  1:368058/1‑148 (MQ=255)
                                                                                                               tCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCg                                     >  2:1720859/1‑148 (MQ=255)
                                                                                                               tCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCg                                     <  1:2362604/148‑1 (MQ=255)
                                                                                                                         tGACATTTTCGGTAAATTCACGCATTACGGCTTTAAAGACACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTGACGTTTGATGTTTTCTTCGAGCGTCAGAATGTCGCTGATGACGGAATAATAATAATAATGATGAGGGGCGACAGt                           >  2:148657/1‑148 (MQ=255)
                                                                                                                         tGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGt                           <  2:1403275/148‑1 (MQ=255)
                                                                                                                          gACATTGGCGGTAAATTCACGCATACCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGAc                              >  1:3020774/1‑144 (MQ=255)
                                                                                                                          gACATTGGCGGTAAATTCACGCATACCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGtt                          >  1:604040/1‑148 (MQ=255)
                                                                                                                          gACATTGGCGGTAAATTCACGCATACCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGtt                          >  1:2855923/1‑148 (MQ=255)
                                                                                                                           aCATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGttt                         <  2:1334705/148‑1 (MQ=255)
                                                                                                                           aCATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGttt                         <  2:2982275/148‑1 (MQ=255)
                                                                                                                            cATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTg                        <  2:1566695/148‑1 (MQ=255)
                                                                                                                                ggCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGa                    >  1:1091831/1‑148 (MQ=255)
                                                                                                                                   ggTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATgg                 <  2:3086873/148‑1 (MQ=255)
                                                                                                                                   ggTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATgg                 <  2:489208/148‑1 (MQ=255)
                                                                                                                                     tAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGt               >  1:3219184/1‑148 (MQ=255)
                                                                                                                                       aaTTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGa             <  1:502866/148‑1 (MQ=255)
                                                                                                                                       aaTTCACGCATACCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGa             <  2:1456728/148‑1 (MQ=255)
                                                                                                                                        aTTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGAt            <  1:293959/148‑1 (MQ=255)
                                                                                                                                          tCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCg          >  1:751612/1‑148 (MQ=255)
                                                                                                                                          tCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCg          >  1:2101230/1‑148 (MQ=255)
                                                                                                                                           cACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCgg         <  1:2170902/148‑1 (MQ=255)
                                                                                                                                             cGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGAt       <  2:2871913/148‑1 (MQ=255)
                                                                                                                                               cATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCt     >  1:1783438/1‑148 (MQ=255)
                                                                                                                                                  tCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGc  >  2:3299272/1‑148 (MQ=255)
                                                                                                                                                  tCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGc  >  2:1488471/1‑148 (MQ=255)
                                                                                                                                                  tCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGc  >  2:727834/1‑148 (MQ=255)
                                                                                                                                                  tCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGAGGGTTACGTATGATGGTTTCGTCGAGCGTCAGAATGACGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTAGCAGATGGGTGATCGGATCTGGc  >  1:243887/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GTGTTGCGTGCGGCTTCTTCGCCCAGCCAGGTGGTGTTGGTGCGTGAACCTTCGGCGCGATATTTTTCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGC  >  CP009273/4605396‑4605689

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: