Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 76,185 G→A 12.1% intergenic (+243/‑891) 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,1850GA12.1% 162.8 / 12.9 72intergenic (+243/‑891)USA300HOU_0072/USA300HOU_0075universal stress protein/lysophospholipase
Reads supporting (aligned to +/- strand):  ref base G (34/29);  new base A (4/5);  total (38/34)
Fisher's exact test for biased strand distribution p-value = 7.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.65e-01

ACTAATATTGCTGCACCTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGG  >  CP000730/76053‑76323
                                                                                                                                    |                                                                                                                                          
aCTAATATTGCTGCACCTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTACAGTATAAttt                                                                                                                                    <  2:269923/141‑1 (MQ=255)
   aaTATTGCTGCACCTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTaa                                                                                                                                 <  2:389144/141‑1 (MQ=255)
   aaTATTGCTGCACCTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTaa                                                                                                                                 <  2:22567/141‑1 (MQ=255)
       ttGCTGCACCTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATGTAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACaaa                                                                                                                             >  2:205133/1‑141 (MQ=255)
        tgcttcACCGTTATCACTTATAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAAGAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACaaaa                                                                                                                            <  1:104934/141‑1 (MQ=39)
        tgctgcACCTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACaaaa                                                                                                                            >  1:42667/1‑141 (MQ=255)
         gctgcACCTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACaaaaa                                                                                                                           >  1:300480/1‑141 (MQ=255)
         gctgcACCTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACaaaaa                                                                                                                           >  1:27326/1‑141 (MQ=255)
             cactggTATCACTTAGAACAGCTAAACAAAGAGTTAGGAATGCAGGGAAGACTAAAACAAGGGCAGTTAATTTAATACGTATAGCAAAATAAAGGTGCTGTGGCAAGGTTGAATCAATAGTATAATTTTAACAGAAAGGAg                                                                                                                       >  2:255948/6‑141 (MQ=255)
             cACCTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAg                                                                                                                       <  1:17802/141‑1 (MQ=255)
                cTGTACCATATAGAGCAAATAACAAAAGTGTTAGTAATCCAGGGAATACTTACAGAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTATGTTGCTAAGTTGAATCTATAGTATCATTTTAACAACAAGAAGtct                                                                                                                    >  2:104935/1‑141 (MQ=255)
                  gTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATCTAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGtcttc                                                                                                                  >  1:308542/1‑141 (MQ=255)
                  gTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTCGAGTATAATTTTAACAAAAAGGAGtcttc                                                                                                                  >  2:170610/1‑141 (MQ=255)
                  atATCACTAAAAACAGCTAACCAAAATTTTAGAAATCCAGGGAATACTAACACAAAGACAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTCGCAAAGTTGAATCTATAGTATAATATTAACAAAAAAGAAtcttc                                                                                                                  <  2:52992/140‑1 (MQ=255)
                    aTCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTg                                                                                                                >  2:151423/1‑141 (MQ=255)
                        tttAGAACAGCTAACCAAAATGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATg                                                                                                            <  1:382727/140‑1 (MQ=255)
                        cTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATg                                                                                                            <  1:43621/141‑1 (MQ=255)
                                   tAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAg                                                                                                 >  1:367886/1‑141 (MQ=255)
                                           gtgtTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAAc                                                                                         >  2:288862/1‑141 (MQ=255)
                                                aGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAACCAAttt                                                                                    >  1:169217/1‑141 (MQ=255)
                                                   aaTCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAAc                                                                                 <  1:19403/141‑1 (MQ=255)
                                                   aaTCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAAc                                                                                 <  2:978/141‑1 (MQ=255)
                                                       cAGGGAATACTAACACAAAGGCAATTAATTTAATAATGATAGCAAAATATAGGGTATGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAAGAAATTCACAAgg                                                                             >  2:237099/1‑141 (MQ=21)
                                                                        aaGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGCCTTCTGTATGAACTATTTCAGCTATAAACAATTTACCAAGGATGTTATCACTGTAGcc                                                            >  2:305220/1‑141 (MQ=255)
                                                                         aGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCg                                                           <  1:14287/141‑1 (MQ=255)
                                                                          ggCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGt                                                          <  1:365392/141‑1 (MQ=255)
                                                                            cAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGAAATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTg                                                        >  2:20760/1‑141 (MQ=255)
                                                                             aattaattTAATAATTATAGCAAAATATAGGTTCTGTTGAAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTgg                                                       <  2:79275/141‑1 (MQ=255)
                                                                              attaattTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGc                                                      >  2:258549/1‑141 (MQ=255)
                                                                                 aatttaatAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGctac                                                   <  2:42667/141‑1 (MQ=255)
                                                                                   tttaatAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGctacta                                                 >  2:33105/1‑141 (MQ=255)
                                                                                   tttaatAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGctacta                                                 <  1:135040/141‑1 (MQ=255)
                                                                                     taataaGTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATc                                               <  1:204806/141‑1 (MQ=255)
                                                                                       ataatTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTactaacta                                             >  2:84955/1‑141 (MQ=255)
                                                                                       ataatTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTa                                             >  1:360595/1‑141 (MQ=255)
                                                                                       ataatTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTa                                             >  1:145166/1‑141 (MQ=255)
                                                                                         aatTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTTTAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAg                                           <  1:192098/141‑1 (MQ=255)
                                                                                             ataGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGat                                         >  1:332479/1‑139 (MQ=255)
                                                                                             ataGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGat                                         <  2:332479/139‑1 (MQ=255)
                                                                                               aGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGc                                     >  1:327654/1‑141 (MQ=255)
                                                                                                 cagttaggaGGTTCTGTTGCAAAGTTGAATTTATAATATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCAt                                   >  1:316209/9‑141 (MQ=11)
                                                                                                 cAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCAt                                   <  2:300480/141‑1 (MQ=255)
                                                                                                 agtaaatgaGGTTCTGTTGCAAAGTTGAATTTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCAt                                   >  1:245167/9‑141 (MQ=1)
                                                                                                    ttaggaGGTTCTGTTGCAAAGTTGAATTTATAATATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGa                                <  2:108980/136‑1 (MQ=2)
                                                                                                    aaattcggTTTTGGTGCAATGTTGAATATATAATATAATTTTAACAAAAAGGAGTCTGCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGa                                <  2:315574/135‑1 (MQ=11)
                                                                                                    aaatgaGGTTCTGTTGCAAAGTTGAATTTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGa                                <  2:128291/136‑1 (MQ=9)
                                                                                                     taggaGGTTCTGTTGCAAAGTTGAATTTATAATATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAg                               >  2:194928/5‑141 (MQ=0)
                                                                                                     atatAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAg                               <  1:367491/141‑1 (MQ=40)
                                                                                                     aatgaGGTTCTGTTGCAAAGTTGAATTTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAg                               <  2:351642/137‑1 (MQ=2)
                                                                                                     aatgaGGTTCTGTTGCAAAGTTGAATTTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAg                               >  2:271317/5‑141 (MQ=2)
                                                                                                      atgaGGTTCTGTTGCAAAGTTGAATTTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGt                              <  1:56602/138‑1 (MQ=2)
                                                                                                      atgaGGTTCTGTTGCAAAGTTGAATTTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGt                              >  2:251580/4‑141 (MQ=2)
                                                                                                      aggaGGTTCTGTTGCAACGTTGAATTTATAATATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAgc                              >  1:312135/4‑140 (MQ=255)
                                                                                                      aggaGGTTCTGTTGCAAAGTTGAATTTATAATATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAgc                              >  2:238224/4‑140 (MQ=0)
                                                                                                        gaGGTTCTGTTGCAAAGTTGAATTTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGtaa                            >  1:71858/2‑139 (MQ=0)
                                                                                                        gaGGTTCTGTTGCAAAGTTGAATTTATAATATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAgcta                            <  2:56276/140‑4 (MQ=2)
                                                                                                           gTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCg                         <  1:63373/141‑1 (MQ=35)
                                                                                                            ttatGTTGCAAAGGTGAATTTATAATATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGt                        <  2:249975/138‑1 (MQ=11)
                                                                                                                gTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTCTTTCAGATATAACCAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGata                    <  2:567/141‑1 (MQ=25)
                                                                                                                 tttcAAAGTTCAATTTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACACGCATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGatat                   <  1:332073/138‑1 (MQ=11)
                                                                                                                  tGCAAAGCTGAACGGATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCGGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGatata                  >  2:232916/1‑141 (MQ=11)
                                                                                                                    cAAAGATGAATCTATAGTATAATTTTAACAAGAAGGAGTCCTCTGTATGAACTATTTCAGATATAAGCAATTTAACAAGGATGTAATCACTGCAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATc                >  1:293818/1‑141 (MQ=11)
                                                                                                                      aaGTTGAATTTATAATATAATTTTAACAAAAAGGAGTCTGCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTg              <  1:269560/141‑1 (MQ=255)
                                                                                                                       aGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCGTCTGTATGAACTAGTGCAGATATAAACGATTTAACAAGGATGTTATCACTGTAGTCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGa             <  1:151423/141‑1 (MQ=11)
                                                                                                                             atttataatataaTTTTAACAAAACGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATAtt       <  2:248915/138‑1 (MQ=2)
                                                                                                                             atttataatataaTTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATAtt       >  1:104961/4‑141 (MQ=0)
                                                                                                                             atttataatataaTTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGACGTTGGCTACTATCTAAGGTATGCATTAAGCTATCGTGAGATATCTGAAATAtt       >  1:19223/4‑141 (MQ=1)
                                                                                                                              tttataatataaTTTTAACAAAAAGGAGTCTTATGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCTTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTa      <  1:267466/139‑1 (MQ=11)
                                                                                                                               ttataatataaTTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGTTGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATAtttc     <  2:103483/140‑3 (MQ=255)
                                                                                                                                tataatataaTTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAg    >  1:85716/1‑141 (MQ=17)
                                                                                                                                tataGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAg    >  1:43093/1‑141 (MQ=17)
                                                                                                                                tataGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAAGATTAAg    >  1:135724/1‑141 (MQ=17)
                                                                                                                                 ataatataatTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAg  <  2:130588/141‑2 (MQ=12)
                                                                                                                                 ataGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAgg   >  2:346943/1‑141 (MQ=17)
                                                                                                                                 ataGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAgg   <  2:355281/141‑1 (MQ=32)
                                                                                                                                  taGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAggg  >  2:125226/1‑141 (MQ=17)
                                                                                                                                  taGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAggg  >  1:313882/1‑141 (MQ=17)
                                                                                                                                  taGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAgga  >  2:82115/1‑140 (MQ=17)
                                                                                                                                  taGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAggg  <  2:367886/141‑1 (MQ=32)
                                                                                                                                    |                                                                                                                                          
ACTAATATTGCTGCACCTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGG  >  CP000730/76053‑76323

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: