Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,190,004 T→C 7.1% G191G (GGT→GGC)  yjaG → DUF416 domain protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,190,0040TC 7.1% 164.1 / 3.1 56G191G (GGT→GGC) yjaGDUF416 domain protein
Reads supporting (aligned to +/- strand):  ref base T (29/23);  new base C (2/2);  total (31/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.17e-01

TCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTAGAGCAAGAATGGGACATTCAGTGGGAAATATTCCGACTTTTAGCTGAGTGCGAAGAACGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGATAACTTATGAACAAGACTCAACTGATTGATGTAATTGCAGAGAAAGCAGAACTGTCCAAAACCCAGGCTAAAG  >  CP009273/4189742‑4190275
                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                               
tCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTAGAGCAAGAATGGGACATTCAGTGGGAAATATTCCGACTTTTAGCTGAGTGCGAAGAACGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTAAg                                                                                                                                                                                                                                                                     >  1:453662/1‑275 (MQ=255)
                                      gTGGAAGTGAGTAAGACCTCAATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCNTAAAGAAANCCCAGCTGTAGATCAAGAATGGGACATTCAGTGGGAAATATTCCGACTTTTAGCTGAGTGCGAAGAACGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAACAAGACCAAAAAACGTGATTTAACGCCTGATTTGt                                                                                                                                                                                                                               <  1:560251/275‑1 (MQ=255)
                                                 taagcacTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACNCAGCTGTANAGCAAGAATGGGTCATTCAGTGGGAAATATTCCGACTTTTAGCTGAGTGCGAAGAACGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAATACGTGATATAACGCCTGATTTGTCGTACCTGGAg                                                                                                                                                                                                                    <  1:540394/269‑1 (MQ=255)
                                                                  ccGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTAGAGCAAGAATGGGACATTCAGTGGGAAATATTCCGACTTTTAGCTGAGTGCGAAGAACGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGNGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAg                                                                                                                                                                                                                                                         >  1:653462/1‑221 (MQ=255)
                                                                                     gACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTAGAGCAAGAATGGGACATTCAGTGGGAAATATTCCGACTTTTAGCTGAGTGCGAAGAACGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCNTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTggg                                                                                                                                                                                >  1:714074/1‑275 (MQ=255)
                                                                                                gTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTAGAGCAAGAATGGGACATTCAGTGGGAAATATTCCGACTTTTAGCTGAGTGCGAAGAACGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGCATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGTAg                                                                                                                                                                                                                    >  2:486020/1‑228 (MQ=255)
                                                                                                    cgAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTAGAGCAAGAATGGGACATTCAGTGGGAAATATTCCGACTTTTAGCTGAGTGCGAAGAACGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGCATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCtt                                                                                                                                                                                                                <  1:486023/228‑1 (MQ=255)
                                                                                                      aaaTGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTAGAGCAAGAATGGGACATTCAGTGGGAAATATTCCGACTTTTAGCTGAGTGCGAAGAACGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAAttt                                                                                                                                                                                                                                                                          >  1:54416/1‑168 (MQ=255)
                                                                                                          gAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTAGAGCAAGAATGGGACATTCAGTGGGAAATATTCCGACTTTTAGCTGAGTGCGAAGAACGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCa                                                                                                                                                                                                                                                                      <  2:54416/168‑1 (MQ=255)
                                                                                                                       tgtgtataagtgacAGCTGTAGAGCAAGCATGGGACATTAAGTGGGAAATATTCCGACTTTTAGATGAGTGCGAAGAACGCNATATCGAGNTGATAAAAGGCCCTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGTTATAATTTTTCAGCAATAAGACCAGAAACCGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGAGGGGCGACAAAAAGTGGCTATCGGTGCGTGTATg                                                                                                                                               <  1:552133/261‑1 (MQ=255)
                                                                                                                       cTTAAAGAAAACCCAGCTGTAGAGCAAGAATGGGACATTCAGTGGGAAATATTCCGACTTTTAGCTGAGTGCGAAGAACGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGt                                                                                                                                                        >  1:416408/1‑265 (MQ=255)
                                                                                                                                     aGCTGTAGAGCAAGAATGGGACATTCAGTGGGAAATATTCCGACTTTTAGCTGAGTGCGCAGACCGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACATGCGTGAGTAGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACAAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCGCTCTGGTCTACATTTGGGGGGCGAAAAAAAGGGGCTATGCTTGGGTGTATGCacaaaatttttt                                                                                                                                  >  2:686584/1‑262 (MQ=255)
                                                                                                                                                 aGAATGGGACATTCAGTGGGAAATATTCCGACTTTTAGCTGAGTGCGAAGAACGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATGGGTGCGTGTATGCAGGAGAGTGCAATTCTGGCATTTCCg                                                                                                                    >  1:137476/1‑275 (MQ=255)
                                                                                                                                                      gggACATTCAGTGGGAAATATTCCGACTTTTAGCTGAGTGAGAAGAACGCGATATCGAGCTGATAAAAGGCCTTAGGGCATACCGTCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCATAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGGTATTCTGGCATTTCCGTCg                                                                                                                 <  1:686590/273‑1 (MQ=255)
                                                                                                                                                        gaatttgagggggAAATATTCCAAATTTTGGGGTAGTGCGAAGAACGCGATATCGAGGTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGNGGGTGAGATCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGTTGCGTGTATGCAGGATAGTGCTATTCTGGCATTTCCGTCGCACt                                                                                                             <  1:623802/265‑1 (MQ=255)
                                                                                                                                                             tCAGTGGGAAATATTCCGACTTTTAGCTGAGTGCGAAGAACGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGGGAGAGCAATATTGGTATAAtt                                                                                                                                                                                                                                                                           >  1:727950/1‑112 (MQ=255)
                                                                                                                                                             tCAGTGGGAAATATTCCGACTTTTAGCTGAGTGCGAAGAACGCGATATCGAGCTGATAAAAGGCCGTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATGTAACGCCTGATTTGTCGTACCTGGAGTCTTCACTTTAGCACCCCGTCTGATCTACATTGGGTGAGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGAGCGATGCTGGCATTTCCggagaactgtatg                                                                                                        >  2:251574/1‑263 (MQ=255)
                                                                                                                                                                 tGGGAAATATTCCGACTTTTAGCTGAGTGCGAAGAACGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAAttt                                                                                                                                                                                                                                                                          >  1:322854/1‑109 (MQ=255)
                                                                                                                                                                    gAAATATTCCGACTTTTAGCTGAGTGCGAAGAACGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTGCCCTTTCGcccc                                                                                                                                                                                                    >  1:648713/1‑176 (MQ=255)
                                                                                                                                                                     aaaTATTCCGACTTTTAGCTGAGTGCGAAGAACGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCa                                                                                                                                                                                                                                                                      <  2:322851/109‑1 (MQ=255)
                                                                                                                                                                        tatTCCGACTTTTAGCTGAGTGCGAAGAACGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGt                                                                                                                                                                                                <  2:648707/176‑1 (MQ=255)
                                                                                                                                                                             cGACTTTTAGCTGAGTGCGAAGAACGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGAtt                                                                                                                                                                                                                                  >  2:353103/1‑137 (MQ=255)
                                                                                                                                                                                 ttttAGCTGAGTGCGAAGAACGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTc                                                                                                                                                                                                                              <  1:353106/137‑1 (MQ=255)
                                                                                                                                                                                          aGTGCGAAGAACGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGcccc                                                                                                                                                                                                    >  2:774403/1‑154 (MQ=255)
                                                                                                                                                                                            tGCGAAGAACGCGATTTCGAGCTGATAAAGGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAANGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGAt                                                                          >  1:576966/1‑274 (MQ=255)
                                                                                                                                                                                              cGAAGAACGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGt                                                                                                                                                                                                <  1:774409/154‑1 (MQ=255)
                                                                                                                                                                                                    aCGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGc                                                                                                                                              >  2:147867/1‑198 (MQ=255)
                                                                                                                                                                                                      gcgATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTa                                                                                                                                                                                                                                            >  1:70421/1‑102 (MQ=255)
                                                                                                                                                                                                      gcgATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGcc                                                                                                                                                                                                                                       >  2:515409/1‑107 (MQ=255)
                                                                                                                                                                                                        gATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGga                                                                                                                                          <  1:147868/198‑1 (MQ=255)
                                                                                                                                                                                                          tatCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGc                                                                                                                                                                                                                                        <  2:70421/102‑1 (MQ=255)
                                                                                                                                                                                                          tatCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGAt                                                                                                                                                                                                                                   <  1:515412/107‑1 (MQ=255)
                                                                                                                                                                                                                gCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGCATAAtt                                                                                                                                                                                                                                                                           >  2:810617/1‑61 (MQ=255)
                                                                                                                                                                                                                    aTAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGCATAATTTTTc                                                                                                                                                                                                                                                                       <  1:810624/61‑1 (MQ=255)
                                                                                                                                                                                                                              ttAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCACTTTCGcctc                                                                                                                                                                                                    >  2:677565/1‑116 (MQ=255)
                                                                                                                                                                                                                                  ggCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGt                                                                                                                                                                                                <  1:677571/118‑1 (MQ=255)
                                                                                                                                                                                                                                            cGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTaa                                                                              >  1:343882/1‑222 (MQ=255)
                                                                                                                                                                                                                                              tGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGaaaaa                                                                                                                                                                      >  2:155199/1‑132 (MQ=255)
                                                                                                                                                                                                                                                aGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGAt                                                                          <  2:343879/222‑1 (MQ=255)
                                                                                                                                                                                                                                                 ggCGGGTGAGAGCAATCTTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGcccc                                                                                                                                                                                                    >  1:205540/1‑99 (MQ=255)
                                                                                                                                                                                                                                                 ggCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTGACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCAGTCg                                                                                                                 >  1:437537/1‑182 (MQ=255)
                                                                                                                                                                                                                                                 ggCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGACATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGATAACTTATTAAc                                                               >  1:515777/1‑232 (MQ=255)
                                                                                                                                                                                                                                                  gCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGt                                                                                                                                                                  <  1:155200/132‑1 (MQ=255)
                                                                                                                                                                                                                                                     ggTGAGAGCAATATTTGTATATTTTTTCATCAATAAAACCAGAAACCGTCATTTAACGCCTTATTTGTCGTACCTGGAGTCTTCCCTTGCGCCCCCCGTCTGGTTTACATGTGCGGAGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACt                                                                                                             <  2:437534/182‑1 (MQ=255)
                                                                                                                                                                                                                                                     ggTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGAc                                                                                                                                                                                                                                                            >  2:385891/1‑39 (MQ=255)
                                                                                                                                                                                                                                                     ggTGAGAGCAAGAGTGGTATAATTATTCAGGAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTGCACTTTCGCCCCCCGt                                                                                                                                                                                                <  2:205537/99‑1 (MQ=255)
                                                                                                                                                                                                                                                         agagCAATATTGGTATAATTTTTCAGCAATAAGACCAGa                                                                                                                                                                                                                                                        <  1:385894/39‑1 (MQ=255)
                                                                                                                                                                                                                                                          gggaggtgttGGTATCATTTTTCAGCAATAAGACCCTACAACGTGATTTATCTCCTGCTTTGTCGTACCTGGAGGCTTCCCTTTCACCCCCCGTGTGGTCTACATTGGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACAGTGTAAGGATAACTTATGAACAAGACTCAACTGATTGATGTAATTGCAGAGAAAGCAGAACTGTCCAAAAcc            <  2:137475/266‑1 (MQ=255)
                                                                                                                                                                                                                                                                tatTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGAt                                                                                                                                                                                                                                   >  1:765992/1‑53 (MQ=255)
                                                                                                                                                                                                                                                                tatTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGc                                                                                                                                              >  1:540456/1‑138 (MQ=255)
                                                                                                                                                                                                                                                                tatTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCCATTCTGGCATTTc                                                                                                                      >  1:693585/1‑162 (MQ=255)
                                                                                                                                                                                                                                                                    tttataATTTTTCAGCAATAAGACCTGTATATGTTATTTAACGCGTGATTTGTCTTACTTGGCGTCTTCCCTTGCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCGTAGCAAGAGATAAACACATTGTAAGGATAACTTATGAACAAGACTCAACTGATTGATGTAATTGCAGAGAAAGCAGAACTGTCCAAAACCCAGGCTAAAg  <  1:448760/272‑1 (MQ=255)
                                                                                                                                                                                                                                                                    ggTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGt                                                                                                                                                                                                                               <  2:765986/53‑1 (MQ=255)
                                                                                                                                                                                                                                                                    ggTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGga                                                                                                                                          <  2:540453/138‑1 (MQ=255)
                                                                                                                                                                                                                                                                    ggTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCCATTCTGGCATTTCCGTc                                                                                                                  <  2:693579/162‑1 (MQ=255)
                                                                                                                                                                                                                                                                     gTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGcccc                                                                                                                                                                                                    >  2:520434/1‑79 (MQ=255)
                                                                                                                                                                                                                                                                     gTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGATAACTTATGAACAAGACTCAACtga                                                  >  2:803628/1‑225 (MQ=255)
                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                               
TCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTAGAGCAAGAATGGGACATTCAGTGGGAAATATTCCGACTTTTAGCTGAGTGCGAAGAACGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGATAACTTATGAACAAGACTCAACTGATTGATGTAATTGCAGAGAAAGCAGAACTGTCCAAAACCCAGGCTAAAG  >  CP009273/4189742‑4190275

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑