Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,809,169 Δ1 bp 28.6% intergenic (‑41/+25) pyrE ← / ← rph orotate phosphoribosyltransferase/ribonuclease PH (defective);enzyme; Degradation of RNA; RNase PH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,809,1690G.28.6% 133.1 / 59.5 56intergenic (‑41/+25)pyrE/rphorotate phosphoribosyltransferase/ribonuclease PH (defective);enzyme; Degradation of RNA; RNase PH
Reads supporting (aligned to +/- strand):  ref base G (18/22);  new base . (10/6);  total (28/28)
Fisher's exact test for biased strand distribution p-value = 3.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACTCAATGCCGGAATCCACCAACGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGATCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGC  >  CP009273/3808938‑3809429
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                     
aaaagaaagccggaagcCACCAACGCGGCAGCGGAAAACCGGACTAACAGTTCCAGAGCGAGCCAGGGATTAAAAAGCCAGGCGTTGAAGAAATAGGGGCGGTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCC                                                                                                                                                                                                                            <  2:617647‑M1/267‑1 (MQ=255)
      aTGCCGGAATCCACCAACGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGATCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGt                                                                                                                                                                                                                                                          >  1:137643/1‑239 (MQ=255)
      aTGCCGGAATCCACCAACGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGATCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGt                                                                                                                                                                                                                                                          <  2:137642/239‑1 (MQ=255)
      aTGCCGGAATCCACCAACGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGATCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCATTTTTCCACCCCCCCCTTCTtct                                                                                                                                                                                                                     >  2:561382‑M1/1‑272 (MQ=255)
               tCCACCAACGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGATCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTCAACTCGACAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTATACAGAC‑AAAAAACAGGCGACTCAACAGACGCCTTAAAAATCAGTTTGCAAGAACCGCCTTCTGCGTccccac                                                                                                                                                                                                              >  2:249454/1‑269 (MQ=255)
               gaaaatacagattcagagtaaaaaNGACCTCACAGTGCCCGATCGCTCACGGTATTAAACAGCCCGGCGTTTAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCGAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCTTTCTCTCCTCATCTTACTTTTCGACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctg                                                                                                                                                                                                             <  2:422990‑M1/264‑10 (MQ=255)
                             aGCGTAAAATCGGCCTAACAGTGCCAGATCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTG                                                                                                                                                                                                                       <  1:164664‑M1/251‑1 (MQ=255)
                                      acttcccaaaaatgccagagcgcgcCCTGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACCTTAACACCTGCTTGCTAAGCGCAAATTCAATAACCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtctt                                                                                                                                                                                      <  1:346089‑M1/263‑33 (MQ=255)
                                                                         acggCCCGGCGTTGCAGAAATAGGCGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTGTAAACCCTGCTTGCTAAGCTCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTGAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGaa                                                                                                                                                  <  1:249457/272‑1 (MQ=255)
                                                                         aCAGCCCTGCGTTGAAGAAATAGGTGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggca                                                                                                                                                   <  1:318400‑M1/275‑68 (MQ=255)
                                                                             cccGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgat                                                                                                                                                                                                 >  1:48334‑M1/1‑204 (MQ=255)
                                                                             cccGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTACTGCGTCGCTACAATGGATTCAATTCCCCTC‑GGTCCAGAGCCACAAAGAAGAGTAACTCTTCATGTGTGAACGGc                                                                                                                                              >  1:317263/1‑274 (MQ=255)
                                                                             cccGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCAAGCGCCGCCTTCTGCGTCTCTACAATGGATCCGATTCCCCTC‑GGGCCAGAGCCAACAAGACGAGTAGCTCTTCATGGGAGAACGGc                                                                                                                                              >  1:226372/1‑274 (MQ=255)
                                                                                 gCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgc                                                                                                                                                                                             <  2:48333‑M1/225‑26 (MQ=255)
                                                                                        agaaATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTAAATAAACTGGCGCTGNTATGGTTTCATGCCTTCGCGCCTCATCTTACTTTTCTACAGAC‑AAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGcc                                                                                                                                  <  1:561389/275‑1 (MQ=255)
                                                                                                                                gCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCT                                                                                                                                                                                                                           >  2:475461‑M1/1‑148 (MQ=255)
                                                                                                                                    aaCTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAATAATCAGTTTGCCAGCGCCGCCTTCTG                                                                                                                                                                                                                       <  1:475467‑M1/148‑1 (MQ=255)
                                                                                                                                               cTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgat                                                                                                                                                                                                 >  1:408123‑M1/1‑138 (MQ=255)
                                                                                                                                                 gcttgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGt                                                                                                                                >  2:137611/1‑220 (MQ=255)
                                                                                                                                                   ttgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgc                                                                                                                                                                                             <  2:408117‑M1/159‑26 (MQ=255)
                                                                                                                                                     gctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTcccc                                                                                                                            <  1:137612/220‑1 (MQ=255)
                                                                                                                                                           cgcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGcc                                                                                                                                                                                                                                     >  1:77869/1‑110 (MQ=255)
                                                                                                                                                           cgcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGTTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGCCCCTGCGGcg                                                                >  1:594462/1‑273 (MQ=255)
                                                                                                                                                               aaTTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAgcg                                                                                                                                                                                                                                 <  2:77868/110‑1 (MQ=255)
                                                                                                                                                               aaTTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCAGCTTAATTTTCTACAGAC‑AAAAAAAAGGc                                                                                                                                                                                                                                                                     <  1:504390/75‑1 (MQ=255)
                                                                                                                                                                             gcgcTGATATGGTTTCTTGCCTTCGATCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaaactgcggctaccatccctttcatcggattggttttcagcttgc                                               <  1:440399‑M1/275‑168 (MQ=255)
                                                                                                                                                                             gcgcTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTC‑GGGCCAGAGCCAACAAGATGa                                                                                                                                                                     >  2:647643/1‑155 (MQ=255)
                                                                                                                                                                                 tGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTC‑GGGCCAGAGCCAACAAGATGAGTAg                                                                                                                                                                 <  1:647651/155‑1 (MQ=255)
                                                                                                                                                                                   atatGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatg                                                                                                                                                                                                >  2:589536‑M1/1‑102 (MQ=255)
                                                                                                                                                                                      tGGTTTCATGCCTTCTCTCCTCATCTTACTTTTCTACAGAA‑AAAAAAAAGTCGACTCATCAGTCGCCTTAAAAATCAGTGTGCCAGCGCCGNCTTCTGCcgtcccctgcacatcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaaactgcggctaccatccctttcatcggattggttttcagcttgccgttttcca                                      <  1:575050‑M1/275‑177 (MQ=255)
                                                                                                                                                                                       ggTTTCATGCCTTCGCTACTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgccc                                                                                                                                                                                           >  1:207600‑M1/1‑98 (MQ=255)
                                                                                                                                                                                             aTGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgc                                                                                                                                                       >  1:327753‑M1/1‑92 (MQ=255)
                                                                                                                                                                                             aTGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGACGCTACAATGGATTCGATTCCCCTCGGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGAt                                          >  1:108224/1‑264 (MQ=255)
                                                                                                                                                                                                 cTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggca                                                                                                                                                   <  2:327748‑M1/155‑68 (MQ=255)
                                                                                                                                                                                                 cTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGACGCTACAATGGATTCGATTCCCCTCGGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCa                                      <  2:108223/264‑1 (MQ=255)
                                                                                                                                                                                                  ttCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCGTCTGCCGTCCCCTCCACTTCACTGATGCGCCCGGCTTCGGTCATCACTACGTTCATGTCGGTCTCTGGGGCACCGCCTTCAACGTCTGCCCGAGCGCACACCGCTTTGcgtc                        >  2:19061/1‑272 (MQ=255)
                                                                                                                                                                                                   tCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTC‑GGGCCAGAGCCAACAAGCTGAGTAGc                                                                                                                                                                >  2:515419/1‑138 (MQ=255)
                                                                                                                                                                                                      ctcctcATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCANTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaaactgcggctaccatccctttca                                                                    <  1:555742‑M1/229‑147 (MQ=255)
                                                                                                                                                                                                       tcctcATCTTAATTTTCTACAGAC‑AAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTC‑GTGCCTTACCCAACAAGATGAGTAGCTCtt                                                                                                                                                            <  1:515425/138‑1 (MQ=255)
                                                                                                                                                                                                              cTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaaactgcggctaccatccctttcatcggattggttttcagcttgccgttttccaccatcttctgtagcgcatctacc               >  1:103570‑M1/1‑75 (MQ=255)
                                                                                                                                                                                                                 aCTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagag                                                                                                                                                >  1:3445‑M1/1‑72 (MQ=255)
                                                                                                                                                                                                                 aCTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaa                                                                                                                                                                                                     >  2:216991‑M1/1‑72 (MQ=255)
                                                                                                                                                                                                                   ttttCTACAGACAAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTC‑ggg                                                                                                                                                                                       >  1:308479/1‑101 (MQ=255)
                                                                                                                                                                                                                     ttCTACAGAC‑AAAAAAAAGGCTAGGCATAAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcgttcctcactacgttcatgtcggtctctgcggcagag                                                                                                                                                >  1:501962‑M1/1‑68 (MQ=255)
                                                                                                                                                                                                                     ttCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtctt                                                                                                                                            <  2:3445‑M1/142‑75 (MQ=255)
                                                                                                                                                                                                                     ttCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgat                                                                                                                                                                                                 <  1:216994‑M1/89‑22 (MQ=255)
                                                                                                                                                                                                                       cTACAGACAAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTC‑GGGCCag                                                                                                                                                                                   <  2:308474/101‑1 (MQ=255)
                                                                                                                                                                                                                        tACAGAC‑AAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCCCTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGt                                                                          >  1:215701/1‑203 (MQ=255)
                                                                                                                                                                                                                         aCAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtctt                                                                                                                                            <  2:501956‑M1/138‑75 (MQ=255)
                                                                                                                                                                                                                          cAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgagcccgtcttcggtcatcactacgttcatgtcggtctatgcggcagagtcttcaacgtattccagatcgaaaaccgcttcgccgttcacaattccgacagaaactgcggctaccagcccttccatcgggttggttttcagcttgccgatttcaaccagcttctgtaacacatctaccagagcaacgcagg  >  1:140880‑M1/1‑63 (MQ=255)
                                                                                                                                                                                                                          cAGAC‑AAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTC‑GGGCCAGTGCCAACAAGATGAGTAGctc                                                                                                                                                              >  1:257301/1‑117 (MQ=255)
                                                                                                                                                                                                                           aGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCG                                                                                                                                                                                                                              >  2:288739‑M1/1‑54 (MQ=255)
                                                                                                                                                                                                                            gAC‑AAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTctct                                                                      <  2:215698/203‑1 (MQ=255)
                                                                                                                                                                                                                                caaaaaaaaGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTC‑GGGCCAGTGCCAACAAGATGAGTAGCTCTTCa                                                                                                                                                          <  2:257297/116‑1 (MQ=255)
                                                                                                                                                                                                                                aaaaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTT                                                                                                                                                                                                                          <  1:288743‑M1/54‑1 (MQ=255)
                                                                                                                                                                                                                                  aaaaaaGGCGACTCAACAGTCGACTTAAAAATCAGTTTGCCAGAGCCGAATTCTGCGTCGCTACAATGGATTCGCTTCACCTC‑GGGCCAGAGCCAACAAGATGAGTAGc                                                                                                                                                                >  1:574571/1‑109 (MQ=255)
                                                                                                                                                                                                                                      aaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCtt                                                                                                                                                            <  2:574564/109‑1 (MQ=255)
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                     
AACTCAATGCCGGAATCCACCAACGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGATCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGC  >  CP009273/3808938‑3809429

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