Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 76,889 G→A 9.1% intergenic (+947/‑187) USA300HOU_0072 → / → USA300HOU_0075 universal stress protein/lysophospholipase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00073076,8890GA 9.1% 194.9 / 5.5 76intergenic (+947/‑187)USA300HOU_0072/USA300HOU_0075universal stress protein/lysophospholipase
Reads supporting (aligned to +/- strand):  ref base G (32/37);  new base A (2/5);  total (34/42)
Fisher's exact test for biased strand distribution p-value = 4.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.66e-01

GGTATCAAAATATCAATACAGCAAAGAATACTTTAAAAGGTATTGAATGTATTTACGCTCTATATAAAAAGAACCGCAGGTCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAA  >  CP000730/76749‑77029
                                                                                                                                            |                                                                                                                                            
ggTATCAAAATATCAATACAGCAAAGAATACTTTAAAAGGCATTGAATGTATTTACGCCCTATATAAAAAGAACCGCAGGTCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAg                                                                                                                                              >  2:932322/1‑141 (MQ=17)
 gTATCCAAATATCAATACAGCAAAGAATCCCTTAAAANGTATTGCATGGATTTACGCTCTATATCAGAAAAATCGTGTGTCCCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTCAGc                                                                                                                                             <  2:105213/141‑1 (MQ=11)
  gaTCAAAATATCAATACAGCAAAGAATACTTTAAAAGGTATTGAATGTATTTACGCTCTATATAAAAAGAACCGCAGGTCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCg                                                                                                                                            <  2:152743/140‑1 (MQ=17)
    tCAAAATATCAATACAGCAAAGAATACTTTAAAAGGTATTGAATGTATTTACGCTCTATATAAAAAGAACCGCAGGTCTCTTCAAATCTACGGATTTTCGCCATGCCACGAAATTAGCCTCATGCTAGCAAATTAAGCGaa                                                                                                                                          >  2:419247/1‑141 (MQ=12)
     aaaaaTAACAATACGGCAAAGAATACTTTAAAAGGTATTGAATGTATTTACGCTCTATATAAAAAGAACCGTAGGTCTCTTCAGATCAACGGATTTTCCCCATGCCACGAAATTACCCTCATGCTACCAATTTAAACGAac                                                                                                                                         <  2:319199/140‑1 (MQ=11)
      aaaaTATCAATACAGCAAAGAATACTTTAAAAGGTATTGAATGTATTTACGCTCTATATAAAAAGAACCGCAGGTCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAaca                                                                                                                                        >  1:578222/1‑141 (MQ=32)
      aaaaTATCAATACAGCAAAGAATACTTTAAAAGGTATTGAATGTATTTACGCTCTATATAAAAAGAACCGCAGGTCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAaca                                                                                                                                        >  1:606114/1‑141 (MQ=32)
      aaaaTATCAATACAGCAAAGAATACTTTAAAAGGTATTGAATGTATTTACGCTCTATATAAAAAGAACCGCAGGTCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAaca                                                                                                                                        >  1:12001/1‑141 (MQ=32)
       aaggataaaTACGGTAAAGAATACTTTAAAAGGTATTGAATGTATTTACTCTCTGTATAAAAAGAACCGTAGGTCTCTTCAGATCTACGGATGTGCCCCATTCCACGAAATTAGCATCAGGCTAGCAAGTTAAACGAacac                                                                                                                                       <  2:155103/134‑1 (MQ=11)
        aaTAGCAACACAGCAAAGAATATTTTAAAAGGTATTGAATGTATTTAAAGCCTTGAGAAAAAGAACCGCAGGTCTCTTCAGAGCTACGGATTGTCGCCATGCCACGAAATTAGCATCAGGCTAGCACGTTAAGCGAACACt                                                                                                                                      <  1:492677/141‑1 (MQ=11)
           ataaaTACAGCAAAGAATACTTTAAAAGGTATTGAGTGTATTTACGCTCTATATAAAAAGAACCGCAGGTCTCTTCAGATCTCCGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGAc                                                                                                                                   <  2:386545/138‑1 (MQ=0)
                    gaaaaaagtACTTTAGAAGGGATGGAATGGATTGACGCTCCATATAAAAAGGACCGTAGGTCTCTTCAGATCTACGGATTTTCCCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGatatat                                                                                                                          <  2:44680/133‑1 (MQ=1)
                         gactACTTTAAAAGGTATTGAGTGTATTTACGCTCTATATAAAAAGAACCGCAGGTCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTg                                                                                                                     <  1:717668/138‑1 (MQ=1)
                               ttgaaaaGGTATTGTGTGTATTTACGCTCTATATAAAAAGAACCGCAGGTCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTACCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGc                                                                                                               <  1:674059/138‑1 (MQ=0)
                                   aacggTATTGAATGTATTTACGCTCTATATAAAAAGAACCGTAGGTCTCTTCAGATCTACGGATTTTCCCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGCTTAGCtata                                                                                                           >  2:428182/4‑141 (MQ=1)
                                        ttttgatggTATTTACGCTCTATATAAAAAGAACCGCAGGTCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATAttttt                                                                                                      <  2:657346/133‑1 (MQ=1)
                                            gagtGTATTTACGCTCTATATAAAAAGAACCGCAGGTCTCTTCCGATCTGCGGGGTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTAc                                                                                                  <  1:444562/138‑1 (MQ=1)
                                                             atataACAAGAACCGCAGGTCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCtat                                                                                 <  2:168359/141‑1 (MQ=14)
                                                             atataAAAAGAACCGCAGGTCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCtat                                                                                 <  1:827319/141‑1 (MQ=31)
                                                               ataAAAGGAACCGCAGGTCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATaa                                                                               <  1:210900/141‑1 (MQ=17)
                                                               ataAAAAGAACCGCAGGTCTCTTCAGTTCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTGCTTTGCAACAGAACCTATaa                                                                               <  1:572275/141‑1 (MQ=17)
                                                                 aaaaaGAACCGCAGGTCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATaaaa                                                                             <  1:380661/141‑1 (MQ=34)
                                                                  aagaGAACCGCAGGTCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTCACACGATAAATTAGTGATTAGCTATA‑TTTTTAACTTTGAAACAGAActattattttt                                                                           >  2:80889/4‑131 (MQ=1)
                                                                  aaaaGAACCGCAGGTCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACAGTATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATaaaaa                                                                            >  1:780214/1‑141 (MQ=25)
                                                                    aggAACCGCAGGGCCCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAAc                                                                          <  2:415489/139‑1 (MQ=255)
                                                                    aaGAGCCGCAGGTCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAAc                                                                          <  1:78946/141‑1 (MQ=37)
                                                                    aaGAACCGCAGGTCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAAc                                                                          <  1:827594/141‑1 (MQ=37)
                                                                      gAACCGCAGGTCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAACCAt                                                                        >  1:20888/1‑141 (MQ=25)
                                                                       aacaacaaatttcgtCATATCTACGTATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAAACTAAAAAAAACATg                                                                       <  2:192393/127‑1 (MQ=11)
                                                                          cGCAGGTCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGc                                                                    <  1:200751/141‑1 (MQ=255)
                                                                          cGCAGGTCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGc                                                                    >  1:224141/1‑141 (MQ=255)
                                                                           gCAGGTCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCt                                                                   <  2:94638/141‑1 (MQ=255)
                                                                             aGGTCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTa                                                                 >  2:895044/1‑141 (MQ=255)
                                                                              ggTCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTaa                                                                <  2:12001/141‑1 (MQ=255)
                                                                              agTCTCTTCAGAACTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTaa                                                                <  2:545140/140‑1 (MQ=255)
                                                                               gtctgtgcataccTACGGATTTTCGCAATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTaaa                                                               <  2:676159/129‑1 (MQ=14)
                                                                               gTCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCCAACGAACCTATAAAAAAACTGGGCTTaaa                                                               >  2:114267/1‑141 (MQ=14)
                                                                               gTCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTtaa                                                               >  2:848779/1‑141 (MQ=255)
                                                                               gTCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTaaa                                                               >  2:335794/1‑141 (MQ=255)
                                                                               gTCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACATATAAAAAACATGTCCTTaaa                                                               >  2:217583/1‑141 (MQ=255)
                                                                               gTCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGAAATGATAAATTCGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTTTAAAAAACATGGGCTTaaa                                                               >  2:372186/1‑141 (MQ=255)
                                                                                 ctctccacaacaaCGGATATTCGTCATGTCACGAAATTAACATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAAAAAATAAAAAAAACATGGGCTTAAAgg                                                             <  1:795567/129‑1 (MQ=14)
                                                                                   ctTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCGGCAGAAGCTATAAAAAAAATGGGCTTAAAGGGt                                                           >  1:404310/1‑141 (MQ=255)
                                                                                   ctTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAAGATGGGCTTAAAGGGt                                                           <  1:700996/141‑1 (MQ=255)
                                                                                   ctTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGt                                                           >  1:460950/1‑141 (MQ=255)
                                                                                    ttccgaaccaCGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGtt                                                          <  1:753335/132‑1 (MQ=255)
                                                                                    ttCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGtt                                                          >  2:783780/1‑141 (MQ=255)
                                                                                     tCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTCCAGGGTTa                                                         >  1:254320/1‑141 (MQ=255)
                                                                                      cagttctacttcttttCGCCATGGCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTaa                                                        <  1:210865/129‑1 (MQ=17)
                                                                                        gATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAAtt                                                      <  1:539324/141‑1 (MQ=255)
                                                                                         aTCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAAttt                                                     >  2:274951/1‑141 (MQ=255)
                                                                                             aCGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAAAATATAAAAAACATGGGCTTAAAGGGTTAATTtatat                                                 >  1:725837/1‑141 (MQ=255)
                                                                                              cGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGATAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTtatata                                                <  1:880668/141‑1 (MQ=255)
                                                                                                 aTTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAGAAAACATGGGCTTAAAGGGTTAATTTATATAGTa                                             <  1:932322/141‑1 (MQ=255)
                                                                                                 aTTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTa                                             >  2:1565/1‑141 (MQ=255)
                                                                                                 aTTTTCGCCATGCCACGAAATCAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATATTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTa                                             >  1:485694/1‑141 (MQ=255)
                                                                                                  ttttCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTAt                                            >  1:603667/1‑141 (MQ=255)
                                                                                                   tttCGCCATGCCACGAAATGAGACTCATGCTAGCACGATAAGCGAACACTGACTTGTCAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTAtt                                           <  1:216839/141‑1 (MQ=18)
                                                                                                     tagCCATGCCACGAAATTAGCAGCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTAc                                         <  2:128943/139‑1 (MQ=255)
                                                                                                       gCCCTGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTGGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGGATAAAGGGGTAAATTATATAGTATTACtt                                       >  2:576982/1‑140 (MQ=17)
                                                                                                       gCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTCCTg                                       >  1:5175/1‑141 (MQ=255)
                                                                                                       gCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTg                                       >  1:623460/1‑141 (MQ=255)
                                                                                                         cATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTCCTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTGCTGCa                                     >  2:3195/1‑141 (MQ=255)
                                                                                                         cATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCa                                     <  1:895044/141‑1 (MQ=255)
                                                                                                             ccACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACa                                 <  1:335794/141‑1 (MQ=255)
                                                                                                                  aaaTTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTCGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAACAATGGGCTTAAAGGGTTCGTTATTATAGTATTAATGCAAACCATGAt                            >  2:838664/1‑141 (MQ=17)
                                                                                                                     ttAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTg                         <  1:232705/141‑1 (MQ=255)
                                                                                                                       aGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTggg                       >  2:754154/1‑141 (MQ=255)
                                                                                                                         catcatGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTc                     <  1:60375/141‑1 (MQ=255)
                                                                                                                          ctcatGCGAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGGTAGATATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGaaaagaaaaa                    <  2:341005/140‑10 (MQ=25)
                                                                                                                          atcatGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCa                    <  1:382806/141‑1 (MQ=255)
                                                                                                                           tcatGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCaa                   <  1:383459/141‑1 (MQ=255)
                                                                                                                           ccatGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTCAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCaa                   <  2:550945/140‑1 (MQ=255)
                                                                                                                                 tAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAAtaa             >  2:718421/1‑139 (MQ=255)
                                                                                                                                      aGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGt        >  1:320767/1‑141 (MQ=255)
                                                                                                                                        ttAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGtgt      >  2:29863/1‑139 (MQ=255)
                                                                                                                                          aaGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTa    >  2:312136/1‑141 (MQ=255)
                                                                                                                                            gCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTaaa  >  1:309257/1‑141 (MQ=255)
                                                                                                                                            |                                                                                                                                            
GGTATCAAAATATCAATACAGCAAAGAATACTTTAAAAGGTATTGAATGTATTTACGCTCTATATAAAAAGAACCGCAGGTCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAA  >  CP000730/76749‑77029

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