Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 35,850 C→T 51.7% intergenic (+81/+27) USA300HOU_0025 → / ← USA300HOU_0026 hypothetical protein/IS431mec transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00073035,8500CT51.7% 34.0 / 78.3 67intergenic (+81/+27)USA300HOU_0025/USA300HOU_0026hypothetical protein/IS431mec transposase
Reads supporting (aligned to +/- strand):  ref base C (18/14);  new base T (14/20);  total (32/35)
Fisher's exact test for biased strand distribution p-value = 3.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

TGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACT‑‑TGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATA  >  CP000730/35715‑35976
                                                                                                                                       |                                                                                                                                
tGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATGAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTa                                                                                                                             >  2:345351/1‑141 (MQ=255)
 gACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTaa                                                                                                                            >  2:199838/1‑141 (MQ=255)
 gACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTaa                                                                                                                            >  1:65327/1‑141 (MQ=255)
            tAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATGACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTAATCATGTc                                                                                                                 >  2:211369/1‑141 (MQ=255)
                   gAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAGAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTc                                                                                                          >  1:131790/1‑141 (MQ=255)
                         aTTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCtcaa                                                                                                    <  1:84301/141‑4 (MQ=255)
                         aTTCAGTATAAAGTTATTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAAAAGGTTCTTTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTaa                                                                                                    <  2:30590/141‑1 (MQ=255)
                            cAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACt‑‑t                                                                                               >  1:79028/1‑141 (MQ=255)
                             aGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACT‑‑Tg                                                                                              >  1:361765/1‑141 (MQ=255)
                               tataAAGTTATTGGCGAAAAATGCTTACTTGCACAAAAGCTACAATTAAATTGTAAAGCAAAGTGAGTTGAAGGTTCTGTTGCAGATTTTGAAAATGTGGCGATGCTCTAATTTATCATGTCAGTGTTCGCTTAACT‑‑TGCt                                                                                            <  2:390835/141‑1 (MQ=12)
                                ataAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAACTAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCAATAATTTATCATGTCAGTGTTCGCTTAAAT‑‑CGCTa                                                                                           >  2:152762/1‑141 (MQ=255)
                                ataAAGTTTTTGGCGAAAAATGATTACTTGAAAAAGAGCTAAAATAAAATTGAGATCTAAAAAAAATAGAAGGTTATGTTGCAAAGTAAAAAAATATAGATAACCACTAATTTATCATGTAAGTGATAGATACACT‑‑TGCTa                                                                                           >  2:152820/1‑141 (MQ=25)
                                agaaaGTTTTTGGCGAAAAATGAGTACGTGAAAAAGAGCTAAAGTAAAATTGTGATCAAATAAAAATAGAAGGTACTGTTGCAAAGAAAAAAAATATAGCTGACCACTAATATAACATGTCAGTGTTCGCATAACtgactc                                                                                             >  2:152768/3‑136 (MQ=255)
                                        tttGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACT‑‑TGCTAGCatgatg                                                                                   >  2:33293/1‑141 (MQ=255)
                                        tttGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACT‑‑TGCTAGCatgatg                                                                                   >  1:51269/1‑141 (MQ=255)
                                           ggCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACT‑‑TGCTAGCATGATGCTa                                                                                <  2:103130/141‑1 (MQ=255)
                                              gAAAAATGATTACTTTAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACT‑‑TGCTAGCATGATGCTAAtt                                                                             >  1:294749/1‑141 (MQ=255)
                                               aaaaaaGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCGAATAACAATAGACGGTTCTGTTGCCAAGTACAAAAATATAGCTAACCACTGATTTATCATGTCAGTGTTCGCTTAACT‑‑TGCTAGCATGATGCTAAttt                                                                            <  1:15612/141‑1 (MQ=21)
                                                        tACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCATAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACT‑‑TGCTAGCATGATGCTAATTTCGTGGCATg                                                                   >  2:220994/1‑141 (MQ=255)
                                                         aCTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATAATGTCAGAGTACGCTTAACT‑‑TGCTAGCATGATGCTAATTTCGTGGCATgg                                                                  >  2:88980/1‑141 (MQ=255)
                                                          cTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTACCCACTAATTTAGCATGTCAGTGTCTGCTTAACT‑‑TGCTAGCATGATGCTAATTTCatggcatggc                                                                 <  1:289646/141‑1 (MQ=21)
                                                            ttaaCAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACT‑‑TGCTAGCATGATGCTAATTTCGTGGCATGGCGa                                                               <  2:50736/139‑1 (MQ=255)
                                                            tGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCTTGTCAGTGTTCGCTTAACT‑‑TGCTAGCATGTTGCTATTTTCGTGGCATGGCGa                                                               <  1:197718/141‑1 (MQ=255)
                                                                           acaaaagagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACT‑‑TGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGaa                                              >  2:363379/15‑141 (MQ=11)
                                                                           acaaaagagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACT‑‑TGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGaa                                              >  2:387443/15‑141 (MQ=11)
                                                                            taaaatTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACT‑‑TGCTAGCATGATGCTAATTTCGTGGCATGGCGAGAATCCGTAGATCTGa                                               >  2:74608/1‑141 (MQ=255)
                                                                            caaaggagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACT‑‑TGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAag                                             >  1:268552/14‑141 (MQ=11)
                                                                             aaaagagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACT‑‑TGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAaga                                            <  2:369002/129‑1 (MQ=11)
                                                                              aaagagtattgtaataaAAATGATGATTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGCTCGATTAGCT‑‑TACTATCATGATGCTAATTTCGTGGCATGACAAAAGTCCGTAGATCTGAagag                                           >  2:337085/12‑141 (MQ=2)
                                                                               aagagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACT‑‑TGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAgaga                                          >  1:311166/11‑141 (MQ=11)
                                                                               aagagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACT‑‑TGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAgaga                                          <  2:188235/131‑1 (MQ=11)
                                                                                agagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCG‑‑TAACTGCCCCGAGCATGATGCTTGAATCTTGGCATGGCGAAAATCCGTAGATTTGAAGAGAc                                         <  1:289441/132‑1 (MQ=2)
                                                                                agagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGGCAGTGTTCGTTTAACT‑‑TGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAc                                         <  1:178887/132‑1 (MQ=11)
                                                                                 gagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACT‑‑TGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAcc                                        <  2:29740/133‑1 (MQ=11)
                                                                                   gtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACT‑‑TGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTg                                      <  2:270921/135‑1 (MQ=11)
                                                                                   gtattgtaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACT‑‑TGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTg                                      <  1:59417/135‑1 (MQ=11)
                                                                                    tgATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACT‑‑TGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAAACAGTAGATCTGCagaaaact                                       >  2:349245/1‑136 (MQ=255)
                                                                                    tattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACT‑‑TGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGc                                     >  2:274299/6‑141 (MQ=11)
                                                                                    gattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACT‑‑TGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGc                                     <  2:93929/136‑1 (MQ=11)
                                                                                     attttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACT‑‑TGCTCGCCTGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCg                                    <  1:391118/137‑1 (MQ=11)
                                                                                     attttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACT‑‑TGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCg                                    <  1:357293/137‑1 (MQ=11)
                                                                                     attttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACT‑‑TGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCg                                    <  1:157762/137‑1 (MQ=11)
                                                                                      ttttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACT‑‑TGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGTCCTGCgg                                   >  2:385487/4‑141 (MQ=11)
                                                                                      ttttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACT‑‑TGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCgg                                   <  2:59277/138‑1 (MQ=11)
                                                                                      ttttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACT‑‑TGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCgg                                   <  2:365927/138‑1 (MQ=11)
                                                                                       tttaataaAAATGATGGTTCTGTTGCAAAGTCAAAAATATAGCTAATCACTAATATATCATGTGAGTGTTCGTTTAACT‑‑TAGTAGCATGATGCTAATTTCGTGGGATGGGAAAAATCGGTAGGTATGAGAAGACCAGCggg                                  >  2:3582/3‑140 (MQ=2)
                                                                                       tttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACT‑‑TGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGt                                  <  2:88204/139‑1 (MQ=11)
                                                                                       tttaataaAAATGATGGTTCTGGTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACT‑‑TGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGt                                  <  2:235082/139‑1 (MQ=11)
                                                                                       tCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACT‑‑TGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCg                                    <  2:131790/141‑1 (MQ=255)
                                                                                         taataaAAATAGAAGGGTCTGTTGCAAAGTAAAAAAATATAGCTGACCCCTAATTTATCATGTCCGTGTTCGCTTAACT‑‑TGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGt                                  <  1:215785/141‑1 (MQ=255)
                                                                                          aataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTACTATATCATGTCAGTGTTCATAAAACT‑‑TAATAGCAAGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGCACTGGAAAGACCTGCCGTCCt                               >  2:238377/1‑141 (MQ=2)
                                                                                          aataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACT‑‑TGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCt                               >  2:246254/1‑141 (MQ=11)
                                                                                          aataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACT‑‑TGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCt                               <  1:213740/141‑1 (MQ=11)
                                                                                           ataaAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACT‑‑TGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTc                                >  1:134379/1‑141 (MQ=38)
                                                                                           ataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACT‑‑TGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCtt                              >  2:166983/1‑141 (MQ=11)
                                                                                           ataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACT‑‑TGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCtt                              >  1:215423/1‑141 (MQ=11)
                                                                                            taaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACT‑‑TGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCttt                             <  2:282856/141‑1 (MQ=11)
                                                                                            taaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACT‑‑TGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCttt                             <  1:311184/141‑1 (MQ=11)
                                                                                            taaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACT‑‑TGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCttt                             >  1:394659/1‑141 (MQ=11)
                                                                                            taaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACT‑‑TGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCttt                             <  1:344232/141‑1 (MQ=11)
                                                                                            taaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACT‑‑TGCTAGCATGATACTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACTTGCGGTTCttt                             >  1:123703/1‑141 (MQ=11)
                                                                                              aaaaTAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACT‑‑TGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGAACTGAAGAGACCTGCGGTTCttt                             <  2:288962/141‑1 (MQ=25)
                                                                                              aaaaGAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACT‑‑TGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCttt                             <  1:313010/141‑1 (MQ=18)
                                                                                                     aaGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACT‑‑TGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTtatata                      <  2:185097/141‑1 (MQ=11)
                                                                                                         ttatGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACT‑‑TGCTAGCATGATGCTAATTTCGTGACATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCg                 <  1:184088/138‑1 (MQ=1)
                                                                                                          tccgTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACT‑‑TGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCg                 <  2:346012/138‑1 (MQ=1)
                                                                                                                         aacaaTATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACT‑‑TGCTAGCATGGTGCTAATTTCGTGGCATGGCGAAAACCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATa  <  2:8910/138‑1 (MQ=0)
                                                                                                                                       |                                                                                                                                
TGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACT‑‑TGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATA  >  CP000730/35715‑35976

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: