Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 35,858 T→A 56.8% intergenic (+89/+19) 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,8580TA56.8% 25.6 / 77.8 67intergenic (+89/+19)USA300HOU_0025/USA300HOU_0026hypothetical protein/IS431mec transposase
Reads supporting (aligned to +/- strand):  ref base T (14/15);  new base A (16/22);  total (30/37)
Fisher's exact test for biased strand distribution p-value = 6.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

TAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACT‑‑TGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTA  >  CP000730/35727‑35988
                                                                                                                                   |                                                                                                                                    
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:387443/15‑141 (MQ=11)
                                                               acaaaagagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACT‑‑TGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGaa                                                          >  2:363379/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                                                      <  2:188235/131‑1 (MQ=11)
                                                                   aagagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACT‑‑TGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAgaga                                                      >  1:311166/11‑141 (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                                          >  1:215423/1‑141 (MQ=11)
                                                                               ataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACT‑‑TGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCtt                                          >  2:166983/1‑141 (MQ=11)
                                                                                taaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACT‑‑TGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCttt                                         <  1:344232/141‑1 (MQ=11)
                                                                                taaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACT‑‑TGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCttt                                         <  1:311184/141‑1 (MQ=11)
                                                                                taaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACT‑‑TGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCttt                                         <  2:282856/141‑1 (MQ=11)
                                                                                taaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACT‑‑TGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCttt                                         >  1:394659/1‑141 (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)
                                                                                                                      gccaatcACTAATATATCATGTCAGTGTTCGTTTAACT‑‑TGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTaaa     <  1:36263/135‑1 (MQ=1)
                                                                                                                         aatcACTAATATATCATGTCAGTGTTCGTTTAACT‑‑TGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTa  >  2:389990/4‑141 (MQ=1)
                                                                                                                         aatcACTAATATATCATGTCAGTGTTCGTTTAACT‑‑TGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTa  <  1:95547/138‑1 (MQ=0)
                                                                                                                                   |                                                                                                                                    
TAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACT‑‑TGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTA  >  CP000730/35727‑35988

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: