Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 110,660 T→A 16.2% I173I (ATT→ATA USA300HOU_0109 → staphylococcal tandem lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730110,6600TA16.2% 177.7 / 24.1 73I173I (ATT→ATAUSA300HOU_0109staphylococcal tandem lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (25/36);  new base A (7/5);  total (32/41)
Fisher's exact test for biased strand distribution p-value = 3.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.33e-01

GGTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGTTCAGGGAAT  >  CP000730/110522‑110798
                                                                                                                                          |                                                                                                                                          
ggTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAAAAATGGCTATTATAAAGACGGAGATAttt                                                                                                                                           <  1:660476/140‑1 (MQ=14)
ggTTAAAATAGTCTATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATAttt                                                                                                                                           <  2:633028/140‑1 (MQ=255)
 gTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTc                                                                                                                                          <  1:627446/140‑1 (MQ=255)
   tAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCat                                                                                                                                        <  2:607879/140‑1 (MQ=255)
     aaaaGGTCAATAATAAAATCATTCCAACAGAATAAATTACAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCatat                                                                                                                                      <  2:371993/140‑1 (MQ=17)
        tGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTAAAAAATTATAAAGACAGAGATATTTAATATAAt                                                                                                                                   >  1:336940/1‑140 (MQ=255)
             gATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATTATTCATTTAATCCaat                                                                                                                              >  2:359454/1‑138 (MQ=255)
                 acaaaaGCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTg                                                                                                                          <  2:256923/138‑1 (MQ=255)
                   aaaaataTTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGcc                                                                                                                        <  2:274670/134‑1 (MQ=255)
                   aaaaTCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGcc                                                                                                                        <  1:267580/140‑1 (MQ=255)
                   aaaaTCACTCCAACAAAAGTAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGcc                                                                                                                        <  1:278493/140‑1 (MQ=255)
                    aaaTTATTCTATCAAAAGAAATCAACGATCAAAAAATCAAAAAAGGATTCGAAAACTTTAAGTTCTTTGTTCAATATGGCGATTTTATAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCg                                                                                                                       <  1:15838/140‑1 (MQ=17)
                    aaaTCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCg                                                                                                                       >  1:124396/1‑140 (MQ=255)
                       tCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGt                                                                                                                    <  1:306927/140‑1 (MQ=255)
                                  aaaGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATATGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTaa                                                                                                         >  1:485955/1‑140 (MQ=255)
                                   aaGAAATTAAAGATGAAGAAATAAAAAAAGAAATCGAAAACCTTAATTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTaaa                                                                                                        <  2:328966/140‑1 (MQ=255)
                                   aaGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTaaa                                                                                                        <  2:601839/140‑1 (MQ=255)
                                   aaGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCAGAGTTATTCGGCTaaa                                                                                                        >  2:182931/1‑140 (MQ=255)
                                     gAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATGATCCAGAGGTGCCGAGTTATTCGGCTAAata                                                                                                      <  1:1754/140‑1 (MQ=255)
                                        aTTAAAGATGACAAAATAAAAAAAAAAATCGAAAACTTTAAGTTCTTTGTTCAATATCGCGACATTAAAAATTTGAAAAATTATAACGACGGAGATATTTCAAATAATGAAGAGTAGCCGAGTTATTCGGCtcactagca                                                                                                   >  2:33225/1‑132 (MQ=17)
                                           aaaGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAAtt                                                                                                >  2:97442/1‑140 (MQ=255)
                                           aaaGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAAtt                                                                                                >  2:460542/1‑140 (MQ=255)
                                           aaaGATAAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAAtt                                                                                                <  2:180169/140‑1 (MQ=255)
                                            aagaataaaaaaTAAAAAAAGAAATCTAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTa                                                                                               <  2:321233/134‑1 (MQ=255)
                                             ataTGAAAAAATAAAAAAAGAAATCGATAACTTTAAGTTCTTTGTTCAATATTGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTGATTTTGCTAAATATCAATTaa                                                                                              <  2:144190/138‑1 (MQ=17)
                                             aGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTCGAAAAATCATCCCGACCGCGCTATTCCATATCCTCCAGCGCTCCCCAGTCACTCGGCTAACTCTCaaccaa                                                                                              >  2:139520/1‑136 (MQ=12)
                                               aTGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACt                                                                                            <  1:106827/140‑1 (MQ=255)
                                                agAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTGTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTCCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTa                                                                                           <  2:168978/139‑1 (MQ=255)
                                                  aaaaaaTAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGAGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAat                                                                                         <  1:550340/140‑1 (MQ=255)
                                                     ataTAAAAAAAGAAATCGAAAACTTTAAGTTCTGTGTTCAATATGGCGTCTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatgat                                                                                      <  1:317793/138‑1 (MQ=255)
                                                      aaTAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatgatg                                                                                     <  2:234261/140‑1 (MQ=255)
                                                       aTAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATTTTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAAtgatga                                                                                    <  1:200981/140‑1 (MQ=255)
                                                          aaaaaaGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATta                                                                                 >  2:580704/1‑140 (MQ=255)
                                                             aaaGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGCAAAATTTCATATAATACAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTaaaa                                                                              >  1:155071/1‑140 (MQ=17)
                                                             aaaGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATaa                                                                              >  1:542904/1‑140 (MQ=255)
                                                              aaGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAAt                                                                             >  1:641631/1‑140 (MQ=255)
                                                               aGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATg                                                                            <  1:643105/140‑1 (MQ=255)
                                                                 aacgcgatgaCTTTAAGTTCTTTTTTCAATATTGGGAGTTTAAAAATTTGAAAAATTATAAAGACGGATATATTTCATATAATCCAGAGGGGCCGAGTTATTCGGCTAAATATAAATTAACTAATGATGATTATAATGTa                                                                          <  1:102284/131‑1 (MQ=17)
                                                                 aaaTCGAAAACTTTCAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTGTAATGTa                                                                          <  1:416035/140‑1 (MQ=255)
                                                                 aaaTCGAAAACTTTAAGTTCTTTGTTCAATACGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTa                                                                          <  2:16744/140‑1 (MQ=255)
                                                                            tttAAGTTCTTTGTTCAATATGGTGACTTTAAAAATTTTAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGTTGCCGAGTTATTCGGCTAAATATCAATTAACTAAGGATGATTATAATGTAAAACAATTACg                                                               <  1:205390/140‑1 (MQ=255)
                                                                                atttCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGATGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCaaa                                                           <  1:78548/138‑1 (MQ=255)
                                                                                        gTTCAATATGGCGGCTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAgatatga                                                   >  1:662986/1‑140 (MQ=255)
                                                                                                 ggCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACcgac                                          >  2:591227/1‑140 (MQ=255)
                                                                                                  gCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACcgacg                                         <  2:382504/140‑1 (MQ=255)
                                                                                                       tttAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACATCTATCCAAAAGATATGATATACCGACGAGTaa                                    <  1:159177/140‑1 (MQ=255)
                                                                                                         tAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATAGcgacgaacaaca                                  >  2:327056/1‑134 (MQ=255)
                                                                                                         tAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAg                                  >  2:228524/1‑140 (MQ=255)
                                                                                                          aaagatttaatgaaTTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGc                                 >  2:199647/13‑140 (MQ=11)
                                                                                                          aaaaaTTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGc                                 <  1:255915/140‑1 (MQ=255)
                                                                                                          aaaaaTTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGc                                 >  2:644096/1‑140 (MQ=255)
                                                                                                           aagatttaatgaaTTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCAACGAGTAAAGcc                                >  1:107342/12‑139 (MQ=11)
                                                                                                           aagatttaatgaaTTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATAAAAGACCGACGAGTAAAGCt                                >  2:601539/12‑140 (MQ=11)
                                                                                                            aaaTTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTc                               >  2:82103/1‑140 (MQ=255)
                                                                                                             gatttaatgaaTTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGAATATAATGTAAAACAATTACGCAAAAGAAAAGATATACCGACGAGAAAAGCTcc                              >  1:318467/10‑140 (MQ=11)
                                                                                                             aaTTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGATTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTcc                              <  1:244433/140‑1 (MQ=255)
                                                                                                             aaTTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATAAACCGACGAGTAAAGCTcc                              >  1:285008/1‑140 (MQ=255)
                                                                                                              aTTTGAAAAATCATAAAGACGTGTGTGCTTTATATAATCCAGAAGCGCCTAATTATTCGGCTAAACATCAATTAACTAATGATGATTACAATGTAAAACAATTACGCAAAAGATATGAGATACCGACGAGTAAAGCTCCa                             <  2:513341/140‑1 (MQ=17)
                                                                                                               tttaatgaaTTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCaa                            >  2:509017/8‑140 (MQ=14)
                                                                                                               tttGAAAATTTATTAAGACGGAGATATTTCATCTAATCCAGTGATGCCGATTTATTCGTCTAAATTTCAATTACCCTATTACGATTATGATGTAAAACAACTCCGCATTATGTATGATATAACGATGAGTAGAGCTCCat                            <  2:207072/140‑2 (MQ=12)
                                                                                                                    tgaaTTATAAAGATGGAGATATATCATATAATCCAGAGGTGGCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTa                       <  1:655658/138‑1 (MQ=255)
                                                                                                                     gaaTTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTAt                      <  1:60097/139‑1 (MQ=255)
                                                                                                                     aaaTTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATAGCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTAt                      >  2:335105/1‑140 (MQ=255)
                                                                                                                        ttATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGt                   <  1:527809/140‑1 (MQ=255)
                                                                                                                        ttATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGt                   >  1:548564/1‑140 (MQ=255)
                                                                                                                          ataaaattgaagaTATAACATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAATTTATTGTTa                 <  2:599165/130‑1 (MQ=17)
                                                                                                                            aaacatccaaatatatCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTACTGATGATTATAATGTAAAACAATTACGCAAAAGATATCATATAATGTTTAGTACAGCTCCAAAGTTATTGTTaaa               <  2:610556/130‑1 (MQ=17)
                                                                                                                             aaGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCTAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTaaaa              <  1:534995/140‑1 (MQ=255)
                                                                                                                              aGACGGAGATATTTCATATAATCCAGATGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAg             <  1:510600/140‑1 (MQ=255)
                                                                                                                              aGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGGGTAAAGCTCCAAAGTTATTGTTAAAAg             >  2:445946/1‑140 (MQ=255)
                                                                                                                                 cGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGtt          <  2:124396/140‑1 (MQ=255)
                                                                                                                                  ggAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGTTc         <  1:92868/140‑1 (MQ=255)
                                                                                                                                   gagaTATATCATATAATCCAGAGGTGCCGAGTTATTTGTCTAAGTCTAAAATGACTAGATATTGGTATGATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTAATTGTTAAAAGGTTCa        >  2:303962/1‑140 (MQ=17)
                                                                                                                                     gaTATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGTTCAgg      >  1:459526/1‑140 (MQ=255)
                                                                                                                                      atatTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGTTCAggg     >  2:330405/1‑140 (MQ=255)
                                                                                                                                       tatTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGTTCAGGGa    >  2:334682/1‑140 (MQ=255)
                                                                                                                                         tttCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATTTACCGACGAGTAAAGCTCCAACGTTATTGTTAATAGGTTGGGGGaag  >  1:594006/1‑139 (MQ=18)
                                                                                                                                          |                                                                                                                                          
GGTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGTTCAGGGAAT  >  CP000730/110522‑110798

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