Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 35,850 C→T 40.6% 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,8500CT40.6% 93.1 / 70.1 78intergenic (+81/+27)USA300HOU_0025/USA300HOU_0026hypothetical protein/IS431mec transposase
Reads supporting (aligned to +/- strand):  ref base C (21/25);  new base T (8/24);  total (29/49)
Fisher's exact test for biased strand distribution p-value = 9.50e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATA  >  CP000730/35719‑35956
                                                                                                                                   |                                                                                                          
aaaGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAAttt                                                                                                   <  2:118183/141‑1 (MQ=255)
aaaGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTGGGCGAAAAATGATTACTTGAACAAGTGCTAAAATAAAATTGTGATCTATTAACAATTGTAGGTTTTGATAAAAAAAAAAAAAATATAGCTAACCACTAAttt                                                                                                   <  2:14841/141‑1 (MQ=255)
aaaGTATTTAATGTAGAAGCCATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGAACTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAAttt                                                                                                   <  2:178454/141‑1 (MQ=255)
 aaGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTa                                                                                                  >  1:177984/1‑141 (MQ=255)
      tttAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATg                                                                                             <  1:323201/141‑1 (MQ=255)
      tttAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATg                                                                                             >  2:134912/1‑141 (MQ=255)
      tttAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATg                                                                                             >  2:289885/1‑141 (MQ=255)
       ttAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGt                                                                                            <  2:44254/141‑1 (MQ=255)
          aTGTACAAGCCACTCAGTATAAAGTTTTTGGCCAAAAATGAGTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGTTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAg                                                                                         <  1:48866/141‑1 (MQ=37)
                   cTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCtt                                                                                <  1:207222/141‑1 (MQ=255)
                   cTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCtt                                                                                <  1:176168/141‑1 (MQ=255)
                        cAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACtt                                                                           >  1:61403/1‑141 (MQ=255)
                          ggataAAGTTTTTGGCGGAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGc                                                                         <  2:23588/139‑1 (MQ=255)
                          ggataAAGTTTATGGAGAAAAATGATTACATGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTATGAAGAAAAAAAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGc                                                                         <  2:182102/139‑1 (MQ=25)
                          ataggaaGTTTTTGGCGGAAAATGATGAATTGCAGTAGGGCCAAATTAAAATTGAAAAATAAAAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGc                                                                         <  2:208359/136‑1 (MQ=17)
                           tataAAGTTTTTGGCGAAAAATGATTACTTGAACACGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAACAATATCGCTAACCACTAATTTTTCATGTCAGTGTTCGCTTAACTTGCt                                                                        <  2:280519/141‑1 (MQ=255)
                           tataAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCt                                                                        <  1:99099/141‑1 (MQ=255)
                               aaGTTTTTGGCGAAACATGCTTACTTGAACAAGAGCTAAAATAAAATTGTGATCTACTAAAACCAGAAGGCTCTGTTGCACAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTCTTCGCTTAACTTGCTaccc                                                                    >  2:349042/1‑138 (MQ=21)
                               aaGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAACAATAGAAGGTTCTGTTGCAAAGTAAAAAAATGTAGCTAACCACTAATTTATCATGTCAGTGTTCGCTGAACTTGCTAGCa                                                                    >  2:168274/1‑141 (MQ=255)
                                    tttGGCGAAAAATGATAACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCatgatg                                                               >  1:195544/1‑141 (MQ=255)
                                    caggggggtagaTGATTACTATAACAAGAGCTAAAATAAAATTGTGATCTAATAACAATAGAAGGTTCTGTTGAAAAGTAAAAAAATATAGCTAACCACTAATGTAGCATGTCAGTGTTCGCTTAACTTGCTAGCatgatg                                                               <  2:94336/130‑1 (MQ=18)
                                     ttGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGc                                                              <  1:343172/141‑1 (MQ=255)
                                     taggaaaaaaaTGACTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTAAAAAAAAAAAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGc                                                              <  2:99470/135‑1 (MQ=21)
                                          gAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGCCAGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTGGCTAGCATGATGCTAAtt                                                         >  1:72633/1‑141 (MQ=255)
                                          gAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGACTTCGCGTACCTTGCTAGCATGATGCTAAtt                                                         >  2:239565/1‑141 (MQ=255)
                                            aaaaTGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTc                                                       >  2:50591/1‑141 (MQ=255)
                                                         gaacaaagagtaaaataaaaTTGTAATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGaa                                          <  2:303541/131‑1 (MQ=255)
                                                              agagCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATcc                                     >  1:320738/1‑141 (MQ=255)
                                                              agagCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGCCACTCGTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATcc                                     >  1:358312/1‑141 (MQ=255)
                                                               gagCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAGCTTGCTAGCATGATGCTAATTTCGTGGCCTGGCGAAAATCCg                                    >  1:2429/1‑141 (MQ=255)
                                                               gagCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCg                                    >  2:259048/1‑141 (MQ=255)
                                                                agCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGt                                   <  2:26462/141‑1 (MQ=255)
                                                                       acaaaagagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGaa                          <  1:3457/127‑1 (MQ=11)
                                                                       acaaaagagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGaa                          >  1:28203/15‑141 (MQ=11)
                                                                        taaaatTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCGTGGCGAAAATCCGTACATCTGa                           >  1:256736/1‑141 (MQ=255)
                                                                        caaaagagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAag                         <  1:54206/128‑1 (MQ=11)
                                                                        caaaagagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATAACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAag                         >  1:230165/14‑141 (MQ=11)
                                                                          aaagagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAagag                       <  1:31791/130‑1 (MQ=11)
                                                                           aagagtattttaataaAAATGATGGTGCTTGTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAgaga                      <  1:362686/131‑1 (MQ=11)
                                                                            agagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATTTATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAc                     <  1:50098/132‑1 (MQ=11)
                                                                            agagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAc                     <  1:236434/132‑1 (MQ=11)
                                                                            agagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAc                     <  1:136768/132‑1 (MQ=11)
                                                                            agagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAc                     <  2:53954/132‑1 (MQ=11)
                                                                            agagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAc                     <  1:60447/132‑1 (MQ=11)
                                                                            agagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGACCTGAAGAGAc                     >  1:288450/10‑141 (MQ=11)
                                                                             gagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAcc                    <  1:361002/133‑1 (MQ=11)
                                                                               ttgATTTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAcc                    <  2:37991/140‑1 (MQ=255)
                                                                               gtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTg                  <  1:243031/135‑1 (MQ=11)
                                                                                tattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTCGATCTGAAGAGACCTGc                 <  2:192795/136‑1 (MQ=11)
                                                                                tattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGc                 <  2:196476/136‑1 (MQ=11)
                                                                                tattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGc                 <  1:366229/136‑1 (MQ=11)
                                                                                 taTCTAATAAAAATAGAAGGGTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATGTATCATGTCAGTGTTCGCTTAACTTGCTAGCATTATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTg                  <  1:364514/140‑1 (MQ=255)
                                                                                 attttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCg                >  2:312142/5‑141 (MQ=11)
                                                                                 attttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCg                >  2:258630/5‑141 (MQ=11)
                                                                                  ttttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCgg               <  1:245200/138‑1 (MQ=11)
                                                                                  ttttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCgg               <  1:81635/138‑1 (MQ=11)
                                                                                  ttttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCGAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCgg               <  2:291873/138‑1 (MQ=11)
                                                                                  aTCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGc                 >  1:216593/1‑141 (MQ=255)
                                                                                   tttaataaAAATTATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTCATATATCATGTCAGTGTTCGTTTAACCTGCTAGTTTGATGCTAATTTCGTGGTATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGt              <  1:174426/139‑1 (MQ=11)
                                                                                   tttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGt              <  2:313926/139‑1 (MQ=11)
                                                                                   tttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGt              >  2:312786/3‑141 (MQ=11)
                                                                                   tttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTCCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGt              <  1:162065/139‑1 (MQ=11)
                                                                                   tttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCAATAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGAGCTGAAGAGACCTGCGGt              >  1:75182/3‑141 (MQ=11)
                                                                                   tttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATACAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGt              <  2:208841/139‑1 (MQ=11)
                                                                                   tCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTAAagtggcgcgcg                >  2:57802/1‑132 (MQ=255)
                                                                                    ttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGtt             <  2:68014/140‑1 (MQ=11)
                                                                                     taataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTc            <  2:244793/141‑1 (MQ=11)
                                                                                      aataaAAATGGAAGTTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTGATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGtt             <  2:182580/141‑1 (MQ=255)
                                                                                       ataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGTAGAGACCTGCGGTTCtt          >  2:257491/1‑141 (MQ=11)
                                                                                        taaAAATAGAAGGTTCTGTTGCAAAGTAAACAAATATAGCTAACCACTAATTTATCACGTCCGTGTTCGCTTAACTTGCTAGCAGGAAGCTAATTTCGTGGCAGGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGGTCt           >  2:39149/1‑141 (MQ=11)
                                                                                        taaAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATTTAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCt           >  2:47803/1‑141 (MQ=38)
                                                                                        taaAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCt           >  1:277386/1‑141 (MQ=38)
                                                                                           aaaTAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGCGTTCGCTCAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCtttt        <  1:182349/141‑1 (MQ=25)
                                                                                             aTAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAGCTGCGGTTCTTTTta      <  1:284926/141‑1 (MQ=18)
                                                                                              ccgAAGTTTCTTTTGCAAAGTAAAAAATATAGCTACTCACTAATATATCATGTCAGTGTTCGTTTAGCTTGCTAGCAGGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTtata    <  2:284489/139‑1 (MQ=11)
                                                                                               aGAAGGTTCTGTGGCAAAGTAAAAAAATATAGCTAACAACTGATTTATCATATCAGTGTTCGATTAACTTACTAGCATGATGCTAATTCCGTGGCATGGCGAAAGTCCGTAGATCTGAAGAGACCTGCGGTGCttttgaga    >  2:95131/1‑137 (MQ=11)
                                                                                                 aaGTTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCCAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTtatata  <  2:256738/141‑1 (MQ=11)
                                                                                                 aaGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGTCCTGCGGTTCTTTTtatata  <  1:379500/141‑1 (MQ=11)
                                                                                                                                   |                                                                                                          
AAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATA  >  CP000730/35719‑35956

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: