Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 76,180 C→T 24.6% intergenic (+238/‑896) 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,1800CT24.6% 178.7 / 42.2 87intergenic (+238/‑896)USA300HOU_0072/USA300HOU_0075universal stress protein/lysophospholipase
Reads supporting (aligned to +/- strand):  ref base C (36/28);  new base T (8/13);  total (44/43)
Fisher's exact test for biased strand distribution p-value = 2.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

ACTAATATTGCTGCACCTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAAT  >  CP000730/76053‑76315
                                                                                                                               |                                                                                                                                       
aCTAATATTGCTGCACCTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAAttt                                                                                                                            >  2:444457/1‑141 (MQ=255)
 cTAATATTGCTGCACCTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAAtttt                                                                                                                           >  1:669563/1‑141 (MQ=255)
   aaaagtgctgcACGTGTATCACAGAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTaa                                                                                                                         <  1:693779/136‑1 (MQ=255)
   aaTATTGCTGCACCTTTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTaa                                                                                                                         <  2:689447/141‑1 (MQ=255)
    atatTGCTGCACCTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTTTAGTATAATTTTAAc                                                                                                                        <  1:384232/141‑1 (MQ=255)
     tatTGCTGCACCTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACa                                                                                                                       <  1:661745/141‑1 (MQ=255)
       ttGCTGCACCTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACaaa                                                                                                                     >  2:249248/1‑141 (MQ=255)
       ttGCTGCACCTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACaaa                                                                                                                     >  2:692221/1‑141 (MQ=255)
        tgctgcACCTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACaaaa                                                                                                                    <  2:11042/141‑1 (MQ=255)
             cACCTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAg                                                                                                               >  1:388397/1‑141 (MQ=255)
               ccTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGtc                                                                                                             >  1:67578/1‑141 (MQ=255)
                cTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGtct                                                                                                            <  1:387190/141‑1 (MQ=255)
                cTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGtct                                                                                                            >  1:84480/1‑141 (MQ=255)
                 tGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACGAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGtctt                                                                                                           >  1:510236/1‑141 (MQ=255)
                  gTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGtcttc                                                                                                          >  1:162603/1‑141 (MQ=255)
                    aTCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTGGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTTTg                                                                                                        >  2:499705/1‑141 (MQ=255)
                    aTCACTTAGAACAGATAAACAAAGTGGTAGTAATTCAGGGAATACTAACACAAAGGCAAGTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTAAAAAAAAAAAAAAAAGGAGTCTTCTg                                                                                                        <  2:317144/141‑1 (MQ=37)
                        cTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATGATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATg                                                                                                    <  2:654539/141‑1 (MQ=255)
                             aaCAGCTAACCAAAGTGTAAGTAATCCAGGGAAGACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTa                                                                                               <  1:327144/141‑1 (MQ=255)
                              aCAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTAt                                                                                              >  2:596054/1‑141 (MQ=255)
                                aGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTAttt                                                                                            >  1:32360/1‑141 (MQ=255)
                                  cTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCa                                                                                          >  2:469406/1‑141 (MQ=255)
                                     aCCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGat                                                                                       >  1:453041/1‑141 (MQ=255)
                                     aCCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGat                                                                                       >  1:312002/1‑141 (MQ=255)
                                     aCCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGat                                                                                       >  1:639875/1‑141 (MQ=255)
                                      ccAAAGTGTTAGTAGTCCAGGGAATACTGACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTGTTTCAGata                                                                                      >  1:501241/1‑141 (MQ=255)
                                      ccAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGaaa                                                                                      >  1:123508/1‑139 (MQ=255)
                                          aGTGTTAGTAATCCAAGAAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATaaa                                                                                  <  2:691903/141‑1 (MQ=255)
                                           gtgtTAGTAATTCTGGGAATACTAACACAAATGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTTCAAAGTTGAATCTATAGTATAATTTTAACAAATAGGAGTCTTCTGTATGAACTATTTCAGATATAAAc                                                                                 <  2:93314/141‑1 (MQ=25)
                                           gtgtTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGTGTCTTCTGTATGAACTATTTTGGATATacac                                                                                 >  1:11387/1‑141 (MQ=255)
                                           gtgtTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAAc                                                                                 <  1:473006/141‑1 (MQ=255)
                                           gtgtTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAAc                                                                                 >  1:471989/1‑141 (MQ=255)
                                            tgtTACTAATCCACGGAATCATAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATCGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACa                                                                                <  1:60322/141‑1 (MQ=21)
                                                aGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACCAAAAGGAGTCTTCTGTATGAACTATGTCAGATATAAACAAttt                                                                            >  1:683555/1‑141 (MQ=255)
                                                aGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAAttt                                                                            >  1:606151/1‑141 (MQ=255)
                                                 gTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAACTTAAAGAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTa                                                                           <  1:337385/141‑1 (MQ=255)
                                                       cAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAgg                                                                     <  2:189957/141‑1 (MQ=255)
                                                          ggAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATg                                                                  <  2:233256/141‑1 (MQ=255)
                                                                cTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTAGAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCa                                                            >  2:654920/1‑141 (MQ=255)
                                                                 tAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTGGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGACTTCTGTATGAACTATTTCACATATAAACATTTTAACCAGGATGTTATAAc                                                           >  1:620034/1‑141 (MQ=21)
                                                                  aaCACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACt                                                          <  1:489786/141‑1 (MQ=255)
                                                                  aaCACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACt                                                          <  1:134185/141‑1 (MQ=255)
                                                                        aaGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGcc                                                    <  1:202165/141‑1 (MQ=255)
                                                                          ggCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGTCGt                                                  >  2:629494/1‑141 (MQ=255)
                                                                            cAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTg                                                <  2:669563/141‑1 (MQ=255)
                                                                              attaattTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGc                                              >  1:304897/1‑141 (MQ=255)
                                                                              attaattTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGc                                              >  1:307930/1‑141 (MQ=255)
                                                                               ttaatttaatAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGct                                             <  1:387796/141‑1 (MQ=255)
                                                                                taatttaatAATTAGAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGcta                                            <  2:453041/141‑1 (MQ=255)
                                                                                 aatttaatAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGctac                                           <  2:639875/141‑1 (MQ=255)
                                                                                      aataatTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCt                                      >  2:361485/1‑141 (MQ=255)
                                                                                            tataGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGata                                >  1:648921/1‑141 (MQ=255)
                                                                                               aGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGc                             >  2:623496/1‑141 (MQ=255)
                                                                                                 cAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATACGctt                           <  2:388397/141‑3 (MQ=255)
                                                                                                 agtaaatgaGGTTCTGTTGCAAAGTTGAATTTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCAt                           <  1:626168/133‑1 (MQ=11)
                                                                                                   taaatgaGGTTCTGTTGCAAAGTTGAATTTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTg                         >  2:24722/7‑141 (MQ=2)
                                                                                                   taaatgaGGTTCTGTTGCAAAGTTGAATTTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTg                         <  1:332793/135‑1 (MQ=255)
                                                                                                    taaggaGGTTCTGTTGCAAAGTTGAATTTATAATAAAAAAAAAAAAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGa                        <  2:706117/136‑1 (MQ=11)
                                                                                                    aaatgaagtacTGTTGCAAAGTTGAATTTATAGTATAATATTAACAAAAAGGAGTATTCTGTATGAACTATTTCAGATACAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACAATATAAGATATGCATTGa                        >  2:69417/11‑141 (MQ=1)
                                                                                                     taggaGGTTCTGTTGCAAAGTTGAATTTATAATATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTTACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAg                       <  2:172971/137‑1 (MQ=255)
                                                                                                     taggaGGTTCTGTTGCAAAGTTGAATTTATAATATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAg                       <  2:663835/137‑1 (MQ=2)
                                                                                                     aatgaGGTTCTGTTGCAAAGTTGAATTTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAg                       <  1:443328/137‑1 (MQ=9)
                                                                                                     aatgaGGTTCTGTTGCAAAGTTGAATTTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAg                       >  1:5409/5‑141 (MQ=2)
                                                                                                      atgaGGTTCTGTTGCAAAGTTGAATTTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGt                      <  2:354534/138‑1 (MQ=2)
                                                                                                      atgaGGTTCTGTTGCAAAGTTGAATTTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGt                      <  1:279820/138‑1 (MQ=9)
                                                                                                      atgaGGTTCTGTTGCAAAGTTGAATTTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGt                      <  2:532222/138‑1 (MQ=9)
                                                                                                      atgaGGTTCTGATGCAAAGTTGAATTTATAGTATAATTCTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGt                      <  1:519122/138‑1 (MQ=11)
                                                                                                      aggaGGTTCTGTTGCAAAGTTGAATTTATAATATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAgc                      <  1:39531/138‑2 (MQ=2)
                                                                                                       tgaGGTTCTGTTGCAAAGTTGAATTTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGtt                     <  1:160034/139‑1 (MQ=9)
                                                                                                       tgaGGTTCTATTGCAAAGTTGAATTTAAAGTATAATTTTAACAAAAAGGCGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACATGGATGTTATCACTGTAGCTTGGGGCTACTATCTAAGATATGCATTGAGtt                     <  1:756837/139‑1 (MQ=11)
                                                                                                       ggaGGTTCTGTTGCAAAGTTGAATTTATAATATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAgct                     >  2:194146/3‑139 (MQ=0)
                                                                                                       ggaGGTTCTGTTGCAAAGATGTATTAATAATATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAgct                     >  2:738968/3‑139 (MQ=11)
                                                                                                       ataGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGtt                     <  2:510236/141‑1 (MQ=38)
                                                                                                        gaGGTTCTGTTGCAAAGTTGAATTTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTa                    >  1:715440/2‑141 (MQ=2)
                                                                                                         aGGTTCTGTTGCAAAGTTGAATTTATAATATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGTATTGAGCTAt                   >  1:79949/1‑141 (MQ=0)
                                                                                                         aGGTTCTGTTGCAAAGTTGAATTTATAATATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTAt                   >  2:140610/1‑141 (MQ=0)
                                                                                                         aGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTGAGATATGCTTTGAGTTAt                   >  1:58374/1‑141 (MQ=25)
                                                                                                          ggTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATc                  <  1:494547/141‑1 (MQ=35)
                                                                                                          ggTTCTGGTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAAAAATTTAACAAGGATATTCTCACAGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATc                  >  2:732202/1‑141 (MQ=11)
                                                                                                           gTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCg                 >  2:613978/1‑141 (MQ=35)
                                                                                                            ttCTGTTGCAAAGTTGACTCTATAGTATAATTTTAACAAACAGGAGTTCTCTTTATGAACTATTTCAATTGCCCTCTTTTTAAAATAACTGATGTCTGTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGt                <  1:754485/141‑1 (MQ=2)
                                                                                                            ttCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATGAATGAAGCCGATGGCTACTATCTAAGATATGCATTGAGTTATCGt                >  1:330090/1‑141 (MQ=255)
                                                                                                              ttggtGCAAAGTTGAATGTATTGTATAATTTTAACAAAAAGGACTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGGTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGa              <  2:30574/137‑1 (MQ=11)
                                                                                                                tttGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGata            <  1:444457/140‑1 (MQ=25)
                                                                                                                    cAAAGTTGAATATATAATATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATc        <  2:660882/141‑1 (MQ=12)
                                                                                                                      aaGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTg      >  2:482196/1‑141 (MQ=35)
                                                                                                                          tGAATCTATAATATAACCTCAGCTAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAGGCTATGCATTCAGCTACCGTGATATAACTGTAat  >  2:545202/1‑141 (MQ=11)
                                                                                                                               |                                                                                                                                       
ACTAATATTGCTGCACCTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAAT  >  CP000730/76053‑76315

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